常微分方程的实际应用

1、常微分方程的实际应用于萍摘要：常微分方程在当代数学中是极为重要的一个分支，它的实用价值很高，应用也很广泛,本文主要介绍常微分方程在几何、机械运动、电磁振荡方面的应用，并举例说明，体会常微分方程对解决实际问题的作用，在解决实际问题过程中通常是建立起实际问题的数学模型，也就是建立反映这个实际问题的微分方程，求解这个微分方程，用所得的数学结果解释实际问题，从而预测到某些物理过程的特定性质，以便达到能动地改造世界，解决实际问题的目的。关键字：常微分方程，几何，机械运动，电磁振荡，应用Abstract:Nomaldifferentialequationisanimportantpartofmathati

2、thasahighpracticalvalue.Thisthesisshowstheuseingeometry,mechaicsandelectrothermalandmakessomeexamples.Also,itsummarizesthenormalmoveofdealingwithpracticalproblemsbythenormaldifferentialequation.Normal,wesetupthemathsmaticmodeloftheproblem,solutethenormaldifferenticalequationmaketheuseoftheresulttoex

3、plainpracticalproblemsandmakeaforecastofsomespecialcharacterofphysicalprocess.Key:Normaldifferetialequationgeometrymechanicselectrothermaluse引言数学分析中所研究的函数，是反映客观现实世界运动过程中量与量之间的一种关系，但在大量的实际问题中遇到稍为复杂的一些运动过程时，反映运动规律的量与量之间的关系（即函数）往往不能直接写出来，却比较容易地建立这些变量和它们的导数（或微分）间的关系式，不同的物理现象可以具有相同的数学模型，这一事实正是现代许多应用数学工作者

4、和工程人员应用模拟方法解决物理或工程问题的理论依据。例如，利用电路来模拟某些力学系统或机械等等在现时已相当普遍。在自然科学和技术科学的其他领域中，例如化学、生物学、自动控制、电力技术等等，都提出了大量的微分方程问题，因此，社会的生产实践是常微分方程理论取之不尽的基本源泉。止匕外，常微分方程与数学的其他分支的关系也是非常密切的。它们往往互相联系、互相促进。例如，几何学、机械运动、电磁振荡就是常微分方程理论的丰富的源泉之一，常微分方程也是解决实际问题不可或缺的武器。一、常微分方程在几何学的应用在几何应用问题中，列的方程常常是含有变限定积分的方程。在求解时要化为相应的微分方程或微分方程初值问题。凡是

5、能用定积分计算的量，一定分布在某个区间(比如b,b】)上，并且对于该区间具有可加性，曲边梯形的面积A与区间a,b】有关，当把b,b1分成n个部分区间时，则所求量A也相应地分成n个部分量AA(i=12，n),而A就等于所有这些部n分之和，即a=£&A，这时我们就称面积A对区间la,b具有可加性，几i1何中的面积、弧长，曲线方程等都具有这种特性。在求解微分方程的应用问题时，列出方程是关键性的一步，一定要逐字逐句地仔细阅读题目，根据题目的要求确定未知函数和自变量，然后利用题设中指出的(或包含的)相等关系列出方程，应用问题常常是初值问题。因而，要从题设中确定未知函数满足的初始条件。常

6、微分方程在解决几何问题的过程中通常采用数形结合，达到简易直观的效果。利用y'表示曲线y=f(x)上(x,y)点处的切线斜率或-2表示曲线dyxy=f(x)上(x,y)点的法线斜率以及ff(t)dt表示由曲线ay=f(x)(f(x)上0),直线x=x,x=a,x轴所围图形的面积等方面的意义，列方程。解方程，在求解过程中一定要对常微分方程的解法熟悉于心，才能得心应手。首先要审视方程，判断方程类型，属于一阶微分方程还是可降阶微分方程或高阶微分方程等等。根据不同类型，确定解题方案。F面就让我们结合具体例题来体会常微分方程在解决几何问题的应例1，、设y=f(x)是第一象限内连接点A(0,1),B

