研究论文

算法927：两点边值问题数值解的MATLAB代码bvptwp.m

作者:

J.R.公司。现金,

D。霍利沃特,

F、。马齐亚,

上午。纳吉作者信息和声明

ACM数学软件交易（TOMS）,体积39,问题2

条款编号：15，页数1-12

https://doi.org/10.1145/2427023.2427032

出版:2013年2月1日出版历史

获取访问权限

摘要

在本文中，我们描述了用于求解两点边值问题的MATLAB代码bvptwp.m。此代码基于著名的Fortran代码twpbvp.f、twpbvpl.f和acdc.f，这些代码采用了基于本地错误估计和这些代码修订版（称为twpbvpc.f、twbbvplc.f和acdcc.f）的网格选择策略，该算法采用基于局部误差估计和表征问题条件的两个参数估计的网格选择策略。代码twpbvp.f/tpbvpc.f使用基于Mono-Implicit Runge-Kutta方法（MIRK）的延迟校正方案；其他代码使用基于Lobatto公式的延迟校正方案。acdc.f/acdcc.f代码实现了自动延续策略。通过执行一些数值测试来检查新求解器的性能和特征，以表明新代码是健壮的，能够解决非常困难的奇异摄动问题。得到的结果表明，bvptwp.m通常能够解决需要严格精度的问题和解决方案变化很大的问题。该代码与现有的边界值代码（如bvp4c.m）相结合，使得MATLAB BVP部分对于一系列非常广泛的问题来说是一个非常强大的部分。

补充材料

PDF文件（a15-cash_appendix.pdf）

证据在电子附录中给出，可在ACM数字图书馆在线获取。

下载
114.26 KB

ZIP文件（927.拉链）

两点边值问题数值解的MATLAB代码bvptwp.m软件

下载
1.36 MB

工具书类

[1]

阿舍尔，U。、克里斯蒂安森，J。和罗素，R。D。1981年。边界值节点的配置软件。ACM事务处理。数学。柔和。7, 2, 209--222.

数字图书馆

[2]

阿舍尔·U·M、马特海伊·R·M·M和罗素·R·D·1995。常微分方程边值问题的数值解法《应用数学经典丛书》，第13卷，宾夕法尼亚州费城SIAM。（1988年原版的更正重印版）

[3]

Bashir-Ali，Z.、Cash，J.R.和Silva，H.H.M.1998年。洛巴托延迟了刚性两点边值问题的修正。计算。数学。申请。36, 10--12, 59--69.

[4]

Bruganano，L.、Mazzia，F.和Trigiante，D.，2011年。五十年的僵硬。在计算和应用数学的最新进展斯普林格，1-22。

[5]

Cash，J.R.1988。关于使用迭代延迟修正的非线性两点边值问题的数值积分。二、。开发和分析高度稳定的延迟修正公式。SIAM J.数字。分析。25, 4, 862--882.

数字图书馆

[6]

Cash，J.R.和Mazzia，F.，2005年。一种新的基于条件的两点边值代码网格选择算法。J.计算。申请。数学。184, 2, 362--381.

数字图书馆

[7]

Cash，J.R.和Mazzia，F.2006a。两点边值问题的求解算法。网址：http://www.ma.ic.ac.uk/~jcash/BVP_software/twpbvp.php。

[8]

Cash，J.R.和Mazzia，F.，2006年b。基于条件的两点边值问题混合网格选择算法。J.数字。分析。Ind.申请。数学。1, 1, 81--90.

[9]

Cash，J.R.和Mazzia，F.，2009年。两点边值问题的条件和混合网格选择算法。可扩展计算。实践。专家。10, 4, 347--361.

[10]

Cash，J.R.和Wright，M.H.，1990年。两点边值问题求解非线性方程组的实现问题。计算45, 1, 17--37.

数字图书馆

[11]

Cash，J.R.和Wright，M.H.，1991年。非线性两点边值问题的延迟修正方法：实现和数值计算。SIAM J.科学。统计师。计算。12, 4, 971--989.

数字图书馆

[12]

Cash，J.R.、Moore，G.和Wright，R.W.，1995年。一种求解奇异摄动线性两点边值问题的自动延拓策略。J.计算。物理学。122, 2, 266--279.

数字图书馆

[13]

Cash，J.R.、Moore，G.和Wright，R.，2001年。一种求解奇摄动非线性边值问题的自动延拓策略。ACM事务处理。数学。柔和。27, 2, 245--266.

