数学模型姜启源微分方程模型.pptx

《数学模型姜启源微分方程模型.pptx》由会员分享，可在线阅读，更多相关《数学模型姜启源微分方程模型.pptx（62页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、5.1 传染病模型传染病模型问题问题 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律，用机理分析方法建立模型第1页/共62页 已感染人数(病人)i(t)每个病人每天有效接触(足以使人致病)人数为 模型模型1 1假设假设若有效接触的是病人，则不能使病人数增加必须区分已感染者(病人)和未感染者(健康人)建模建模?第2页/共62页模型模型2 2区分已感染者(病人)和未感染者(健康人)假设假设1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为,且使接触的健康人致病建模建模 日接触率SI 模型模型第3页/

2、共62页模型模型21/2tmii010ttm传染病高潮到来时刻 (日接触率)tm Logistic 模型病人可以治愈！?t=tm,di/dt 最大第4页/共62页模型模型3传染病无免疫性病人治愈成为健康人，健康人可再次被感染增加假设增加假设SIS 模型模型3）病人每天治愈的比例为 日治愈率建模建模 日接触率1/感染期 一个感染期内每个病人的有效接触人数，称为接触数。第5页/共62页模型模型3i0i0接触数 =1 阈值感染期内有效接触感染的健康者人数不超过病人数1-1/i0模型2(SI模型)如何看作模型3(SIS模型)的特例idi/dt01 10ti 11-1/i0t 1di/dt 1/i(t)

3、先升后降至0P2:s01/i(t)单调降至01/阈值P3P4P2S0第10页/共62页模型模型4SIR模型模型预防传染病蔓延的手段预防传染病蔓延的手段 (日接触率)卫生水平 (日治愈率)医疗水平 传染病不蔓延的条件s01/的估计 降低 s0提高 r0 提高阈值 1/降低 (=/),群体免疫第11页/共62页模型模型4SIR模型模型被传染人数的估计被传染人数的估计记被传染人数比例x 03)经济增长的条件第19页/共62页劳动力增长率小于初始投资增长率每个劳动力的产值 Z(t)=Q(t)/L(t)增长dZ/dt03)经济增长的条件第20页/共62页5.3 正规战与游击战正规战与游击战战争分类：正规

4、战争，游击战争，混合战争只考虑双方兵力多少和战斗力强弱兵力因战斗及非战斗减员而减少，因增援而增加战斗力与射击次数及命中率有关建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例第一次世界大战Lanchester提出预测战役结局的模型第21页/共62页一般模型一般模型 每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为u(t),v(t)f,g 取决于战争类型x(t)甲方兵力，y(t)乙方兵力模型模型假设假设模型模型第22页/共62页正规战争模型正规战争模型 甲方战斗减员率只取决于乙方的兵力和战斗力双方均以正规部队作战 忽略非战斗减员 假设没

5、有增援f(x,y)=ay,a 乙方每个士兵的杀伤率a=ry py,ry 射击率，py 命中率第23页/共62页0正规战争模型正规战争模型为判断战争的结局，不求x(t),y(t)而在相平面上讨论 x 与 y 的关系平方律 模型乙方胜第24页/共62页游击战争模型游击战争模型双方都用游击部队作战 甲方战斗减员率还随着甲方兵力的增加而增加 忽略非战斗减员 假设没有增援f(x,y)=cxy,c 乙方每个士兵的杀伤率c=ry pyry射击率py 命中率py=sry/sxsx 甲方活动面积sry 乙方射击有效面积第25页/共62页0游击战争模型游击战争模型线性律 模型第26页/共62页0混合战争模型混合战

6、争模型甲方为游击部队，乙方为正规部队乙方必须10倍于甲方的兵力设 x0=100,rx/ry=1/2,px=0.1,sx=1(km2),sry=1(m2)第27页/共62页5.4 药物在体内的分布与排除药物在体内的分布与排除 药物进入机体形成血药浓度(单位体积血液的药物量)血药浓度需保持在一定范围内给药方案设计 药物在体内吸收、分布和排除过程 药物动力学 建立房室模型药物动力学的基本步骤 房室机体的一部分，药物在一个房室内均匀分布(血药浓度为常数)，在房室间按一定规律转移 本节讨论二室模型中心室(心、肺、肾等)和周边室(四肢、肌肉等)第28页/共62页 中心室周边室给药排除模型假设模型假设 中心

