中国传媒大学《MATLAB基础与应用》课件-第9章微分方程问题的求解.pptVIP

将改进的欧拉法推广： ) , ( ) , ( ] [ 1 2 1 2 2 1 1 1 phK y ph x f K y x f K K K h y y i i i i i i + + = = + + = + l l 考察(xi , xi+p)内任一点：xi+p=xi+ph 首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开 Step 2: 将 K2 代入第1式，得到 Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 若取 则得到另一个龙格——库塔法 的计算公式： ) , ( ) , ( 1 2 1 2 1 K y x f K y x f K K h y y i i+ i i i i + = = + = + 1 2 h 2 此即中点方式。此公式称作变形的欧拉公式。 为了进一步提高精度，设除xi+p外再考虑一点： xi+q=xi+qh (p≤q≤1) 用xi , xi+p xi+q点的斜率值K1、K2、K3 加权平均作为平均斜率的近似值，这时计算格式具有形式： yi+1= yi+h(λ1K1+λ2K2+λ3K3 ) 其中K1、K2仍用前述所取的形式。 用K1、K2的加权平均作为(xi,xi+q)上的平均斜率，得到y(xi+q)的预报值，再调用f 产生斜率值K3： yi+q= yi+qh(rK1+sK2) K3= f(xi+q,yi+q)=f ( xi+q , yi+qh(rK1+sK2) ) 仿二阶情形的推导，只要把K1,K2,K3代入yn+1的表达式中在( xn,yn)处作泰勒展开, 再与y(xn+1)在xn处的泰勒展开式比较,欲使公式的局部截断误差=O(h4),可得含参数的方程组,类似知有无穷多解；得到的公式称为三阶龙格——库塔公式。 一个比较简单而又重要的三阶龙格——库塔公式，称为库塔公式，其为： yi+1= yi+ (K1+4K2+K3 ) ) , ( ) , ( 1 2 1 K y x f K y x f K i i+ i i + = = 1 2 h 2 K3 = f ( xi+1 , yi + h(-K1+2K2) ) 继续上述这个过程,可以进一步讨论四阶龙格——库塔公式。 ? 最常用的是4阶经典龙格-库塔法： h 6 9.2 线性多步法 用若干节点处的 y 及 y’ 值的线性组合来近似y(xi+1)。 ) ... ( ... 1 1 0 1 1 1 1 0 1 k i k i i i k i k i i i f f f f h y y y y - - + - - - + + + + + + + + + = b b b b a a a 其通式可写为： ? 基于数值积分的构造法 将 在 上积分，得到 只要近似地算出右边的积分 ，则可通过 近似y(xi+1) 。而选用不同近似式 Ik，可得到不同的计算公式。 ? 亚当姆斯显式公式 利用k+1 个节点上的被积函数值 构造 k 阶牛顿后插多项式 , 有 /* 显式计算公式 */ 局部截断误差为： 例：k=1 时有 ? 亚当姆斯隐式公式 利用k+1 个节点上的被积函数值 fi+1 , fi , …, fi?k+1 构造 k 阶牛顿前插多项式。与显式多项式完全类似地可得到一系列隐式公式，并有 ，其中 与 fi+1 , fi , …, fi?k+1 的系数亦可查表得到。 ~ ~ 1 0 1 2 3 k fi+1 f