数值计算课件数值分析第6章 常微分方程数值解法教案.pptVIP

将(6―45)、(6―46)代入(6―42)得到近似差分方程为 (6―47) 其中 pi=p(xi)，qi=q(xi)，fi=f(xi) 将式(4―47)整理后，得 其中 (6―49) 这个方程组的系数矩阵是三对角的，可以利用追赶法求解。 根据(6―45)、(6―46)式知，对于y(xi)-yi的误差方程只要把fi取成O(h2)即可。现剩下两个问题：其一是线性方程组(6―48)［即差分方程(6―47) 解的存在唯一性;其二是(6―48)的解的收敛性，也即当h→0时，解是否收敛于微分方程(6―42)的准确解。 下面对差分方程(6―47)给出结论，不加以证明。 定理若ai＞0，ci＞0，-bi≥ai+ci，i=1，2，…，n-1，则二阶差分方程 aiyi-1+biyi+ciyi+1=di y0=α，yn=β，i=1，2，…，n-1 例5 用差分法解边值问题 y″-y=x y(0)=0，y(1)=1，0≤x≤1，h=01 解 这里步长h=1/10，则基点 表 6―5 3.2 龙格库塔法 前面已经知道，初值问题(6―1)等价于 龙格库塔法的基本思想是：用f(x，y)在几个不同点的 数值加权平均来代替f(xi+θh，y(xi+θh))的值，而使截断 误差的阶尽可能高。 1. 二阶龙格库塔公式 将预估校正公式(6―8)改写成更一般的形式 (6―16) 适当选取%α、β、ω1、ω2%的值，使截断误差 y(x i+1)-y i+1的阶数尽可能高。这里仍假定yi=y(xi)，显然 2. 四阶龙格库塔公式 二阶龙格库塔公式是由使用在两个不同点上的函数值的线性组合而得到的。同样，我们用四个不同点上的函数值的线性组合就可得到四阶龙格库塔公式。设 yi+1=yi+h(ω1k1+ω2k2+ω3k3+ω4k4) (6―20) 这里k1、k2、k3、k4为四个不同点上的函数值，分别设其为 k1=f(xi，yi) k2=f(xi+λ1h，yi+μ11k1h) k3=f(xi+λ2h，yi+μ21k1h+μ22k2h) k4=f(xi+λ3h，yi+μ31k1h+μ32k2h+μ33k3h) 其中ω1、ω2、ω3、ω4、λ1、λ2、λ3、μ11、μ21、μ22、μ31、μ32、μ33均为待定系数。 类似于前面的讨论，把k2、k3、k4分别在xi点展成h的幂级数，代入线性组合式(6―20)中，将得到的公式与y(xi+1)在xi点上的泰勒展开式比较，使其两式右端直到h4的系数相等，经过较复杂的运算便可得到关于ωi，λi，μij的一组特解 λ1=λ2=μ11=μ22=1/2 μ21=μ31=μ32=0 λ3=μ33=1 ω1=ω4=1/6 ω2=ω3=1/3 (6―22) 从而得到常用的标准四阶龙格库塔公式： (6―23) 图 6.5 表 6―3 §4 阿达姆斯方法 我们已经知道，初值问题(6―1)等价于积分方程(6―9)，即 对积分式分别采用矩形公式和梯形公式可得到欧拉公式和改进的欧拉公式，截断误差分别为O(h2)和O(h3)。为此，我们自然可以想到，若用更高次的插值多项式来代替f(x，y)，则所得公式的精度会更高。这就是线性多步法的起源思想。 本章前面介绍的方法称为单步法，因为在计算yi+1时，只用到前面yi的值。而对于线性多步法是要利用前面已经算出的若干个值yi-k，…，yi-1，yi来求yi+1。