数学模型第四版第四章--微分方程模型.ppt

动态模型,描述对象特征随时间,(,空间,),的演变过程,.,分析对象特征的变化规律,.,预报对象特征的未来性态,.,研究控制对象特征的手段,.,根据函数及其变化率之间的关系确定函数,.,微分方程建模,根据建模目的和问题分析作出简化假设,.,按照内在规律或用类比法建立微分方程,.,5.1 传染病模型,描述传染病的传播过程,.,分析受感染人数的变化规律,.,预报传染病高潮到来的时刻,.,预防传染病蔓延的手段,.,不是从医学角度分析各种传染病的特殊机理,而是按照传播过程的一般规律建立数学模型.,背景 与,问题,传染病的极大危害(艾滋病、SARS、,),基本方法,已感染人数(病人),i,(,t,),每个病人每天有效接触,(,足以使人致病,),人数为,模型1,假设,若有效接触的是病人,则不能使病人数增加,必须区分已感染者(病人)和未感染者(健康人),建模,?,模型,3,传染病无免疫性病人治愈成为健康人，健康人可再次被感染.,增加假设,SIS,模型,3）病人每天治愈的比例为,日,治愈率,建模,日接触率,1/,感染期,一个感染期内,每个病人的有效接触人数，称为,接触数,.,m,l,s,/,=,模型,3,i,0,i,0,接触数,=1,阈值,感染期内,有效接触使健康者感染的人数不超过原有的病人数,1-1/,i,0,模型2(SI模型)如何看作模型3(SIS模型)的特例,i,di/dt,0,1,1,0,t,i,1,1-1/,i,0,t,1,di,/,dt,1,i,0,1/,i,(,t,),先升后降至0,P,2,:,s,0,1/,i,(,t,),单调降至0,1/,阈值,P,3,P,4,P,2,S,0,模型,4,SIR,模型,预防传染病蔓延的手段,(日接触率),卫生水平,(日,治愈率),医疗水平,传染病不蔓延的条件,s,0,1/,的估计,降低,s,0,提高,r,0,提高阈值,1/,降低,(=,/,),群体免疫,忽略,i,0,模型,4,预防传染病蔓延的手段,降低日接触率,提高日,治愈率,提高移出比例,r,0,以最终未感染比例,s,和病人比例最大值,i,m,为度量指标.,1,/,s,0,i,0,s,i,1,0.3,0.3,0.98,0.02,0.0398,0.3449,0.6,0.3,0.5,0.98,0.02,0.1965,0.1635,0.5,0.5,1.0,0.98,0.02,0.8122,0.0200,0.4,0.5,1.25,0.98,0.02,0.9172,0.0200,1,0.3,0.3,0.70,0.02,0.0840,0.1685,0.6,0.3,0.5,0.70,0.02,0.3056,0.0518,0.5,0.5,1.0,0.70,0.02,0.6528,0.0200,0.4,0.5,1.25,0.70,0.02,0.6755,0.0200,s,0,(,r,0,),s,i,m,s,i,m,模型,4,SIR,模型,被传染人数的估计,记被传染人数比例,x,0,3)经济增长的条件,劳动力相对增长率,每个劳动力的产值,Z,(,t,)=,Q,(,t,)/,L,(,t,)增长,dZ/dt,0,3)经济增长的条件,劳动力增长率小于初始投资增长率,5.3 正规战与游击战,战争分类：正规战争，游击战争，混合战争.,只考虑双方兵力多少和战斗力强弱.,兵力因战斗及非战斗减员而减少，因增援而增加.,战斗力与射击次数及命中率有关.,建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例.,第一次世界大战,Lanchester,提出预测战役结局的模型.,一般模型,每方战斗减员率取决于双方的兵力和战斗力,.,每方非战斗减员率与本方兵力成正比,.,甲乙双方的增援率为,u,(,t,),v,(,t,).,f,g,取决于战争类型,x,(,t,)甲方兵力，,y,(,t,)乙方兵力,模型假设,模型,正规战争模型,甲方战斗减员率只取决于乙方的兵力和战斗力,双方均以正规部队作战,忽略非战斗减员,假设没有增援,f,(,x,y,)=,ay,a,乙方每个士兵的杀伤率,a,=,r,y,p,y,r,y,射击率，,p,y,命中率,0,正规战争模型,为判断战争的结局，不求,x,(,t,),y,(,t,)而在相平面上讨论,x,与,y,的关系.,平方律 模型,乙方胜,游击战争模型,双方都用游击部队作战,甲方战斗减员率还随着甲方兵力的增加而增加,忽略非战斗减员,假设没有增援,f,(,x,y,)=,cxy,c,乙方每个士兵的杀伤率,c,=,r,y,p,y,r,y,射击率,p,y,命中率,p,y,=s,ry,/,s,x,s,x,甲方活动面积,s,ry,乙方射击有效面积,0,游击战争模型,线性律 模型,0,混合战争模型,甲方为游击部队，乙方为正规部队,乙方必须10倍于甲方的兵力!,设,x,0,=100,r,x,/,r,y,=1/2,p,x,=0.1,s,x,=1(km,2,),s,ry,=1(m,2,),5.4,药物在体内的分布与排除,药物进入机体形成,血药浓度,(,单位体积血液的药物量,).,血药浓度需保持在一定范围内,给药方案设计,.,药物在体内吸收、分布和排除过程,药物动力学,.,建立,房室模型,药物动力学的基本步骤,.,房室,机体的一部分，药物在一个房室内均匀分布,(,血药浓度为常数,),，在房室间按一定规律转移,.,本节讨论,二室模型,中心室,(,心、肺、肾等,),和周边室,(,四肢、肌肉等,).,模型假设,中心室,(1),和周边室,(2),容积不变,.,药物在房室间转移速率及向体外排除速率与该室血药浓度成正比,.,药物从体外进入中心室，在二室间相互转移,从中心室排出体外,.,模型建立,中心室,周边室,给药,排除,c,1,(,t,),x,1,(,t,),V,1,c,2t,),x,2,(,t,),V,2,转移,线性常系数非齐次方程,对应齐次方程通解,模型建立,几种常见的给药方式,1,.快速静脉注射,t,=0,瞬时注射剂量,D,0,的药物进入中心室,血药浓度立即为,D,0,/,V,1,给药速率,f,0,(,t,),和初始条件,2,.恒速静脉滴注,t,T,c,1,(,t,)和,c,2,(,t,)按指数规律趋于零,0,t,T,药物以速率,k,0,进入中心室,3,.口服或肌肉注射,相当于药物(剂量,D,0,)先进入吸收室，吸收后进入中心室.,吸收室药量,x,0,(,t,),吸收室,中心室,D,0,参数估计,各种给药方式下的,c,1,(,t,),c,2,(,t,)取决于参数,k,12,k,21,k,13,V,1,V,2,t,=0,快速静脉注射,D,0,在,t,i,(,i,=1,2,n,),测得,c,1,(,t,i,),由较大的 用最小二乘法定,A,由较小的 用最小二乘法定,B,参数估计,进入中心室的药物全部排除,建立,房室模型,研究体内,血药浓度,变化过程,确定,转移速率、排除速率等参数,为制订给药方案提供依据,.,机理分析确定模型形式，测试分析估计模型参数,.,药物在体内的分布与排除,房室模型：,一室模型,二室模型,多室模型,非线性(一室)模型,c,1,较小时近似于线性,一级,排除过程,如,c,1,较大时近似于常数,零级,排除过程,过滤嘴的作用与它的材料和长度有什么关系,?,人体吸入的毒物量与哪些因素有关，其中什么因素影响大，什么因素影响小,?,模型分析,分析吸烟时毒物进入人体的过程，建立吸烟过程的数学模型,.,设想一个,“,机器人,”,在典型环境下吸烟，吸烟方式和外部环境在整个过程中不变,.,问题,5.5,香烟过滤嘴的作用,模型假设,定性分析,1）,l,1,烟草长，,l,2,过滤嘴长，,l=l,1,+,l,2,，,毒物量,M,均匀分布，密度,w,0,=,M,/,l,1,.,2）点燃处毒物随烟雾进入空气和沿香烟穿行的数量比是,a,:,a,a,+,a,=1.,3）未点燃的烟草和过滤嘴对随烟雾穿行的毒物的(单位时间)吸收率分别是,b,和,.4）烟雾沿香烟穿行速度是常数,v,，香烟燃烧速度是常数,u,，,v u.,Q,吸一支烟毒物进入人体总量,模型建立,0,t,=0,x,=0，点燃香烟,q,(,x,t,)毒物流量,w,(,x,t,)毒物密度,1)求,q,(,x,0)=,q,(,x,),流量守恒,t,时刻，香烟燃至,x=ut,1)求,q,(,x,0)=,q,(,x,),2)求,q,(,l,t,),3)求,w,(,ut,t,),考察,t,内毒物密度的增量,(单位长度烟雾毒物被吸收部分),4)计算,Q,Q,吸一支烟毒物进入人体总量,结果分析,烟草,为什么有作用?,1）,Q,与,a,M,成正比，,aM,是毒物集中在,x,=,l,处的吸入量,2)过滤嘴因素，,l,2,负指数,作用,是毒物集中在,x,=,l,1,处的吸入量,3）,(,r,)烟草的吸收作用,b,l,1,线性,作用,带过滤嘴,不带过滤嘴,结果分析,4)与另一支不带过滤嘴的香烟比较，,w,0,b,a,v,l,均相同，吸至,x=l,1,扔掉.,提高,-,b,与加长,l,2,，效果相同.,香烟过滤嘴的作用,在基本合理的简化假设下，用精确的数学工具解决一个看来不易下手的实际问题,.,引入两个基本函数：流量,q,(,x,t,),和密度,w,(,x,t,),，运用物理学的守恒定律建立微分方程，构造动态模型,.,对求解结果进行定性和定量分析，得到合乎实际的结论,.,背景,年 1625 1830 1930 1960 1974 1987 1999,人口(亿)5 10 20 30 40 50 60,世界人口增长概况,中国人口增长概况,年 1908 1933 1953 1964 1982 1990 1995 2000,人口(亿)3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0,研究人口变化规律,控制人口过快增长,5.6 人口预测和控制,做出较准确的预报,建立人口数学模型,指数增长模型马尔萨斯,1798年,提出,常用的计算公式,x,(,t,),时刻,t,的,人口,基本假设,:人口(相对)增长率,r,是常数,今年人口,x,0,年增长率,r,k,年后人口,随着时间增加，人口按指数规律无限增长.,与常用公式的一致,rt,e,x,t,x,0,),(,=,?,指数增长模型的应用及局限性,与,19,世纪以前欧洲一些地区人口统计数据吻合,.,适用于,19,世纪后迁往加拿大的欧洲移民后代,.,可用于短期人口增长预测,.,不符合,19,世纪后多数地区人口增长规律,.,不能预测较长期的人口增长过程,.,19世纪后人口数据,人口增长率,r,不是常数(逐渐下降),阻滞增长模型逻辑斯蒂(,Logistic,),模型,人口增长到一定数量后，增长率下降的原因：,资源、环境等因素对人口增长的阻滞作用,且阻滞作用随人口数量增加而变大,假设,r,固有增长率(,x,很小时),x,m人口容量（资源、环境能容纳的最大数量）,r,是,x,的减函数,dx,/,dt,x,0,x,m,x,m,/2,t,x,0,x,增加先快后慢,x,m,x,0,x,m,/2,阻滞增长模型(,Logistic,模型),指数增长模型,Logistic,模型的应用,经济领域中的增长规律,(,耐用消费品的售量,).,种群数量模型,(,鱼塘中的鱼群,森林中的树木,).,S,形曲线,参数估计,用指数增长模型或阻滞增长模型作人口,预报，必须先估计模型参数,r,或,r,x,m,.,模型的参数估计、检验和预报,指数增长模型,阻滞增长模型,由统计数据用,线性最小二乘法,作参数估计,例：美国人口数据(百万),1860 1870 1880 1960 1970 1980 1990 2000,31.4 38.6 50.2 179.3 204.0 226.5 251.4 281.4,t,x,r,=0.2022/10,年，,x,0,=6.0450,模型的参数估计、检验和预报,指数增长模型,阻滞增长模型,r,=0.2557/10,年，,x,m,=392.0886,年,实际人口,计算人口,(指数增长模型),计算人口,(阻滞增长模型),1790,3.9,6.0,3.9,1800,5.3,7.4,5.0,1960,179.3,188.0,171.3,1970,204.0,230.1,196.2,1980,226.5,281.7,221.2,1990,251.4,344.8,245.3,2000,422.1,指数增长模型,阻滞增长模型,模型检验,用模型计算2000年美国人口,误差不到3%,与实际数据(2000年,281.4,)比较,=,274.5,模型的参数估计、检验和预报,为做模型检验在参数估计时未用,2000,年实际数据,模型应用,加入2000年人口数据后重新估计模型参数,r,=0.2490,x,m,=434.0,x,(2010)=306.0,预报美国,2010,年的人口,考虑年龄结构和生育模式的人口模型,年龄分布对于人口预测的重要性,.,只考虑自然出生与死亡，不计迁移,.,人口发展方程,F,(,r,t,)人口分布函数(年龄,r,的人口),p,(,r,t,)人口密度函数,N,(,t,)人口总数,r,m,(,),最高年龄,人口发展方程,一阶偏微分方程,人口发展方程,0,t,r,定解条件,已知函数(人口调查),生育率(控制手段),生育率,f,(,t,),的分解,总和生育率,h,生育模式,0,人口控制系统,总和生育率控制生育的多少,生育模式控制生育的早晚和疏密,正反馈系统,滞后作用很大,输入,输入,输出,反馈,人口指数,1）人口总数,2）平均年龄,3）平均寿命,t,时刻出生的人，死亡率按,(,r,t,),计算的平均存活时间,4）老龄化指数,控制生育率,控制,N,(,t,),不过大,控制,(,t,),不过高5.7,烟雾的扩散与消失,现象和,问题,炮弹在空中爆炸，烟雾向四周扩散，形,成圆形不透光区域,.,不透光区域不断扩大，然后区域边界逐,渐明亮，区域缩小，最后烟雾消失,.,建立模型描述烟雾扩散和消失过程，分,析消失时间与各因素的关系,.,问题分析,无穷空间由瞬时点源导致的扩散过程，,用二阶偏微分方程描述烟雾浓度的变化,.,观察到的烟雾消失与烟雾对光线的吸收、,以及仪器对明暗的灵敏程度有关,.,模型假设,1）烟雾在无穷空间扩散，不受大地和风,的影响；扩散服从扩散定律.,2）光线穿过烟雾时光强的相对减少与烟雾,浓度成正比；无烟雾的大气不影响光强.,3）穿过烟雾进入仪器的光线只有明暗之分，,明暗界限由仪器灵敏度决定.,模型建立,1）烟雾浓度 的变化规律,扩散定律：,单位时间通过单位法向面积的流量,q,与浓度,C,的,梯度成正比.,曲面积分奥氏公式,1）烟雾浓度 的变化规律,的微分形式，并利用积分中值定理,初始条件,Q,炮弹释放的烟雾总量,单位强度的点源函数,对任意,t,C,的等值面是球面,x,2,+,y,2,+,z,2,=,R,2,;,R,C,仅当,t,对任意点,(,x,y,z,),C,0,1）烟雾浓度 的变化规律,2）光强穿过烟雾时的变化规律,假设2）光强的相对减少与烟雾浓度成正比.,I,(,l,)沿,l,方向的光强，,C,(,l,)沿,l,方向的烟雾强度,记未进入烟雾(,l,l,0,),时光强为,I,(,l,0,)=,I,0,3）仪器灵敏度与烟雾明暗界限,烟雾浓度连续变化,烟雾中光强连续变化,仪器,z,-,设光源在z=-,仪器在z=,则观测到的,明暗界限为,不透光区域有扩大、缩小、消失的过程,穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定.,不透光区域边界,4）不透光区域边界的变化规律,对任意,t,不透光区域边界是圆周,不透光区域边界半径,r,(,t,),r,m,0,t,1,t,2,t,结果分析,观测到不透光区域边界达到最大的时刻,t,1,，可以预报烟雾消失的时刻,t,2,5.8 万有引力定律的发现,背景,航海业发展,天文观测精确,“地心说”动摇,哥白尼：“日心说”,伽利略：落体运动,开普勒：行星运动三定律,变速运动的计算方法,牛顿：一切运动有力学原因,牛顿运动三定律,牛顿：研究变速运动，发明微积分（流数法）,开普勒三定律,牛顿运动第二定律,万有引力定律,自然科学之数学原理(1687),模型假设,极坐标系(,r,),1.行星轨道,a,长半轴,b,短半轴,e,离心率,3.行星运行周期,T,行星位置：向径,2.单位时间,扫过面积为常数,A,m,行星质量,绝对常数,4.行星运行受力,太阳(0,0),O,r,P,行星,模型建立,O(太阳),P,(行星),r,向径 的基向量,模型建立,只需证明,4,A,2,/,p=kM,（,A,2,/,p,与哪一颗行星无关）,A,单位时间,扫过面积,与万有引力定律,比较,T,运行周期,可以证明(习题10),l,p,/,/,2,2,=,p,A,4,A,2,/,p=kM,在正确假设基础上用数学建模方法得到万有引力定律,