常用算法--几数字积分法.docVIP

几种常用的数字积分方法（微分方程的数字解） 2－5数字积分法 欧拉法（折线法） t0 t0 t1 t y(t) h y y(t0) y(t1) 图2-5-1 t0 t1 t f(t) f f(t0) 图2-5-2 h 设一阶微分方程 由图可知，过（t0, y0）点的斜率为 如果 离 很近，即 很小，曲线y(t)可用切线来近似，其切线方程 其微分方程在t=t1 时，可近似表示为 重复上述近似过程，当时， 则有一般近似公式 如果令，称为计算步矩，则 (1) 这就是欧拉法数字积分的递推计算公式。 由公式可看出，只要我们给出方程的初值（t0, y0）以及相应的步距，逐步进行递推就可获得微分方程的近似数字解。 欧拉法的计算是十分简单的，其计算误差正比于，由此，要获得高精度解，必须减小步距，但这使得计算次数增加，又由于计算机的字长有限，h减小得过小，将引入舍入误差，所以此方法的精度提高有限，实际应用中较少采用。 梯形法（预报――校正法） 欧拉法精度低，却给我们一些启发，对微分方程 可改写成 当 时，则 从此式可以看出，要求得 的值，等式右边中含有未知函数，所以不能得到的值，但如果我们用已知的函数值来代替，用不变取代变化的函数，即 实际上右边是一个矩形面积 则 递推公式为 用此矩形的面积