[理学]微分方程模型.pptVIP

第一讲 微分方程模型 一、人口增长模型（种群增长模型） 二、传染病模型 三、捕食系统的 Volterra 模型 四、用Matlab软件求常微分方程的数值解 三、捕食系统的Volterra方程 3. 能否对Logistic模型做进一步修改？如何修改？ 二、传染病模型 问题： 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律，用机理分析方法建立模型 已感染人数 (infective) i(t) 每个病人每天有效接触(足以使人致病)人数为? 模型1 假设： 若有效接触的是病人，则不能使病人数增加 必须区分已感染者(病人)和未感染者(健康人) 建模： ? 模型2 区分已感染者(infective)和易感染者(susceptible) 假设： 1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为?, 且使接触的健康人致病 建模： ? ~ 日 接触率 SI 模型 病人的数的增加率 模型2 1/2 tm i i0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ? t=tm, di/dt 最大 卫生水平 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染 增加假设： SIS 模型 3）病人每天治愈的比例为? ? ~日治愈率 建模： ? ~ 日接触率 1/? ~感染期 ? ~ 一个感染期内每个病人的有效接触人数，称为接触数。 模型3 i0 i0 接触数? =1 ~ 阈值 感染期内有效接触感染的健康者人数不超过病人数 1-1/? i0 模型2(SI模型)如何看作模型3(SIS模型)的特例 i di/dt 0 1 ? 1 0 t i ? 1 1-1/? i 0 t ? ?1 di/dt 0 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者(Removed) SIR模型 假设： 1）总人数N不变，病人、健康人和移出者的比例分别为 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模： 需建立 的两个方程 对移出者 模型4 SIR模型 无法求出 的解析解 在相平面 上 研究解的性质 设 作数值计算，i(t)~s(t) 和i~s的图形如下： i(t) s(t) i~s function ill_ex ts=0:50; x0=[0.02,0.98]; [t,x]=ode45(ill,ts,x0);// [t,x]=ode45(@ill,ts,x0); [t,x] figure(1) plot(t,x(:,1),t,x(:,2)),grid figure(2) plot(x(:,2),x(:,1)),grid function y=ill(t,x) a=1; b=0.3; y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]; 模型4 消去dt SIR模型 相轨线 的定义域 相轨线 1 1 s i 0 D 在D内作相轨线 的图形，进行分析 在相平面i~s上研究解的性质 s i 1 0 1 D 模型4 SIR模型 相轨线 及其分析 传染病蔓延 传染病不蔓延 s(t)单调减?相轨线的方向 P1 s0 im P1: s01/ ? i(t)先升后降至0 P2: s01/ ? i(t)单调降至0 1/ ~阈值 P3 P4 P2 S0 模型4 SIR模型 预防传染病蔓延的手段 ? (日接触率)? ? 卫生水平? ?(日治愈率)? ? 医疗水平? 传染病不蔓延的条件——s01/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 群体免疫 模型4 SIR模型 被传染人数的估计 记被传染人数比例 xs0 i 0 P1 ? i0 ?0, s0 ?1 ? 小, s0 ? ?1 提高阈值1/? ?降低被传染人数比例 x s0 - 1/? = ? 意大利生物学家D’Ancona曾致力于鱼类种群相互制约关系的研究，在研究过程中他无意中发现了一些第一次世界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很理想的鱼类占总渔获量的百分比。在 1914~1923年期间，意大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加. 他知道，捕获的各种鱼的比例近似地反映了