7、(1,0)的一段连续曲线，M(x,y)为该曲线上任意一点，点C为M在x轴上的投影。O为坐标3原点，若梯形OCMA的面积与曲边三角形CBM的面积之和为X+1,求63f(x)的表达式。解：根据题意有:f(0)=1,f=0且x1+f(x)+f12x将上式两边对x求导数,xf(x)-f(t)dt=当0<xE1时，可化为一阶线性微分方程:一1一1f(x)-f(x)=x-xx方程两边同除x,即得盘=1-2xx积分可得必=x1cxx于是，方程通解为f(x)=x21cx把f(1)=0代入通解，可确定常数c=-2故所求函数f(x)的表达式为：-2_2一f(x)=x1-2x=(x-1),0.x.1例2、在上

8、半平面求一条向上凹的曲线，其任一点p(x,y)处的曲率等于此曲线在该点的法线段PQ长度的倒数，(Q是法线与x轴的交点),且曲线在点(1,1)处切线与x轴平行。解：见图，所求曲线为y=f(x),于是其在p(x,y)点处的曲率为：(1y2)2(1y2)2(.曲线为凹的，.二y">0)曲线y=f(x)在p(x,y)点处的法线方程：Y-y=-1(X-x)(y-二0)y它与x轴的交点Q的坐标Q(x+yy：0),1于是PQ=4(yy)2+y2=y(1+y'2)2,由题设k=1PQ即=11(1y2)2y(1y2)2=yy'=1+y'2这是不显含x的方程初始条件为，y|

9、x1,y'|xT0令y'=p,y"=pdp,于是方程变为dydp2pdyyp1pdp二dy1py12x.二-ln(1+p)=lny+c1,代入y'|x=。，得c1=022-：2.pp=y-1=p=±qy-1,积分得ln(y-y2-1)=(x-1)c2代入y|xg=1,得C2=0故所求曲线为：2/(xJ)1xx)yqy-1=e，即y=2(ee)例3、已知曲线过(1,1)点，如果把曲线上任一点P处的切线与y轴的交点记作Q,则以PQ为直径所做的圆都经过点F(1,0),求此曲线方程。解：见图所求曲线设为y=f(x)于是切线方程为Y-y=y(X-x)切线PQ与

10、y轴的交点Q的坐标为Q(0,y-xy)设M点为切线段PQ的中点，坐标为'xxy'y122)圆经过点F(1,0)MQ=MF,=1于是得方程yy-xy.y|xj1令y2=z,则方程1/2、1=2(y)、y2.12o21_:z=_z-'2'_dz2,dx=Inz=2lnxIncz=cx2(2)令z=c(x)x2为的解,代入并整理，得c(x)x2=-22c(x)x=c(x)=-x故的通解为:cx2=2x-1+dx即方程的通解为22y=2x-1cx,代入初值y|x/=1,得=0x故所求曲线为y2=2x-1例4、在制造探照灯的反射镜面时，总是要求将点光源射出的光线平行地反射

11、出去，以保证探照灯有良好的方向性，试求反射镜面的几何形状。解：取光源所在处为坐标原点，而x轴平行于光的反射方向，(见图)。设所求曲面由曲线V=f(x),z=0绕x轴旋转而成，则求反射镜面问题归结为求xy平面上的曲线y=f(x)的问题。过曲线y=f(x)上任一点M(x,y)作切线NT则由反射定律：入射角等于反射角，容易推知=“2从而OM=on注意到dy二tg：2=MPdxNP及5P=x,MP=y,OM=Jx2y2就得到函数y=f(x)所满足的微分方程式吆=一这是齐次方程。dxxxy设二y,将它化为变量分离方程求解x得y2=c(c+2x)c为任意常数故反射镜面的形状为旋转抛物面y2+z2=c(c+

12、2x)二、常微分方程在机械振动中的应用常微分方程与物理联系甚为广泛，下面我们就一起来看一下常微分方程在机械振动中的应用，常微分方程解决力学问题需要：建立坐标系，对所研究物体进行受力分析；根据牛顿第二定律F=ma,列方程；解方程。下面，让我们从实例中体会常微分方程在力学中的作用。例1：一个质量为m的船以速度v0行驶，在t=0时，动力关闭，假设水的阻力正比于vn,其中n为一常数，v为瞬时速度，求速度与滑行距离的函数关系。解：船所受的净力=向前推力-水的阻力=0-kvn,加速度=速度对时间的导数，即a=dv,dt于是，由题设有dvnm=-kvdt、V|t4=v0现在要求的不是速度与时间的关系,而是速

