20113564 胡骐薪工商1112微分方程的基本应用微分方程是数学的重要分支, 用微分方程来刻画许多自然科学、经济科学甚至社会科学领域中的一些规律,这是微分方程应用的重要领域,也是其发展的动力.在这里我重点介绍了几个利用微分方程常来解决的问题的例子,从中我们可以了解到微分方程用的广泛性以及解决具体问题时常采用的一般步骤.微分方程是与微积分一起形成发展起来的重要数学分支,已有悠久的历史,早在17~18世纪,牛顿、莱布尼兹、贝努里和拉格朗日等人在研究力学和几何学中就提出了微分方程【1,2】.随着科学的发展,它在力学、电学、天文学和其他数学物理领域内的应用不断获得成功,有力地推动了这些学科的发展,已成为研究自然科学和社会科学的一个强有力工具.如今,微分方程仍继续保持着进一步发展的活力,其主要原因是它的根源深扎在各种实际问题之中,许多实际问题可以通过建立微分方程模型得以解决.常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的. 数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具.微分方程可以精确地表述事物变化所遵循的基本规律. 随着微分方程的理论的逐步完善，只要列出相应的微分方程并找到解方程的方法, 微分方程也就成了最有生命力的数学分支. 事实上，大部分的常微分方程求不出十分精确的解，而只能得到近似解. 当然，这个近似解的精确程度是比较高的.现在，常微分方程在很多学科领域内有着重要的应用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等. 这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题. 应该说，应用常微分方程理论已经取得了很大的成就. 解常微分方程大致有分离变量法、变量替换法、常数变易法以及积分因子法等等，其中，积分因子法尤为重要，本论文主要讨论积分因子存在条件及其解法，通过积分因子使常微分方程化为全微分方程形式来求解.微分方程在科学技术和实际生活中都有着广泛的应用。

