第二讲 常微分方程数值解法.ppt

初始条件则相应的化为 9.3 刚性方程组 在求解方程组 时，经常出现解的分量数量级差别很大的情形，这给数值求解带来很大困难，这种问题称为刚性问题。 若线性系统 ，其中， 中A的特征值 满足条件0(j=1,…,N), 且 则称 系 统为刚性方程，称s为刚性比 §7 习题与总结 1数值例题 2数值练习 3总结 1、数值例题 我们已经学习了很多数值算法，他们的效果到底如何呢？ 下面我们来分析一道例题，看看那些方法，就这个问题，最能接近真实值 求初值问题 的数值解，取h=0.1并同精确解 比较 （1）用欧拉法来计算这个初值问题 （2）用各阶的Runge—Kutta方法来计算这个初值问题 （3）用四阶的Adams定步长，Adams定步长预测校正方 法来 计算这个初值问题。然后将数值结果列成表格： 自变量 解析解 欧拉法 后退欧拉法 梯形法 改进的欧拉法 0.1 1.0048 1.0000 1.0091 1.0048 1.0050 0.2 1.0187 1.0100 1.0264 1.0186 1.0190 0.3 1.0408 1.0290 1.0513 1.0406 1.0412 0.4 1.0703 1.0561 1.0830 1.0701 1.0708 0.5 1.1065 1.0905 1.1209 1.1063 1.1071 0.6 1.1488 1.1314 1.1645 1.1485 1.1494 0.7 1.1966 1.1783 1.2132 1.1963 1.1972 0.8 1.2493 1.2305 1.2665 1.2490 1.2500 0.9 1.3066 1.2874 1.3241 1.3063 1.3072 1.0 1.3679 1.3488 1.3855 1.3676 1.3685 自变量 解析解 4阶显RK 4阶Adams显法 预测校正法 0.1 1.0048 1.0048 1.0048 1.0048 0.2 1.0187 1.0187 1.0187 1.0187 0.3 1.0408 1.0408 1.0408 1.0408 0.4 1.0703 1.0703 1.0703 1.0703 0.5 1.1065 1.1065 1.1065 1.1065 0.6 1.1488 1.1488 1.1488 1.1488 0.7 1.1966 1.1966 1.1966 1.1966 0.8 1.2493 1.2493 1.2493 1.2493 0.9 1.3066 1.3066 1.3066 1.3066 1.0 1.3679 1.3679 1.3679 1.3679 2、数值练习 计算实习题： 1.使用改进的Euler算法求初值问题 的数值解，取步长h=0.1 ,并同精确解 比较 2. 利用4阶Runge-Kutta算法求初值问题 的数值解，取步长：(1)h=0.1;(2)h=0.01.求数值解及精确解 3.使用4阶定步长Adams预测－校正算法求初值问题 的数值解，取h=0.1,并同精确解 比较 4.使用4阶定步长Adams预测－校正算法求初值问题 的数值解，使用h=0.5 5．综合题 请用已学的数值方法来计算下面初值问题，看看那种精度最好， 的数值解。 总 结 本章研究的是常微分方程初值问题的数值解法，我们从单步法、多步法两大类来介绍经典算法，并且对相容性、收敛性、稳定性这些理论进行系统的、全面介绍，为同学们学习偏微分方程数值解法作了较好的铺垫。 梯形法、R-K方法、变步长法以及Adams方法，都是一些高精度的方法，但是它们也有缺点，四阶R-K方法必须要求f(x,y)有较好的光滑性，若光滑性差，还不如欧拉法或改进欧拉公式。 它的另一个缺点是计