传染病模型优秀课件.ppt

传染病模型第1页，本讲稿共39页1 1 微分方程简介微分方程简介 第2页，本讲稿共39页第3页，本讲稿共39页第4页，本讲稿共39页第5页，本讲稿共39页第6页，本讲稿共39页第7页，本讲稿共39页第8页，本讲稿共39页第9页，本讲稿共39页第10页，本讲稿共39页第11页，本讲稿共39页第12页，本讲稿共39页第13页，本讲稿共39页2 2 传染病模型传染病模型问题问题 描述传染病的传播过程描述传染病的传播过程 分析受感染人数的变化规律分析受感染人数的变化规律 预报传染病高潮到来的时刻预报传染病高潮到来的时刻 预防传染病蔓延的手段预防传染病蔓延的手段 按照传播过程的一般规律，按照传播过程的一般规律，用机理分析方法建立模型用机理分析方法建立模型第14页，本讲稿共39页 已感染人数已感染人数(病人病人)i(t)每个病人每天有效接触每个病人每天有效接触(足以使人致病足以使人致病)人数为人数为 模型模型1 1假设假设若有效接触的是病人，若有效接触的是病人，则不能使病人数增加则不能使病人数增加必须区分已感染者必须区分已感染者(病病人人)和未感染者和未感染者(健康人健康人)建模建模?第15页，本讲稿共39页模型模型2 2区分已感染者区分已感染者(病人病人)和未感染者和未感染者(健康人健康人)假设假设1）总人数）总人数N不变，病人和健康不变，病人和健康 人的人的 比例分别为比例分别为 2）每个病人每天有效接触人数）每个病人每天有效接触人数为为,且且使接触的健康人致病使接触的健康人致病建模建模 日日接触率接触率SI 模型模型第16页，本讲稿共39页模型模型21/2tmii010ttm传染病高潮到来时刻传染病高潮到来时刻 (日接触率日接触率)tm Logistic 模型病人可以治愈！病人可以治愈！?t=tm,di/dt 最最大大第17页，本讲稿共39页模型模型3传染病无免疫性传染病无免疫性病人治愈成病人治愈成为健康人，健康人可再次被感染为健康人，健康人可再次被感染增加假设增加假设SIS 模型模型3）病人每天治愈的比例为）病人每天治愈的比例为 日日治愈率治愈率建模建模 日接触率日接触率1/感染期感染期 一个感染期内一个感染期内每个病人的每个病人的有效接触人数，称为有效接触人数，称为接触数接触数。第18页，本讲稿共39页模型模型3i0i0接触数接触数 =1 阈值阈值感染期内感染期内有效接触感染的健康者有效接触感染的健康者人数不超过病人数人数不超过病人数1-1/i0模型模型2(SI模型模型)如何看作模型如何看作模型3(SIS模型模型)的特例的特例idi/dt01 10ti 11-1/i0t 1di/dt 0第19页，本讲稿共39页模型模型4传染病有免疫性传染病有免疫性病人治愈病人治愈后即移出感染系统，称后即移出感染系统，称移出者移出者SIR模型模型假设假设1）总人数）总人数N不变，病人、健康人和移不变，病人、健康人和移出者的比例分别为出者的比例分别为2）病人的日接触率）病人的日接触率 ,日日治愈率治愈率,接触数接触数 =/建模建模需建立需建立 的两个方程的两个方程第20页，本讲稿共39页模型模型4SIR模型模型无法求出无法求出 的解析解的解析解在相平面在相平面 上上研究解的性质研究解的性质第21页，本讲稿共39页模型模型4消去消去dtSIR模型模型相轨线相轨线 的定义域的定义域相轨线相轨线11si0D在在D内作相轨线内作相轨线 的图形，进行分析的图形，进行分析第22页，本讲稿共39页si101D模型模型4SIR模型模型相轨线相轨线 及其分析及其分析传染病蔓延传染病蔓延传染病不蔓延传染病不蔓延s(t)单调减单调减相轨线的方向相轨线的方向P1s0imP1:s01/i(t)先升后降至先升后降至0P2:s01/i(t)单调降至单调降至01/阈阈值值P3P4P2S0第23页，本讲稿共39页模型模型4SIR模型模型预防传染病蔓延的手段预防传染病蔓延的手段 (日接触率日接触率)卫生水平卫生水平 (日日治愈率治愈率)医疗水平医疗水平 传染病不蔓延的条件传染病不蔓延的条件s01/的估计的估计 降低降低 s0提高提高 r0 提高阈值提高阈值 1/降低降低 (=/),群体免疫群体免疫第24页，本讲稿共39页模型模型4SIR模型模型被传染人数的估计被传染人数的估计记被传染人数比例记被传染人数比例xs0i0P1i0 0,s0 1 小小,s0 1提高阈值提高阈值1/降低降低被被传染人数比例传染人数比例 xs0-1/=第25页，本讲稿共39页战争分类：正规战争，游击战争，混合战争战争分类：正规战争，游击战争，混合战争只考虑双方兵力多少和战斗力强弱只考虑双方兵力多少和战斗力强弱兵力因战斗及非战斗减员而减少，因增援而增加兵力因战斗及非战斗减员而减少，因增援而增加战斗力与射击次数及命中率有关战斗力与射击次数及命中率有关建模思路和方法为用数学模型讨论社会建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例领域的实际问题提供了可借鉴的示例第一次世界大战第一次世界大战Lanchester提出预测战役结局的模型提出预测战役结局的模型3 3 战争模型战争模型 第26页，本讲稿共39页一般模型一般模型 每方战斗减员率取决于双方的兵力和战斗力每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为甲乙双方的增援率为u(t),v(t)f,g 取决于战争类型取决于战争类型x(t)甲方兵力，甲方兵力，y(t)乙方兵力乙方兵力模型模型假设假设模型模型第27页，本讲稿共39页正规战争模型正规战争模型 甲方战斗减员率只取决于乙方的兵力和战斗力甲方战斗减员率只取决于乙方的兵力和战斗力双方均以正规部队作战双方均以正规部队作战 忽略非战斗减员忽略非战斗减员 假设没有增援假设没有增援f(x,y)=ay,a 乙方每个士兵的杀伤乙方每个士兵的杀伤率率a=ry py,ry 射击率，射击率，py 命中命中率率第28页，本讲稿共39页0正规战争模型正规战争模型为判断战争的结局，不求为判断战争的结局，不求x(t),y(t)而在相平面上讨论而在相平面上讨论 x 与与 y 的关系的关系平方律平方律 模型模型乙方胜乙方胜第29页，本讲稿共39页游击战争模型游击战争模型双方都用游击部队作战双方都用游击部队作战 甲方战斗减员率还随着甲方兵力的增加而增加甲方战斗减员率还随着甲方兵力的增加而增加 忽略非战斗减员忽略非战斗减员 假设没有增援假设没有增援f(x,y)=cxy,c 乙方每个士兵的杀伤率乙方每个士兵的杀伤率c=ry pyry射击率射击率py 命中率命中率py=sry/sxsx 甲方活动面积甲方活动面积sry 乙方射击有效面积乙方射击有效面积第30页，本讲稿共39页0游击战争模型游击战争模型线线性性律律 模模型型第31页，本讲稿共39页0混合战争模型混合战争模型甲方为游击部队，乙方为正规部队甲方为游击部队，乙方为正规部队乙方必须乙方必须10倍于甲方的兵力倍于甲方的兵力设设 x0=100,rx/ry=1/2,px=0.1,sx=1(km2),sry=1(m2)第32页，本讲稿共39页 再生资源（渔业、林业等）与再生资源（渔业、林业等）与非再生资源（矿业等）非再生资源（矿业等）再生资源应适度开发再生资源应适度开发在持续稳在持续稳产前提下实现最大产量或最佳效益。产前提下实现最大产量或最佳效益。问题问题及及 分析分析 在在捕捞量稳定捕捞量稳定的条件下，如何控的条件下，如何控制捕捞使产量最大或效益最佳。制捕捞使产量最大或效益最佳。如果使捕捞量等于自然增长量，如果使捕捞量等于自然增长量，渔渔场鱼量将保持不变场鱼量将保持不变，则捕捞量稳定。，则捕捞量稳定。背景背景4 最优捕鱼问题 第33页，本讲稿共39页产量模型产量模型假设假设 无捕捞时鱼的自然增长服从无捕捞时鱼的自然增长服从 Logistic规律规律 单位时间捕捞量与渔场鱼量成正比单位时间捕捞量与渔场鱼量成正比建模建模 捕捞情况下捕捞情况下渔场鱼量满足渔场鱼量满足 不需要求解不需要求解x(t),只需知道只需知道x(t)稳定的条件稳定的条件r固有增长率固有增长率,N最大鱼量最大鱼量h(x)=Ex,E捕捞强度捕捞强度x(t)渔场鱼量渔场鱼量第34页，本讲稿共39页一阶微分方程的平衡点及其稳定性一阶微分方程的平衡点及其稳定性一阶非线性（自治）方程一阶非线性（自治）方程F(x)=0的根的根x0 微分方程的微分方程的平衡点平衡点设设x(t)是方程的解，若从是方程的解，若从x0 某邻域的任一初值出发，某邻域的任一初值出发，都有都有称称x0是方程是方程(1)的的稳定平衡点稳定平衡点不求不求x(t),判断判断x0稳定性的方法稳定性的方法直接法直接法(1)的近似线性方程的近似线性方程第35页，本讲稿共39页产量模型产量模型平衡点平衡点稳定性判断稳定性判断x0 稳定稳定,可得到稳定产量可得到稳定产量x1 稳定稳定,渔场干枯渔场干枯E捕捞强度捕捞强度r固有增长率固有增长率第36页，本讲稿共39页产量模型产量模型在捕捞量稳定的条件下，控制在捕捞量稳定的条件下，控制捕捞强度使产量最大捕捞强度使产量最大图解法图解法P的横坐标的横坐标 x0平衡点平衡点P的纵坐标的纵坐标 h产量产量产量产量最大最大f 与与h交点交点P控制渔场鱼量为最大鱼量的一半控制渔场鱼量为最大鱼量的一半y=rxhPx0hmx0*=N/2P*y=E*xy0y=h(x)=ExxNy=f(x)第37页，本讲稿共39页效益模型效益模型假设假设 鱼销售价格鱼销售价格p 单位捕捞强度费用单位捕捞强度费用c 单位时间利润单位时间利润在捕捞量稳定的条件下，控制捕捞在捕捞量稳定的条件下，控制捕捞强度使效益最大强度使效益最大.稳定平衡点稳定平衡点求求E使使R(E)最大最大渔场渔场鱼量鱼量收入收入 T=ph(x)=pEx支出支出 S=cE第38页，本讲稿共39页EsS(E)T(E)0rE捕捞捕捞过度过度 封闭式捕捞封闭式捕捞追求利润追求利润R(E)最大最大 开放式捕捞开放式捕捞只求利润只求利润R(E)0R(E)=0时的捕捞强度时的捕捞强度(临界强度临界强度)Es=2ER临界强度下的渔场鱼量临界强度下的渔场鱼量捕捞过度捕捞过度ERE*令令=0第39页，本讲稿共39页