常微分方程教程课件

常微分方程教程

丁同仁、李承治编常微分方程教程

丁同仁、李承治编1主要参考书：东北师范大学数学系编写的高等学校教材《常微分方程》复旦大学数学系金福临等编写的《常微分方程》（上海科技出版社第二版）；南京大学数学系叶严谦等编写的《常微分方程讲义》；中山大学数学希望高雄等编写的《常微分方程》（高教第二版）.

主要参考书：东北师范大学数学系编写的高等学校教材《常微分方程2第一章基本概念§1.1微分方程模型第一章基本概念§1.1微分方程模型3例1求平面上过点(1,3)且每点切线斜率为横坐标2倍的曲线方程.解:设所求的曲线方程为由导数的几何意义,应有即又由条件:曲线过(1,3),即于是得故所求的曲线方程为:例1求平面上过点(1,3)且每点切线斜率为横坐标2倍的曲4例2物理冷却过程的数学模型

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

试决定此物体的温度

和时间

的关系,并计算20分钟后物体的温度.这里假设空气的温度保持在例2物理冷却过程的数学模型将某物体放置于空气中,5解:Newton冷却定律:1.热量总是从温度高的物体向温度低的物体传导;2.在一定的温度范围内,一个物体的温度变化速度与这一物体的温度与其所在的介质的温度之差成正比.

设物体在时刻

的温度为

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

由Newton冷却定律,得到

其中

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

和

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

与

之间的关系式,以后再介绍.解:Newton冷却定律:设物体在时刻6例3R-L-C电路

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

例3R-L-C电路如图所示的R-L-C电路.它7解:电路的Kirchhoff第二定律:

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

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

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

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

解:电路的Kirchhoff第二定律:设8例4

数学摆

数学摆是系于一根长度为

的线上而质量为

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

解:Newton第二定律:

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

的正方向.则由Newton第二定律,得到摆的运动方程为附注1:如果研究摆的微小振动,即当

比较小时,可以取的近似值

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

例4数学摆数学摆是系于一根长度为9附注2:假设摆是在一个有粘性的介质中作摆动,如果阻力系数为

则摆的运动方程为:附注3:假设摆还沿着摆的运动方向受到一个外力F(t)的作用,则摆的运动方程为:附注2:假设摆是在一个有粘性的介质中作摆动,如果阻力系10例5追击问题

假设敌舰从原点出发以速度a沿y轴正向行驶，同时有鱼雷从点出发以速度b追击敌舰，求鱼雷的运动路线（称为追线）。

解:设时间t，敌舰到达R(0,at)，鱼雷到达P(x,y),追赶时,鱼雷总是向敌舰所在的位置方向追赶，即dt时刻内，鱼雷的运动方向为，设在时刻t+dt，鱼雷到达则P,P’,R三点一线，如右图所示。例5追击问题解:设时间t，敌舰到达R(0,at)，鱼雷11从而即此为带有微分的方程.从而即此为带有微分的方程.12战争分类：正规战争，游击战争，混合战争只考虑双方兵力多少和战斗力强弱兵力因战斗及非战斗减员而减少，因增援而增加战斗力与射击次数及命中率有关第一次世界大战Lanchester提出预测战役结局的模型例6正规战与游击战战争分类：正规战争，游击战争，混合战争只考虑双方兵力多少和战13每方战斗减员率取决于双方的兵力和战斗力每方非战斗减员率与本方兵力成正比甲乙双方的增援率为u(t),v(t)f,g

取决于战争类型x(t)~甲方兵力，y(t)~乙方兵力模型假设:每方战斗减员率取决于双方的兵力和战斗力每方非战斗减员率与14忽略非战斗减员假设没有增援f(x,y)=ay,a~乙方每个士兵的杀伤率a=rypy,ry~射击率，