7、室(1)和周边室(2),容积不变 药物在房室间转移速率及向体外排除速率，与该室血药浓度成正比 药物从体外进入中心室，在二室间相互转移,从中心室排出体外模型建立模型建立第29页/共62页线性常系数非齐次方程对应齐次方程通解模型建立模型建立第30页/共62页几种常见的给药方式几种常见的给药方式1.快速静脉注射t=0 瞬时注射剂量D0的药物进入中心室,血药浓度立即为D0/V1给药速率 f0(t)和初始条件第31页/共62页2.恒速静脉滴注t T,c1(t)和 c2(t)按指数规律趋于零药物以速率k0进入中心室0Tt第32页/共62页吸收室中心室3.口服或肌肉注射相当于药物(剂量D0)先进入吸收室，吸

8、收后进入中心室吸收室药量x0(t)第33页/共62页参数估计参数估计各种给药方式下的各种给药方式下的 c1(t),c2(t)取决于参数取决于参数k12,k21,k13,V1,V2t=0快速静脉注射D0,在ti(i=1,2,n)测得c1(ti)由较大的 用最小二乘法定A A,由较小的 用最小二乘法定B,第34页/共62页参数估计参数估计进入中心室的药物全部排除第35页/共62页 过滤嘴的作用与它的材料和长度有什么关系过滤嘴的作用与它的材料和长度有什么关系 人体吸入的毒物量与哪些因素有关，其中哪人体吸入的毒物量与哪些因素有关，其中哪些因素影响大，哪些因素影响小。些因素影响大，哪些因素影响小。模型模

9、型分析分析 分析吸烟时毒物进入人体的过程，建立吸分析吸烟时毒物进入人体的过程，建立吸烟过程的数学模型。烟过程的数学模型。设想一个设想一个“机器人机器人”在典型环境下吸烟，在典型环境下吸烟，吸烟方式和外部环境认为是不变的。吸烟方式和外部环境认为是不变的。问题问题5.5 香烟过滤嘴的作用香烟过滤嘴的作用第36页/共62页模型模型假设假设定性分析定性分析1）l1烟草长，烟草长，l2过滤嘴长，过滤嘴长，l=l1+l2，毒物量毒物量M均匀分布，密度均匀分布，密度w0=M/l12）点燃处毒物随烟雾进入空气和沿香烟穿）点燃处毒物随烟雾进入空气和沿香烟穿行的数量比是行的数量比是a:a,a+a=13）未点燃的烟

10、草和过滤嘴对随烟雾穿行的）未点燃的烟草和过滤嘴对随烟雾穿行的毒物的毒物的(单位时间单位时间)吸收率分别是吸收率分别是b和和 4）烟雾沿香烟穿行速度是常数）烟雾沿香烟穿行速度是常数v，香烟燃，香烟燃烧速度是常数烧速度是常数u，v uQ 吸一支烟毒物进入人体总量第37页/共62页模模型型建建立立0t=0,x=0，点燃香烟q(x,t)毒物流量w(x,t)毒物密度1)求q(x,0)=q(x)第38页/共62页t时刻，香烟燃至 x=ut1)求q(x,0)=q(x)2)求q(l,t)第39页/共62页3)求w(ut,t)第40页/共62页4)计算 Q第41页/共62页结果结果分析分析烟草为什么有作用?1）

11、Q与a,M成正比，aM是毒物集中在x=l 处的吸入量2)过滤嘴因素，,l2 负指数作用是毒物集中在x=l1 处的吸入量3）(r)烟草的吸收作用b,l1 线性作用第42页/共62页带过滤嘴不带过滤嘴结果结果分析分析4)与另一支不带过滤嘴的香烟比较，w0,b,a,v,l 均相同，吸至 x=l1扔掉提高 -b 与加长l2，效果相同第43页/共62页5.6 人口预测和控制人口预测和控制 年龄分布对于人口预测的重要性 只考虑自然出生与死亡，不计迁移人口人口发展发展方程方程第44页/共62页人口发展方程人口发展方程一阶偏微分方程第45页/共62页人口发展方程人口发展方程已知函数（人口调查）生育率（控制人口

