梯形法、龙格库塔法求解微分方程 数值解 笔记

首此发表于misskissC博客：http://blog.csdn.net/misskissc/article/details/8913941

•       根据积分表达式,微分方程的数值解关键在于微分方程的初值及计算微分方程式在tm(上一时刻)与tm+1（下一时刻）与坐标轴围成面积，若这个面积计算得越准确则得到的数值解也就越精确。微分表达式中与坐标轴围成的面积可表示如下，再实施算法的时候可以结合这个图更加直观点：

 

图1 微分方程中微分表达式中tm 与tm+1与坐标中组成阴影面积图

 

1.梯形法求数值解

       显然阴影部分的面积的近似计算中，用梯形法近似比用矩形近似更为准确，梯形的斜边用tm时刻的f的值代替，即原函数在tm时刻的导数。即用梯形法的数值解的表达式变为y(tm+1) = y(tm) + ( f（ y(tm), tm) + f( y1(tm+1), tm+1) ) / 2 * (tm+1 - tm)

从这个表达式来看，如果要用梯形法求得y(tm+1)的值则需要先知道y1(tm+1)的值，这个值可以用欧拉法求得