第一章一阶微分方程的应用优秀课件.ppt

《第一章一阶微分方程的应用优秀课件.ppt》由会员分享，可在线阅读，更多相关《第一章一阶微分方程的应用优秀课件.ppt（43页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第一章一阶微分方程的应用第1页，本讲稿共43页1：应用大意应用大意适应范围适应范围与变化率有关的各种实际问题应用三步曲应用三步曲(1)建立模型建立模型 Modelling(2)模型求解模型求解 Solving(3)模型应用模型应用 Application建议:模型要详略得当第2页，本讲稿共43页应用一：曲线族的等角轨线应用一：曲线族的等角轨线设给定一个平面上以设给定一个平面上以C C为参数的曲线族为参数的曲线族（*）我们设法求出另一个以我们设法求出另一个以k k为参数的曲线族为参数的曲线族（*）使得曲线族（使得曲线族（*）中的任一条曲线与曲线族）中的任一条曲线与曲线族样的曲线族（样的曲线族（*

2、）是已知曲线族（）是已知曲线族（*）的）的（*）中的每一条曲线相交时成定角）中的每一条曲线相交时成定角 则称这则称这等角轨线族。等角轨线族。第3页，本讲稿共43页当当 时，称曲线族（时，称曲线族（*）是（）是（*）的）的正交轨线族。正交轨线族。例如：曲线族例如：曲线族 是曲线族是曲线族 的正交轨线族。的正交轨线族。第4页，本讲稿共43页第5页，本讲稿共43页设设 y=y(x)为为(C)中任一条曲线，于是存在相应的中任一条曲线，于是存在相应的C，使得，使得 因为要求因为要求x,y,y 的关系，将上式对的关系，将上式对x求导数，得求导数，得 (1.84)这样，将上两式联立，即由这样，将上两式联立，

3、即由上述关系式成为曲线族满足的微分方程上述关系式成为曲线族满足的微分方程第6页，本讲稿共43页例例1 1 求抛物线族求抛物线族 的正交轨线族。的正交轨线族。解：解：对方程两边关于对方程两边关于x x求导得求导得由由 解出解出C C代入上式得曲线族代入上式得曲线族 在点在点 处切线斜率为处切线斜率为 第7页，本讲稿共43页由于所求曲线族的曲线与由于所求曲线族的曲线与 中的曲线在中的曲线在 正交，故满足方程正交，故满足方程 这是一个变量可分离方程求解得这是一个变量可分离方程求解得 的正交的正交曲线族为曲线族为这是一个椭圆，如右图这是一个椭圆，如右图放大此图 图图2.162.16第8页，本讲稿共43

4、页第9页，本讲稿共43页应用二：应用二：雨滴的下落雨滴的下落考虑雨滴在高空形成后下落的过程中速度的变化三种不同的假设(1)自由落体运动自由落体运动(2)小阻力的情况(3)大阻力的情况大阻力的情况第10页，本讲稿共43页(1)自由落体运动下落过程中没有任何阻力下落过程中没有任何阻力第11页，本讲稿共43页(2)小阻力的情况下落过程中阻力与速度和半径的乘积成比例下落过程中阻力与速度和半径的乘积成比例 第12页，本讲稿共43页(3)大阻力的情况下落过程中阻力与速度和半径的乘积平方成比例下落过程中阻力与速度和半径的乘积平方成比例第13页，本讲稿共43页三、药物设计医生给病人开处方是必须注意两点：医生给

5、病人开处方是必须注意两点：服药的剂量和服药的时间间隔。服药的剂量和服药的时间间隔。超剂量的药物会对患者产生严重不良后果，超剂量的药物会对患者产生严重不良后果，甚至死亡；剂量不足，则不能达到治疗的效果。甚至死亡；剂量不足，则不能达到治疗的效果。第14页，本讲稿共43页一次给药的药时曲线一次给药的药时曲线血血药药浓浓度度mg/l时间时间残留期残留期持续期持续期药峰时间药峰时间潜伏期潜伏期药峰浓度药峰浓度最低中毒浓度最低中毒浓度最低有效浓度最低有效浓度安安全全范范围围转化排泄过程转化排泄过程第15页，本讲稿共43页多次给药的药时曲线多次给药的药时曲线0 01 11 12 22 23 34 45 56

