第八章_常微分方程的数值解法.ppt

四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ e(h)=O(h5) 4阶精度 ?优缺点： 优点是：(1)都是一步法，因此只要给定一个初始值就 可以一直计算下去； (2)精度相对较高，如经典R-K法为四阶精度 缺点是：(特别是三阶和四阶法)计算量较大。 ?其他问题讨论： R-K法推导基于Taylor展式，因而要求y(x)有较好的光滑性（即有高阶导数）。 最常用的是四阶公式，它适用于一般的问题，准确、稳定、易于编程。 步长h减小，局部误差O(h5)减小 但步数增加，舍入误差积累增加 h要适当， 总误差才最小 四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ 例：取h=0.2，用四阶龙格－库塔法求解初值问题 举例 解： 这里 ，经典的四阶龙格-库塔公式为 举例 表中列出了计算结果，同时列出了相应的精确解.比较本章第一个例子的计算结果，显然龙格-库塔方法的精度高. (Matlab) 习题八 3，7 第八章 作 业 龙格－库塔法 （Runge－Kutta Method ） * 提问学生：数值解法的步进式特点，让学生回答这个特点是怎么回事，然后让学生按照递归的特点，得出计算表达式的特点。 * 第八章 常微分方程数值解法(Numerical Methods for Ordinary Differential Equations) 引言 引言 ? 考虑一阶常微分方程的初值问题 (Initial-Value Problem ): 只要 f (x, y) 在x∈ [a, b] 上连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述常微分方程存在唯一解。 引言 数值解法就是要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 y0 y1 … yn 节点间距 称为步长，通常采用等距节点，即取 hi = h (常数)。 解决： 数值解法的一个基本特点是“步进式”，即求解时顺着节点排列的次序一步步地向前推进。 单步： yk-1 yk 多步： yk-p … yk-2 ，yk-1 yk 注意：与“迭代法”区别 第八章 常微分方程数值解法 欧拉方法（ Euler’s Method ） 显式Euler公式 由两点公式求导数，在[xj , xj+1]子区间上有： 其中?j ? [xj , xj+1] 代入方程 有 显式Euler公式 x0 x1 y(x1) y(x0) h 显式Euler公式的误差 局部截断误差 显式Euler公式 定义 　　　在假设 yj = y(xj)，即第 i 步计算是精确的前提下，考虑的截断误差 ej+1 = y(xj+1) ? yj+1 称为局部截断误差 /* local truncation error */。 （P230）例1：取h=0.1，分别用显式Euler法、显隐结合的预测校正系统求解初值问题 解： 用显式Euler法求解，有： …依次下去计算结果见P231 隐式Euler公式 由两点公式求导数，在[xj , xj+1]子区间上有： 其中?j ? [xj , xj+1] 代入方程 有 局部截断误差 隐式Euler公式 是一个关于yj+1的方程, 要从中解出yj+1 x0 x1 显隐结合的预测校正系统——避免求解方程（predictor-corrector method ） Step 1: 先用显式欧拉公式作预测，算出预测值 Step 2: 再用隐式欧拉公式作校正，得到校正值 写成一个公式为： y0 ＝y(a)－ y1－ y1－ y2－ y2 … － yn－ yn 计算顺序： （P230）例1：取h=0.1，分别用显式Euler法、显隐结合的预测校正系统求解初值问题 ……计算过程见P231 解： 用预测校正系统求解，有： y0=1 p阶精度 　　若某算法的局部截断误差e(h)满足： e(h)= O(hp+1)，即有： e(h)/ hp+1=c(常数）， 则称该算法有p 阶精度。 定义 ? 欧拉法的局部截断误差： 欧拉法具有 1 阶精度。 显式： 隐式： 梯形公式 （trapezoid form