数值分析课件14.pptVIP

* * 7.3.2 二阶龙格-库塔方法 原来改进欧拉公式k2为xi+1点的近似平均值, 现对k2重新设置取 内任一点 使用xi和xi+l两点斜率值k1,k2加权平均作为 , 即 其中 为待定系数 k1仍为x1上的斜率值，即 问题是如何预测xi+l处的斜率k2. 类似于改进的欧拉方法，先利用欧拉公式得到斜率值 所以计算公式如下: 其中 为调节参数, 使其具有2阶精度 考虑截断误差 对 做二阶泰勒展开 再使用yi+l, 通过f(x,y)形式, 得到斜率值 对 作为二元函数1阶泰勒展开 其中 ，余项是关于 的高阶无穷小. 对节点 即 ，此时 对 ① 利用复合函数求导：对 其中 有 可对①式关于x求导： ② 将①和②带入整理得： 要使 和 为一族系数与l有关． 当 时, 时, 即为改进欧拉公式 ; 当 该公式为变形欧拉公式, 即 特殊地 ①在区间 上取两点，预测相应斜率值 ②对 作加权得平均斜率近似值 ③对截断误差公式，要达到 ，则 的系数须为0, 从而建立方程组求系数. 总结二阶龙格-库塔构造方程: 7.3.3 高阶龙格-库塔公式： 为进一步提高精度（如三阶精度），在[xi, xi+l]上除xi和xi+l 外,可再增加一点 以xi, xi+l, xi+m三点处的斜率值k1,k2,k3加权平均, 计算公式： 其中k1,k2仍为前面定义，即 需要求K3, 可在 上利用二阶龙格-库塔公式 得到预测值后，再利用f(x,y)得到斜率值 整个计算公式为： 需求 使其具有3阶精度. 可利用截断误差公式求解，具体略, 可得条件: 该方程组得到的是一族解， 常用的三阶龙格库塔公式有： 同样地，精度还可以提高到4阶， 具体略。 在计算yi+1之前，已经求出y0、y1、……yi. 是否可以用前面已算好 的yi，得到高精度的yi+1的值？ 7.4线性多步法 该方法为多步法，（前面方法为单步法, 只涉及到yi、yi+1） 线性r步计算公式可写成： ① yi+1=yi+hf(xi,yi), yi+1=yi+ [f(xi,y(xi))+f(xi+1,y(xi+1))]， 其中, 例如欧拉公式： 其中, 梯形公式： 式①中，αj, βj为常数。 当β-1=0时为显式公式, β-1≠0时为隐式公式. Ri+1=y(xi+1)- 常微分方程数值解常用方法： 利用数值积分方法（另外还有微分中值定理） y(xi+1)=y(xi)+ 对 提高精度（原来只使用一点、两点插值）. 7.4.1阿当姆斯内插公式（Admas） 使用xi,xi+1作插值节点，节点少，精度不高， 可增加节点插值（但节点数不能太多） 误差余项: , 采用更好的节点插值方法构造， ∴需增加节点, 若先取［xi,xi+1］内部点，函数值未知（可通过龙格-库塔方法 计算，但相对复杂）可选取［xi,xi+1］外部的已知节点. 阿当姆斯内插公式： 选取节点为xi-2, xi-1, xi, xi+1来计算y(xi+1)。 构造三次插值多项式代替 ，近似计算 ∴被积函数表示为 ∈(xi-2,xi+1) ∴ （利用积分第一中值定理: x∈[xi,xi+1]内 不变号） 多项式函数可求出积分值（可通过换元法） ∴前两项可作为近似计算公式，最后一项作为误差余项。 具体化简为： ∴ 内插公式为: 是关于yi+1的隐式公式，截断误差为4阶. 需计算xi,xi-1,xi-2三点函数值，初值点y0,y1,y2称为三步法。 （可通过欧拉方法或龙格-库塔方法计算） 7.4.2 阿当姆斯外推公式 将隐式转为显式公式，即将节点xi+1转为xi-3, 取xi-3，xi-2，xi-1，xi为插值节点，求y(xi+1). 采用插值多项式，类似地得到 其误差为: ∴外推公式： 该公式是显式公式，截断误差为4阶. 需要节点值 即初值需 ，可通过欧拉 方法或龙格—库塔方法求得，∴称为四步法。 比较内插和外推公式： 内插公式：隐式（计算复杂），但误差相对较小（余项系数较小）， 且计算公式（三步法）中的4个系数相对较小，产生的 误差影响小。 外推公式：显式（计算相对简单），但误差相对较大（余项系数 较大），且计算（四步法）公式中的4个系数相对较大， 误差影响大。 ∴可将两公式联合使用，构成阿当姆斯校正公式： 即先通过外推公式求得 预测值, 代入到内插公式中： 再将 7.5 一阶方程组和高阶方程7.5.1 一阶方程组 原来一阶方