6、 6血血药药浓浓度度时间时间CmaxCmaxCminCmin治疗窗口治疗窗口第16页，本讲稿共43页药物消除类型药物消除类型 1 一级动力学消除（恒比消除）：一级动力学消除（恒比消除）：单位时间内按血药浓度的恒比进行消除。单位时间内按血药浓度的恒比进行消除。消除速度与血药浓度成正比。消除速度与血药浓度成正比。若以血药浓度（若以血药浓度（C）的对数与时间）的对数与时间（t）作图，为）作图，为一直线。一直线。第17页，本讲稿共43页2.零级动力学消除（恒量消除）：零级动力学消除（恒量消除）：单位时间内始终以一个恒定的数量进行消除。单位时间内始终以一个恒定的数量进行消除。消除速度与血药浓度无关。消除

7、速度与血药浓度无关。第18页，本讲稿共43页是指包括零级和一级动力学消除在内的混合型消除方式。如当药物剂量急剧增加或患者有某些疾病，血浓达饱和时，消除方式则可从一级动力学消除转变为零级动力学消除。如乙醇血浓0.05 mg/ml时，则可转成按零级动力学消除。3米氏消除动力学（米氏消除动力学（混合型消除混合型消除）：第19页，本讲稿共43页模型及其数值实现模型及其数值实现第20页，本讲稿共43页阅读材料:服药问题医生给病人开处方时必须注明两点医生给病人开处方时必须注明两点:服药的剂服药的剂量和服药的时间间隔量和服药的时间间隔.超剂量的药品会对身体产生严超剂量的药品会对身体产生严重不良后果重不良后果

8、,甚至死亡甚至死亡,而剂量不足而剂量不足,则不能达到治病则不能达到治病的目的的目的.已知患者服药后已知患者服药后,随时间推移随时间推移,药品在体内逐药品在体内逐渐被吸收渐被吸收,发生生化反应发生生化反应,也就是体内药品的浓度逐渐也就是体内药品的浓度逐渐降低降低.药品浓度降低的速率与体内当时药品的浓度成药品浓度降低的速率与体内当时药品的浓度成正比正比.当服药量为当服药量为A A、服药间隔为、服药间隔为T,T,试分析体内药的试分析体内药的浓度随时间的变化规律浓度随时间的变化规律.第21页，本讲稿共43页第22页，本讲稿共43页第23页，本讲稿共43页体内药的浓度随时间的变化规律第24页，本讲稿共4

9、3页Model 3:Population dynamicsIn this section we examine equations of the form y=f(y),called autonomous equations,where the independentvariable t does not appear explicitly.The main purpose of this section is to learn how geometricmethods can be used to obtain qualitative informationdirectly from dif

10、ferential equation without solving it.Simplest model:population growth rate is proportional tocurrent size of the population:Solution:exponential growth):第25页，本讲稿共43页Model 3:Population dynamicsLogistic Growth An exponential model y=ry,with solution y=ert,predictsunlimited growth,with rate r 0 indepe

11、ndent of population.Assuming instead that growth rate depends on populationsize,replace r by a function h(y)to obtain y=h(y)y.We want to choose growth rate h(y)so that h(y)r when y is small,h(y)decreases as y grows larger,and h(y)0.Our differential equation then becomes This equation is known as the

12、 Verhulst,or logistic,equation.第26页，本讲稿共43页 The logistic equation from the previous slide is This equation is often rewritten in the equivalent form where K=r/a.The constant r is called the intrinsic growth rate,and as we will see,K represents the carrying capacity of the population.A direction fiel

13、d for the logisticequation with r=1 and K=10is given here.第27页，本讲稿共43页Equilibrium solutions of the logistic equation Our logistic equation is Two equilibrium solutions are clearly present:In direction field below,with r=1,K=10,note behavior ofsolutions near equilibrium solutions:y=0 is unstable,y=K=

14、10 is asymptotically stable.第28页，本讲稿共43页Qualitative analysis of the logistic equation To better understand the nature of solutions to autonomousequations y=f(y),we start by graphing f(y)vs.y.In the case of logistic growth,that means graphing thefollowing function and analyzing its graph using calcul

