第7章节 常微分方程1.pptVIP

第7章 常微分方程的数值解法 §7.1 引言 一阶常微分方程初值问题是 其中f(x,y)是给定的已知函数。y=y(x)是未知函数， y(x0)=y0是 初值条件。 Th7.1 设f是在区域D={(x,y)|a≤x≤b,y∈R}上的连续函数，其中，a,b∈R，f对y满足Lipschitz条件，即存在常数L0，使得对任意x∈[a,b]及y1 ,y2 ∈R ，有 | f(x,y1) - f(x,y2) |≤L| y1 - y2 | 则对于任意的x0∈[a,b]， y0 ∈R，初值问题（*）在[a,b]存在唯一的连续解。 数值解法: 计算精确解y=y(x)在区间 [a,b]上一系列离散节点 a≤x0≤ x1 ≤…≤ xn ≤b 处的近似值 y0 ,y1 ,…,yn 一般的对区间[a,b]做等距离分割，即取等距离节点xi=a+ih(i=0,1,…,n)，h=(b-a)/n。 求解常微分方程初值问题采用的是“步进式”: 从已知的y0出发（ 由初始条件给定），通过某一计算公式算出y1 ，再算出y2，依此类推。 若在计算yn时仅用到数据yn-1 ，称这种方法为单步法； 若在计算yn时，不仅用到数据yn-1 ，还用到数据yn-2, yn-3 ,…, yn-k+1 (k≥2），称这种方法为多步法。 §7.2 简单的数值方法 7.2.1 欧拉方法 /* Euler’s Method */ ? 显式Euler公式： 向前差商近似导数 记为 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ ? 隐式Euler公式： 记为 向后差商近似导数 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。 ) 1 , ... , 0 ( ) , ( 1 1 1 - = + = + + + n i y x f h y y i i i i y1 y0 y2 x2 x1 x0 一般的，单步法可以表示为 7.2.2 单步法的局部截断误差和阶 其中Φ称为增量函数。 当Φ和yi+1无关时，（*）式称为显示单步法。 当Φ和yi+1 相关时，（*）式称为隐式单步法。 整体截断误差 ei=y(xi)-yi (计算复杂) 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ti = y(xi+1) ? yi+1 称为局部截断误差 /* local truncation error */。 定义 　　　若某算法的局部截断误差O(hp+1)，即 Ti = y(xi+1) ? yi+1= O(hp+1) 则称该算法有p 阶精度。 定义 欧拉法的局部截断误差： Ri 的主项 /* leading term */ 欧拉法具有 1 阶精度。 ? 隐式欧拉法的局部截断误差： 即隐式欧拉公式具有 1 阶精度。 7.2.3 梯形公式 /* trapezoid formula */ 思想 ：显、隐式两种算法的平均值。 故梯形公式是2阶隐式单步法。 梯形公式的局部截断误差 7.2.4 改进Euler法 /* modified Euler’s method */ Step 1: 先用显式Euler公式作预测，算出 ) , ( 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 故梯形公式是2阶隐式单步法。 改进的Euler公式的局部截断误差 y1 y0 y2 x2 x1 x0 例1 取步长h=0.1，分别用Euler法，改进的Euler法求初值问题 解：Euler法的公式为 改进的Euler公式为 精确解为 xn Euler法yn 改进欧拉法yn 精确解y(xn) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1 1.1 1.1918 1.2774 1.3582 1.4351 1.509 1.5803 1.6498 1.7178 1.7848 1 1.0959 1.1841 1.2662 1.3434 1.4164 1.486 1.5525 1.6165 1.6782 1.7279 1 1.0954