数值分析-第7章常微分方程数值解法.pptVIP

第7章 常微分方程数值解法 内容概要 常微分方程的初值问题 欧拉方法 龙格-库塔方法 线性多步法 7．1 常微分方程的初值问题 7．2 欧拉方法 1. 欧拉公式 1. 欧拉公式 欧拉公式的精度 欧拉公式的精度 2. 欧拉公式的MATLAB实现 7．2．2 后退的欧拉公式 后退欧拉公式的求解 后退欧拉公式的精度 7．2．3 梯形公式 梯形方法的几何意义 梯形方法的精度 梯形方法的求解 7．3 龙格-库塔方法 7．3．1 泰勒级数法 泰勒级数法的MATLAB实现（四阶泰勒法） 7．3．2 龙格-库塔方法的基本思想 平均斜率 阶龙格-库塔法 7．3．3 二阶龙格-库塔方法 7．3．3 二阶龙格-库塔方法 7．3．3 二阶龙格-库塔方法 7．3．4 四阶龙格-库塔法 四阶龙格-库塔法的MATLAB实现 程序7-3 function R=rk4(f, a, b, ya, N) %y=f(x,y) , a, b左右端点。 %N为迭代步数, h为步长。 %ya为初值。 h=(b-a)/N; T=zeros(1, N+1); Y=zeros(1, N+1); T=a: h: b; Y(1)=ya; for j=1:N k1=h*feval(f, T(j), Y(j)); k2=h*feval(f, T(j)+h/2,Y(j)+k1/2); k3=h*feval(f, T(j)+h/2,Y(j)+k2/2); k4=h*feval(f, T(j)+h,Y(j)+k3); Y(j+1)=Y(j)+(k1+2*k2+2*k3+k4)/6; end R=[T Y] 7．4 线性多步法 线性多步法的基本思想就是充分利用前面多步的信息来预测，使算法尽可能获得更高的精度。 欧拉公式、后退欧拉公式、梯形公式和四阶龙格-库塔法都是线性一步法。 7．4． 1 用待定系数法构造线性多步法 7．4． 1 用待定系数法构造线性多步法 7．4． 1 用待定系数法构造线性多步法 7．4．2 用数值积分法构造线性多步法公式 The End c1+c2 =1 a2c2 =1/2 阶 随意考察[ xn , xn+1 ]内一点: xn+a2 = xn+a2h K1 = f (xn , yn)= yn’ 0a2≤1 K2 = f (xn+a2h , yn+a2hK1 ) f (xn+a2h , y(xn ) +a2hy’(xn ) ) 二元函数泰勒展开 = f (xn , y(xn)) + a2h f ’(xn , y(xn)) + O(h2) = y’(xn) + a2h y’’(xn) + O(h2) yn+1 = yn + h(c1K1+c2K2 ) = yn + (c1+c2 )h yn’ + a2c2h2 yn’’ + O(h3) 设 yn = y( xn ) 而 y(xn+1 ) 的泰勒展开式为: y(xn+1 ) = y(xn ) + hy’(xn ) + y’’(xn ) + O(h3) h2 2 yn+1 = yn + h(c1K1+c2K2 ) = yn + (c1+c2 )h yn’ + a2c2h2 yn’’ + O(h3) 设 yn = y( xn ) 比较系数得: c1+c2 =1 a2c2 =1/2 ( 7-17 ) ( 7-18 ) ( 要求 yn+k ) ( 7-21 ) ( 7-19 ) yn+k = 设αk = 1 * 数值分析——第 7章 《数值分析》 冶金工业出版社 数值分析——第 7章 知识要点 : 掌握解常微分方程初值问题的Euler方法的三种构造手段(Taylor级数法、数值积分法和数值微分法)； 会用Euler方法及改进的Euler方法和经典的Runge-Kutta方法解常微分方程初值问题。 ( 7-1 ) 实际上, 由常微分方程解的存在惟一性定理知, 只要函数 f ( x, y )满足利普希兹条件: | f ( x, y1 ) – f ( x, y2 )|≤ L | y1 - y2| 且 f ( x, y )满足一定的光滑度, 理论上就可以保证初值问题: dx = f ( x, y ) y( x0 ) = y0 ( 7-1 ) 的解 y = y( x ) 存在并且惟一。 dy 顶点Pn和Pn+1的坐标有下列关系： yn+1 - yn xn+1 - x