资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第五章 微分方程模型,5.1,传染病模型,5.2,经济增长模型,5.3,正规战与游击战,5.4,药物在体内的分布与排除,5.5,香烟过滤嘴的作用,5.6,人口预测和控制,5.7,烟雾的扩散与消失,5.8,万有引力定律的发现,动态模型,描述对象特征随时间,(,空间,),的演变过程,分析对象特征的变化规律,预报对象特征的未来性态,研究控制对象特征的手段,根据函数及其变化率之间的关系确定函数,微分方程建模,根据建模目的和问题分析作出简化假设,按照内在规律或用类比法建立微分方程,5.1,传染病模型,问题,描述传染病的传播过程,分析受感染人数的变化规律,预报传染病高潮到来的时刻,预防传染病蔓延的手段,按照传播过程的一般规律，用机理分析方法建立模型,已感染人数,(,病人,),i,(,t,),每个病人每天有效接触,(,足以使人致病,),人数为,模型,1,假设,若有效接触的是病人，则不能使病人数增加,必须区分已感染者,(,病人,),和未感染者,(,健康人,),建模,?,模型,2,区分已感染者,(,病人,),和未感染者,(,健康人,),假设,1,）总人数,N,不变，病人和健康 人的 比例分别为,2,）每个病人每天有效接触人数为,且,使接触的健康人致病,建模,日,接触率,SI,模型,模型,2,1/2,t,m,i,i,0,1,0,t,t,m,传染病高潮到来时刻,(,日接触率,) ,t,m,Logistic,模型,病人可以治愈！,?,t=t,m,di,/,dt,最大,模型,3,传染病无免疫性,病人治愈成为健康人，健康人可再次被感染,增加假设,SIS,模型,3,）病人每天治愈的比例为,日,治愈率,建模,日接触率,1/,感染期,一个感染期内,每个病人的有效接触人数，称为,接触数,。,模型,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,(,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/,降低,(=,/,),群体免疫,模型,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),容积不变,药物在房室间转移速率及向体外排除速率，与该室血药浓度成正比,药物从体外进入中心室，在二室间相互转移,从中心室排出体外,模型建立,线性常系数非齐次方程,对应齐次方程通解,模型建立,几种常见的给药方式,1,.,快速静脉注射,t,=0,瞬时,注射剂量,D,0,的药物进入中心室,血药浓度立即为,D,0,/,V,1,给药速率,f,0,(,t,),和初始条件,2,.,恒速静脉滴注,t,T,c,1,(,t,),和,c,2,(,t,),按,指数规律趋于零,药物以速率,k,0,进入中心室,0,T,t,吸收室,中心室,3,.,口服或肌肉注射,相当于药物,(,剂量,D,0,),先进入吸收室，吸收后进入中心室,吸收室药量,x,0,(,t,),参数估计,各种给药方式下的,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,参数估计,进入中心室的药物全部排除,过滤嘴的作用与它的材料和长度有什么关系,人体吸入的毒物量与哪些因素有关，其中哪些因素影响大，哪些因素影响小。,模型分析,分析吸烟时毒物进入人体的过程，建立吸烟过程的数学模型。,设想一个,“,机器人,”,在典型环境下吸烟，吸烟方式和外部环境认为是不变的。,问题,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,),4),计算,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,，,效果相同,5.6,人口预测和控制,年龄分布对于人口预测的重要性,只考虑自然出生与死亡，不计迁移,人口发展方程,人口发展方程,一阶偏微分方程,人口发展方程,已知函数（人口调查）,生育率（控制人口手段）,0,t,r,生育率的分解,总和生育率,h,生育模式,0,人口发展方程和生育率,总和生育率,控制生育的多少,生育模式,控制生育的早晚和疏密,正反馈系统,滞后作用很大,人口指数,1,）人口总数,2,）平均年龄,3,）平均寿命,t,时刻出生的人，死亡率按,(,r,t,),计算的平均存活时间,4,）老龄化指数,控制生育率,控制,N,(,t,),不过大,控制,(,t,),不过高,5.7,烟雾的扩散与消失,现象和,问题,炮弹在空中爆炸，烟雾向四周扩散，形成圆形不透光区域。,不透光区域不断扩大，然后区域边界逐渐明亮，区域缩小，最后烟雾消失。,建立模型描述烟雾扩散和消失过程，分析消失时间与各因素的关系。,问题分析,无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述烟雾浓度的变化。,观察的烟雾消失与烟雾对光线的吸收，以及仪器对明暗的灵敏程度有关。,模型假设,1,）烟雾在无穷空间扩散，不受大地和风的影响；扩散服从热传导定律。,2,）光线穿过烟雾时光强的减少与烟雾浓度成正比；无烟雾的大气不影响光强。,3,）穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定。,模型建立,1,）烟雾浓度 的变化规律,热传导定律：单位时间通过单位法向面积的流量与浓度梯度成正比,曲面积分的奥氏公式,1,）烟雾浓度 的变化规律,初始条件,Q,炮弹释放的烟雾总量,单位强度的点源函数,对任意,t, C,的等值面是球面,x,2,+,y,2,+,z,2,=,R,2,;,R,C,仅当,t,对任意点,(,x,y,z,),C,0,1,）烟雾浓度 的变化规律,2,）穿过烟雾光强的变化规律,光强的减少与烟雾浓度成正比,3,）仪器灵敏度与烟雾明暗界限,烟雾浓度连续变化,烟雾中光强连续变化,仪器,z,-,设光源在,z,=-,仪器在,z,=,则观测到的,明暗界限为,不透光区域有扩大、缩小、消失的过程,穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定。,不透光区域边界,4,）不透光区域边界的变化规律,对任意,t,不透光区域边界是圆周,不透光区域边界半径,r,(,t,),r,m,0,t,1,t,2,t,结果分析,观测到不透光区域边界达到最大的时刻,t,1,，,可以预报烟雾消失的时刻,t,2,5.8,万有引力定律的发现,背景,航海业发展,天文观测精确,“,地心说”动摇,哥白尼：“日心说”,伽里,略：落体运动,开普,勒：行星运动三定律,变速运动的计算方法,牛顿：一切运动有力学原因,牛顿运动三定律,牛顿：研究变速运动，发明微积分（流数法）,开普,勒三定律,牛顿运动第二定律,万有引力定律,自然科学之数学原理,(1687),模型假设,极,坐标系,(,r,),太阳,(0,0,),1.,行星轨道,a,长半轴,b,短半轴,e,离心率,3.,行星运行周期,T,行星位置：向径,O (,太阳,),P (,行星,),r,2.,单位时间,扫,过,面积为常数,A,m,行星质量,绝对常数,4.,行星运行受力,模型建立,O (,太阳,),P (,行星,),r,向径 的基向量,模型建立,万有引力定律,需,证明,4,A,2,/,p =,kM,（与,哪一颗行星无关）,A,单位时间,扫,过,面积,O (,太阳,),P (,行星,),r,(,习题,),l,p,/,/,2,2,=,p,A,
展开阅读全文