matlab微分方程实例,Matlab常微分方程的解法

【实例简介】

和Matlab应用有关的，具体介绍常微分方程的使用和解法，原理性介绍，帮助理解。

局部截断误差指的是,按()式计算由到这一步的计算值与精确值

之差

+。为了估计它,由

展开得到的精确值

是

()、()两式相减(注意到=

)得

即局部截断误差是阶的,而数值算法的精度定义为:

若一种算法的局部截断误差为

,则称该算法具有阶精度。

显然葳大,方法的精度葳高。式()说明,向前

方法是一阶方法,因此

它的精度不高。

改进的方法

梯形公式

利用数值积分方法将微分方程离散化时,若用梯形公式计算式中之右端积分,

即

并用代替

,则得计算公式

这就是求解初值问题()的梯形公式。

直观上容易看出,用梯形公式计算数值积分要比矩形公式好。梯形公式为二阶方法

梯形公式也是隐式格式,一般需用迭代法求解,迭代公式为

由于函数

关于满足

条件,容易看出

其中为

常数。因此,当

如果实际计算时精度要求不太高,用公式()求解时,每步可以只迭代一次,由此导

出一种新的方法一改进

改进

法

按式()计算问题()的数值解时,如果每步只迭代一次,相当于将公式

与梯形公式结合使用:先用公式求的个初步近似值,称为预测值,然

后用梯形公式校正求得近似值,即

预测

校正

式()称为由公式和梯形公式得到的预测一校正系统,也叫改进法。

为便于编制程序上机,式()常改写成

改进