微分方程及其应用

1、.第九章 微分方程及其应用§9.1 微分方程及其相关概念所谓微分方程，就是含有自变量、自变量的未知函数以及未知函数的导数（或微分）的方程。例如，以下各式都是微分方程： . . .只含一个自变量的微分方程,称为常微分方程,自变量多于一个的称为偏微分方程。本章只研究常微分方程，因而以后各节提到微分方程时均指常微分方程。微分方程中所含有的未知函数最高阶导数的阶数，称为该微分方程的阶。例如，、为一阶方程，、为二阶方程，而为n阶方程。微分方程中可以不含有自变量或未知函数，但不能不含有导数，否则就不成为微分方程。微分方程与普通代数方程有着很大的差别，建立微分方程的目的是寻找未知函数本身。如果P1

2、96有一个函数满足微分方程，即把它代入微分方程后，使方程变成（对自变量的）恒等式，这个函数就叫做微分方程的解。例如显然是的解，因为。若方程解中含有独立的任意常数的个数等于微分方程的阶数，则称此解为微分方程的通解，例如就是的通解。从通解中取定任意常数的一组值所得到的解，称为微分方程的特解。例如就是的一个特解。用来确定通解中任意常数值的条件称为定解条件，当自变量取某个值时，给出未知函数及其导数的相应值的条件称为初始条件。在本章中，我们遇到的用来确定任意常数值的条件一般为初始条件。例如，如果的初始条件为，则在代入到通解后，可以求得，从而得到特解。一般的，因为阶微分方程的通解中含有个独立的任意常数。需

3、要有个（一组）定解条件，所以阶方程的初始条件为：其中为个给定常数。微分方程的解所对应的几何图形叫做微分方程的积分曲线。通解的几何图形是一族积分曲线，特解所对应的几何图形是一族积分曲线中的某一条。例如，方程的积分曲线族如图9所示。其中就是满足初始条件的特解。§9.2 微分方程的经典案例例1 自由落体运动的规律自由落体运动是指物体在仅受到地球引力的作用下，初速度为零的运动。根据经典力学的牛顿第二定律：物体动量变化的大小与它所受到的外力成正比，其方向与外力的方向一致。当物体的运动速度的绝对值不大（与光速=3km/s相比较）时，其质量可以是一恒量。于是这一运动定律能表达成或 （1）其中表示物

4、体所受外力的合力。对于仅受到地球引力作用的自由落体的运动，则有：这里表示重力加速度，其大小一般取为：；表示自由落体运动的路程，其大小以表示之。注意到的方向与的方向一致，将代入式后得到自由落体运动立场大小变化的规律： 或 （2）运动规律式表示一个微分方程问题。等式（2）的左端是路程大小的二次微商它的右端是常数。这里和之间不是普通的函数关系，而是二微商的关系。例2 单摆运动单摆又称为钟摆或数学摆。所谓单摆运动是指一质量为>0的小球，用长度为的柔软细绳拴住，细绳的一端固定在某点O处。小球在铅垂平面内运动，略去空气的阻力和细绳在O点处的摩擦力。并且认为细绳的长度不变，仅考虑地球的引力和细绳对小球

5、的拉力（见图92）。在铅垂平面内引进以O为坐标原点的极坐标系统，由于细绳长度不变且细绳总是直的，所以小球的位置用一个坐标就能表示。这里表示细绳和铅垂方向之间的夹角。铅垂方向即是小球的平衡方向，它对应的为零。作用在小球上的地球引力的大小为，其方向铅垂向下。重力沿细绳方向的分力的大小为，其方向沿细绳指向外。这个力与小球运动所需要的向心力刚好平衡。所以小球沿细绳方向没有运动。重力在垂直于细绳方向的分力的大小为，它的方向与角增加的方向相反。根据牛顿第二定律得到单摆运动的规律为： （3）根据圆周运动规律有：于是从式（3）得出： （4）关系式（4）是包含及其二接微商的方程，并且不是线性而是非线性地出现在方

