培训上--微分方程.pptVIP

模型求解 σ1时，t →+∞则 i(t) →1-1/σ。 画出解的图象为 ： σ1，t →+∞时i(t) →0. σ=λ∕μ 1-1/σ i t i0 i0 模型结果分析 i i0 0 t σ1，t →+∞时i(t) →0. σ=λ∕μ 1、假设：这里的假设类似于模型B，只是引入R类人群。分别记s(t)、i(t)、r(t)为病人、易感人群、移出者在总人口中所占的比例。s(t)+ i(t)+ r(t) = 1。另外，日接触率λ，日治愈率μ。 3、SIR模型（免疫模型） 根据假设，模型被修正为 初值条件为i(0) = i0，r(0) = r0，s(0) = s0。 注意：此方程组无法求解析解。 可以求数值解 模型求解 将方程组转化成下面的形式： 其中s≥0，i≥0且s+i≤1。 此方程是可以求解析解 求解得到： 考察随着时间的推移，s(t)、i(t)、r(t)的变化规律。 首先，t →+∞时，分别以s ? , i ? , r?记各自的极限，这些极限都存在。 模型分析 i ? = 0 ?（用反证法） 假设i ? ? 0 ，那么必然有 i ? = ? 0。 根据极限的定义，对于充分大的t，都应该有i(t)ε/2，把这个结论代入方程组。 模型分析 dr/dt=μi με/2 这会导致r(t)→+∞，这跟上面r(t)的极限也存在的结论有矛盾。 所以只能有： i ? = 0 ，即传染病最终将消失。 其次，考虑随着t的变化，i-s平面上解的轨线变化情况。大概的走势图为： 模型分析 i 1 0 1/σ s σ=λ∕μ i 1 0 1/σ s s0 1/?时，i(t) 先升后降至0 传染病蔓延 s0 1/?时，i(t) 单调降至0 传染病不会蔓延开来 1/σ是一个边界点，为了让传染病不蔓延，需要调整s0和1/σ。具体的方法：一是降低s0，如接种疫苗，使S类人群直接变成R类； 二是提高1/σ使之大于s0，σ=λ/μ,也就是降低λ而提高μ，强化卫生教育和隔离病人，同时提高医疗水平。 模型分析 对参数σ的估计： 令解两端同时取t→+∞，因为 i ? = 0 ，得到 参数估计 根据历史数据和此公式就可以得到σ的估计值。 关于传染病模型，我们还可以进一步考虑更复杂的情形，如考虑出生率、死亡率、防疫措施的作用、潜伏期等。 其他类型的传染病模型 SIES模型——健康－染病－潜伏期－健康不免疫 SIER模型——健康－染病－潜伏期－移出系统 SIRS模型——健康－染病－短时免疫－健康(易感) 考虑抵抗能力 考虑地域传播 考虑传播途径(接触、空气、昆虫、水源等) 四、正规战与游击战 战争分类：正规战争，游击战争，混合战争 只考虑双方兵力和战斗力 兵力因战斗及非战斗减员而减少，因增援而增加 战斗力与射击次数及命中率有关 建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例 第一次世界大战Lanchester提出预测战役结局的模型 一般模型 每方战斗减员率取决于双方的兵力和战斗力 每方非战斗减员率与本方兵力成正比 甲乙双方的增援率为u(t), v(t) f, g 取决于战争类型 x(t) ~甲方兵力，y(t) ~乙方兵力 模型假设 模型 正规战争模型 甲方战斗减员率只取决于乙方的兵力和战斗力 双方均以正规部队作战 忽略非战斗减员 假设没有增援 f(x, y)=?ay, a ~ 乙方每个士兵的杀伤率 a=ry py, ry ~射击率， py ~命中率 0 正规战争模型 为判断战争的结局，不求x(t), y(t)而在相平面上讨论 x 与 y 的关系 平方律 模型 乙方胜 游击战争模型 双方都用游击部队作战 甲方战斗减员率还随着甲方兵力的增加而增加 忽略非战斗减员 假设没有增援 f(x, y)=?cxy, c~ 乙方每个士兵的杀伤率 c = ry py ry~射击率 py ~命中率 py=sry /sx sx ~ 甲方活动面积 sry ~ 乙方射击有效面积 0 游击战争模型 线性律 模型 0 混合战争模型 甲方为游击部队，乙方为正规部队 乙方必须10倍于甲方的兵力 设 x0=100, rx/ry=1/2, px=0.1, sx=1(km2), sry=1(m2) 五、减肥计划 减肥有很多办法：加强运动、节食、服用药物、手术、采用专用仪器等。 不过于健康无害的减肥方式却只有两种：运动和节食。 减肥的一个重要原则是不能影响健