数值分析第八章常微分方程数值解.ppt

《数值分析第八章常微分方程数值解.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第八章常微分方程数值解数值分析华长生制作18.1引言(基本求解公式)在工程和科学技术的实际问题中，常需要求解微分方程只有简单的和典型的微分方程可以求出解析解而在实际问题中的微分方程往往无法求出解析解在高等数学中我们见过以下常微分方程：-----------(1)-----------(2)华长生制作2-----------(3)(1),(2)式称为初值问题,(3)式称为边值问题-----------(4)另外,在实际应用中还经常需要求解常微分方程组:本课程主要研究问题(1)的数值解法,对(2)~(4)只作

2、简单介绍我们首先介绍初值问题(1)的解存在的条件华长生制作3定理1.对于问题(1),要求它的数值解华长生制作4-----------(1)从(1)的表达式可以看出,求它的数值解的关键在于求解微分方程的数值方法数值积分数值微分华长生制作5而数值积分问题我们已经学习过,下考虑数值微分方法微积分中，关于导数的定义如下：自然而又简单的方法就是，取极限的近似值，即差商!华长生制作6向前差商由Taylor展开因此，有误差华长生制作7向后差商误差:中心差商华长生制作8常微分方程数值解的基本思想华长生制作9华长生制作10

3、华长生制作11华长生制作12华长生制作138.2欧拉(Euler)方法华长生制作14以上三种方法推导出同一个数值求解公式:这个数值公式称为欧拉(尤拉)(Euler)公式P299:例9.2.1华长生制作15华长生制作164.改进的欧拉(Euler)公式（显式方法）华长生制作17局部截断误差/方法的阶评价一个微分方程求解公式的标准当然是其精度而在求解公式中误差项称为在处的整体截断误差华长生制作18定义1(a).定义1(b).---------(*)上两定义本质是一样的,前者意义直观,后者用于计算推导较方便!华

4、长生制作19定义3.华长生制作20Euler公式的局部截断误差为具有1阶精度后退Euler公式的局部截断误差为也具有1阶精度具有2阶精度华长生制作21显然一个求解公式的精度越高,计算解的精确性也就越好从前面的分析可知,Euler法的精度并不算高因此有必要找寻精度更高的求解公式华长生制作22Don’tgoaway!华长生制作23