ch9-常微分方程初值问题数值解法 1.pptVIP

这种预报-校正格式是四步法，它在计算yi+1时不但用到前一步的信息 ，而且要用到再前面三步的信息 ，因此它不能自行启动。在实际计算时，可借助于某种单步法，譬如四阶龙格—库塔法提供开始值 。 例9.5 取步长h=0.1,用亚当姆斯预报-校正公式求解 初值问题 的数值解。 解: 用四阶龙格-库塔公式求出发值 ,计算得： 表中的 ,yi和y(xi)分别为预报值、校正值和准确解( ),以比较计算结果的精度。 再使用亚当姆斯预报-校正公式(9.25), §9.6 一阶常微分方程组与高阶方程 我们已介绍了一阶常微分方程初值问题的各种数值解法，对于一阶常微分方程组，可类似得到各种解法，而高阶常微分方程可转化为一阶常微分方程组来求解。? 9.6.1 一阶常微分方程组 对于一阶常微分方程组的初值问题 （9.31） 可以把单个方程 中的f 和y看作向量来处理，这样就可把前面介绍的各种差分算法推广到求一阶方程组初值问题中来。 设 为节点上的近似解， 则有改进的Euler格式为 预报： 校正： （9.32） 又，相应的四阶龙格—库塔格式（经典格式）为 （