常微分方程数值解59514.pptVIP

例：用欧拉法解初值问题 取步长 h = 0.2.计算过程保留4位小数. 解：f(x,y)=－y－xy2 ， h = 0.2，由欧拉公式得 故y(0.2)?y1=0.2×1(4－0×1)＝0.800 0 y(0.4)?y2=0.2×0.8×(4－0.2×0.8)＝0.614 4 y(0.6)?y3=0.2×0.614 4×(4－0.4×0.4613)＝0.800 0 例1：用预报－校正公式求解初值问题 取步长 h = 0.2,计算 y(1.2), y(1.4)的近似值，计算过程保 留5位小数. 解： * 常微分方程数值解 计算出解函数y(x)在一系列节点a=x0x1… xn=b处的近似值yi=y(xi)(i=1,2…n),即数值解。 §1 欧拉方法（ Euler’s Method ） ? 欧拉公式： 向前差商近似导数 x0 x1 )) ( , ( ) ( 0 0 0 1 x y x f h y x y + ? ) 1 , ... , 0 ( ) , ( 1 - = + = + n i y x f h y y i i i i 求解过程顺着节点排列次序一步步向前推进，即按递推公式由已知的y0,y1,…yi，求出yi+1。 y o x0 y(x) p0 x1 p1 p2 x2 x xi Pi+1 Pi xi+1 定义 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差 定义 　　　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 可见: (1)每一步都会产生误差. (2) 与前面每一步计算产生的误差都有关，具有整体性，所以分析和确定它是很复杂的。 ? 欧拉公式求解基本思路 ： ? 欧拉公式的改进： ? 隐式欧拉法 /* implicit Euler method */ 向后差商近似导数 x0 x1 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? ) 1 , ... , 0 ( ) , ( 1 1 1 - = + = + + + n i y x f h y y i i i i 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。 一般先用显式计算一个初值，再迭代求解。 ? 隐式欧拉法的局部截断误差： 即隐式欧拉公式具有 1 阶精度。 ? 梯形公式 /* trapezoid formula */ — 显、隐式两种算法的平均 注：的确有局部截断误差 ， 即梯形公式具有2 阶精度，比欧拉方法有了进步。但注意到该公式是隐式公式，计算时不得不用到迭代法，其迭代收敛性与欧拉公式相似。 ? 中点欧拉公式 /* midpoint formula */ 中心差商近似导数 x0 x2 x1 假设 ，则可以导出 即中点公式具有 2 阶精度。 需要2个初值 y0和 y1来启动递推 过程，这样的算法称为双步法 /* double-step method */，而前面的三种算法都是单步法 /* single-step method */。 方 法 ? ? 显式欧拉 隐式欧拉 梯形公式 中点公式 简单 精度低 稳定性最好 精度低, 计算量大 精度提高 计算量大 精度提高, 显式 多一个初值, 可能影响精度 ? 改进欧拉法 /* modified Euler’s method */ Step 1: 先用显式欧拉公式作预测，算出 ) , ( 1 i i i i y x f h y y + = + Step 2: 再将 代入隐式梯形公式的右边作校正，得到 1 + i y )] , ( ) , ( [ 2 1 1 1 + + + + + = i i i i i i y x f y x f h y y 注：此法亦称为预测-校正法 /* predictor-corrector method */。可以证明该算法具有 2 阶精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。后面将看到，它的稳定性高于显式欧拉法。 §2 龙格 - 库塔法 /* Runge-Kutta Method */ 建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发，以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点