《第一章一阶微分方程的应用PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第一章一阶微分方程的应用PPT讲稿.ppt（43页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第一章一阶微分方程的应用第1页，共43页，编辑于2022年，星期一1：应用大意应用大意适应范围适应范围与变化率有关的各种实际问题应用三步曲应用三步曲(1)建立模型建立模型 Modelling(2)模型求解模型求解 Solving(3)模型应用模型应用 Application建议:模型要详略得当第2页，共43页，编辑于2022年，星期一应用一：曲线族的等角轨线应用一：曲线族的等角轨线设给定一个平面上以设给定一个平面上以C C为参数的曲线族为参数的曲线族（*）我们设法求出另一个以我们设法求出另一个以k k为参数的曲线族为参数的曲线族（*）使得曲线族（使得曲线族（*）中的任一条曲线与曲线族）中的任一
2、条曲线与曲线族样的曲线族（样的曲线族（*）是已知曲线族（）是已知曲线族（*）的）的（*）中的每一条曲线相交时成定角）中的每一条曲线相交时成定角 则称这则称这等角轨线族。等角轨线族。第3页，共43页，编辑于2022年，星期一当当 时，称曲线族（时，称曲线族（*）是（）是（*）的）的正交轨线族。正交轨线族。例如：曲线族例如：曲线族 是曲线族是曲线族 的正交轨线族。的正交轨线族。第4页，共43页，编辑于2022年，星期一第5页，共43页，编辑于2022年，星期一设设 y=y(x)为为(C)中任一条曲线，于是存在相应的中任一条曲线，于是存在相应的C，使得，使得 因为要求因为要求x,y,y 的关系，将上
3、式对的关系，将上式对x求导数，得求导数，得 (1.84)这样，将上两式联立，即由这样，将上两式联立，即由上述关系式成为曲线族满足的微分方程上述关系式成为曲线族满足的微分方程第6页，共43页，编辑于2022年，星期一例例1 1 求抛物线族求抛物线族 的正交轨线族。的正交轨线族。解：解：对方程两边关于对方程两边关于x x求导得求导得由由 解出解出C C代入上式得曲线族代入上式得曲线族 在点在点 处切线斜率为处切线斜率为 第7页，共43页，编辑于2022年，星期一由于所求曲线族的曲线与由于所求曲线族的曲线与 中的曲线在中的曲线在 正交，故满足方程正交，故满足方程 这是一个变量可分离方程求解得这是一个
4、变量可分离方程求解得 的正交的正交曲线族为曲线族为这是一个椭圆，如右图这是一个椭圆，如右图放大此图 图图2.162.16第8页，共43页，编辑于2022年，星期一第9页，共43页，编辑于2022年，星期一应用二：应用二：雨滴的下落雨滴的下落考虑雨滴在高空形成后下落的过程中速度的变化三种不同的假设(1)自由落体运动(2)小阻力的情况(3)大阻力的情况第10页，共43页，编辑于2022年，星期一(1)自由落体运动下落过程中没有任何阻力下落过程中没有任何阻力第11页，共43页，编辑于2022年，星期一(2)小阻力的情况下落过程中阻力与速度和半径的乘积成比例下落过程中阻力与速度和半径的乘积成比例 第1
5、2页，共43页，编辑于2022年，星期一(3)大阻力的情况下落过程中阻力与速度和半径的乘积平方成比例下落过程中阻力与速度和半径的乘积平方成比例第13页，共43页，编辑于2022年，星期一三、药物设计医生给病人开处方是必须注意两点：医生给病人开处方是必须注意两点：服药的剂量和服药的时间间隔。服药的剂量和服药的时间间隔。超剂量的药物会对患者产生严重不良后果，超剂量的药物会对患者产生严重不良后果，甚至死亡；剂量不足，则不能达到治疗的效果。甚至死亡；剂量不足，则不能达到治疗的效果。第14页，共43页，编辑于2022年，星期一一次给药的药时曲线一次给药的药时曲线血血药药浓浓度度mg/l时间时间残留期残留
6、期持续期持续期药峰时间药峰时间潜伏期潜伏期药峰浓度药峰浓度最低中毒浓度最低中毒浓度最低有效浓度最低有效浓度安安全全范范围围转化排泄过程转化排泄过程第15页，共43页，编辑于2022年，星期一多次给药的药时曲线多次给药的药时曲线0 01 11 12 22 23 34 45 56 6血血药药浓浓度度时间时间CmaxCmaxCminCmin治疗窗口治疗窗口第16页，共43页，编辑于2022年，星期一药物消除类型药物消除类型 1 一级动力学消除（恒比消除）：一级动力学消除（恒比消除）：单位时间内按血药浓度的恒比进行消除。单位时间内按血药浓度的恒比进行消除。消消除速度与血药浓度成正比。除速度与血药浓度成
7、正比。若以血药浓度（若以血药浓度（C）的对数与时间）的对数与时间（t）作图，）作图，为一直线。为一直线。第17页，共43页，编辑于2022年，星期一2.零级动力学消除（恒量消除）：零级动力学消除（恒量消除）：单位时间内始终以一个恒定的数量进行消单位时间内始终以一个恒定的数量进行消除。除。消除速度与血药浓度无关。消除速度与血药浓度无关。第18页，共43页，编辑于2022年，星期一是指包括零级和一级动力学消除在内的混合型消除方式。如当药物剂量急剧增加或患者有某些疾病，血浓达饱和时，消除方式则可从一级动力学消除转变为零级动力学消除。如乙醇血浓0.05 mg/ml时，则可转成按零级动力学消除。3米氏消
8、除动力学（米氏消除动力学（混合型消除混合型消除）：第19页，共43页，编辑于2022年，星期一模型及其数值实现模型及其数值实现第20页，共43页，编辑于2022年，星期一阅读材料:服药问题医生给病人开处方时必须注明两点医生给病人开处方时必须注明两点:服药的剂服药的剂量和服药的时间间隔量和服药的时间间隔.超剂量的药品会对身体产生严超剂量的药品会对身体产生严重不良后果重不良后果,甚至死亡甚至死亡,而剂量不足而剂量不足,则不能达到治则不能达到治病的目的病的目的.已知患者服药后已知患者服药后,随时间推移随时间推移,药品在体内药品在体内逐渐被吸收逐渐被吸收,发生生化反应发生生化反应,也就是体内药品的浓也
9、就是体内药品的浓度逐渐降低度逐渐降低.药品浓度降低的速率与体内当时药品的药品浓度降低的速率与体内当时药品的浓度成正比浓度成正比.当服药量为当服药量为A A、服药间隔为、服药间隔为T,T,试分析体内试分析体内药的浓度随时间的变化规律药的浓度随时间的变化规律.第21页，共43页，编辑于2022年，星期一第22页，共43页，编辑于2022年，星期一第23页，共43页，编辑于2022年，星期一体内药的浓度随时间的变化规律第24页，共43页，编辑于2022年，星期一Model 3:Population dynamicsIn this section we examine equations of the
10、 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 differential equation without solving it.Simplest model:population growth r
11、ate is proportional tocurrent size of the population:Solution:exponential growth):第25页，共43页，编辑于2022年，星期一Model 3:Population dynamicsLogistic Growth An exponential model y=ry,with solution y=ert,predictsunlimited growth,with rate r 0 independent of population.Assuming instead that growth rate depends
12、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 Verhulst,or logistic,equation.第26页，共43页，编辑于2022年，星期一 The logi
13、stic 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 field for the logisticequation with r=1 and K=10is given
14、 here.第27页，共43页，编辑于2022年，星期一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=10 is asymptotically stable.第28页，共43页，编辑于2
15、022年，星期一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 calculus.第29页，共43页，编辑于2022年，星期一Qualita
16、tive 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页，编辑于2022年，星期一Qualitative analysis,increasing/decreasing Note dy/dt 0 for 0 y K,so y is an increasing functio
17、n 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页，编辑于2022年，星期一Qualitative analysis,concavity Next,to examine concavity of y(t),we find y:Thus the graph of y is concave up when
18、 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页，编辑于2022年，星期一Qualitative analysis,curve sketching Combining the information on the previous sli
19、des,we have: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页，编辑于2022年，星期一Qua
20、litative analysis,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 t
21、o the line y=K,for 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页，编辑于2022年，星期一Model 3:Population dynamics Exac
22、t solution:separating variables 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页，编辑于2022年，星期一Exact solution:resolving for explicit solutionWe have:If 0 y
23、0 K,then 0 y K and henceRewriting,using properties of logs:第36页，共43页，编辑于2022年，星期一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页，编辑于2022年，星期一模型预测的动态行为与大量的实验和观测数模型预测的动态行为与大量的
24、实验和观测数据吻合据吻合第38页，共43页，编辑于2022年，星期一 水的流出问题水的流出问题一横截面积为一横截面积为 A，高为，高为 H 的水池内盛满了水，有池底一的水池内盛满了水，有池底一横截面积为横截面积为 B 的小孔放水。设水从小孔流出的速度为的小孔放水。设水从小孔流出的速度为 ，求在任意时刻的水面高度和将水放空所需的时间，求在任意时刻的水面高度和将水放空所需的时间AB水面1水面2第39页，共43页，编辑于2022年，星期一问题分析问题分析 从水面从水面1降到水面降到水面2所失去的水量等于从小孔流出的水量所失去的水量等于从小孔流出的水量 容器内水的体积为零容器内水的体积为零 即为容器内
25、水的高度为零即为容器内水的高度为零建模与求解建模与求解1）从水面）从水面1将到水面将到水面2所失去的体积为所失去的体积为-在时间内在时间内 ，实际损失的体积是，实际损失的体积是2）在同样时间内，水从小孔流出的体积为在同样时间内，水从小孔流出的体积为-是水在是水在 时间内从小孔流出保持水平前进时所经过时间内从小孔流出保持水平前进时所经过的距离的距离第40页，共43页，编辑于2022年，星期一则则两端同除以两端同除以 ，并令并令 取极限得取极限得由于由于可得一阶方程：可得一阶方程：解为解为第41页，共43页，编辑于2022年，星期一下面求将水放空的时间下面求将水放空的时间 t*t*令令h=0 代入上式得代入上式得例如例如,设设 A=0.54 m2,B=0.001 m2,H=4.9m,g=9.8m/s2.则将则将水放空的时间为水放空的时间为第42页，共43页，编辑于2022年，星期一作业nP612，3，6第43页，共43页，编辑于2022年，星期一