数字图书馆

[14]

Cash，J.R.、Moore，D.R.、Sumarti，N.和Daele，M.V.，2003年。非线性两点边值问题系统的一种带插值的高稳定延迟校正格式。J.计算。申请。数学。155, 2, 339--358.

数字图书馆

[15]

Díaz，J.C.、Fairweather，G.和Keast，P.1983。算法603。COLROW和ARCECO:Fortran包，用于通过修改交替行和列消除来求解某些几乎块对角线性系统。ACM事务处理。数学。柔和。9, 3, 376--380.

数字图书馆

[16]

Enright，W.E.和Muir，P.H.，1986年。两点边值问题的有效Runge-Kutta方法类。计算37, 4, 315--334.

数字图书馆

[17]

Enright，W.E.和Muir，P.H.，1996年。Runge-Kutta软件，用于边界值节点的缺陷控制。SIAM J.科学。计算。17, 479--497.

数字图书馆

[18]

福克斯，L.1957。常微分方程两点边值问题的数值解牛津大学出版社，纽约。

[19]

Kierzenka，J.和Shampine，L.F.2001。一种基于残差控制和MATLAB pse的BVP求解器。ACM事务处理。数学。柔和。27, 3, 299--316.

数字图书馆

[20]

Kierzenka，J.和Shampine，L.F.2008年。控制残差和误差的BVP解算器。J.数字。分析。Ind.申请。数学。三, 1--2, 27--41.

[21]

林德伯格，B.1980。离散化算法的误差估计和迭代改进。位20, 4, 486--500.

数字图书馆

[22]

Mathworks，T.2012。Matlab版本2012a。http://www.mathworks.com/。

[23]

Mazzia，F.和Magerini，C.2008年。初始值问题解决程序测试集，2.4版。巴里大学数学系和巴里研究所INdAM。网址：http://www.dm.uniba.it/~测试集。

[24]

Shampine，L.F.、Gladwell，I.和Thompson，S.，2003年。用MATLAB求解常微分方程剑桥大学出版社，英国剑桥。

数字图书馆

[25]

Shampine，L.F.、Muir，P.H.和Xu，H.，2006年。用户友好的Fortran BVP解算器。J.数字。分析。Ind.申请。数学。1, 2, 201--217.

[26]

Skeel，R.D.1982年。证明延迟修正准确度结果的理论框架。SIAM J.数字。分析。19, 1, 171--196.

数字图书馆

[27]

Wright，R.、Cash，J.R.和Moore，G.1994年。刚性两点边值问题的网格选择。数字。阿尔戈。7, 2--4, 205--224.

引用人

佐丹奴D阿莫迪奥P伊韦纳罗F马齐亚F瓦恩·P苏黎世M(2024)范德瓦尔斯气体自引力等温球的流体静力学流体物理学10.1063/5.020633436:5在线发布日期：2024年5月23日
https://doi.org/10.1063/5.0206334
达拉南G埃利斯A(2024)非对称拍卖：扰动、ε-平衡和平衡游戏与经济行为2016年10月10日/j.geb.2024.06.003147(1-18)在线发布日期：2024年9月
https://doi.org/10.1016/j.geb.2024.06.003
Settanni G公司(2024)HOFiD_bvp码在求解不同类型二阶边值问题中的潜力应用数值数学2016年10月10日/j.apnum.2023.08.008200(379-388)网上发布日期：2024年6月
https://doi.org/10.1016/j.apnum.2023.08.008
显示更多引用者

索引术语

算法927：两点边值问题数值解的MATLAB代码bvptwp.m
1. 计算数学
  1. 数学分析
    1. 微分方程
      1. 常微分方程

建议

具有全局误差和缺陷控制的Runge-Kutta BVODE求解器

边界值常微分方程（BVODEs）是在两个或多个不同的点上施加边界条件的常微分方程组。BVODE数值解的全局误差（GE）是数值解的差异量。。。
基于条件的两点边值码网格选择新算法

我们提出了一种混合网格选择策略，用于两点边值问题的数值求解。这种新的网格策略基于两个参数的估计，这两个参数表征了连续网格的调节。。。
奇摄动非线性边值问题解的自动延拓策略

在最近的一篇论文中，作者导出了求解线性奇异摄动问题的自动延拓算法。该算法被合并到求解边值问题的两个通用代码中，并被证明是有效的。。。

