[工学]2010数模讲座--微分方程模型.ppt

模型求解: 1、分别用m-文件shier.m和shier.m定义上述两个方程 2、建立主程序shark.m, 求解两个方程，并画出两种情况下鲨鱼数在鱼类总数中所占比例 x2(t)/[x1(t)+x2(t)] 实线为战前的鲨鱼比例，“*”线为战争中的鲨鱼比例 结论：战争中鲨鱼的比例比战前高！ 思考题： 讨论资金积累、国民收入、与人口增长的关系. （1）若国民平均收入x与按人口平均资金积累y成正比，说明仅当总资金积累的相对增长率k大于人口的相对增长率r时，国民平均收入才是增长的. （2）作出k(x)和r(x)的示意图，分析人口激增会引起什么后果. 要求： 1、上网收集有关数据； 2、建立数学模型； 3、模型求解(理论分析、算法设计、结果分析等)； 4、结论和建议。 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者(Removed) SIR模型 假设： 1）总人数N不变，病人、健康人和移出者的比例分别为 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模： 需建立 的两个方程 对移出者 模型4 SIR模型 无法求出 的解析解 在相平面 上 研究解的性质 模型4 消去dt SIR模型 相轨线 的定义域 相轨线 1 1 s i 0 D 在D内作相轨线 的图形，进行分析 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/? = ? 三、捕食系统的Volterra方程 意大利生物学家D’Ancona曾致力于鱼类种群相互制约关系的研究，在研究过程中他无意中发现了一些第一次世界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很理想的鱼类占总渔获量的百分比。在 1914~1923年期间，意大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加. 他知道，捕获的各种鱼的比例近似地反映了地中海里各种鱼类的比例。战争期间捕鱼量大幅下降，但捕获量的下降为什么会导致鲨鱼、鳐鱼等食肉鱼比例的上升，即对捕食者有利而不是对食饵有利呢？他百思不得其解，无法解释这一现象，就去求教当时著名的意大利数学家V.Volterra，希望他能建立一个数学模型研究这一问题。 年代 1914 1915 1916 1917 1918 百分 比 11.9 21.4 22.1 21.2 36.4 年代 1919 1920 1921 1922 1923 百分比 27.3 16.0 15.9 14.8 10.7 Volterra将鱼划分为两类。一类为食用鱼（食饵），数量记为x1(t)，另一类为食肉鱼（捕食者），数量记为x2(t)，并建立双房室系统模型。 1、模型建立 大海中有食用鱼生存的足够资源，可假设食用鱼独立生存将按增长率为r1的指数律增长（Malthus模型），即设： 对于食饵（Prey）系统 ： 由于捕食者的存在，食用鱼数量因而减少，设减少的速率与两者数量的乘积成正比（竞争项的统计筹算律），即： λ1反映了捕食者掠取食饵的能力 对于捕食者（Predator）系统 ： 捕食者设其离开食饵独立存在时的死亡率为r2，即： 但食饵提供了食物，使生命得以延续。这一结果也要通过竞争来实现，再次利用统计筹算律，得到： 综合以上分析，建立P-P模型（Volterra方程） 的方程组： （１） 2、模型分析 方程组（１）是非线性的，不易直接求解。容易看出，该方程组共有两个平衡点，即： Po(0,0)是平凡平衡点且明显是不稳定，没必要研究 方程组还有两组平凡解： 和 和 当x