第7章 常微分方程（组）的数值解法-天津大学理学院数学系.ppt

7: ODE * 求此问题数值解的标准四阶Runge-Kutta公式为 7: ODE * 7.5.3 高阶常微分方程组初值问题的数值解法 对于高阶常微分方程组的初值问题, 考察下面的二阶常微分方程组初值问题: 令 将此问题转化为一阶常微分方程组初值问题 7: ODE * 等价的一阶常微分方程组初值问题 7: ODE * 本章小结 基本知识 建立初值问题数值解的计算格式的方法 数值微分方法、Taylor 展开法和数值积分方法 截断误差与方法的精度 Runge-Kutta方法 基本思想 一阶、二阶、四阶Runge-Kutta方法 一阶常微分方程组数值解法 高阶方程（组）的数值解法 * * * 当前的整体截断误差=上一步的整体截断误差+上一步整体截断误差引起的导函数f(x,y)的误差+当前的局部截断误差。 * 7: ODE * 7.1.2 截断误差与方法的精度 定义 7.1.1 称误差 en+1 = y ( xn+1 ) - yn+1 为数值方法在点 xn+1 的截断误差, 又称整体截断误差. 设 yk= y ( xk ) (k = 0,1,. . . , n)，则 为数值方法在点 xn+1 的局部截断误差. 7: ODE * 整体截断误差 en+1 是在没有引进舍入误差的情况下, 纯粹因为不准确的计算格式造成的, 故又称为方法误差.它不仅与 x = xn+1 这一步的计算有关, 而且和 xn , xn-1 ,. . . , x1 这几步的计算都有关系. 局部截断误差是假设 xn 之前各数值解没有误差, 仅由 xn 到 xn+1 这一步计算由计算格式引起的误差. 7: ODE * 如Euler公式 在点 xn+1 的整体截断误差 en+1 = y (xn+1)- yn+1 局部截断误差 当前的整体截断误差=上一步的整体截断误差+上一步整体截断误差引起的导函数 f (x, y) 的误差+当前的局部截断误差。 7: ODE * 定义7.1.2 若某数值方法的局部截断误差为 则称该方法具有 P 阶精度, 或称其为 P 阶方法. 可以证明: Euler 方法的局部截断误差 其具有一阶精度. 梯形方法的局部截断误差 其具有二阶精度. 改进的 Euler 方法的局部截断误差 具有二阶精度. 7: ODE * 7.2 Runge-Kutta 方法 继续讨论前面的 Taylor 展开法。 设 y∈C2[a , b], 由 Taylor 公式有 由 故上式即为 略去余项, 并以 yn , yn+1 代替 y (xn) 及 y (xn+1), 得到Euler 公式. 7: ODE * 进一步假设 y∈Cp+1[a, b],由 Taylor 公式有 其中 由 故式(7.2.1)即为 7: ODE * 7.2 Runge-Kutta 方法 略去余项, 并以数值解 yn , yn+1 代替 (7.2.3) 中的 解析解y (xn) 及 y (xn+1), 可得到一个差分方程，即 其中余项可写成 注：这里 7: ODE * 在(7.2.3)中略去余项,用yn, yn+1代替y(xn)及y(xn+1) 其中 称 (7.2.4) 式为求解常微分方程初值问题数值解Taylor的格式 . 7: ODE * 由于局部截断误差 可知它是一个 p 阶方法。当p=1时，上式正是Euler 公式。但当 p 2 时，需要计算f (x, y(x) ) 的高阶导数，特别是对于复杂函数 f (x, y(x)) 的求导，这无疑是大大增加计算量，这是它最大的缺点。因此高阶的Taylor方法是不实用的。 德国数学家C.Runge及M.W.Kutta提出了一种改进策略，得到了至今还作为高精度单步法被广泛使用龙格-库塔法(Runge-Kutta method)。 7: ODE * 7.2.1 Runge-Kutta方法的基本思想 Runge-Kutta方法是利用 f 在某些点处函数值的线性组合替代（7.2.4）步长 h 后面括号中的因子来构造差分方程, 从而避免了高阶导数的计算, 这就是 Runge-Kutta 方法的基本思想. 用f 在某些点处函数值的线性组合替代这一部分 7: ODE * 其一般形式为: 其中 r 是上式中调用 f 的个数, r 称为级数，