【实例简介】
和Matlab应用有关的,具体介绍常微分方程的使用和解法,原理性介绍,帮助理解。
局部截断误差指的是,按()式计算由到这一步的计算值与精确值
之差
+。为了估计它,由
展开得到的精确值
是
()、()两式相减(注意到=
)得
即局部截断误差是阶的,而数值算法的精度定义为:
若一种算法的局部截断误差为
,则称该算法具有阶精度。
显然葳大,方法的精度葳高。式()说明,向前
方法是一阶方法,因此
它的精度不高。
改进的方法
梯形公式
利用数值积分方法将微分方程离散化时,若用梯形公式计算式中之右端积分,
即
并用代替
,则得计算公式
这就是求解初值问题()的梯形公式。
直观上容易看出,用梯形公式计算数值积分要比矩形公式好。梯形公式为二阶方法
梯形公式也是隐式格式,一般需用迭代法求解,迭代公式为
由于函数
关于满足
条件,容易看出
其中为
常数。因此,当
如果实际计算时精度要求不太高,用公式()求解时,每步可以只迭代一次,由此导
出一种新的方法一改进
改进
法
按式()计算问题()的数值解时,如果每步只迭代一次,相当于将公式
与梯形公式结合使用:先用公式求的个初步近似值,称为预测值,然
后用梯形公式校正求得近似值,即
预测
校正
式()称为由公式和梯形公式得到的预测一校正系统,也叫改进法。
为便于编制程序上机,式()常改写成
改进