第6章_常微分方程数值解法教材教学课件.pptVIP

第6章 常微分方程的数值解法;只要 f (x, y) 在[a, b] ? R1 上连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述问题存在唯一解。; 微分方程的数值方法。求的是在一系列离散点列上，未知函数y(x)在这些点上的值的近似。;6.1 欧拉Euler 法与改进欧拉法;用向后差商公式代替导数项;例6.1：考虑初值问题： 取步长h = 0.1，并与准确解 比较。 解：因为xi = 0.1i，而f(x, y) = -y + x，故Euler格式 yi+1=yi+hf(xi, yi) = (1-h)yi + hxi (i=0,1,…) 隐式Euler格式为 yi+1=yi+hf(xi+1, yi+1) = yi – hyi+1+hxi+1 (i=0,1,…) 整理得 yi+1= (yi+hxi+1 ) / (1+h) (i=0,1,…) 用Matlab计算：clear;x=0,y=0; for i=1:5 y=(1-h)*y+h*x x=x+h; end;6.1.2 欧拉格式的改进;类似，可以算出梯形格式的误差估计式：;? 改进欧拉格式 /* modified Euler’s Formula */;注：此法亦称为预测-校正法 /* predictor-corrector method */。可以证明该算法具有 2 阶精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。它的稳定性高于显式欧拉法。;称之为中点公式,或可写为;其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 ; 四阶经典R-K公式; 解 四阶标准R-K公式为;n; 求解初值问题; 定义6.2 设y(x)是初值问题(6.1)的解,yn是单步法(*)产生的近似解.如果对任意固定的点xn,均有;用Euler法得近似解表达式; 一般地，有如下定理;利用Lipschitz条件得; 设?(x,y)连续且关于y满足Lipschitz条件,对于Euler方法,由于?(x,y,h)=?(x,y),故Euler方法是收敛的.;6.3.2 稳定性; 将Euler方法应用于方程y?=?y, 得到 ;类似前面分析,可知绝对稳定区域为; 解 因y0=1,计算得y10=1024,而y(1)=9.357623?10-14.; 单步显式方法的稳定性与步长密切相关, 在一种步长下是稳定的差分公式,取大一点步长就可能是不稳定的.;6.4 线性多步方法（Adams格式）; 例 选取参数?,?0,?1,?2,使三步方法 ; ?=1, ?0+?1+?2=1, ?1+2?2=-1/2 , ?1+4?2=1/3;设pr(x)是函数?(x,y(x))的某个r次插值多项式,则有; ???已求得精确解y(x)在步长为h的等距节点xn-r,…,xn上的近似值yn-r ,…,yn , 记?k=?(xk,yk) ,利用r+1个数据(xn-r,?n-r),…,(xn,?n)构造r次Lagrange插值多项式;下面列出几个带有局部截断误差主项的Adams显式公式 ; 2. Adams隐式格式 ; 3. Adams预估-校正格式 ; 解 用四阶R-K公式提供起始值,计算结果如下;6.5 一阶微分方程组和高阶微分方程;以两个方程为例;经典四阶龙格-库塔格式;Adams外插格式;6.5.2 高阶常微分方程;以二阶方程为例;练习题