应用微分方程解决实际问题,其实就是建立微分方程数学模型,通过建立微分方程、确定定解条件、求解及对解的分析可以揭示许多自然界和科学技术中的规律.应用微分方程解决具体问题的主要步骤:(1)分析问题,将实际问题抽象,设出未知函数，建立微分方程,并给出合理的定解条件;(2)求解微分方程的通解及满足定解条件的特解,或由方程讨论解的性质;(3)由所求得的解或解的性质,回到实际问题,解释该实际问题,得出客观规律. 微分方程的应用举例几何问题1.等角轨线我们来求这样的曲线或曲线族,使得它与某已知曲线族的每一条曲线相交成给定的角度.这样的曲线轨线已知曲线的等角轨线.当所给定的角为直角时,等角轨线就轨线正交轨线.等角轨线在很多学科（如天文,气象等）中都有应用.下面就来介绍等角轨线的方法.首先把问题进一步提明确一些.设在（x,y ）平面上,给定一个单参数曲线族（C ）:()0,,=c y x ϕ求这样的曲线l ,使得l 与(C)中每一条曲线的交角都是定角α .设l 的方程为1y =)(1x y .为了求)(1x y ,我们先来求出)(1x y 所对应满足的微分方程,也就是要求先求得x , 1y ,'1y 的关系式.条件告诉我们l 与（C ）的曲线相交成定角α,于是,可以想象,1y 和'1y 必然应当与（C ）中的曲线y =)(x y 及其切线的斜率'y 有一个关系.事实上,当α≠2π时,有 k y y y y ==+-αtan 1'1'''1 或1'1'1'+-=ky k y y当α=2π时,有 '1'1y y -=又因为在交点处,)(x y =)(1x y ,于是,如果我们能求得x , 1y ,'1y 的关系()0,,'=y y x F采用分析法.设y =)(x y 为（C ）中任一条曲线,于是存在相应的C,使得()()0,,≡c x y x ϕ因为要求x ,y, '1y 的关系,将上式对x 求导,得()()()()()0,,,,'''≡+x y c x y x c x y x y xϕϕ 这样,将上两式联立,即由()()()⎩⎨⎧=+=0,,,,0,,'''y c y x c y x c y x y x ϕϕϕ 消去C,就得到()()x y x y x ',,所应当满足的关系()0,,'=y y x F这个关系称为曲线族（C ）的微分方程.于是,等角轨线（α≠2π）的微分方程就是01,,'1'11=⎥⎦⎤⎢⎣⎡+-ky k y y x F 而正交轨线的微分方程为01,,'11=⎥⎦⎤⎢⎣⎡-y y x F为了避免符号的繁琐,以上两个方程可以不用1y ,而仍用y,只要我们明确它是所求的等角轨线的方程就行了.为了求得等角轨线或正交轨线,我们只需求上述两个方程即可. 例1 求直线束cx y =的等角轨线和正交轨线.解 首先求直线束cx y =的微分方程.将cx y =对x 求导,得'y =C,由⎩⎨⎧==cy cxy '消去C,就得到cx y =的微分方程xy dx dy = 当α≠2π时,由（2.16）知道,等角轨线的微分方程为 x y dxdy kkdx dy =+-1 或kydxxdy ydy xdx -=+ 及22221y x ydxxdy k yx ydy xdx +-⋅=++ 即22211⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛⋅=++x y x y d k y x ydy xdx 积分后得到()c xyk y x ln arctan 1ln 2122+=+ 或xycey x arctan 2122=+如果α=2π,由（2.17）可知,正交轨线的微分方程为 x y dxdy =-1 即yx dx dy -= 或 0=+ydy xdx 故正交轨线为同心圆族222c y x =+.例2 抛物线的光学问题在中学平面解析几何中已经指出,汽车前灯和探照灯的反射镜面都取为旋转抛物面,就是将抛物线绕对称轴旋转一周所形成的曲面.将光源安置在抛物线的焦点处,光线经镜面反射,就成为平行光线了.这个问题在平面解析几何中已经作了证明,现在来说明具有前述性质的曲线只有抛物线,由于对称性,只有考虑在过旋转轴的一个平面上的轮廓线l ,如图,以旋转轴为Ox 轴,光源放在原点O(0,0).设l 的方程为y=y(x,y).由O 点发出的光线经镜面反射后平行于Ox 轴.设M(x,y)为l 上任一点,光线OM 经反射后为MR.MT 为l 在M 点的切线,MN 为l 在M 点的法线,根据光线的反射定律,有∠OMN=∠NMR从而tan ∠OMN=tan ∠NMR因为MT 的斜率为'y ,MN 的斜率为-'1y ,所以由正切公式,有 tan ∠OMN='1'1xy y x y y ---, tan ∠NMR='1y从而'1y =-yxy yy x -+''即得到微分方程2'yy +2x 'y -y=0由这方程中解出'y ,得到齐次方程'y =-1)(2+±yxy x 令xy=u,即y=xu,有 dx dy =u+dxdu x 代入上式得到dx du x =uu u 221)1(+±+- 分离变量后得=+±+221)1(u u udu xdx-令1+22t u =上式变为xdxt dt -=±1.积分后得 ln xC t ln1=+ 或112±=+xcu .两端平方得2211⎪⎭⎫⎝⎛+=+x c u化简后得x c xc u 2222+=以222c cx y xyu +==代入，得.这是一族以原点为焦点的抛物线. 2．动力学问题动力学是微分方程最早期的源泉之一.我们都知道动力学的基本定律是牛顿第二定律ma f =这也是用微分方程来解决动力学的基本关系式.它的右端明显地含有加速度a,a 是位移对时间的二阶导数.列出微分方程的关键就在于找到外力f 和位移对时间的导数－速度的关系.只要找到这个关系,就可以由ma f =列出微分方程了. 在求解动力学问题时,要特别注意力学问题中的定解条件,如初值条件等. 例：物体由高空下落,除受重力作用外,还受到空气阻力的作用,在速度不太大的情况下,空气阻力可看做与速度的平方成正比试证明在这种情况下,落体存在极限速度1v .解 设物体质量为m,空气阻力系数为k,又设在时刻t 物体下落的速度为v,于是在时刻t 物体所受的合外力为2kv mg f -=（重力-空气阻力）从而,根据牛顿第二定律可得出微分方程2kv mg dtdv m -= 因为是自由落体,所以有()00=v⎰⎰=-t vdt kvmg mdv 002 积分得t kvmg kv mg mg m=-+ln 21 或mkg tkvm g kv m g 2ln=-+ 解出v,得⎪⎪⎭⎫⎝⎛+⎪⎪⎭⎫⎝⎛-=1122m kg t m kg t e k e m g v 当∞→t 时,有1lim v kmgv t ==+∞→ 据测定,s k αρ=,其中 为与物体形状有关的常数,为介质密度,s 为物体在地面上的投影面积.人们正是根据公式1lim v kmgv t ==+∞→ ,来为跳伞者设计保证安全的降落伞的直径大小的.在落地速度1v ,m, α,与 一定时,可定出s 来. 3.流体混合问题中学代数中有这样一类问题：某容器中装有浓度为1c 的含某种物质A 的液体V 升,从其中取出1v 升后,加入浓度为2c 的液体2v 升,要求混合后的液体的浓度以及物质A 的含量.这类问题用初等代数就可以解决了. 但是,容器内装有含物质A 的流体.设时刻t=0时,流体的体积为0V ,物质A 的质量为0x .今以速度2v （单位时间的流量）放出流体,而同时又以速度1v 注入浓度为1c 的流体,试求时刻t 时容器中物质A 的质量及流体的浓度. 这类问题称为流体混合问题.它是不能用初等数学解决的,必须用微分方程来计算.首先,我们用微元发来列方程.设在时刻t,容器内物质A 的质量为x=x(t),浓度为2c ,经过时间dt 后,容器内物质A 的质量增加了dx.于是,有关系式()dt v c v c dt v c dt v c dx 22112211-=-=因为()tv v v xc 2102-+=代入上式有()dt t v v v xv v c dx ⎥⎦⎤⎢⎣⎡-+-=210211或()tv v v xv v c dt dx210211-+-= 这是一个线性方程.求物质A 在时刻t 的质量的问题就归结为求方程满足初始条件x(0)= 0x 的解的问题.例: 某厂房容积为45m ×15m ×6m,经测定,空气中含有0.2﹪的2CO .开通通风设备,以360s m 3的速度输入含有0.05﹪的2CO 的新鲜空气,同时又排出同等数量的室内空气.问30min 后室内所含2CO 的百分比.解 设在时刻t,车间内2CO 的百分比为x(t) ﹪,当时间经过dt 后,室内2CO 的该变量为45×15×6×dx ﹪=360×0.05﹪×dt-360×x ﹪×dt于是有关系式4050dx=360(0.05-x)dt或()dt x dx -=05.0454初值条件为x(0)=0.2.将方程分离变量并积分,初值解满足dt x dxt x ⎰⎰=-02.045405.0求出x,有X=0.05+0.15t e454-以t=30min=1800s 代入,得x ≈0.05.即开动通风设备30min 后,室内的2CO 含量接近0.05﹪,基本上已是新鲜空气了. 4.变化率问题若某未知函数的变化率的表达式为已知,那么据此列出的方程常常是一阶微分方程.例：在某一个人群中推广技术是通过其中已掌握新技术的人进行的,设该人群的总人数为N,在t=0时刻已掌握新技术的人数为0x ,在任意时刻t 已掌握新技术的人数为x(t)(将x(t)视为连续可微变量),其变化率与已掌握新技术人数和未掌握新技术人数之积成正比,比例系数k ＞0,求x(t).解 由题意立即有()()00,x x x N kx dt dx=-= 按分离变量法解之,()kdt x N x dx=-,即kNdt dx x N x =⎪⎭⎫ ⎝⎛-+11 积分并化简的通解kNtkNtceNce x +=1 由初值条件得特解kNtkNtex x N e Nx x 000+-= 总结：通过以上几个简单的例子，我们发现用微分方程解决一些实际问题其实很方便，也很普遍，所以在以后的学习中，除了学习必须的理论与方法外，更应该加强理论与实际的联系，将学习的知识更好的用于解决实际问题中.参考文献[1] 欧阳瑞,孙要伟.常微分方程在数学建模中的应用[J].宿州教育学院学报,2008,11(2)[2] 朱思铭,王高雄.常微分方程[M],第三版.高等教育出版社.2006.[3] 东北师范大学微分方程教研室.常微分方程,第二版.高等教育出版社.2005 [4] 蔡燧林,常微分方程,第二版.武汉大学出版社.2003.。