西安电子科技大学数学建模讲义第三讲.pptVIP

数学建模讲义 最小二乘法 稳定性问题 捕食系统的Volterra方程 考察两阶微分方程组： （3.29） 令 ，作一坐标平移，不妨仍用x记x’，则平衡点xo的稳定性讨论转化为原点的稳定性讨论了。将f(x1,x2)、g(x1,x2)在原点展开，（3.29）又可写成: 考察（3.29）的线性近似方程组: （3.30） 其中： 记 λ1、λ2为A的特征值则λ1、λ2是方程： det（A-λI）=λ2- (a+b) λ+ (ad – bc )=0的根 令p=a+d, q=ad-bc=|A|，则 ，记 。 讨论特征值与零点稳定的关系 （1）若△0，可能出现以下情形： ① 若q0，λ1λ20。 当p0时，零点不稳定； 当p0时，零点稳定 若q0，λ1λ20 零点为不稳定的鞍点 ③ q=0，此时λ1=p，λ2=0，零点不稳定。 （2） △=0，则λ1=λ2: 当p0时，零点不 稳定 当p0时，零点稳定 （2） △0，此时 若a0，零点稳定, a 0,零点不稳定 若a=0，有零点为中心的周期解 综上所述：仅当p0且q0时， （3.30）零点才是渐近稳定的；当p=0且q0时（3.30）有周期解，零点是稳定的中心（非渐近稳定）；在其他情况下，零点均为不稳定的。 非线性方程组（3.29）平衡点稳定性讨论可以证明有下面定理成立: 定理2 若（3.30）的零点是渐近稳定的，则（3.29）的平衡点 也是渐近稳定的；若（3.30）的零点是不稳定的，则（3.29） 的平衡点也是不稳定的。 问题背景： 意大利生物学家D’Ancona曾致力于鱼类种群相互制约关系的研究，在研究过程中他无意中发现了一些第一次世界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很理想的鱼类占总渔获量的百分比。在 1914~1923年期间，意大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加： 年代 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 他知道，捕获的各种鱼的比例近似地反映了地中海里各种鱼类的比例。战争期间捕鱼量大幅下降，但捕获量的下降为什么会导致鲨鱼、鳐鱼等食肉鱼比例的上升，即对捕食者有利而不是对食饵有利呢？他百思不得其解，无法解释这一现象，就去求教当时著名的意大利数学家V.Volterra，希望他能建立一个数学模型研究这一问题。 1-1/? i di/dt 0 1 ? 1 i0 i0 0 t i ? 1 1-1/? i0 i 0 t ? ?1 di/dt 0 接触数? =1 ~ 阈值 如果感染期内有效接触感染的健康者人数不超过病人数,患者就会全部治愈。 模型4 传染病有免疫性——病人治愈后即移出感染系统，称移出者 SIR模型 假设 1）总人数N不变，病人、健康人和移出者的比例分别为 2）病人的日接触率? , 日治愈率?, 接触数 ? = ? / ? 建模 需建立 的两个方程 模型4 SIR模型 无法求出 的解析解 在相平面 上 研究解的性质 模型4 消去dt SIR模型 相轨线 的定义域 相轨线 1 1 s i 0 D 在D内作相轨线 的图形，进行分析 s i 1 0 1 D 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/