章常微分方程数值方法.docVIP

第六章 常微分方程数值方法 连续问题的离散处理——寻求微分方程的解在某些离散点上的值——在处的近似值； 记号：——所求函数在处的（准确）函数值， ——计算得到的处的函数（近似）值， ——节点， ——步长； 等步长：；以下若不另作说明，一般总记为等步长。 § 6.1 初值问题的数值方法 考察微分方程初值问题： 由微分方程理论可知：若函数关于满足条件，即存在与无关的常数，使 初值问题的解存在且唯一； 6.1.1 法及其变形 1、法 由Taylor 展开式： 作局部化假设： ，并略去项，便有 将此式右端作为的近似，便得到公式 Euler公式 两者的差（即略去的项）称为局部截断误差，记作： Euler公式也可由其他方法导出，例如由第四章数值导数公式，可有： ； 解出，并由 替换，便可得Euler公式。 又如，根据Newton-leibniz公式： 将以“0次插值多项式”替换，即以 代入积分，得到数值积分（左矩形）公式，及其误差： 又得到与Taylor 展开式相同的表达式，从而又导出Euler公式。 几何意义——折线法 若一个公式的局部截断误差为，则称该公式的精度为阶，或该公式为阶公式。 Euler公式是1阶公式。 注意，以上的截断误差是在局部化假设的前提下得到的，即认定。倘若在每一步都按局部化假设，我们