微分方程方法.pptVIP

四 .SARS的传播问题 模型2 1/2 tm i i0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ? t=tm, di/dt 最大 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染. 增加假设 SIS 模型 3）病人每天治愈的比例为? ? ~日治愈率 建模 ? ~ 日接触率 1/? ~感染期 ? ~ 一个感染期内每个病人的有效接触人数，称为接触数. 四 .SARS的传播问题 i0 i0 接触数? =1 ~ 阈值 感染期内有效接触感染的健康者人数不超过病人数 i0 模型2(SI模型)如何看作模型3(SIS模型)的特例 0 t i ? 1 1-1/? i 0 t ? ?1 di/dt 0 四 .SARS的传播问题 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者. SIR模型 假设 1）总人数N不变，病人、健康人和移出者的比例分别为 . 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模 需建立 的两个方程. 四 .SARS的传播问题 无法求出 的解析解 先做数值计算, 再在相平面上研究解析解性质 四 .SARS的传播问题 SIR模型的数值解 i(t)从初值增长到最大; t??, i?0 s(t)单调减; t??, s?0.04 设?=1, ?=0.3, i0=0.02, s0=0.98, 用MATLAB计算作图i(t), s(t)及i(s) s i 相轨线i(s) 四 .SARS的传播问题 模型4 消去dt SIR模型的相轨线分析 相轨线 的定义域 相轨线 四 .SARS的传播问题 在D内作相轨线 的图形，进行分析 相轨线 及其分析 s i 1 0 1 D P1 s0 im P3 P4 P2 S0 s(t)单调减?相轨线的方向 传染病蔓延 传染病不蔓延 P1: s01/?? i(t)先升后降至0 P2: s01/?? i(t)单调降至0 1/?~阈值 四 .SARS的传播问题 ? (日接触率)? ? 卫生水平? ?(日治愈率)? ? 医疗水平? 传染病不蔓延的条件——s01/? ? 的估计 降低 s0 提高 r0 提高阈值 1/? 降低 ?(=?/?) ? ?, ? ? 群体免疫 忽略i0 四 .SARS的传播问题 降低日接触率? 提高日治愈率? 提高移出比例r0 以最终未感染比例s?和病人比例最大值im为度量指标. ? ? 1/? s0 i0 s? i? 1 0.3 0.3 0.98 0.02 0.0398 0.3449 0.6 0.3 0.5 0.98 0.02 0.1965 0.1635 0.5 0.5 1.0 0.98 0.02 0.8122 0.0200 0.4 0.5 1.25 0.98 0.02 0.9172 0.0200 1 0.3 0.3 0.70 0.02 0.0840 0.1685 0.6 0.3 0.5 0.70 0.02 0.3056 0.0518 0.5 0.5 1.0 0.70 0.02 0.6528 0.0200 0.4 0.5 1.25 0.70 0.02 0.6755 0.0200 ? ?, ? ? s0? (r0 ?) s? ?, im ? s? ?, im ? 四 .SARS的传播问题 传染病模型 模型1 模型2 模型3 (SIS) 模型4 (SIR) 区分病人和健康人 考虑治愈 模型3, 4: 描述传播过程, 分析变化规律, 预报高潮时刻, 预防蔓延手段. 模型4: 数值计算与理论分析相结合. 四 .SARS的传播问题 * * 一、微分方程的一般理论 3. 微分方程的稳定性问题 （1）有限区间的稳定性 （2）无限区间的稳定性 （3）渐近稳定性 （4）经常扰动下的稳定性 第三章 微分方程方法 实际中，对于很多问题的微分方程模型并不需要求其一般解，而是需要求其某种理想状态下的解，这种解称为平衡点。 * 二 .微分方程的平衡点及其稳定性 １．平衡点的概念 * 二 .微分方程的平衡点及其稳定性 １．平衡点的概念 问题：如何来断别平衡点的稳定性呢？ * 二 .微分方程的平衡点及其稳定性 １．平衡点的概念 * 二 .微分方程的平衡点及其稳定性 2. 一阶方程的平衡点及稳定性 * 二 .微分方程的平衡点及其稳定性