常微分方程数值ppt课件.ppt

《常微分方程数值ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第七章常微分方程数值解数值分析1第七章常微分方程数值解7.1引言(基本求解公式)7.2Runge-Kutta法7.3微分方程组和高阶方程解法简介2本章要点：本章主要研究基于微积分数值解法的常微分方程数值解，主要方法有线性单步法中的Euler方法、Euler方法的改进方法、Runge-Kutta方法高阶微分方程和微分方程组的数值解法37.1引言(基本求解公式)在工程和科学技术的实际问题中，常需要求解微分方程只有简单的和典型的微分方程可以求出解析解而在实际问题中的微分方程往往无法求出解析解在高等数学中我们见过以下常微分方程：----

2、-------(1)-----------(2)4-----------(3)(1),(2)式称为初值问题,(3)式称为边值问题本课程主要研究问题(1)的数值解法我们首先介绍初值问题(1)的解存在的条件5定理1.对于问题(1),要求它的数值解6-----------(1)从(1)的表达式可以看出,求它的数值解的关键在于而数值微分或数值积分问题我们都已经学习过7一、基于数值微分的常微分方程数值解法-----------(1)对于初值问题(1)为了讨论方便,假设以下节点为等距节点8(一)Euler公式将在按泰勒公式展开，有9----

3、----(2)--------(3)记其中(2)和(3)式称为求解初值问题(1)的(前进)Euler公式和误差项在按泰勒公式展开10同理，将在按泰勒公式展开，有11记其中--------(4)--------(5)(4)和(5)式称为求解初值问题(1)的后退Euler公式和误差项12从(2)或(4)式不难看出,这种类型的方法称为单步格式或单步法Euler方法的几何体现:前进Euler公式后退Euler公式13Euler1.m例1.解:由前进Euler公式14得依此类推,有01.00000.10001.10000.20001.19

4、180.30001.27740.40001.35820.50001.43510.60001.50900.70001.58030.80001.64980.90001.71781.00001.784815由于后退Euler公式是隐形公式,计算例1将很麻烦事实上大多数情况下用后退Euler公式都较困难就可得到新的Euler公式--------(6)此方法称为预测—校正系统16用Euler公式的预测——校正系统求解例1.例2.解:由(6)式,有Euler1.m17依此类推,得01.00000.10001.09180.20001.1763

5、0.30001.25460.40001.32780.50001.39640.60001.46090.70001.52160.80001.57860.90001.63211.00001.6819比较不同的结果18(二)常微分方程数值解的截断误差评价一个微分方程求解公式的标准当然是其精度而在求解公式中误差项19定义1.因为一般情况下,求解公式的每一步都存在误差,因此有定义2.定义3.2021Euler公式的局部截断误差为具有1阶精度后退Euler公式的局部截断误差为也具有1阶精度显然一个求解公式的精度越高,计算解的精确性也就越好从前

6、面的分析可知,Euler法的精度并不算高因此有必要找寻精度更高的求解公式22二、基于数值积分的常微分方程数值解法-----------(1)对于初值问题-----------(7)23矩形求积公式梯形求积公式,误差为将以上求积公式代入(7)式,并加以处理就可得到相对应的求解公式Simpson求积公式,误差为24(一)矩形求解公式由可得令-----------(8)(8)式称为矩形公式(矩形法)实际上就是Euler求解公式25(二)梯形求解公式由可得令------(9)称(9)式为梯形求解公式(梯形法)注意:(9)式是隐形公式26

7、则梯形公式第k步的截断误差为显然梯形法具有二阶精度由于梯形公式为隐形公式,一般情况下不易显化27------(10)以上公式称为改进的Euler求解公式(改进Euler法),即------(11)28例3.用Euler公式、梯形公式和改进Euler公式求解初值问题,并比较结果的精度解:(1)Euler公式29(2)梯形公式3031(3)改进Euler公式xy01.00000.10.90500.20.81900.30.74120.40.67080.50.6071使用MATLAB软件Euler2.m结果为320.90500.8190

8、0.74120.67080.6071Euler公式梯形公式改进Euler公式结果比较Euler法的精度不如梯形公式3334