13、度与距离的关系，设距离为X,于是，上述方程可化为:dvdvdxdvmmmv-dtdxdtdx-kvn=mv1"dv=-kdx格)当n#2时，两边积分，得2-nmv2-n-kxc把v|y=V0,x|t白=0代入上式，得2-nmV。2-n故v2'k(2-n)2Hxv0m当n=2时，（户mv/dv=-kdx,kx积分得v=cem,将初值代入，得c=v。kx故v-v0em例2、两个质量相同的重物挂于弹簧下端，其中一个坠落，求另一个重物的运动规律，已知弹簧挂一个重物伸长为a。解：如图所示，建立坐标系设弹簧自由状态时长度为l,取l十a处(即挂一重物时弹簧的长度)为坐标原点，取x轴铅直向下

14、，设在t时刻，重物在x处,由虎克定律知，此时弹性恢复力为-kx,k为弹性系数，负号“一”是因为弹性恢复力与位移反向，由牛顿第二定律有:212m=-ktdtx(0)=ax(0)=0挂两重物时，弹簧伸长2a,由虎克定律有:mg2mg=k2a=k=a2方程=-9x,dt2a于是方程通解为x=gcos.gtc2sin.gtaa把初始条件x(0)=a,x'(0)=0代入以上两式其特征方程：，2=_g=.gaa得c1=a,c2=0所求重物的运动规律为x=acos.gt数学摆是系于一根长度为l的线上而质量为m的质点M在重力作用下，它在垂直于地面的平面上沿圆周动运。如图所示，试确定摆的运动方程。点M向

15、着平衡位置A的方向运动，即当角中为正时，向减小中的方向运动,当角邛为负时，向增大中的方向运动，所以MP的数值等于-mgsin%因此，摆的运动方程是mdv=-mgsin%即曹=.9而中。dtdt2l(1)如果只研究摆的微小振动，即当中比较小时的情况，我们可以取sin平的近似值中代入上式，这样就得到微小振动时摆的运动方程:d2;dt2(2)如果我们假设摆是在一个粘性的介质中摆动，那么，沿摆的运动方向就存在一个与速度v成比例的阻力，若阻力系数为N,则摆动方程为d2dt2mdt(3)如果沿摆的运动方向恒有一个外力F(t)作用于它，这时摆的运动称为强迫微小振动，其方程为：d4+匕当+9*lF(t)。dt

16、2mdtlml当要确定摆的某一个特定运动时，我们应给出摆的初始状态：当t=0d;日寸，甲=%=Wo。dt这里中0代表摆的初始位置，Wo代表初始角速度。例4：生产实践中很多机械问题都归结为弹性振动问题，下面便是一个弹簧振动的典型例子。设有弹性系数c而自然长度为l的弹簧竖着悬挂着。它的上端固定，下端悬挂，一个质量为m的物体，物体受到垂直干扰力f=f4)，求物体的运动规律所满足的微分方程。解：如图所示，取通过悬挂点的直线为x轴,向下记为正方向，原点取在系统平衡位置，为确定物体运动规律，先分析它的位置，x=x处的受力情况。(1)弹簧弹性力f0,依虎克定律f0=-c(6+x),其中6为弹簧在物体重力作用

17、下的伸长量。(2)物体所受重力p=mg(3)介质阻力R与物体运动速度成正比，与运动方向相反，dxR-v二一dt其中R为常数，称为阻尼系数。(4)重力干扰力f=fi(t)因此，这时物体所受合外力F=f0PRf=_cE，xmg-；Xflm/=-c、xmg-史dtdtfi(t)由于系统的平衡位置处，弹性力f0=-c6与重力p=mg平衡，故有再由牛二定律，得方程:-cmg=0于是上述方程写成2dx,dx1小m-cx=f1(t)dt2dt若记一=2n,=k2mmfi(t)mf(t)则可写成d2x八dx72nkx=f(t)dt2dt这就是该物体在外力f(t)作用下运动规律。x=x(t)所满足的微分方程若物