15、us.第29页，本讲稿共43页Qualitative analysis,critical points The intercepts of f occur at y=0 and y=K,correspondingto the critical points of logistic equation.The vertex of the parabola is(K/2,rK/4),as shown below.第30页，本讲稿共43页Qualitative analysis,increasing/decreasing Note dy/dt 0 for 0 y K,so y is an increa

16、sing function oft there(indicate with right arrows along y-axis on 0 y K(indicatewith left arrows along y-axis on y K).In this context the y-axis is often called the phase line.第31页，本讲稿共43页Qualitative analysis,concavity Next,to examine concavity of y(t),we find y:Thus the graph of y is concave up wh

17、en f and f have samesign,which occurs when 0 y K.The graph of y is concave down when f and f have oppositesigns,which occurs when K/2 y K.Inflection point occurs at intersection of y and line y=K/2.第32页，本讲稿共43页Qualitative analysis,curve sketching Combining the information on the previous slides,we h

18、ave:Graph of y increasing when 0 y K.Slope of y approximately zero when y 0 or y K.Graph of y concave up when 0 y K.Graph of y concave down when K/2 y K.Inflection point when y=K/2.Using this information,we cansketch solution curves y fordifferent initial conditions.第33页，本讲稿共43页Qualitative analysis,

19、curve sketching Using only the information present in the differential equationand without solving it,we obtained qualitative informationabout the solution y.For example,we know where the graph of y is the steepest,and hence where y changes most rapidly.Also,y tendsasymptotically to the line y=K,for

20、 large t.The value of K is known as the carrying capacity,orsaturation level,for the species.Note how solution behavior differsfrom that of exponential equation,and thus the decisive effect ofnonlinear term in logistic equation.第34页，本讲稿共43页Model 3:Population dynamics Exact solution:separating variab

21、les Provided y 0 and y K,we can rewrite the logistic ODE:Expanding the left side using partial fractions,Thus the logistic equation can be rewritten as Integrating the above result,we obtain第35页，本讲稿共43页Exact solution:resolving for explicit solutionWe have:If 0 y0 K,then 0 y K and henceRewriting,usin

22、g properties of logs:第36页，本讲稿共43页Exact solution:resolving for explicit solutionWe have:for 0 y0 K.Also,this solution contains equilibrium solutions y=0 and y=K.Hence solution to logistic equation is第37页，本讲稿共43页模型预测的动态行为与大量的实验和观测数模型预测的动态行为与大量的实验和观测数据吻合据吻合第38页，本讲稿共43页 水的流出问题水的流出问题一横截面积为一横截面积为 A，高为，高为

23、H 的水池内盛满了水，有池底一的水池内盛满了水，有池底一横截面积为横截面积为 B 的小孔放水。设水从小孔流出的速度为的小孔放水。设水从小孔流出的速度为 ，求在任意时刻的水面高度和将水放空所需的时间，求在任意时刻的水面高度和将水放空所需的时间AB水面1水面2第39页，本讲稿共43页问题分析问题分析 从水面从水面1降到水面降到水面2所失去的水量等于从小孔流出的所失去的水量等于从小孔流出的水量水量 容器内水的体积为零容器内水的体积为零 即为容器内水的高度为零即为容器内水的高度为零建模与求解建模与求解1）从水面）从水面1将到水面将到水面2所失去的体积为所失去的体积为-在时间内在时间内 ，实际损失的体积

24、是，实际损失的体积是2）在同样时间内，水从小孔流出的体积为在同样时间内，水从小孔流出的体积为-是水在是水在 时间内从小孔流出保持水平前进时所经过时间内从小孔流出保持水平前进时所经过的距离的距离第40页，本讲稿共43页则则两端同除以两端同除以 ，并令并令 取极限得取极限得由于由于可得一阶方程：可得一阶方程：解为解为第41页，本讲稿共43页下面求将水放空的时间下面求将水放空的时间 t*t*令令h=0 代入上式得代入上式得例如例如,设设 A=0.54 m2,B=0.001 m2,H=4.9m,g=9.8m/s2.则则将水放空的时间为将水放空的时间为第42页，本讲稿共43页作业nP612，3，6第43页，本讲稿共43页