常微分方程数值解法简介

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

1、科学计算与数学建模中南大学数学科学与计算技术学院第7章常微分方程数值解法简介简单的数值方法与基本概念2数值线性多步法3非线性高阶单步法——Runge-Kutte法45一阶方程组和高阶方程的初值问题6常微分方程边值问题的数值解法1实际问题的微分方程模型第七章　常微分方程数值解法简介微分方程在科学和工程技术中有很广泛的应用。许多实际问题的数学模型都可以用微分方程来描述，归结为常微分方程的定解问题。很多偏微分方程问题，也可以化为常微分方程问题来近似求解，但是求出所需的解绝非易事。实际上，除了极特殊情形外，人们不可能求出微分方程的解析解，只能用各种近似方法得到满足一定精度的近似解。在常微分方程中已经熟

2、悉了级数解法和Picard逐步逼近法，这些方法可以给出解的近似表达式，称为近似解析方法。另一类方法只给出解在一些离散点上的值，称为数值方法。数值方法应用范围更广，特别适合用计算机计算，本章主要介绍常用的常微分方程数值解法。函数是事物的内部联系在数量方面的反映,如何寻找变量之间的函数关系,在实际应用中具有重要意义.在许多实际问题中，往往不能直接找出变量之间的函数关系,但是有时却容易找出变量的改变量之间的关系，从而建立描述问题的微分方程模型.§1实际问题的微分方程模型例7.1将初始温度的一碗汤放置于环境温度保持在的桌上，10分钟后测得汤的温度为1000C。如果汤的温度低于550C才可以喝，问再过2

3、0分钟后这碗汤能喝了吗？设物体在时刻的温度为，从，温度从，注意到热量总是从温度高的物体向温度低的物体传导，因而，所以温度差恒正，又因物体将随时间而逐渐冷却，则温度的改变量为两边除以，并令得温度变化速度为解：为了解决这一问题，需要了解有关热力学的一些基本规律.热量总是从温度高的物体向温度低的物体传导的；在一定的温度范围内，一个物体的温度变化速度与这个物体的温度和其所在介质温度的差值成正比。其中是比例常数.从而得出描述物体冷却过程的微分方程模型为容易求出这个一阶微分方程初值问题的解为根据问题所给的条件知,当时，，得到将，代入，得(7.1.1)(7.1.2)从而得到这碗汤的温度随时间变化的函数

4、关系为于是，将代入计算得到再过20min汤的温度，这说明再过20min后这碗汤能喝了.不过，并不是所有的微分方程模型都可求出解析解。例如，看似简单的微分方程，自德国数学家WilhelmvonLeibniz提出100多年后才被法国数学家JosephLiouville证明它没有解析解，只能借助于数值的方法求数值解.(7.1.3)例7.2某地区发现一种有免疫性的传染病，为了控制疫情扩散对该地人群进行隔离处理.为了分析受感染人数的变化规律，需要建立描述传染病传播过程的数学模型.解设该地区的总人数为常数，任意时刻病人、健康人和病人治愈后移出感染系统的移出者的比例分别为，病人的日接触率，日治愈率，则容

5、易得出从时刻，病人和健康人的改变量为每个方程两边除以，并令，化简后得(7.1.4)其中（对任意的t）.式（7.1.4）就是描述病人和健康人的比例和随时间变化的微分方程模型，这是一个微分方程组的初值问题.但是，这一初值问题的解析解是无法求出的，因此不能直接利用和的解析式来分析和解决问题。在数学建模课程中学到的大量数学模型都是用微分方程形式给出的，各类微分方程本身和它们的解所具有的特性在常微分方程及数学物理方程中有所解释.虽然，求解微分方程有许多解析方法，但解析方法只能够求解一些特殊类型的方程，在实际应用中人们更关心的是某些特定的自变量在某一个定义范围内的一系列离散点上的近似值.这样一组

6、近似解称为微分方程在该范围内的数值解，寻找微分方程数值解的过程称为微分方程的数值解法。§2简单的数值方法与基本概念设在区域上连续，求满足其中是已知常数,这就是一阶常微分方程的初值问题.为使问题（7.2.1）的解存在、唯一且连续依赖初值，即初值问题（7.2.1）适定，还必须对右端项加以适当限制，通常要求关于是已知函数,且满足Lipschitz条件，即存在常数L，使7.2.1常微分方程初值问题(7.2.1)(7.2.2)对所有及成立。本章总假定满足条件（7.2.2）。1.Euler方法的导出与几何意义最简单的数值解法是Euler法。将区间作N等分，小区间的长度称为步长，点列称为节点，。

7、由已知初值，可算出在的导数。7.2.2Euler法及改进的Euler法其中，并略去二阶小量，得下面用3种方法导出Euler法.本章用表示函数在点的精确值，表示的近似值。就是的近似值。利用可算出，如此下去可算出在所有节点上的值，它的一般递推公式为(7.2.3)1）幂级数展开法利用Taylor展式(7.2.4)这就是Euler法。实际上，初值问题（7.2.1）的解是xy平面上过点的