第七章 常微分方程数值解法汇.ppt

本章主要内容： 7.1 欧拉法和改进的欧拉法 7.2 龙格-库塔法 7.3 线性多步法 引言： 引言： 引言： 引言： 引言： 7.1.1 欧拉法及其截断误差 7.1.1 欧拉法及其截断误差 7.1.1 欧拉法及其截断误差 7.1.1 欧拉法及其截断误差 7.1.1 欧拉法及其截断误差 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 改进的欧拉法及预测-校正公式 7.1.2 向后 (隐式)欧拉公式 7.1.2 两步欧拉公式 引言：（回顾） 7.1 欧拉法和改进的欧拉法 7.1 欧拉法和改进的欧拉法 引言： 引言： 引言： 7.2.1 二阶龙格-库塔公式 7.2.1 二阶龙格-库塔公式 7.2.1 二阶龙格-库塔公式 7.2.1 二阶龙格-库塔公式 7.2.1 二阶龙格-库塔公式 7.2.2 四阶龙格-库塔公式 7.2.2 四阶龙格-库塔公式 7.2.2 四阶龙格-库塔公式 7.3 线性多步法 7.3.0 多步法的概念 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.1 四阶阿达姆斯外插公式 7.3.2 四阶阿达姆斯内插公式 7.3.2 四阶阿达姆斯内插公式 7.3.3 初始出发值的计算 7.3.4 阿达姆斯预测-校正公式 7.3.4 阿达姆斯预测-校正公式 7.3.4 阿达姆斯预测-校正公式 由T(h)的泰勒公式 为使T(h)=O(h3),令 解出 ，得 整理得 利用 可以推出 取x=xi并略去O(h3)便得到二阶龙格－库塔公式 或 可以推出 取x=xi并略去O(h3)便得到二阶龙格－库塔公式 或 仿照上述的讨论，可导出四阶龙格－库塔公式: 例 取步长h=0.2，用四阶龙格－库塔公式求下面初值问题的数值解。 解 由公式得 数值解yi与准确解y(xi)的对照见表 准确解是 0 0.2 0.4 0.6 0.8 1.0 1 1.18323 1.34167 1.48328 1.61251 1.73214 1 1.18322 1.34164 1.48324 1.61245 1.73205 xi yi y(xi) 7.3.1 四阶阿达姆斯(Adams)外插公式 7.3.2 四阶阿达姆斯(Adams)内插公式 7.3.0 多步法的概念 7.3.3 初始出发值的计算 7.3.4 阿达姆斯预测-校正公式 单步法 ----计算yi+1时只使用yi的值。 多步法 ----计算yi+1时使用前面的k个y值，即由 yi-k+1 , yi-k+2 , …,yi-1, yi计算yi+1。(k=1,2,……) 线性多步法 ----计算yi+1的公式由yi-k+1 , yi-k+2 , …,yi-1, yi 的线性组合表达。 设想用yi-3 ,yi-2 ,yi-1 ,yi 的值计算yi+1。为方便讨论由 出发计算y(x+h), 由初值问题的方程y′=f(x,y(x))两边从x到x+h积分, 可得到等价的积分方程 设想运用数值积分方法，取 x-3h, x-2h, x- h, x 为插值基点，做 f(s,y(s)) 的三次拉格朗日插值，用它近似计算上式的积分 。 这样得到的数值积分公式是f(s,y(s))在4个插值基点处的函数值的线性组合。 由于 f(x-ih,y(x-ih))= y(x-ih)，所得到的计算 y(x+h) 的近似公式形为： 为达到四阶精度，希望确定参数b0 ,b1 ,b2 ,b3使满足 运用在 h=0 处的泰勒公式得 为达到四阶精度，希望确定参数b0 ,b1 ,b2 ,b3使满足 运用在 h=0 处的泰勒公式得 为使误差等于O(h5)，令h, h2,h3,h4 的系数为0， 得方程组： 求得 代入前面的公式得 令x=xi 并记 略去余项，得到四阶阿达姆斯外插公式 ： 这是显式公式。公式的截断误差为O(h5)。 把四阶阿达姆斯外插公式中使用的yi-3 ,yi-2 ,yi-1 ,yi 改为 ： yi-2 ,yi-1 ,yi ,yi+1 , 经类似的推导可得近似公式： 确定待定系数的方程组为 解得 得到四阶阿达姆斯内插公式 这是一个隐式公式，截断误差也是O(h5