(整理)微分方程及其应用

1、精品文档精品文档第九章微分方程及其应用 9.19.1 微分方程及其相关概念所谓微分方程，就是含有自变量、自变量的未知函数以及未知函数的导数（或微分）的方程。 例如，以下各式都是微分方程：型+P(x)y =Q(x).dxF(x, y, y y(n)=0. .只含一个自变量的微分方程，称为常微分方程，自变量多于一个的称为偏微分方程。本章只研究常微分方程, 因而以后各节提到微分方程时均指常微分方程。微分方程中所含有的未知函数最高阶导数的阶数，称为该微分方程的阶。例如，、为一阶方程，、为二阶方程，而为n n 阶方程。微分方程中可以不含有自变量或未知函数，但不能不含有导数，否则就不成为微分方程。微分方程

2、与普通代数方程有着很大的差别，建立微分方程的目的是寻找未知函数本身。如果P196P196有一个函数满足微分方程，即把它代入微分方程后，使方程变成（对自变量的）恒等式，这个函数就叫做若方程解中含有独立的任意常数的个数等于微分方程的阶数，则称此解为微分方程的通解，例如y =-x3二就是的通解。3从通解中取定任意常数的一组值所得到的解，称为微分方程的特解。例如y=-x3二就是的一个特解。3用来确定通解中任意常数值的条件称为定解条件，当自变量取某个值时，给出未知函数及其导数的相应值 的条件称为初始条件。在本章中，我们遇到的用来确定任意常数值的条件一般为初始条件。例如，如果的初J1313始条件为y 0

3、1=寫，则在代入到通解yx3c后，可以求得c=庶，从而得到特解y x3*。33一般的，因为n阶微分方程的通解中含有n个独立的任意常数。需要有n个（一组）定解条件，所以n阶方程的初始条件为：y(x)= y, y (x )= %, y (x)= y?,，y(n，)=y.其中yo, y1, y2,为n个给定常数。微分方程的解所对应的几何图形叫做微分方程的积分曲线。通解的几何图形是一族积分曲线，特解所对应 的几何图形是一族积分曲线中的某一条。1例如，方程的积分曲线族如图9 91所示。其中y x3二就是满足初始条件y 0 -二的特解。3dy2xdxd2xm 2dt2hx$ kx =dtf (t).竽dt

4、29sin“0. .dt l微分方程的解。例如-x3显然是的解，因为313、d(-x)dx精品文档精品文档9.29.2微分方程的经典案例例 1 1 自由落体运动的规律自由落体运动是指物体在仅受到地球引力的作用下，初速度为零的运动。根据经典力学的牛顿第二定律：物体动量变化的大小与它所受到的外力成正比，其方向与外力的方向一致。当物体的运动速度V的绝对值不大（与光速=3=3105km/skm/s 相比较）时，其质量m可以是一恒量。于是这一运动定律能表达成ddvmv二F,或m F（1 1）dtdt其中F表示物体所受外力的合力。对于仅受到地球引力作用的自由落体的运动，则有：-1- dS -2F =mg,

5、v这里g表示重力加速度，其大小一般取为：g =9.8m/ s2；dtS表示自由落体运动的路程，其大小以S表示之。-_- _ dS注意到S的方向与g的方向一致，将F =mg,v代入式后得到自由落体运动立场dt大小变化的规律：mS二mg或d,g(2 2)dt2dt2运动规律式表示一个微分方程问题。等式（2 2）的左端是路程大小S的二次微商它的右端是常数g。这里S和g之间不是普通的函数关系，而是二微商的关系。例 2 2 单摆运动单摆又称为钟摆或数学摆。所谓单摆运动是指一质量为m00 的小球，用长度为I的柔软细绳拴住，细绳的一端固定在某点o o 处。小球在铅垂平面内运动，略去空气的阻力和细绳在 O O

6、 点处的摩擦力。并且认为细绳的长度I不变，仅考虑地球的引力和细绳对小球的拉力（见图 9 9 2 2）。在铅垂平面内引进以 O O 为坐标原点的极坐标系统，由于细绳长度不变且细绳总是直的，所以小球的位置用一个坐标 二t就能表示。这里 二表示细绳I和铅垂方向之间的夹角。铅垂方向即是小球的平衡方向，它对应的二为零。作用在小球上的地球引力的大小f为mg,其方向铅垂向下。重力沿细绳方向的分力的大小为mgcosr， 其方向沿细绳指向外。 这个力与小球运动所需要的向心力刚好平衡。所以mg sin v，它的方向与精品文档精品文档小球沿细绳方向没有运动。重力在垂直于细绳方向的分力的大小为角二增加的方向相反。根据

7、牛顿第二定律得到单摆运动的规律为:mv一-mg s乔dta(3)精品文档精品文档关系式（4 4）是包含及其二接微商的方程，并且 二不是线性而是非线性地出现在方程 中（以si nr这种非线性形式）。从方程（4 4）来求出二随着时间变化规律的分析表达式是困难的。当|州比较小时，对微分方程（4 4）能够进行线性化出处理,用二来近似si nr。这样得到式（4 4）的线性化微分方程:近似规律。通常将式5”写成如下的规范形式:其中k2=g。l例 3 3 真空中的抛射体运动在真空中运动的抛射体，它的运动规律十分复杂。这里仅考虑在真空中抛射体的运动 规律。即忽略抛射体所受的空气阻力，而仅考虑质量为m的抛射体受

8、地球引力作用而引起的运动。取一直角坐标系Oxyz，Ox轴沿水平方向；Oy轴垂直于Ox轴；Oz轴垂直于xOy平面，并与Ox轴、Oy轴一起组成右手坐标系。依牛顿第二定律，抛射体的运动规律为:d2x dt2d2zm2dt2d2ydt2抛射体的初始状态取为:x 0二y0二z0 =0;dx dy c dz 二=vcosa,：= 0,二 =sinaL_dtdtdt根据圆周运动规律有:于是从式（3 3）得出:l jdtd2.-gsin J即用 v 代替sinv，或者说,dt2在相同初始条件下服从微分方程5”求得的二随时间t变化的规律二t是单摆运动的dSdt2kS - o(5(5)-mg(6)精品文档精品文档

9、其中Vo是抛射体的初始速度，位于xOy平面内，vo表示vo的大小；表示vo与水平方向 （即Ox轴）之间的夹角（见图 9-49-4）。例 4 4 深水炸弹的水下运动一质量为m的深水炸弹，从高为h m处自由下落到海中。 这里不考虑深水炸弹在水平 方向的运动，而仅考虑它在铅直方向的运动。由经典力学知：物体由高为至海平面时，其铅垂方向的速度v0为：v0二2gh这里g为重力加速度。按如下方式取定坐标系：坐标原点O取在海平面上某处，Ox轴沿铅垂向下，（见图9-59-5 ）。深水炸弹m自高度为h m处自由下落至海平面的时间为t0。于是深水炸弹的初始状态为:深水炸弹在海中运动时，我们不考虑海水对它的浮力，这时

10、炸弹受到两个力的作用， 是地球引力mg，其方向铅垂向下；另一个是海水对炸弹的摩擦力。这个摩擦力是很复杂的，它和炸弹的形状、速度等因素有关，这里近似的认为摩擦力的大小和炸弹的速度v成正比，比例系数即摩擦系数u为常数。摩擦力的方向与炸弹的速度方向相反，因而是铅垂向上h m处自由下落的。于是摩擦力f能表示为:dx=u v = -udt精品文档精品文档根据牛顿第二定律知深水炸弹在水下运动的规律为：m史u v mg dtd2x u dx2dt m dt例 5 5 放射性元素的衰变放射性元素铀由于不断的有原子放射出微粒子而变成其他元素，铀的含量就不断减少，这种现象叫做衰变。由原子物理学知道，铀的衰变速度与

11、当时未衰变的原子的含量M成正比。已知t =0时铀的含量为Mo，求在衰变过程中铀含量M t随时间t变化的规律。解铀的衰变速度就是M t对时间t的导数如，由于铀的衰变速度与其含量成正dt比，故得微分方程dM dt =其中 o是常数，叫做衰变系数。前置符号是由于当t增加时M单调减少，即0的缘故。dt按题意，初始条件为M |tz0二M0例 6 6例 7 7指数增长模型（马尔萨斯人口模型）英国人口学家马尔萨斯（MalthusMalthus，1766-18341766-1834）根据百余年的人口统计资料，于17981798年提出了著名的人口指数增长模型。这个模型的基本假设是：人口的增长率是常数，或者 说，

12、单位时间内人口的增长量与当时的人口成正比。记时刻t的人口为x t，当考察一个国家或一个很大地区的人口时，x t是很大的整数。为了利用微积分这一数学工具，将x t视为连续、可微函数。记初始时刻t=0的人口为X。，人口增长率为 r r，r r 是单位时间内x t的增量与x t的比例系数。于是，x t满足如下的微工dxrx分方程：dtx 0 =xe表明人口将按指数规律无限增长r 0。例 8 8例 9 9 阻滞增长模型（LogisticLogistic 模型）(7(7)精品文档精品文档例 6 6 中的指数增长模型在 1919 世纪前比较符合人口增长情况，但从1919 世纪以后，就与人口事实上的增长情况产生了较大的差异。产生上述现象的主要原因是，随着人口的增加，自然资源，环境条件等因素对人口继 续增长的阻滞作用越来越显著。如果当人口较少时（相对于资源而言）人口增长率还可以 看作常数的话，那么当人口增加到一定数量后，增长率就会随着人口的继续增加而逐渐减 少。为了使人口预报特别是长期预报更好的符合实际情况，必须修改指数增长模型关于人 口增长率是常数这个基本假设。