第四篇 常微分方程数值解法.pptVIP

四　常微分方程数值解法 常微分方程数值解法 引言（常微分方程数值解法概述） 显式欧拉法、隐式欧拉法、二步欧拉法 局部截断误差与精度 改进的欧拉方法 龙格-库塔方法 收敛性与稳定性简述 一阶常微分方程组与高阶常微分方程 引言 一阶常微分方程初值问题： 引言（续） 实际生产与科研中，除少数简单情况能获得初值问题的初等解（用初等函数表示的解）外，绝大多数情况下是求不出初等解的。 有些初值问题即便有初等解，也往往由于形式过于复杂而不便处理。 实用的方法是在计算机上进行数值求解：即不直接求 y(x) 的显式解，而是在解所存在的区间上，求得一系列点 xn (n = 0, 1, 2, …) 上解的近似值。 欧拉（Euler）方法 方法一　化导数为差商的方法 代入初值问题表达式可得： 欧拉方法（续） 方法二　数值积分法 欧拉方法的几何意义 隐式欧拉法 在数值积分法推导中，积分的近似值取为积分区间宽度与右端点处的函数值乘积，即： 二步欧拉法 在数值积分法推导中，积分区间宽度选为两步步长，即积分区间为：[xn-1, xn+1]，则： 梯形公式欧拉法 在数值积分法中，如果用梯形公式近似计算 f (x, y) 在区间 [xn, xn+1] 上的积分，即： 例 用显式欧拉法、隐式欧拉法、梯形法求解初值问题： 例：（续） 计算结果： 局部截断误差 为了简化分析某常微分方程数值算法的误差，现假设 yn =