第五章微分方程模型 5.1传染病模型 5.2经济增长模型 5.3正规战与游击战 54药物在体内的分布与排除 5.5香烟过滤嘴的作用 5.6人口预测和控制 5.7烟雾的扩散与消失 58万有引力定律的发现
第五章 微分方程模型 5.1 传染病模型 5.2 经济增长模型 5.3 正规战与游击战 5.4 药物在体内的分布与排除 5.5 香烟过滤嘴的作用 5.6 人口预测和控制 5.7 烟雾的扩散与消失 5.8 万有引力定律的发现
动态 描述对象特征随时间(空间)的演变过程 模型 分析对象特征的变化规律 预报对象特征的未来性态 研究控制对象特征的手段 微分·根据函数及其变化率之间的关系确定函数 方程·根据建模目的和问题分析作出简化假设 建模,按照内在规律或用类比法建立微分方程
动态 模型 • 描述对象特征随时间(空间)的演变过程 • 分析对象特征的变化规律 • 预报对象特征的未来性态 • 研究控制对象特征的手段 微分 • 根据函数及其变化率之间的关系确定函数 方程 建模 • 根据建模目的和问题分析作出简化假设 • 按照内在规律或用类比法建立微分方程
5.1传染病模型 问题·描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律, 用机理分析方法建立模型
5.1 传染病模型 问题 • 描述传染病的传播过程 • 分析受感染人数的变化规律 • 预报传染病高潮到来的时刻 • 预防传染病蔓延的手段 • 按照传播过程的一般规律， 用机理分析方法建立模型
模型1已感染人数(病人)i(0 假设·每个病人每天有效接触 (足以使人致病)数为x 建模(t+△t)-()=(t)△t dt ni d i(t)=ioe (O)=nt→→1→? 若有效接触的是病人, 必须区分已感染者(病 则不能使病人数增加 人和未感染者(健康人
已感染人数 (病人) i(t) • 每个病人每天有效接触 (足以使人致病)人数为 模型1 假设 i(t + t) − i(t) = i(t)t 若有效接触的是病人， 则不能使病人数增加 必须区分已感染者(病 人)和未感染者(健康人) 建模 0 i(0) i i dt di = = t → i → t i t i e 0 ( ) = ?
模型2区分已感染者病人)和未感染者(健康人) 假设1)总人数N不变,病人和健康 人的比例分别为i(t),s(t) SI模型 2)每个病人每天有效接触人数~日 为,且使接触的健康人致病接触率 建模N(t+△)-i()=[s(t)N(t)At = asi =Ai(1-i) S(t)+i(t)=1 i(o)=i
si dt di = s(t) + i(t) =1 模型2 区分已感染者(病人)和未感染者(健康人) 假设 1）总人数N不变，病人和健康 人的 比例分别为 i(t),s(t) 2）每个病人每天有效接触人数 为, 且使接触的健康人致病 建模 N[i(t + t) −i(t)] = [s(t)]Ni(t)t = = − 0 (0) (1 ) i i i i dt di ~ 日 接触率 SI 模型
模型2 1(1-1) LogistIc模型 (O)= i(t) 1+ =tn,lid最大 m~传染病高潮到来时刻t→>→i→1? (日接触率八→tn个病人可以治愈!
t e i i t − + − = 1 1 1 1 ( ) 0 = = − 0 (0) (1 ) i i i i dt di 模型2 1/2 tm i i0 1 0 t = − − 1 1 ln 0 1 i t m tm~传染病高潮到来时刻 (日接触率) → tm t → i →1 Logistic 模型 病人可以治愈！ ? t=tm, di/dt 最大
模型3传染病无免疫性—病人治愈成 为健康人,健康人可再次被感染SS模型 增加假设3)病人每天治愈的比例为~日治愈率 建模M[i(t+△)-i()]=ANs(0)i(t)△t-NVi(t)△t di_ni( -l-llL 元~日接触率 dt i(o)=i 1/~感染期 σ=/a~一个感染期内每个病人的 有效接触人数,称为接触数
模型3 传染病无免疫性——病人治愈成 为健康人，健康人可再次被感染 增加假设 SIS 模型 3）病人每天治愈的比例为 ~日治愈率 建模 N[i(t + t) − i(t)] = Ns(t)i(t)t − Ni(t)t = / ~ 日接触率 1/ ~感染期 ~ 一个感染期内每个病人的 有效接触人数，称为接触数。 = = − − 0 (0) (1 ) i i i i i dt di
模型3 =i(1-1)-1na=M/ai didt ni-(1--) o>1 小→()按形曲线增长感染期内有效接触感染的 健康者人数不超过病人数 模型2(S模型)如何看作模型3(SIS模型)的特例
− = 0, 1 , 1 1 1 ( ) i )] 1 [ (1 = −i i − − dt 模型 di 3 i0 i0 接触数 =1 ~ 阈值 = / 1 i(t) i(t)按S形曲线增长 感染期内有效接触感染的 0小 健康者人数不超过病人数 1 i 1-1/ i0 i i i dt di = (1− ) − 模型2(SI模型)如何看作模型3(SIS模型)的特例 i di/dt 0 1 >1 0 t i >1 1-1/ i 0 t 1 di/dt < 0
模型4传染病有免疫性病人治愈 后即移出感染系统,称移出者 SIR模型 假设1)总人数N不变,病人、健康人和移 出者的比例分别为i(t),S(t),r(t) 2)病人的日接触率,日治愈率 接触数σ=4/ 建模S()+i(1)+r(t)=1 需建立i(t),S(t),r(t)的两个方程
模型4 传染病有免疫性——病人治愈 后即移出感染系统，称移出者 SIR模型 假设 1）总人数N不变，病人、健康人和移 出者的比例分别为 i(t),s(t),r(t) 2）病人的日接触率 , 日治愈率, 接触数 = / 建模 s(t) + i(t) + r(t) =1 需建立 i(t),s(t),r(t) 的两个方程
模型4 SIR模型 M(t+△)-i(t)=N(t)i(t)△t-1Ni(t)△t Ms(t+△)-S()=-AN(0)()△t ddd asi- ui 无法求出i(t),S(t) asi 的解析解 i(0)=0,s(0)=So 在相平面S~i上 研究解的性质 +Sn≈1(通常r(O)=r很小)
N[i(t + t) − i(t)] = Ns(t)i(t)t − Ni(t)t 模型4 SIR模型 i 0 + s0 1(通常r(0) = r0 很小） 无法求出 的解析解 i(t),s(t) 在相平面 上 研究解的性质 s ~ i N[s(t + t) − s(t)] = −Ns(t)i(t)t = = = − = − 0 0 i(0) i , s(0) s si dt ds si i dt di