求解刚性常微分方程的两类函数拟合Rosenbrock方法

发布时间：2020-11-05 19:22

　　 随着对刚性常微分方程初值问题的不断深入研究,国内外学者给出了许多较为有效的特殊Runge-Kutta方法,主要包括对角隐式Runge-Kutta方法和Rosenbrock方法。函数拟合方法是一类在局部区间上用一个有理函数来近似地表示刚性常微分方程的解的方法,考虑在其解区间上构造指数拟合的函数,来使得其近似逼近原方程的解曲线,其中比较有效和精确的算法是将Runge-Kutta方法与指数拟合相结合来处理刚性问题。同时,对于一类解可由集合{eiωx,e-iωx}或由集合{cos(ωx),sin(ωx)}(ω0)线性组合的一阶常微分方程初值问题,可将三角拟合思想应用于Runge-Kutta方法上,来得到一类相比传统方法来说更具优势的新方法。学者们将对角隐式Runge-Kutta方法结合指数拟合和三角拟合的研究已经做了较多的工作,而未见采用指数拟合和三角拟合Rosenbrock方法,且Rosenbrock方法相对于对角隐式Runge-Kutta方法来说,具有更小的计算量,故本文将针对常微分方程初值问题这一模型,利用Rosenbrock方法结合指数拟合和三角拟合思想来得到一类在误差和精度上具有更好的表现的方法。在第一章,介绍了关于刚性常微分方程初值问题的研究背景和国内外现状,且给出了早期学者们对Rosenbrock方法的发展与改进。在第二章和第三章,构造了一类二级2阶指数拟合和三角拟合Rosenbrock方法,得到了定系数的二级2阶指数拟合和三角拟合Rosenbrock方法的具体格式,并证明了不存在此类三级3阶指数拟合和三角拟合Rosenbrock方法。最后验证了定系数的二级2阶指数拟合和三角拟合Rosenbrock方法是A-稳定的。在第四章,利用数值实验验证了指数拟合Rosenbrock方法的有效性,主要应用了几组不同的数值实验来比较方法的收敛性和计算时间,通过实验结果得到与理论基本一致的结论。在最后一章,主要是对本文做一个总结,并针对本文未涉及到的其它模型和方法给出了笔者后期的一些想法与计划。
【学位单位】：华中科技大学
【学位级别】：硕士
【学位年份】：2018
【中图分类】：O241.81
【文章目录】：
摘要
Abstract
1 绪论
    1.1 刚性问题的研究背景
    1.2 指数拟合方法的国内外研究现状
    1.3 Rosenbrock方法的国内外研究现状
    1.4 本文研究内容
2 指数拟合Rosenbrock方法
    2.1 指数拟合Rosenbrock方法的构造
    2.2 指数拟合Rosenbrock方法
    2.3 稳定性分析
3 三角拟合Rosenbrock方法
    3.1 三角拟合Rosenbrock方法的构造
    3.2 三角拟合Rosenbrock方法及其稳定性
4 数值算例
5 总结与展望
致谢
参考文献

【参考文献】

相关期刊论文 前7条

1 张静静;吴旭;;一类精确指数拟合的Runge-Kutta方法[J];河南理工大学学报(自然科学版);2012年05期

2 陈全发;肖爱国;;一类A-稳定对角隐式Runge-Kutta法的指数拟合[J];吉林大学学报(理学版);2008年06期

3 杨逢建,张爱国,陈新明;怎样利用函数exp(q)的含参有理逼近构造高阶指数拟合方法[J];应用数学和力学;1999年09期

4 阮保庚;指数拟合的Runge-Kutta公式[J];数值计算与计算机应用;1999年02期

5 阮保庚;高阶双参量指数拟合方法[J];数学物理学报;1996年02期

6 阮保庚;高阶指数拟合公式及其稳定性[J];数值计算与计算机应用;1995年04期

7 向开理，张光裕;两类求解刚性常微分方程的指数拟合法[J];数学杂志;1995年03期

相关硕士学位论文 前3条

1 陈蓉;求解刚性振荡问题的Rosenbrock方法[D];湘潭大学;2011年

2 李聪颖;对角隐式龙格—库塔法及其指数拟合[D];湘潭大学;2007年

3 文志武;几类刚性微分方程Rosenbrock方法的定量误差分析[D];湘潭大学;2005年

本文编号：2872072