12、手段）0tr第46页/共62页生育率的分解 总和生育率总和生育率h生育模生育模式式0第47页/共62页人口发展方程和生育率人口发展方程和生育率总和生育率总和生育率控制生育的多少控制生育的多少生育模式生育模式控制生育的早晚和疏密控制生育的早晚和疏密 正反馈系统 滞后作用很大第48页/共62页人口指数人口指数1）人口总数2）平均年龄3）平均寿命t时刻出生的人，死亡率按 (r,t)计算的平均存活时间4）老龄化指数控制生育率控制 N(t)不过大控制 (t)不过高第49页/共62页5.7 烟雾的扩散与消失烟雾的扩散与消失现象现象和和问题问题炮弹在空中爆炸，烟雾向四周扩散，形成圆形不透光区域。不透光区域不

13、断扩大，然后区域边界逐渐明亮，区域缩小，最后烟雾消失。建立模型描述烟雾扩散和消失过程，分析消失时间与各因素的关系。问题问题分析分析无穷空间由瞬时点源导致的扩散过程，用二阶偏微分方程描述烟雾浓度的变化。观察的烟雾消失与烟雾对光线的吸收，以及仪器对明暗的灵敏程度有关。第50页/共62页模型模型假设假设1）烟雾在无穷空间扩散，不受大地和风的影响；扩散服从热传导定律。2）光线穿过烟雾时光强的减少与烟雾浓度成正比；无烟雾的大气不影响光强。3）穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏度决定。模型模型建立建立1）烟雾浓度 的变化规律热传导定律：单位时间通过单位法热传导定律：单位时间通过单位法向

14、面积的流量与浓度梯度成正比向面积的流量与浓度梯度成正比 第51页/共62页曲面积分的奥氏公式1）烟雾浓度 的变化规律第52页/共62页 初始条件Q炮弹释放的烟雾总量 单位强度的点源函数 对任意t,C的等值面是球面 x2+y2+z2=R2;RC 仅当 t,对任意点(x,y,z),C01）烟雾浓度 的变化规律第53页/共62页2）穿过烟雾光强的变化规律光强的减少与烟雾浓度成正比第54页/共62页3）仪器灵敏度与烟雾明暗界限烟雾浓度连续变化烟雾中光强连续变化仪器z-设光源在z=-,仪器在z=,则观测到的明暗界限为不透光区域有扩大、缩小、消失的过程穿过烟雾进入仪器的光线只有明暗之分，明暗界限由仪器灵敏

15、度决定。不透光区域边界第55页/共62页4）不透光区域边界的变化规律对任意t,不透光区域边界是圆周不透光区域边界半径第56页/共62页r(t)rm0t1t2t结果分析结果分析观测到不透光区域边界达到最大的时刻t1，可以预报烟雾消失的时刻t2第57页/共62页5.8 万有引力定律的发现万有引力定律的发现背景背景航海业发展天文观测精确“地心说”动摇哥白尼：“日心说”伽里略：落体运动开普勒：行星运动三定律变速运动的计算方法牛顿：一切运动有力学原因牛顿运动三定律牛顿：研究变速运动，发明微积分（流数法）开普勒三定律牛顿运动第二定律万有引力定律自然科学之数学原理(1687)第58页/共62页模型假设模型假设极坐标系(r,)太阳(0,0)1.行星轨道a长半轴,b短半轴,e离心率3.行星运行周期 T行星位置：向径O(太阳)P(行星)r2.单位时间 扫过面积为常数 Am 行星质量 绝对常数4.行星运行受力 第59页/共62页模型建立模型建立O(太阳)P(行星)r向径 的基向量第60页/共62页模型建立模型建立万有引力定律需证明 4A2/p=kM（与哪一颗行星无关）A单位时间 扫过面积O(太阳)P(行星)r(习题)lp/22=pA第61页/共62页感谢您的观看。第62页/共62页