常微分方程数值解.pptVIP

第九章常微分方程数值解计算方法

第九章常微分方程数值解9.1引言9.3Runge-Kutta法9.2简单数值方法

本章要点：本章作业本章主要研究基于微积分数值解法的常微分方程数值解，主要方法有Euler方法、改进Euler方法Runge-Kutta方法P316.1-4选二；5.(1)(2)选一；6,7选一.

9.1引言在工程和科学技术的实际问题中，常需要求解微分方程只有简单的和典型的微分方程可以求出解析解而在实际问题中的微分方程往往无法求出解析解在高等数学中见过以下常微分方程：-----------(1)-----------(2)

-----------(3)(1),(2)式称为初值问题,(3)式称为边值问题-----------(4)另外,在实际应用中还经常需要求解常微分方程组:这里主要研究问题(1)的数值解法。我们首先介绍初值问题(1)的解存在的条件

定理1.对于问题(1),要求它的数值解

-----------(1)从(1)的表达式可以看出,求它的数值解的关键在于

-----------(1)对于初值问题(1)在下列子区间上分别应用两点数值微分公式9.2简单的数值方法9.2.1欧拉法与后退欧拉法假设节点为等距节点，步长

--------(6)记(6)式称为求解初值问题(1)的(前进)Euler公式

（1）右端的积分用左矩形公式（2）右端的积分用右矩形公式---显式的---隐式的

或由也可得由于后退Euler公式是隐形公式,计算例1将很麻烦事实上大多数情况下用后退Euler公式都较困难

用迭代法逐步显化，得到新的Euler公式--------(10)由Lipstz条件因此，只要就有

--------(7)此方法称为预测—校正系统

这种类型的方法称为单步格式或单步法Euler方法的几何体现:前进Euler公式后退Euler公式

例1.解:由前进Euler公式Euler1.m

得依此类推,有01.00000.10001.10000.20001.19180.30001.27740.40001.35820.50001.43510.60001.50900.70001.58030.80001.64980.90001.71781.00001.7848

用Euler公式的预测——校正系统求解例1.例2.解:由(7)式,有Euler1.m

依此类推,得01.00000.10001.09180.20001.17630.30001.25460.40001.32780.50001.39640.60001.46090.70001.52160.80001.57860.90001.63211.00001.6819比较不同的结果

9.2.2梯形法改进的Euler法--------(10)得梯形方法平均化形式

9.2.4单步法的局部截断误差与阶单步法一般形式显式单步法（*）实际上是假设，则有为单步法（*）的局部截断误差。定义1称(**)

Euler公式的局部截断误差为主项称单步法（*）具有阶精度。定义2若存在最大整数，使单步法（*）的截断误差满足对Euler公式

对后退Euler公式主项为

对梯形法主项为

定义2.

例3.用Euler公式、梯形公式和改进Euler公式求解初值问题,并比较结果的精度解:(1)Euler公式

(2)梯形公式

(3)改进Euler公式xy01.00000.10.90500.20.81900.30.74120.40.67080.50.6071使用MATLAB软件Euler2.m结果为

0.9050