常微分方程的应用与基本概念-课件

常微分方程1ppt课件新课导入1.课程介绍：2.最终成绩评定方式：最终成绩=平时成绩(30分)+期末考试成绩（70分）平时成绩=平时作业（5×3=15分）+出勤（5×3=15分）2ppt课件

1.什么是微分方程？含有自变量，未知函数及其导数的关系式称为微分方程。

2.什么是常微分方程？联系自变量、未知函数及未知函数导数（或微分）的关系式称为微分方程；常微分方程是只含有一个自变量的微分方程；

《常微分方程》则是研究含有一个自变量的微分方程的求解方法及解的性质的科学。

研究对象：只含一个自变量的微分方程。

基本任务：求解问题及解的性质。什么是常微分方程？3ppt课件3.学习《常微分方程》这门课的意义？

在自然界中，很多时候为了刻划客观对象的运动规律或变化规律，经常需要描述变量之间的函数关系。但针对实际问题，我们通常很难直接找到这种函数关系，却容易建立起变量所满足的微分方程。如果方程可求解，则可以得到描述客观对象运动规律或变化规律的函数关系。这就是说微分方程有着深刻而生动的实际背景。举例说明：4ppt课件1.《社会学》中的Malthus人口模型2.《种群生态学》中的Logistic虫口模型3.《医学》中的传染病动力学模型应用举例5.《电学》中的R-L-C电路模型4.《力学》中的变力作用模型5ppt课件10.《化学》中的流体浓度问题模型9.《光学》中的光学仪器成象原理模型8.《气象学》中的长期天气预报不可能问题7.《原子物理学》中的裂变问题模型6.《热学》中的物体冷却模型应用举例赝品的鉴定6ppt课件基本概念1.常微分方程与偏微分方程：3.线性常微分方程与非线性常微分方程：2.常微分方程的阶：5.常微分方程的显式解和隐式解：6.常微分方程的特解（特积分）与通解（通积分）：7.常微分方程定解问题的提法：4.齐次微分方程和非齐次微分方程：8.初等积分法求解常微分方程：7ppt课件作业与思考练习1：列车在平直的线路上以20米/秒的速度行驶,当制动时列车获得加速度—0.4米/秒,试建立该问题的微分方程模型，并求列车开始制动后多少时间列车才能停住？以及列车在这段时间内行驶了多少路程？

8ppt课件练习2：数学摆模型

数学摆是系于一根长度为

的线上而质量为

的质点M.在重力作用下,它在垂直于地面的平面上沿圆周运动.如图所示.试确定摆的运动方程.

9ppt课件

一定质量的镭，随着时间的变化，它的质量会减少。已知裂变速度与它的存余量成正比，假设时刻t镭的质量R。

分析：很难直接得到R与t函数关系，但却容易建立变量之间的一个微分方程：

这里为比例常数。借助于微积分知识，可对方程求解。原方程变形为例：镭的裂变模型10ppt课件两端积分得：亦即注意到初始条件：即得于是这里的比例常数可通过实验获得。例1的解决过程就是数学建模的过程。返回11ppt课件解设制动后t秒钟行驶s米，s=s(t),则注意到在t=0时，练习1：列车在平直的线路上以20米/秒的速度行驶,当制动时列车获得加速度—0.4米/秒,问开始制动后多少时间列车才能停住？以及列车在这段时间内行驶了多少路程？

12ppt课件于是代入条件后知：从而故从开始制动到列车完全停住共需秒。13ppt课件练习2：

求平面上过点(1,3)且每点切线斜率为横坐标2倍的曲线方程.解:设所求的曲线方程为由导数的几何意义,应有即又由条件:曲线过(1,3),即于是得故所求的曲线方程为:14ppt课件

例2.长期天气预报不可能问题为了刻划空气气流对流的规律，Lorenz得到了如下的微分方程组这里，x正比于对流运动强度，y正比于水平方向温度的变化，z正比于竖直方向温度的变化，参数为正常数。15ppt课件

通过对Lorenz方程组稳定性的讨论，可得到长期天气预报的不可能性。另外，人口模型、传染病模型、水波运动模型、光纤通讯中脉冲波胞的传输模型等等，均是微分方程模型。返回16ppt课件Kermack-McKendrick的SIR仓室模型所谓仓室模型就是针对某类传染病的传播特征和环境情况将某地区的人群（或某一种群）分成三类：

易感者类：其数量记为S(t),表示t时刻尚未染病但有可能被该类病菌或病毒感染的个体数；

