第五章常微方程初值问题2材料.docVIP

第五章 常微方程初值问题 数值微分 我们讨论的导数问题。实际计算及讨论函数在具体某点的导数值。 我们自然想到微积分中的导数定义： 用作为的近似。 理论上这种近似，h越小越好，但实际计算相近的数相减会损失有效数字，实际计算中要注意这些问题。 二阶导数的三点公式： 精度高一个量级 往往利用插值方法先求出插值函数，在对函数求导，用之代替原函数的导数。 数值积分我们也曾这样做。 本章我们主要解决初值问题的数值解法。即： 数值解就是要求出在 一系列点上的近似值 。 假设在所讨论问题中连续。由Leibniz公式（或对方程积分）有： 依据上式可得到不同的数值解法。 §1 欧拉法和改进欧拉法、预估-校正法 由 积分可用 左矩形公式，得： 显式 用 右矩形公式，得： 隐式 梯形公式，得： 改进欧拉公式 精度比上两个高一个量级。 隐式和改进欧拉法都是一个方程。如果方程解时比较复杂，则用迭代发近似求解。 先用显式欧拉公式计算初值，再用改进欧拉公式迭代求解。 （数学证明当在所讨论的范围有界在h取得很小时，这种迭代收敛）。 公式为: 程序结构： F(x,y)=Y*Y WRITE(*,*) INPUT X0,Y0 READ(*,*) X0,Y0 WRITE(*,*) INPUT [x0，b