计算方法4_常微分方程数值解法.ppt

对象 思路 计算过程： 4.1 Euler 公式 ?思想: 用向前差商近似代替微商. 几何意义 欧拉法（续） ?思想: 用向后差商近似代替微商. 欧拉法（续） 欧拉法（续） 截断误差 Def4.1 设y(xn) 是(4.1)式的精确解，yn是(4.2)式欧拉法得到的近似解，称y(xn)-yn为欧拉法的整体截断误差. 二步欧拉法的局部截断误差 Def4.5 假设yn=y(xn) , yn－1=y(xn－1),称Rn+1=y(xn+1)-yn+1为二步欧拉法的局部截断误差. 数值积分法 ? 梯形公式:将显示欧拉公式,隐式欧拉公式平均可得 从而y(xn+1)=yn+1+h fy(xn+1,η)[y(xn+1)-yn+1] /2 +O(h3) ∴y(xn+1)-yn+1 = h fy(xn+1,η)[y(xn+1)-yn+1] /2 +O(h3) ∴y(xn+1)-yn+1=O(h3)/[1-hfy(xn+1,η)/2]=O(h3) ∴梯形公式的截断误差为O(h3)，其精度是2阶。 欧拉公式的比较 4.2 改进的Euler法 4.3龙格—库塔方法 为了构造函数φ使得(4.5)式成为高阶方法，Taylor 令 Runge-Kutta法的思想：在［xn,xn+1］内多预报几个点的ki值并用其加权平均值作为k*近似值。而构造出具有更高精度的计算公式。 二阶龙格—库塔方法 上述方程组有四