《微分方程-微分方程应用模型举例.pdf》由会员分享，可在线阅读，更多相关《微分方程-微分方程应用模型举例.pdf（63页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第九节微分方程应用模型举例第十二章第十二章一、主要内容一、主要内容二、典型例题二、典型例题三、同步练习三、同步练习四、同步练习解答四、同步练习解答1 建模建模(常微分方程模型常微分方程模型)1.基本步骤基本步骤2 求解求解(精确解或近似解精确解或近似解)或对解或对解作定性分析作定性分析(研究解的性态研究解的性态)3 解的实际意义解的实际意义(解释与预测解释与预测)注注 数学模型数学模型:对实际问题的对实际问题的简化简化而而本质本质的数学描述的数学描述.一、主要内容一、主要内容(一一)利用微分方程解决实际问题利用微分方程解决实际问题的基本步骤与方法的基本步骤与方法2.基本方法基本方法建立常微分方
2、程模型的建立常微分方程模型的主要方法主要方法:(1)根据规律列方程根据规律列方程;(2)微元分析法微元分析法;(3)模拟近似法模拟近似法.(物理物理,力学等力学等)(生物生物,经济经济,医学等医学等)(二二)几何应用几何应用(1)根据几何关系列方程根据几何关系列方程;(2)确定定解条件确定定解条件;(3)求通解求通解,并根据定解条件确定特解并根据定解条件确定特解.1.解微分方程解微分方程几何应用题几何应用题的方法和步骤：的方法和步骤：2.常见的几何公式常见的几何公式(1)平面曲线上一点的切线斜率：平面曲线上一点的切线斜率：xyydd=(2)平面曲线弧段平面曲线弧段 L的弧长：的弧长：,),(:
3、,d12baxxfyLxysba=+=+=,),(:,d22 =+=+=rrLrrs(3)平面曲线的曲率：平面曲线的曲率：(4)平面图形的面积：平面图形的面积：xyO ab)(xfy=baxxfAd)(23)1(2yyK+=(3)立体的体积立体的体积=DyxyxfVdd),(),(yxfz=xyO zD=baxxxfVd)(2 xyO ab)(xfy=vVd),(绕绕x轴旋转一周而成的旋转体体积：轴旋转一周而成的旋转体体积：曲顶柱体的体积：曲顶柱体的体积：一般立体的体积：一般立体的体积：yxfz=xyO z(三三)物理应用物理应用1.解微分方程解微分方程物理应用题物理应用题的方法和步骤：的方法
4、和步骤：(1)根据物理定律及实验规律列方程根据物理定律及实验规律列方程;(2)确定初始条件确定初始条件;(3)求通解求通解,并根据初始条件确定特解并根据初始条件确定特解.2.常用的物理定律常用的物理定律(1)牛顿运动定律；牛顿运动定律；(2)虎克定律；虎克定律；(3)万有引力定律；万有引力定律；(4)基尔霍夫电流及电压定律基尔霍夫电流及电压定律.(三三)其他应用其他应用;1)(00)0()0()()1(=+=+=xexffxxfy和满足条件函数假设：和满足条件函数假设：二、典型例题二、典型例题例例1；和分别相交与点和；和分别相交与点和与曲线与曲线轴的动直线平行于 轴的动直线平行于211)()2
5、(PPeyxfyMNyx=.)(),()3(21的表达式长度，求函数的恒等于轴所围封闭图形的面积与直线曲线的表达式长度，求函数的恒等于轴所围封闭图形的面积与直线曲线xfyPPxMNxfy=)(xfy=MN2P 1P 1=xeyxyO=xxxfS0d)(点的坐标分别为点的坐标分别为21,PP),(,(1xfxP，)1,(2 xexP封闭图形的面积：由题设条件的方程为设动直线封闭图形的面积：由题设条件的方程为设动直线)3(,xxMN解解=)1,(2 xexPMN1=xeyxyO)(,(1xfxP)(xfy=故由题设，有故由题设，有)(1d)(0 xfexxfxx=两两端求导，得端求导，得)()(x
6、fexfx=)(1d)(0 xfexxfxx=即即xexfxf=+=+)()(它的通解为它的通解为()Cxeeexfxxx+=+=d)(dd()Cxeeexxx+=+=dxxCee+=+=21一阶非齐次线性方程一阶非齐次线性方程,0)0(=f由由()xxeexf=21)(因此所求函数为因此所求函数为,21=C得得xxCeexf+=+=21)(例例2)轴与直线轴与直线曲曲线上连续，若由，在设函数线上连续，若由，在设函数xttxxxfyxf)1(,1),(1)(=+)1()(3)(2ftfttV=.92)(2的解件并求该微分方程满足条所满足的微分方程，试求=的解件并求该微分方程满足条所满足的微分方
7、程，试求=xyxfy旋转体体积为轴旋转一周所成的所围成的平面图形绕旋转体体积为轴旋转一周所成的所围成的平面图形绕xxyO)(xfy=1t解解公式有依题意及旋转体的体积 公式有依题意及旋转体的体积=txxftV12d)()(即即)1()(d)(3212ftftxxft=求导，得求导，得两边对两边对t)()(2)(322tftttftf+=+=将上式改写为将上式改写为xyyyx2322=)1()(32ftft=,xyu=令=令即即)1(3dd=uuxux分离变量分离变量xxuuud3)1(d=xyxyxy23dd2=为齐次方程为齐次方程，得，得两边同乘以两边同乘以21xuuxuxu23dd2=+则
8、有=+则有代入上式代入上式，将条件将条件922=xy，得得1=C从而所求解为从而所求解为1,3=xyxxy.1,13+=xxxy或或代回有将代回有将xyu=yCxxy3=即即31Cxuu=积分得积分得Cxuulnln3ln)1ln(+=例例3满足微分方程设函数满足微分方程设函数)(xfy=,xxeyyy223=+处的切线与曲线且其图形在点=+处的切线与曲线且其图形在点)1,0(在该点的切线重合，在该点的切线重合，12+=+=xxy).(xf求函数求函数分析分析件件，本题变相给出了初始条本题变相给出了初始条由于图形由于图形1)0(,1)0(=yy可得出初始条件可得出初始条件,12有公共点+=有公
9、共点+=xxy及在该点与过点及在该点与过点)1,0(解解特征方程为特征方程为0232=+rr特征方程为特征方程为0232=+rr特征根为特征根为2,121=rr解对应齐次线性方程的通 解对应齐次线性方程的通xxeCeCY221+=+=故可设为特征根故可设为特征根由于由于,1=xAxey=代入原方程得求代入原方程得求 yy,2=A故原方程的通解为故原方程的通解为xxxxeeCeCy2221+=+=,1)1,0(2有公共切线处与有公共切线处与由由于曲线在点+=于曲线在点+=xxy由由此得此得121=+CC12221=+CC解得解得0,121=CC故所求曲线为故所求曲线为.)21(xexy=.1)0
10、(,1)0(=yy故有故有xxxxeeCeCy2221+=原方程通解：+=原方程通解：设木板的的速度穿出以设木板的的速度穿出以的木板后的木板后垂直打入厚子弹一初速度垂直打入厚子弹一初速度./80,10/20010smvcmsmv=例例4.,所需的时间求子弹穿过木板成正比所需的时间求子弹穿过木板成正比解解阻力与子弹的速度平方阻力与子弹的速度平方分析分析,律的本题是讨论物体运动规在物体律的本题是讨论物体运动规在物体注意注意,即木板的阻力即木板的阻力由牛顿第二定律得由牛顿第二定律得222ddkvFtsm=,个力运动方向上子弹只受一 个力运动方向上子弹只受一)(2为特定的比例系数其中 为特定的比例系数
11、其中kkv.200,000=ttvs其初始条件为其初始条件为由导数的物理意义知由导数的物理意义知22ddddtstv=积分得积分得2ddkvtvm=分离变量分离变量tmkvvdd2=代入方程得代入方程得11Ctmkv=即即11Ctmkv+=+=,2000=tv代入初始条件代入初始条件20011=C得对得对两边积分得两边积分得，得代入，得代入00=ts200ln2kmC=2)2001ln(Ctmkkms+=+=)2(20011dd+=+=tmkvts)1(所需的时间所需的时间，为了求出子弹穿过木板为了求出子弹穿过木板，时刻为时刻为不妨设子弹穿出木板的不妨设子弹穿出木板的0t,/8010smvvt
12、t=时，时，由已知条件由已知条件，当，当有故由式有故由式)1(20011800+=+=tmk)3(,10110,mcmtt=路径为时子弹在木板中走过的当路径为时子弹在木板中走过的当另一方面另一方面式有从而由式有从而由)2(200ln2001ln1010kmtmkkm+=+=200ln801lnkmkm+=+=80200lnkm=25lnkm=，所以所以25ln10=mk得代入式得代入式)3(25ln10140030=t)(12001s ox钉子需多少时间钉子需多少时间米，试问整个链条滑过米，试问整个链条滑过下下垂米，另一边的一边下垂上，运动开始时，链条一无摩擦的钉子一质量均匀的链条挂在垂米，另
13、一边的一边下垂上，运动开始时，链条一无摩擦的钉子一质量均匀的链条挂在108解解,米链条下滑了经过时间设链条的线密度为米链条下滑了经过时间设链条的线密度为xt 则由牛顿第二定律得则由牛顿第二定律得,)8()10(dd1822gxgxtx+=+=.0)0(,0)0(,99=xxgxgx即即例例5m8m10)(8mx)(10mx+xgxf)10(1+=gxf)8(2=.0)0(,0)0(,99=xxgxgx特征方程特征方程092=gr特征根特征根32,1gr=通解为对应的齐次线性方程的 通解为对应的齐次线性方程的tgtgeCeCtx312311)(+=+=由观察可知：非齐次线性方程有特解由观察可知：
14、非齐次线性方程有特解 x=-1.非齐次线性方程的通解为非齐次线性方程的通解为,8,=x即整个链条滑过钉子即整个链条滑过钉子代入上式得代入上式得)()809ln(3秒+=秒+=gt1)(312311+=+=t gt geCeCtx，得，得由由0)0(,0)0(=xx.2121=CC,1)(21)(3131+=+=tgt geetx例例6解解体积的四分之一,问其余部分在多长时内融化完?体积的四分之一,问其余部分在多长时内融化完?假定一个雪球是半径为假定一个雪球是半径为r 的球,其融化时体积的变化率与雪球的表的球,其融化时体积的变化率与雪球的表雪球融化问题雪球融化问题面积成正比面积成正比,比例常数为
15、比例常数为 k 0 (k 与空气温度等有关与空气温度等有关),已知两小时内融化了其已知两小时内融化了其由于雪球体积的变化率正比与其表面积由于雪球体积的变化率正比与其表面积24ddrktV=代入上式，得将代入上式，得将234rV=24ddrktV=)0dd,(xxy上任意一点过曲线上任意一点过曲线)(,1)0(,0)(xyyyxy=上述两直线轴的垂线作曲线的切线及上述两直线轴的垂线作曲线的切线及,),(xyxP 上以区间为轴所围的三角形面积记上以区间为轴所围的三角形面积记与与xSx,0,1并设记为为曲边的曲边梯形面积并设记为为曲边的曲边梯形面积,2)(Sxyy=.)(,1221的方程求此曲线恒为
16、的方程求此曲线恒为xyySS=1)(xyy=S1xxyO),(yxPS2)1,0(,11)(2且在曲线上点且在曲线上点意一点处的曲率为意一点处的曲率为，其上任是一向上凸的连续曲线设，其上任是一向上凸的连续曲线设yxyy+=+=3.求该曲线的方程，处的切线方程为求该曲线的方程，处的切线方程为,1+=xy.)(的极值并求函数的极值并求函数xyy=目标的跟踪目标的跟踪 设位于坐标原点的设位于坐标原点的甲甲舰向位于舰向位于x 轴上点轴上点A(1,0)处的处的乙乙舰发射制导导弹，导弹舰发射制导导弹，导弹4.,5)0vy速度为轴的直线航行，导弹的沿平行于数速度为轴的直线航行，导弹的沿平行于数.求导弹运行的
17、曲线方程又问求导弹运行的曲线方程又问乙乙舰行驶多远时舰行驶多远时,头始终指向头始终指向乙乙舰.如果舰.如果乙乙舰以最大速度舰以最大速度v0(v0是常是常它将被导弹击中它将被导弹击中?设弹簧上端固定设弹簧上端固定,有两个相同的重物有两个相同的重物(质量为质量为m)挂于弹簧下端挂于弹簧下端,使弹簧伸长使弹簧伸长2a,今突然取走了一个重物，使弹簧由静止开始振动，求所挂重物的运动规律今突然取走了一个重物，使弹簧由静止开始振动，求所挂重物的运动规律.5.6.放射性元素铀由于不断地有原子放射出微粒子而变成其他元素，铀的含量就不断减少放射性元素铀由于不断地有原子放射出微粒子而变成其他元素，铀的含量就不断减少
18、,这种现象叫做衰变这种现象叫做衰变.由原子物理学知道，见到的衰变速度与当时未衰变的原子的含量由原子物理学知道，见到的衰变速度与当时未衰变的原子的含量M成正比成正比.已知已知t=0 时铀的含量为时铀的含量为M0，求在衰变过程中铀含量，求在衰变过程中铀含量M(t)随时间随时间t 变化的规律变化的规律.盐溶液的浓度盐溶液的浓度一容器内盛有一容器内盛有100L盐水盐水,其中含盐其中含盐10kg.今用每分钟今用每分钟2L的速度把净水注入容器的速度把净水注入容器(假定净水与盐水立即调和假定净水与盐水立即调和),又以同样速度使盐水流出又以同样速度使盐水流出.试求容器内盐量随时间变化的规律试求容器内盐量随时间
19、变化的规律.7.四、同步练习解答四、同步练习解答1.,)0,2(,),(),(0上一定点为上任意一点上一定点为上任意一点为为的极坐标方程为设曲线的极坐标方程为设曲线LMLrMrrL =所围成的曲边扇形与曲线若极径所围成的曲边扇形与曲线若极径LOMOM,0解解面面积可用定积分所围成的曲边扇形与积可用定积分所围成的曲边扇形与由由极径极径LOMOM10,rrd2102,表示表示两点间弧长上两点间弧长上而而曲线曲线MML,0两点间弧长值的一半，上两点间弧长值的一半，上面面积值等于积值等于MML,0.的方程求曲线的方程求曲线 L故有故有rrrd21d2102202+=+=得求导 得求导两两边对边对,22
20、2rrr+=+=12=rrr即即分离变量分离变量rrrd1d2=rrd022+=+=MMss0d 两边积分两边积分=22111d1drrrrrr1arcsin=C+=得代入上式将条件得代入上式将条件,2)0(=r,6C=的方程为故所求曲线的方程为故所求曲线Lr1)6sin(=m即即).6csc(m=r2.二阶可导，且设函数二阶可导，且设函数)0()(xxy上任意一点过曲线上任意一点过曲线)(,1)0(,0)(xyyyxy=上述两直线轴的垂线作曲线的切线及上述两直线轴的垂线作曲线的切线及,),(xyxP 上以区间为轴所围的三角形面积记上以区间为轴所围的三角形面积记与与xSx,0,1并设记为为曲边
21、的曲边梯形面积并设记为为曲边的曲边梯形面积,2)(Sxyy=.)(,1221的方程求此曲线恒为的方程求此曲线恒为xyySS=1)(xyy=S1xxyO),(yxPS2解解处的切线方程为在点处的切线方程为在点曲线曲线),()(yxPxyy=)(xXxyyY =),0,)(xyyxx轴的交点为轴的交点为它与它与1)(xyy=S1xxyO),(yxPS2,1)0(,0)(=yxy由于由于时时，即当即当0 x为单调增加函数，为单调增加函数，)(xy于是从而于是从而,0)(xy)(211xyyxxyS=yy=22yyS=221又曲边梯形的面积又曲边梯形的面积=xttyS02d)(知知由条件由条件1221
22、=SS1d)(02=xttyyy，得，得由由1)0(=y.1)0(=y1d)(02=xttyyy上式两边对上式两边对x求导，得求导，得0)()(2222=yyyyyy整理，得整理，得2)(yyy=型型属于属于),(yyfy=，则令，则令yppyypydd)(=代入上述方程代入上述方程2ddpypyp=pypy=dd)0(=pyQpypy=dd分离变量、积分得分离变量、积分得1lnlnlnCyp+=yCp1=即即=y分离变量、积分得分离变量、积分得，得，得由由1)0(,1)0(=yy.11=Cyy=于是于是2lnCxy+=得得由由,1)0(=y.02=C.)(xexy=且在曲线为其上任意一点处的
23、曲率，是一向上凸的连续曲线设且在曲线为其上任意一点处的曲率，是一向上凸的连续曲线设,11)(2yxyy+=3.求该曲线的处的切线方程为上点求该曲线的处的切线方程为上点,1)1,0(+=xy解解,12yyK+=有曲率的定义知曲率由体设得=有曲率的定义知曲率由体设得因曲线是上凸的因曲线是上凸的，所以，所以,0 y23211)1(yyy+=+=+.)(的极值方程，并求函数的极值方程，并求函数xyy=,112=+=+yy则则,pypy=令=令)1(2pp+=+=分离变量得分离变量得xppd)1(d2=+=+积分得积分得xCp=1arctan,1)1,0()(+=xyxyy处的切线方成为在处的切线方成为
24、在因为因为100=xxpy所以所以21yy=即即从而上方程化为从而上方程化为代入上式得代入上式得41C=故故xp=4arctan 所以所以)4tan(xpy=积分得积分得2)4cos(lnCxy+=+=),1,0(因因为曲线过点代入上式得为曲线过点代入上式得2ln2112+=+=C,10=xy所以所以故所求曲线方程为故所求曲线方程为,)(的极值的极值下下面求面求xyy=,1)4cos(x 因因为为,4时所以当时所以当x=,1)4cos(,4=xx 时且当时且当2ln211+=+=y函函数取得极大值数取得极大值回回顾求解过程，顾求解过程，,01)4(=，这种现象叫做衰变，这种现象叫做衰变.由原子
25、物理学知道，见到的衰变速度与当时未衰变的原子的含量由原子物理学知道，见到的衰变速度与当时未衰变的原子的含量M成正比成正比.已知已知t=0时铀的含量为时铀的含量为M0，求在衰变过程中铀含量，求在衰变过程中铀含量M(t)随时间随时间t变化的规律变化的规律.MMt00|=分离变量后得分离变量后得.ddtMM=.dd=tMM M0,得得,lnlnCtM+=以初始条件代入上式，得以初始条件代入上式，得.e tCM=,e00CCM=.e0 tMM=盐溶液的浓度盐溶液的浓度一容器内盛有一容器内盛有100L盐水盐水,其中含盐其中含盐10kg.今用每分钟今用每分钟2L的速度把净水注入容器的速度把净水注入容器(假
26、定净水与盐水立即调和假定净水与盐水立即调和),又以同样速度又以同样速度7.使盐水流出使盐水流出.试求容器内盐量随时间变化的规律试求容器内盐量随时间变化的规律.解解).0d(d,d,+QQQQttt降至盐量由溶液内含内在微小时间间隔设降至盐量由溶液内含内在微小时间间隔设设在设在t t时刻溶液内的含盐量时刻溶液内的含盐量Q（t）为，现）为，现利用微分元素法建立未知函数利用微分元素法建立未知函数Q（t）满足）满足的微分方程的微分方程),:(d2Lt 单位器内流出的溶液量为在这一段时间内，从容单位器内流出的溶液量为在这一段时间内，从容):(100)(LkgtQt单位时刻浓度盐水浓度视为单位时刻浓度盐水浓度视为)61.12(0d2100d=tQQ应满足的微分方程此即未知函数应满足的微分方程此即未知函数)(tQQ=含盐量的改变量含盐量的改变量因此因此,的的通解通解为为方程方程)61.12(50tCeQ=10|0=tQ由由变化规律为故容器内含盐量随时间 变化规律为故容器内含盐量随时间5010)(tetQ=,10kgt 时溶液内含盐量为时溶液内含盐量为由于由于10|0=tQ初始条件为那么 初始条件为那么,10=C得得