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

4．R-K-Fehhlberg 方法 R-K-Fehhlberg方法是在R-K方法的基础上引进误差和步长 控制的办法。即利用5阶R-K法 估计4阶R-K的局部误差，其中 注：R-K-Fehhlberg比4阶R-K方法具有更大的优越性，他是计算稳定，高精度的方法，他的显著优点是，每一步仅需计算f的6个值；若用4阶R-K-L与5阶R-K-L在一起使用，则每步需要计算f的10个值。推荐使用！ 5. 隐式R-K方法 类似于显式R-K公式，稍加改变，就得到隐式R-K方法。 它与显式R-K公式的区别在于：显式公式中对系数求和的上限是i-1，从而构成的矩阵是一个严格下三角阵。而在隐式公式中对系数求和的上限是L，从而构成的矩阵是方阵，需要用迭代法求出Ki，。推导隐式公式的思路和方法与显式R-K法类似。通常，同级的隐式公式获得比显式公式更高的阶。 通常，同级的隐式公式获得比显式公式更高的阶。常用的隐式 R-K法有: 1级2阶中点公式 : 2级2阶梯形公式： 2级4阶R-K公式： 6.变步长方法 在单步法中每一积分步步长实际上是相互独立的，步长的 选择具有了灵活性。根据合理地选择每一积分步的步长， 既保证精度的要求，又可以减少计算量，从而减小舍入误 差。其方便的控制手段是基于误差的事后估计式。 对于给定的精度 ,如果 ,反复将步长折半进行计算，直至 为止,这时取最终得到的作为结果； 如果 ，将反复将步长加倍，直到 为止，这时再将步长折半一次，就得到所要的结果。 这种通过加倍或折半处理步长的计算方法称为变步长方法。 注：推荐使用精度好计算量低的变步长方法。 用四阶显式R-K方法做变步长方法是实践中较好的方法！ 例 分别用改进的欧拉格式和四阶龙格—库塔格式解初值问题（取步长h=0.2）： 节点 改进欧拉法 四阶龙格—库塔法 准确解 0 1 1 1 0.2 1.186667 1.183229 1.183216 0.4 1.348312 1.341667 1.341641 0.6 1.493704 1.483281 1.483240 0.8 1.627861 1.612514 1.612452 1 1.754205 1.732142 1.732051 表 （注：已指出过准确解 ） 利用 Euler方法求初值问题 解 此时的 Euler公式为 例 的数值解.此问题的精确解是 分别取步长h=0.2 ,0.1 ,0.05,计算结果如下 h xn yn y(xn) y(xn)-yn h=0.2 0.00 0.40 0.80 1.20 1.60 2.00 0.00000 0.37631 0.54228 0.52709 0.46632 0.40682 0.00000 0.34483 0.48780 0.49180 0.44944 0.40000 0.00000 -0.03148 -0.05448 -0.03529 -0.01689 -0.00682 h=0.1 0.00 0.40 0.80 1.20 1.60 2.00 0.00000 0.36085 0.51371 0.50961 0.45872 0.40419 0.00000 0.34483 0.48780 0.49180 0.44944 0.40000 0.00000 -0.01603 -0.02590 -0.01781 -0.00928 -0.00419 h=0.05 0.00 0.40 0.80 1.20 1.60 2.00 0.00000 0.35287 0.50049 0.50073 0.45425 0.40227 0.00000 0.34483 0.48780 0.49180 0.44944 0.40000 0.00000 -0.00804 -0.01268 -0.00892 -0.00481 -0.00227 * * 《计算方法》 7 微分方程数值解 飞机设计教研室 白巍 提纲挈领 §1、引言 §2、初值问题的数值解法--单步法 §3、龙格-库塔方法 §4、收敛性与稳定性 §5、初值问题的数值解法―多步法 微分方程的一般概念 定义 联系自变量和未知函数及其导数的等式。 分类 按自变量的个数：常微分方程，偏微