18、体振动过程中，未受外力干扰，即f(t)=0,则微分方程,2.dxcdx52nkx=0dt2dt三、常微分方程在电磁振荡中的应用建立起实际问题的数学模型一般是比较困难的，因为这需要对与问题有关的自然规律有一个清晰的了解，如前面所求的力学问题就要对牛二定律有清楚的认识，同时也需要有一定的数学知识，为了要建立起实际问题的数学模型，一定要学习有关的自然科学和工程技术的专业知识，微分方程往往可以看作是各种不同物理现象的数学模型，我们在建立微分方程的时候，只能考虑影响这个物理现象的一些主要因素，而把其它一些次要因素忽略掉，如果的确考虑到了那些最主要的因素，那么，我们所得到的微分方程，它的解和所考虑的物理现

19、象就是比较接近的，这时，我们得到的数学模型是有用的，否则，我们还应考虑其它一些因素，以便建立起更为合理的数学模型，为了解决热电学问题，需要了解其中的一些基本规律，如下面将用到牛顿冷却定律，其内容为热量总是从物体中温度高的向温度低的物体传导；在一定温度范围内，一个物体的温度变化速度与这一物体的温度和其所在介质温度差值成比例，等等，我们将在实例中一一解答。常微分方程解决电磁振荡问题通常建立起电热学问题的数学模型，也就是反映这个实际问题的微分方程。求解这个微分方程。用所得的数学结果解释实际问题，从而预测到某些物理过程的特定性质，以便达到能动地改造世界，解决实际问题的目的。接下来，就让我们从实例中体会

20、常微分方程在电热方面的应用。例1.R-L电路，如图，它包含电感L,电阻R和电源E,设t=0时，电路中没有电流，我们要求建立：当开关k闭合后，电流I应该满足的微分方程，假设R,L,E都是常数。k解：为了建立电路的微分方程，我三E们引用关于电路的基尔霍夫第二定律：ILL在闭合回路中，所有支路上的电压的代数和等于零。注意到经过电阻R的电压降是RI,而经过电感L的电压降是LdIdt由基尔霍夫第二定律得到E-L*R-。dlR,E即一I=dtLL求出的I=1应满足条件:当t=。时，I=。，如果假定在t=t。时，I=I。，电源E突然短路，因而E变为零，此后亦保持为零，那么电流I满足方程。dI+RI=。，及条

21、件t=t。时，I=I。dtL。，。例2R-L-C电路，如图所示，它包括电感L,电阻R和电容C,设RL,C均为常数，电源e(t)是时间t的已知函数，我们要求建立：当开关k闭合后，电流I应满足的微分方程。解：注意到经过电感L,电阻R和电容C的电压降分别为L且，RI和Q,dtC其中Q为电量，因此由基尔霍夫第二定律得到e(t)=L'RIQdtC1噜,微分上式得到2_d2IRdII1de(t)dt2LdtLCLdt这就是电流I应满足的微分方程，如果e(t)=常数，得到d2Idt2RdlI十LdtLC如果又有R=0,则得到d2Idt21二0LC例3".电容器的充电和放电，如图所示R-C电

22、路，开始时电容C上没有电荷，电容两端电压为零，我们把开关k闭合“1”后，电池E就对电容C充电，电容C两端电压Uc逐渐升高，经过相当时间后，电容充电完毕，我们再把开关k合上“2”，这时电容就开始放电过程，现在要求找出充、放电过程中，电容C两端的电压Uc随时间t的变化规律。解：对于充电过程，由闭合回路的基尔霍夫第二定律有UC+RI=E对电容C充电时，电容上的电量Q逐渐增多，根据Q=CUc得到：I=dQ=d(CUc)=Cdtdtdt将代入，得Uc满足的微分方程:dUcRCcUC=EdtIIC这里R,C,E都是常数，方程属于变量分离方程，将变量分离得到dUc二dtUc-ERCc两边积分，得到lnUc-E=cRCCi_ltLtRCRC=C2e这里C2=土eCl为任意常数。将初始条件：t=0时，Uc=0代入得至口C2=-E这就是R-C电路充电过程电容C两端的电压变化规律，由知道，电压Uc从零开始逐渐增大，且当tTy时，Uc=E,在电工学中，通常称RC为时间常数，当t=3i时，Uc=0.95E,就是说，经过3的时间后，电容C上的电压已达到外加电压的95%,实际上，通常认为这时电容C的充电过程已基本结束，易见充电结果Uc=E,对于放电过程，可以类似地进行。例41.将某物体放置于空气中，在时刻t=0