传染病模型(微分方程).docVIP

传染病模型(微分方程) 传染病模型(微分方程) PAGE / NUMPAGES 传染病模型(微分方程) . 微分方程建模（传染病模型）的求解。 1、模型 1： SI 模型。 假设： （1） t 时刻人群分为易感者（占总人数比例的 s(t) ）和已感染者（占总人数比例的 y(t ) ） （2）每个病人每天有效接触的平均人数是常数 ， 称为日接触率，当健康者与病人接触 时，健康者受感染成为病人。 分析：根据假设，每个患者每天可以使 s(t) 个健康者变为病人，因为病人数为 Ny (t ) ，所 以每天共有 Ns(t ) y(t) 个健康者变为病人。即： N dy Nsy ，且 s(t ) y(t ) 1，设初始时刻病人比例为 b ，则： dt dy y(1 y) ，用 MATLAB解此微分方程： dt y(0) b syms a b f=dsolve(Dy=a*y*(1-y),y(0)=b,t) f = 1/(1-exp(-a*t)*(-1+b)/b) %y(t) 1 1 b 1 e t 1 ( 1 1)e t 1 b b 当 b 0.09, 0.1 时 ， 分 别 在 坐 标 系 oty 中 作 出 y(t) 的 图 像 ， 坐 标 系 oyy 中 作 出 y y(1 y) 的图像， 1/(1+91/9 exp(-1/10 t)) a=0.1; 1 0.9