Ch9 常微分方程初值问题数值解法.pptVIP

定理1表明p≥1时单步法收敛, 并且当y(x)是初值问题(1.1),(1.2)的解, (4.1)具有p阶精度时, 则有展开式 所以p≥1的充分必要条件是 ，而 ，于是可给出如下定义： * 定义4 若单步法(4.1)的增量函数?满足 以上讨论表明p阶方法(4.1)当p≥1时与(1.1), (1.2)相容，反之相容方法至少是1阶的. 于是由定理1可知方法(4.1)收敛的充分必要条件是此方法是相容的. 则称单步法(4.1)与初值问题(1.1),(1.2)相容. * 9.4.2 绝对稳定性与绝对稳定域 前面关于收敛性的讨论有个前提，必须假定数值方法本身的计算是准确的. 实际情形并不是这样，差分方程的求解还会有计算误差. 譬如由于数字舍入而引起的小扰动. 这类小扰动在传播过程中会不会恶性增长，以至于“淹没”了差分方程的“真解”呢？这就是差分方程的稳定性问题. 在实际计算时，我们希望某一步产生的扰动值，在后面的计算中能够被控制，甚至是逐步衰减的. * 定义5 若一种数值方法在节点值yn上大小为δ的扰动，于以后各节点值ym(mn)上产生的偏差均不超过δ，则称该方法是稳定的. 下面以欧拉法为例考察计算稳定性. 例4 用欧拉公式求解初值问题 解 用欧拉法解方程y=-100y 得 其准确解 是一个按指数曲线衰减很快的函数. * 若取步长h=0.025，则欧拉公式的具体形式为 节点xn 欧拉方法yn 后退欧拉方法yn 0.025 0.050 0.075 0.100 -1.5 2.25 -3.375 5.0625 0.2857 0.0816 0.0233 0.0067 计算结果见表, 明显计算过程不稳定, 但取h=0.005, yn+1=0.5yn, 则计算过程稳定. 　　对后退的欧拉公式，取h=0.025时，则计算公式为yn+1=(1/3.5)yn .计算结果见表, 这时计算过程是稳定的. * 例题表明稳定性不但与方法有关，也与步长h有关，当然与方程中的f(x, y)有关. 为了只考察数值方法本身，通常只检验数值方法用于解模型方程的稳定性，模型方程为 其中λ为复数，这个方程分析较简单，对一般方程可以通过局部线性化化为这种形式，例如在(?x, ?y)的邻域，可展开为 * 略去高阶项，再做变换即可得到的形式 u?=λu. 对于m个方程的方程组, 可线性化为y?=Ay, 这里A为m×m雅可比矩阵(?fi/?yj)，若A有m个特征值λ1,λ2,?,λm，其中λi可能是复数，所以，为了使模型方程结果能推广到方程组，方程(4.8)中λ为复数. 为保证微分方程本身的稳定性，还应假定Re(λ)0. 　　下面先研究欧拉方法的稳定性. 模型方程y?=λy的欧拉公式为 * 设在节点yn上有一扰动值εn，它的传播使节点值yn+1产生大小为的扰动值εn+1，假设用yn*=yn+εn，按欧拉公式得出yn+1*=yn+1+εn+1的计算过程不再有新的误差，则扰动值满足 可见扰动值满足原来的差分方程(4.9). 这样，如果差分方程的解是不增长的，即有 则它就是稳定的. 这一论断对于下面将要研究的其它方法同样适用. * 　　显然，为要保证差分方程(4.9)的解是不增长的，只要选取h充分小，使 　　在μ=hλ的复平面上，这是以(-1,0)为圆心，1为半径的单位圆. 称为欧拉法的绝对稳定域，一般情形可由下面定义. 定义6 单步法(4.1)用于解模型方程y?=λy，若得到的解yn+1=E(hλ)yn，满足|E(hλ)|1，则称方法(4.1)是绝对稳定的. 在μ=hλ的平面上, 使|E(hλ)|1的变量围成的区域，称为绝对稳定区域，它与实轴的交称为绝对稳定区间. * 0 - 1 - 2 Re Img 　　对欧拉法E(hλ)=1+hλ，其绝对稳定域为|1+hλ|1,绝对稳定区间为-2λ0,在例5中λ=-100,-2-100h0,即0h2/100=0.02为稳定区间，在例4中取h=0.025，故它是不稳定的，当取h=0.005时它是稳定的. 对二阶R-K方法，解模型方程(4.1)可得到 故 绝对稳定域由|E(hλ)|1得到，于是可得绝对稳定区间为-2hλ0，即0h2/λ. * 这样建立的预测—校正系统通常称为改进的欧拉公式: 或表为下列平均化形式 (2.13) 预测 校正 *