py~命中率如正规战争：双方均以正规部队作战甲方战斗减员率只取决于乙方的兵力和战斗力b~甲方每个士兵的杀伤率忽略非战斗减员假设没有增援f(x,y)=ay,a15事实上，很多问题的处理都归结为一个(或几个）带有微分的方程，如：传染病模型经济增长模型正规战与游击战药物在体内的分布与排除香烟过滤嘴的作用人口预测和控制烟雾的扩散与消失万有引力定律的发现事实上，很多问题的处理都归结为一个(或几个）带有微分的方程，16

把联系自变量、未知函数及未知函数导数（或微分）的关系式称为微分方程.

§1.2基本概念如果其中的未知函数只与一个自变量有关，则称为常微分方程；如果未知函数是两个或两个以上自变量的函数，并且在方程中出现偏导数，则称为偏微分方程.本书主要介绍常微分方程.有时就简称微分方程或方程.未知函数最高阶导数的阶数，称为方程的阶.

定义1.1:把联系自变量、未知函数及未知函数导数§1.2基17例如下面的方程都是微分方程一阶常微分方程一阶常微分方程二阶常微分方程四阶常微分方程二阶偏微分方程例如下面的方程都是微分方程一阶常微分方程一阶常微分方程二阶常18(6)式是n阶常微分方程的一般形式（n阶隐式方程）。

n阶显式方程的一般形式为

n阶常微分方程在方程(6)中，如果左端函数F对未知函数y和它的各阶导数y′,y″,…,的全体而言是一次的，则称为线性常微分方程，否则称它为非线性常微分方程.这样，一个以y为未知函数，以x为自变量的n阶线性微分方程具有如下形式：

(6)式是n阶常微分方程的一般形式（n阶隐式方程）。19通解与特解微分方程的解就是满足方程的函数，可定义如下:定义1.2在区间

上连续，且有直到n阶的代入方程(6)，得到在区间

则称为方程(6)在区间

上的一个解.导数.如果把

上关于

的恒等式,即设函数通解与特解微分方程的解就是满足方程的函数，可定义如下:定义120从定义1.2可以直接验证：

1.函数是方程在区间(-∞，+∞)上的解，其中C是任意的常数.2.函数是方程和是独立的在区间(-∞，+∞)上的解，其中任意常数.

从定义1.2可以直接验证：1.函数213.函数是方程在区间(-∞，+∞)上的解，其中是任意的常数.从上面的例子中，可以看到一个重要事实，那就是微分方程的解中可以包含任意常数，其中任意常数的个数可以多到与方程的阶数相等（也可以不含任意常数）.3.函数22把n阶常微分方程定义1.3的含有n个独立的任意常数的解称为该方程的通解，

如果方程的解不包含任意常数，则称它为特解.

由隐式表出的通解称为通积分

把n阶常微分方程定义1.3的含有n个独立的任意常数的解称为该23例如函数是方程在区间(-∞，+∞)上的通解，其中是任意的常数.而

y=1而是方程的一个特解。

例如函数24例自由落体

设质量为m的物体，在时间t=0时，在距地面高度为H处以初始速度

垂直地面下落，求此物体下落时距离与时间的关系.解:如图建立坐标系.设y=y(t)为t时刻物体的位置坐标.则易得物体下落所满足的方程为y’’=-g（*）其中g是重力加速度.初值问题例自由落体

设质量为m的物体，在时间t=0时，在距25这表明方程(*)有无数个解，原因是未考虑初始状态。为了确定相应的运动，考虑初始条件：于是，得到所要求的自由落体的运动方程为得是通解，其中是两个任意常数。容易验证这表明方程(*)有无数个解，原因是未考虑初始于是，得到所要求26综上所述，自由落体的问题可归结为求如下初值问题（柯西问题（cauchy)）的解：对于一般的n阶方程，初值问题（柯西问题）的一般提法是：

综上所述，自由落体的问题可归结为求如下初值问题（对于一般的n27例求方程

的满足初值条件的解.解:方程通解为求导数后得将初值条件代入，得到方程组

例求方程

的满足初值条件的28解出

故所求特解为解出故所求特解为29§1.2微分方程及其解的几何解释为了便于研究方程解的性质，我们常常考虑解的图象.一阶方程的一个特解的图象是平面上的一条曲线，称为方程的积分曲线，而通解的图象是平面上的一族曲线，称为积分曲线族.

线素，线素场，方向场，等斜线

§1.2微分方程及其解的几何解释为了便于研究方程解的性质30常微分方程教程

丁同仁、李承治编常微分方程教程

丁同仁、李承治编31主要参考书：东北师范大学数学系编写的高等学校教材《常微分方程》复旦大学数学系金福临等编写的《常微分方程》（上海科技出版社第二版）；南京大学数学系叶严谦等编写的《常微分方程讲义》；中山大学数学希望高雄等编写的《常微分方程》（高教第二版）.

主要参考书：东北师范大学数学系编写的高等学校教材《常微分方程32第一章基本概念§1.1微分方程模型第一章基本概念§1.1微分方程模型33例1求平面上过点(1,3)且每点切线斜率为横坐标2倍的曲线方程.解:设所求的曲线方程为由导数的几何意义,应有即又由条件:曲线过(1,3),即于是得故所求的曲线方程为:例1求平面上过点(1,3)且每点切线斜率为横坐标2倍的曲34例2物理冷却过程的数学模型

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

试决定此物体的温度

和时间

的关系,并计算20分钟后物体的温度.这里假设空气的温度保持在例2物理冷却过程的数学模型将某物体放置于空气中,35解:Newton冷却定律:1.热量总是从温度高的物体向温度低的物体传导;2.在一定的温度范围内,一个物体的温度变化速度与这一物体的温度与其所在的介质的温度之差成正比.

设物体在时刻

的温度为

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

由Newton冷却定律,得到

其中

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

和

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

与

之间的关系式,以后再介绍.解:Newton冷却定律:设物体在时刻36例3R-L-C电路

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

例3R-L-C电路如图所示的R-L-C电路.它37解:电路的Kirchhoff第二定律:

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

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

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

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

解:电路的Kirchhoff第二定律:设38例4

数学摆

数学摆是系于一根长度为

的线上而质量为

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

解:Newton第二定律:

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

的正方向.则由Newton第二定律,得到摆的运动方程为附注1:如果研究摆的微小振动,即当

比较小时,可以取的近似值

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

例4数学摆数学摆是系于一根长度为39附注2:假设摆是在一个有粘性的介质中作摆动,如果阻力系数为

则摆的运动方程为:附注3:假设摆还沿着摆的运动方向受到一个外力F(t)的作用,则摆的运动方程为:附注2:假设摆是在一个有粘性的介质中作摆动,如果阻力系40例5追击问题

假设敌舰从原点出发以速度a沿y轴正向行驶，同时有鱼雷从点出发以速度b追击敌舰，求鱼雷的运动路线（称为追线）。

解:设时间t，敌舰到达R(0,at)，鱼雷到达P(x,y),追赶时,鱼雷总是向敌舰所在的位置方向追赶，即dt时刻内，鱼雷的运动方向为，设在时刻t+dt，鱼雷到达则P,P’,R三点一线，如右图所示。例5追击问题解:设时间t，敌舰到达R(0,at)，鱼雷41从而即此为带有微分的方程.从而即此为带有微分的方程.42战争分类：正规战争，游击战争，混合战争只考虑双方兵力多少和战斗力强弱兵力因战斗及非战斗减员而减少，因增援而增加战斗力与射击次数及命中率有关第一次世界大战Lanchester提出预测战役结局的模型例6正规战与游击战战争分类：正规战争，游击战争，混合战争只考虑双方兵力多少和战43每方战斗减员率取决于双方的兵力和战斗力每方非战斗减员率与本方兵力成正比甲乙双方的增援率为u(t),v(t)f,g

取决于战争类型x(t)~甲方兵力，y(t)~乙方兵力模型假设:每方战斗减员率取决于双方的兵力和战斗力每方非战斗减员率与44忽略非战斗减员假设没有增援f(x,y)=ay,a~乙方每个士兵的杀伤率a=rypy,ry~射击率，

py~命中率如正规战争：双方均以正规部队作战甲方战斗减员率只取决于乙方的兵力和战斗力b~甲方每个士兵的杀伤率忽略非战斗减员假设没有增援f(x,y)=ay,a45事实上，很多问题的处理都归结为一个(或几个）带有微分的方程，如：传染病模型经济增长模型正规战与游击战药物在体内的分布与排除香烟过滤嘴的作用人口预测和控制烟雾的扩散与消失万有引力定律的发现事实上，很多问题的处理都归结为一个(或几个）带有微分的方程，46

把联系自变量、未知函数及未知函数导数（或微分）的关系式称为微分方程.

§1.2基本概念如果其中的未知函数只与一个自变量有关，则称为常微分方程；如果未知函数是两个或两个以上自变量的函数，并且在方程中出现偏导数，则称为偏微分方程.本书主要介绍常微分方程.有时就简称微分方程或方程.未知函数最高阶导数的阶数，称为方程的阶.

定义1.1:把联系自变量、未知函数及未知函数导数§1.2基47例如下面的方程都是微分方程一阶常微分方程一阶常微分方程二阶常微分方程四阶常微分方程二阶偏微分方程例如下面的方程都是微分方程一阶常微分方程一阶常微分方程二阶常48(6)式是n阶常微分方程的一般形式（n阶隐式方程）。

n阶显式方程的一般形式为

n阶常微分方程在方程(6)中，如果左端函数F对未知函数y和它的各阶导数y′,y″,…,的全体而言是一次的，则称为线性常微分方程，否则称它为非线性常微分方程.这样，一个以y为未知函数，以x为自变量的n阶线性微分方程具有如下形式：

(6)式是n阶常微分方程的一般形式（n阶隐式方程）。49通解与特解微分方程的解就是满足方程的函数，可定义如下:定义1.2在区间

上连续，且有直到n阶的代入方程(6)，得到在区间

则称为方程(6)在区间

上的一个解.导数.如果把

上关于

的恒等式,即设函数通解与特解微分方程的解就是满足方程的函数，可定义如下:定义150从定义1.2可以直接验证：

1.函数是方程在区间(-∞，+∞)上的解，其中C是任意的常数.2.函数是方程和是独立的在区间(-∞，+∞)上的解，其中任意常数.

从定义1.2可以直接验证：1.函数513.函数是方程在区间(-∞，+∞)上的解，其中是任意的常数.从上面的例子中，可以看到一个重要事实，那就是微分方程的解中可以包含任意常数，其中任意常数的个数可以多到与方程的阶数相等（也可以不含任意常数）.3.函数52把n阶常微分方程定义1.3的含有n个独立的任意常数的解称为该方程