6、程中（以这种非线性形式）。从方程（4）来求出随着时间变化规律的分析表达式是困难的。当比较小时，对微分方程（4）能够进行线性化出处理，即用代替，或者说，用来近似。这样得到式（4）的线性化微分方程： 在相同初始条件下服从微分方程求得的随时间变化的规律是单摆运动的近似规律。通常将式写成如下的规范形式： （5）其中 。例3 真空中的抛射体运动在真空中运动的抛射体，它的运动规律十分复杂。这里仅考虑在真空中抛射体的运动规律。即忽略抛射体所受的空气阻力，而仅考虑质量为的抛射体受地球引力作用而引起的运动。取一直角坐标系，轴沿水平方向；轴垂直于轴；轴垂直于平面，并与轴、轴一起组成右手坐标系。依牛顿第二定律，抛射

7、体的运动规律为： （6）抛射体的初始状态取为： 其中是抛射体的初始速度，位于平面内，表示的大小；表示与水平方向（即轴）之间的夹角（见图9-4）。例4 深水炸弹的水下运动一质量为的深水炸弹，从高为处自由下落到海中。这里不考虑深水炸弹在水平方向的运动，而仅考虑它在铅直方向的运动。由经典力学知：物体由高为处自由下落至海平面时，其铅垂方向的速度为：这里为重力加速度。按如下方式取定坐标系：坐标原点取在海平面上某处， 轴沿铅垂向下，（见图9-5）。深水炸弹自高度为处自由下落至海平面的时间为。于是深水炸弹的初始状态为： 深水炸弹在海中运动时，我们不考虑海水对它的浮力，这时炸弹受到两个力的作用，：一是地球引力

8、，其方向铅垂向下；另一个是海水对炸弹的摩擦力。这个摩擦力是很复杂的，它和炸弹的形状、速度等因素有关，这里近似的认为摩擦力的大小和炸弹的速度成正比，比例系数即摩擦系数为常数。摩擦力的方向与炸弹的速度方向相反，因而是铅垂向上的。于是摩擦力能表示为： 根据牛顿第二定律知深水炸弹在水下运动的规律为：或 （7）例5 放射性元素的衰变放射性元素铀由于不断的有原子放射出微粒子而变成其他元素，铀的含量就不断减少，这种现象叫做衰变。由原子物理学知道，铀的衰变速度与当时未衰变的原子的含量成正比。已知时铀的含量为，求在衰变过程中铀含量随时间变化的规律。解 铀的衰变速度就是对时间的导数，由于铀的衰变速度与其含量成正比

9、，故得微分方程 其中是常数，叫做衰变系数。前置符号是由于当增加时单调减少，即的缘故。按题意，初始条件为 例6 指数增长模型（马尔萨斯人口模型）英国人口学家马尔萨斯（Malthus，1766-1834）根据百余年的人口统计资料，于1798年提出了著名的人口指数增长模型。这个模型的基本假设是：人口的增长率是常数，或者说，单位时间内人口的增长量与当时的人口成正比。记时刻的人口为，当考察一个国家或一个很大地区的人口时，是很大的整数。为了利用微积分这一数学工具，将视为连续、可微函数。记初始时刻的人口为，人口增长率为r，r是单位时间内的增量与的比例系数。于是，满足如下的微分方程： 表明人口将按指数规律无限增长。例7 阻滞增长模型（Logistic模型）例6中的指数增长模型在19世纪前比较符合人口增长情况，但从19世纪以后，就与人口事实上的增长情况产生了较大的差异。产生上述现象的主要原因是，随着人口的增加，自然资源，环境条件等因素对人口继续增长的阻滞作用越来越显著。如果当人口较少时（相对于资源而言）人口增长率还可以看作常数的话，那么当人口增加到一定数量后，增长率就会随着人口的继续增加而逐渐减少。为了使人口预报特别是长期预报更好的符合实际情况，必须修改指数增长模型关于人口增长率