[理学]数值分析第9章 常微分方程初值问题数值解法.pptVIP

和四阶隐式亚当姆斯格式 （9.24） 9.5.2 亚当姆斯预报-校正格式 参照改进的欧拉格式的构造方法，以四阶亚当姆斯为例，将显式（9.22）和隐式（9.24）相结合，用显式公式做预报，再用隐式公式做校正，可构成亚当姆斯预报-校正格式 （9.25） 预报： 校正： 这种预报-校正格式是四步法，它在计算yi+1时不但用到前一步的信息 ，而且要用到再前面三步的信息 ，因此它不能自行启动。在实际计算时，可借助于某种单步法，譬如四阶龙格—库塔法提供开始值 。 例9.5 取步长h=0.1,用亚当姆斯预报-校正公式求解 初值问题 的数值解。 解: 用四阶龙格-库塔公式求出发值 ,计算得： 表中的 ,yi和y(xi)分别为预报值、校正值和准确解( ),以比较计算结果的精度。 再使用亚当姆斯预报-校正公式(9.25), §9.6 一阶常微分方程组与高阶方程 我们已介绍了一阶常微分方程初值问题的各种数值解法，对于一阶常微分方程组，可类似得到各种解法，而高阶常微分方程可转化为一阶常微分方程组来求解。? 9.6.1 一阶常微分方程组 对于一阶常微分方程组的初值问题 （9.31） 可以把单个方程 中的f 和y看作向量来处理，这样就可把前面介绍的各种差分算法推广到求一阶方程组初值问题中来。 设 为节点上的近似解， 则有改进的Euler格式为 预报： 校正： （9.32） 又，相应的四阶龙格—库塔格式（经典格式）为 （9.33） 式中 （9.34） 把节点xi上的yi和zi值代入式(7.34), 依次算出 , 然后把它们代入式(7.33), 算出节点xi+1上的yi+1 和zi+1值。 对于具有三个或三个以上方程的方程组的初值问题,也可用类似方法处理,只是更复杂一些而已。此外,多步方法也同样可以应用于求解方程组初值问题。 例7.6 用改进的Euler法求解初值问题 取步长h=0.1，保留六位小数。 解: 改进的Euler法公式为 预报： 校正： 将 及h=0.1代入上式,得 由初值 ,计算得 9.6.2 高阶方程组 高阶微分方程(或方程组)的初值问题,原则上都 可以归结为一阶方程组来求解。例如,有二阶微分方 程的初值问题 （9.35） 在引入新的变量 后,即化为一阶方程组初值问题: （9.36） 式（9.36）为一个一阶方程组的初值问题，对此可用前面中介绍的方法来求解。例如应用四阶龙格-库塔公式得 （9.37） （9.38） 消去 ，上式简化为： （9.39） （9.40） 上述方法同样可以用来处理三阶或更高阶的微分方程（或方程组）的初值问题 例9.7 求解下列二阶微分方程的初值问题 取步长h=0.1 解:先作变换：令 ，代入上式，得一阶方程组 用四阶龙格-库塔方法求解,按式(9.37)及(9.38)进行计算： 取步长 , , , 时 然后计算 时的 y2和z2；依此类推，直到i=9时的y10和z10，即可得到其数值解。 本章小结 本章介绍了常微分方程初值问题的基本数值解法。包括单步法和多步法。单步法主要有欧拉法、改进欧拉法和龙格—库塔方法。多步法是亚当姆斯法。它们都是基于把一个连续的定解问题离散化为一个差分方程来求解，是一种步进式的方法。用多步法求常微分方程的数值解可获得较高的精度。 实际应用时，选择合适的算法有一定的难度，既要考虑算法的简易性和计算量，又要考虑截断误差和收敛性、稳定性。 龙格-库塔法较为常用，适用于多步方法中作初值计算和函数f(x,y)较为简单的场合。四阶标准龙格—库塔法精度高，程序简单，易于改变步长，比较稳定，也是一个常用的方法，但计算量较大。当函数f(x,y)较为复杂，可用显式亚当姆斯方法或亚当姆斯预测—校正方法，不仅计算量较小，稳定性也比较好，但不易改变步长。 一般采用龙格—库塔法提供初值y1, y2, y3，然后用亚当姆斯外推公式求得预测值 ，再由亚当姆斯内插值求得校正值yi+1，如此求得的值近似程度好且节省计算量，是一种较好的方法。 作业 习