微分方程数值解法

1、微分方程数值解法论文常微分方程初值问题数值解法的讨论与比较一、一阶常微分方程的初值问题科学技术中常常要求解常微分方程的定解问题，这类问题最简单的形式是一阶方程的初值问题我们知道，只要函数适当光滑，譬如关于满足Lipschitz条件 （1.3）理论上就可以保证初值问题(1.1),(1.2)的解存在并且唯一。虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法。二、欧拉法我们知道，在平面上，微分方程(1.1)的解称作它的积分曲线。积分曲线上一点的切线斜率等于函数的值，如果按函数在平面上建立一个方向场，那么，积分曲线上每一点的切

2、线方向均与方向场在该点的方向相一致，基于上述几何解释，从初始点出发，先依方向场在该点的方向推进到上一点，然后再从依方向场的方向推进到上一点，循此前进推出一条折线，一般地，设已做出该折线的顶点，过依方向场的方向再推进到，显然两个顶点，的坐标有关系即(2.1)这就是著名的欧拉(Euler)公式。若初值已知，则依公式(2.1)可逐步算出例1 求解初值问题解 欧拉公式的具体形式为取步长，计算结果如下表：表1 计算结果对比0.10.20.30.40.51.10001.19181.27741.35821.43511.09541.18321.26491.34161.41420.60.70.80.91.01.

3、50901.58031.64981.71781.78481.48321.54921.61251.67331.7321初值问题(2.2)有解，按这个解析式子算出的准确值同近似值一起列在表1， 两者相比较可以看出欧拉方法的精度很差。三、改进欧拉法为得到比欧拉法精度高的计算公式，在等式如果对方程(1.1)从到积分,得(3.1)右端积分中若用梯形求积公式近似，并用代替，代替，则得 (3.2)称为改进欧拉法改进欧拉方法是隐式单步法，可用迭代法求解用欧拉方法提供迭代初值，则改进欧拉法的迭代公式为(3.3)为了分析迭代过程的收敛性，将(3.1)式与(3.2)相减，得，于是有，式中为对满足Lipschitz常

4、数，如果选取充分小，使得，则当时有，这说明迭代过程(3.3)是收敛的例2用改进的欧拉方法求解初值问题(1.1)解改进的欧拉公式为仍取，计算结果见下表同例中欧拉法的计算结果比较，改进欧拉法明显改善了精度表2 计算结果对比0.10.20.30.40.51.09591.18411.26621.34341.41641.09541.18321.26491.34161.41420.60.70.80.91.01.48601.55251.61531.67821.73791.48321.54921.61251.67331.7321四、 线性多步法在逐步推进的求解过程中，计算之前事实上已经求出了一系列的近似值，如

5、果充分利用前面多步的信息来预测，则可以期望会获得较高的精度这就是构造所谓线性多步法的基本思想构造多步法的主要途径是基于数值积分方法和基于泰勒展开方法，前者可直接由方程(1.1)两端积分后利用插值求积公式得到一般的线性多步法公式可表示为， (4.1)其中为的近似，为常数，及不全为零，则称(4.1)为线性步法，计算时需先给出前面个近似值，再由（4.1)逐次求出如果，称(4.1)为显式步尖，这时可直接由(4.1)算出；如果，则(4.1)称为隐式步法，求解时与改进欧拉法相同，要用迭代法方可算出，(4.1)中系数及可根据方法的局部截断误差及阶确定，其定义为：设是初值问题(1.1),(1.2)的准确解。阿

6、当姆斯显式与隐式公式考虑形如，(4.2)的步法，称为阿当姆斯(Admas)方法为显式方法，亦称Adams-Bashforth公式；为隐式方法，亦称Adams-Monlton公式，直接由方程(1.1)两端从到积分求得。例3积分法导出三步显式阿当姆斯方法解对(1.1)两端从到积分考虑以为节点的二次代数多项式插值，则从而得到同理得再由积分中值定理也不难得到，即的阿当姆斯显式公式为，(4.3)其局部截断误差为例4用四阶阿当姆斯显式和隐式方法解初值问题取步长解本题从四阶阿当姆斯显式公式得到对于四阶阿当姆斯隐式公式得到由此可直接解出而不用迭代，得到计算结果如表5，其中显式方法中的及隐式方法中的均用准确解计

7、算得到，对一般方程，可用四阶R-K方法计算初始近似表5例4计算结果阿当姆斯显式公式阿当姆斯隐式格式0.30.40.50.60.70.80.91.01.040818221.070320051.106530661.148811641.196585301.249328961.306569661.367879441.0700322921.106535481.148818411.196593401.249338161.306579621.367889962.87e-0064.82e-0066.77 e-0068.10 e-0069.20 e-0069.96 e-0061.05 e-0061.040818011.070319661.106530141.148811011.196584591.249328191.306568841.367878592.1 e-0073.9 e-0075.2 e-0076.3 e-