评论

信息和贡献者

问询处

发布于

数学软件上的封面图像ACM事务

ACM数学软件汇刊第39卷第2期

2013年2月

151页

国际标准编号：0098-3500

EISSN公司：1557-7295

内政部：10.1145/2427023

版权所有©2013 ACM。

如果复制品不是为了盈利或商业利益而制作或分发的，并且复制品的第一页载有本通知和完整引文，则允许免费制作本作品的全部或部分数字或硬拷贝以供个人或课堂使用。必须尊重ACM以外的其他人对本作品组成部分的版权。允许用信用证进行摘要。要以其他方式复制或重新发布，在服务器上发布或重新发布到列表，需要事先获得特定许可和/或付费。从请求权限[电子邮件保护]

出版商

计算机协会

美国纽约州纽约市

出版历史

出版：2013年2月1日

认可的：2012年5月1日

修订过的：2012年5月1日

收到：2011年4月1日

在TOMS中发布体积39,问题2

权限

请求对此文章的权限。

检查更新

徽章

作者标记

限定符

研究文章
研究
推荐

贡献者

其他指标

查看文章指标

文献计量学和引文

文献计量学

文章指标

19
引文总数
查看引文
962
总下载次数

下载次数（过去12个月）23
下载次数（最近6周）0

反映截至2024年9月20日的下载量

其他指标

查看作者指标

引文

引用人

佐丹奴D阿莫迪奥P伊韦纳罗F马齐亚F瓦恩·P苏黎世M(2024)范德瓦尔斯气体自引力等温球的流体静力学流体物理学10.1063/5.020633436:5在线发布日期：2024年5月23日
https://doi.org/10.1063/5.0206334
达拉南G埃利斯A(2024)非对称拍卖：扰动、ε-平衡和平衡游戏与经济行为2016年10月10日/j.geb.2024.06.003147(1-18)在线发布日期：2024年9月
https://doi.org/10.1016/j.geb.2024.06.003
Settanni G公司(2024)HOFiD_bvp码在求解不同类型二阶边值问题中的潜力应用数值数学2016年10月10日/j.apnum.2023.08.008200(379-388)网上发布日期：2024年6月
https://doi.org/10.1016/j.apnum.2023.08.008
拉莫斯H辛格G(2022)用三阶导数混合块积分器精确求解二阶两点边值问题应用数学与计算2016年10月10日/j.amc.2022.126960421(126960)在线发布日期：2022年5月
https://doi.org/10.1016/j.ac.2022.126960
马齐亚F(2022)边值问题TOM代码的网格选择策略费拉拉大学ANNALI DELL10.1007/s11565-022-00440-768:2(555-574)在线发布日期：2022年8月30日
https://doi.org/10.1007/s11565-022-00440-7
德马里尼斯A伊韦纳罗F马齐亚F(2022)无人机最小时间避障路径规划算法数值算法10.1007/s11075-021-01167-w89:4(1639-1661)在线发布日期：2022年4月1日
https://dl.acm.org/doi/10.1007/s11075-021-01167-w
普卡拉泽五世罗杰斯S(2021)内部点质量驱动滚球机器人的数值模拟数值代数、控制与优化10.3934/naco.202002111:2(143)在线发布日期：2021年
https://doi.org/10.3934/naco.2020021
马齐亚FSettanni G公司(2021)求解最优控制问题的BVP码数学10.3390/天92026189:20(2618)在线发布日期：2021年10月17日
https://doi.org/10.3390/math9202618
纳赛尔N凯尼西M罗氏J(2020)自适应非重叠区域分解法求解非线性微分边值问题适用分析10.1080/00036811.2020.1800649101:6(2044-2065)在线发布日期：2020年8月6日
https://doi.org/101080/00036811.2020.1800649
武拉诺维奇Nhan T公司(2020)奇异摄动对流扩散问题的高阶鲁棒混合格式应用数学与计算2016年10月10日/j.amc.2020.125495386(125495)在线发布日期：2020年12月
https://doi.org/10.1016/j.amc.2020.125495
显示更多引用者

视图选项

获取访问权限

登录选项

检查您是否可以通过登录凭据或您的机构访问本文。

完全访问权限

获取此文章

查看选项

PDF格式

以PDF文件的形式查看或下载。

电子阅读器

使用联机查看电子阅读器.

电子阅读器

媒体

数字

其他

桌子

查看问题目录