染病者类：其数量记为I(t),表示t时刻已经染病并且有感染力的个体数；移出者类：其数量记为R(t),表示t时刻从染病中康复后移出的个体数；17ppt课件SIR模型的三个基本假设：（1）不考虑人口的出生与死亡，环境封闭（没有流入和流出），从而成员的总数始终保持一常数K,即：（2）一个染病者一旦与易感者接触就必然有一定的感染力。假设t时刻单位时间内已染病者传染易感者的数目与此时此刻易感者的数量S(t)成正比，比例系数为；（3）假设t时刻单位时间内从染病者中康复后移出的成员数与此时此刻已感染者的数量I(t)成正比，比例系数为，并且假设康复者具有永久免疫力，康复后不会再次被此疾病感染；18ppt课件S(易感者）I（染病者）R（移出者）SIR模型传染病示意图

对每一仓室的成员变化率建立平衡方程式，变得到了下面经典的SIR模型：返回19ppt课件

例3.物理冷却过程的数学模型

将某物体放置于空气中,在时刻时,测得它的温度为10分钟后测量得温度为

试决定此物体的温度

和时间

的关系,并计算20分钟后物体的温度.这里假设空气的温度保持在20ppt课件解:Newton冷却定律:1.热量总是从温度高的物体向温度低的物体传导;2.在一定的温度范围内,一个物体的温度变化速度与这一物体的温度与其所在的介质的温度之差成正比.

设物体在时刻

的温度为

根据导数的物理意义,则温度的变化速度为

由Newton冷却定律,得到

其中

为比例系数.此数学关系式就是物体冷却过程的数学模型.注意:此式子并不是直接给出

和

之间的函数关系,而只是给出了未知函数的导数与未知函数之间的关系式.如何由此式子求得

与

之间的关系式,以后再介绍.返回21ppt课件例4.R-L-C电路模型

如图所示的R-L-C电路.它包含电感L,电阻R,电容C及电源e(t).设L,R,C均为常数,e(t)是时间t的已知函数.试求当开关K合上后,电路中电流强度I与时间t之间的关系.

22ppt课件解:电路的Kirchhoff第二定律:

在闭合回路中,所有支路上的电压的代数和为零.

设当开关K合上后,电路中在时刻t的电流强度为I(t),则电流经过电感L,电阻R和电容的电压降分别为

其中Q为电量,于是由Kirchhoff第二定律,得到

因为

于是得到这就是电流强度I与时间t所满足的数学关系式.

返回23ppt课件例5.

数学摆模型

数学摆是系于一根长度为

的线上而质量为

的质点M.在重力作用下,它在垂直于地面的平面上沿圆周运动.如图所示.试确定摆的运动方程.

24ppt课件解:Newton第二定律:

取反时针运动方向为计量摆与铅垂线所成的角

的正方向.则由Newton第二定律,得到摆的运动方程为

附注1:如果研究摆的微小振动,即当

比较小时,可以取

的近似值

代入上式,这样就得到微小振动时摆的运动方程:

25ppt课件

为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。一般生态系统的分析可以通过一些简单模型的复合来研究，大家若有兴趣可以根据生态系统的特征自行建立相应的模型。

美丽的大自然

种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量，甚至允许它为可微变量，由此引起的误差将是十分微小的。

离散化为连续，方便研究Malthus模型与Logistic模型26ppt课件例.马尔萨斯（Malthus）模型

马尔萨斯在分析人口出生与死亡情况的资料后发现，人口净增长率r基本上是一常数，（r=b-d,b为出生率，d为死亡率），即：

或

（3.5）

（3.6）

（3.1）的解为：其中N0=N(t0)为初始时刻t0时的种群数。

马尔萨斯模型的一个显著特点：种群数量翻一番所需的时间是固定的。令种群数量翻一番所需的时间为T，则有：故27ppt课件模型检验

比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如，1961年世界人口数为30.6（即3.06×109），人口增长率约为2%，人口数大约每35年增加一倍。检查1700年至1961的260年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每34.6年增加一倍，两者也几乎相同。模型预测

假如人口数真能保持每34.6年增加一倍，那么人口数将以几何级数的方式增长。例如，到2510年，人口达2×1014个，即使海洋全部变成陆地，每人也只有9.3平方英尺的活动范围，而到2670年，人口达36×1015个，只好一个人站在另一人的肩上排成二层了。故马尔萨斯模型是不完善的。几何级数的增长Malthus模型实际上只有在群体总数不太大时才合理，到总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。所以Malthus模型假设的人口净增长率不可能始终保持常数，它应当与人口数量有关。返回28ppt课件例.Logistic模型

人口净增长率应当与人口数量有关，即：r=r(N)

从而有：（3.7）r(N)是未知函数，但根据实际背景，它无法用拟合方法来求。为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）对马尔萨斯模型引入一次项（竞争项），令r(N)=r-aN

此时得到微分方程：

或（3.8）

（3.8）被称为Logistic模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（Verhulst）首先提出的。一次项系数是负的，因为当种群数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。（3.8）可改写成：

（3.9）

(3.9)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为K（近似地将K看成常数），N表示当前的种群数量，K-N恰为环境还能供养的种群数量，（3.9）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（3.9）也被称为统计筹算律的原因。29ppt课件图3-5对（3.9）分离变量：两边积分并整理得：令N(0)=N0，求得：故（3.9）的满足初始条件N(0)=N0的解为：（3.10）易见：

N(0)=N0

，N(t)的图形请看图3.530ppt课件模型检验

用Logistic模型来描述种群增长的规律效果如何呢？1945年克朗皮克（Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（E·F·Gauss）也做