姜启源数学建模课件第六章稳定性模型幻灯片.ppt

6.5 种群的弱肉强食(食饵-捕食者模型) 种群甲靠丰富的天然资源生存，种群乙靠捕食甲为生，形成食饵-捕食者系统，如食用鱼和鲨鱼，美洲兔和山猫，害虫和益虫。 模型的历史背景——一次世界大战期间地中海渔业的捕捞量下降(食用鱼和鲨鱼同时捕捞)，但是其中鲨鱼的比例却增加，为什么？ 食饵（甲）数量 x(t), 捕食者（乙）数量 y(t) 甲独立生存的增长率 r 乙使甲的增长率减小，减小量与 y成正比 乙独立生存的死亡率 d 甲使乙的死亡率减小，减小量与 x成正比 方程(1),(2) 无解析解 食饵-捕食者模型(Volterra) a ~捕食者掠取食饵能力 b ~食饵供养捕食者能力 Volterra模型的平衡点及其稳定性 平衡点 稳定性分析 P点稳定性不能用近似线性方程分析 p =0, q 0 P: 临界状态 q 0 P′ 不稳定 t x(t) y(t) 0 20.0000 4.0000 0.1000 21.2406 3.9651 0.2000 22.5649 3.9405 0.3000 23.9763 3.9269 … … … 5.1000 9.6162 16.7235 5.2000 9.0173 16.2064 … … … 9.5000 18.4750 4.0447 9.6000 19.6136 3.9968 9.7000 20.8311 3.9587 用数学软件MATLAB求微分方程数值解 x~y 平面上的相轨线 计算结果（数值，图形） x(t), y(t)是周期函数，相图(x,y)是封闭曲线 观察，猜测 x(t), y(t)的周期约为9.6 xmax? 65.5, xmin ? 6, ymax ? 20.5, ymin ? 3.9 用数值积分可算出 x(t), y(t)一周期的平均值： x(t)的平均值约为25, y(t)的平均值约为10。 食饵-捕食者模型(Volterra) 消去dt 用相轨线分析 点稳定性 c 由初始条件确定 取指数 x0 fm f(x) x 0 g(y) gm y0 y 0 在相平面上讨论相轨线的图形 用相轨线分析 点稳定性 相轨线 时无相轨线 以下设 y2 y1 x Q3 Q4 q y1 y2 x1 x2 p y y0 x x0 P 0 x1 x2 Q1 Q2 Q1(x1,y0),Q2(x2,y0) Q3(x,y1), Q4(x,y2) 相轨线退化为P点 存在x1x0x2, 使f(x1)=f(x2)=p 存在y1y0y2,使g(y1)=g(y2)=q 相轨线是封闭曲线族 x Q3 Q4 f(x) x x0 fm 0 g(y) gm y0 y 0 相轨线 P~中心 相轨线是封闭曲线 x(t), y(t)是周期函数(周期记 T) 求x(t), y(t) 在一周期的平均值 轨线中心 用相轨线分析 点稳定性 ? T2 T3 T4 T1 P T1 T2 T3 T4 x(t) 的“相位”领先 y(t) 模型解释 初值 相轨线的方向 模型解释 r ~食饵增长率 d ~捕食者死亡率 b ~食饵供养捕食者能力 捕食者 数量 食饵数量 P r/a d/b a ~捕食者掠取食饵能力 捕食者数量与r成正比, 与a成反比 食饵数量与d成正比, 与b成反比 模型解释 一次大战期间地中海渔业的捕捞量下降，但是其中鲨鱼的比例却在增加，为什么？ r?r-?1, d?d+?1 捕捞 战时捕捞 r?r-?2, d?d+?2 , ?2 ?1 ? ? ? x y 食饵(鱼)减少， 捕食者(鲨鱼)增加 自然环境 还表明：对害虫(食饵)—益虫(捕食者)系统，使用灭两种虫的杀虫剂, 会使害虫增加，益虫减少。 食饵-捕食者模型(Volterra)的缺点与改进 Volterra模型 改写 多数食饵—捕食者系统观察不到周期震荡,而是趋向某个平衡状态,即存在稳定平衡点 加Logistic项 有稳定平衡点 相轨线是封闭曲线，结构不稳定——一旦离开某一条闭轨线，就进入另一条闭轨线，不恢复原状。 自然界存在的周期性平衡生态系统是结构稳定的，即偏离周期轨道后，内部制约使系统恢复原状。 食饵-捕食者模型(Volterra)的缺点与改进 r1=1, N1=20, ?1=0.1, w=0.2, r2=0.5, ?2=0.18 相轨线趋向极限环 结构稳定 两种群模型的几种形式 相互竞争 相互依存 弱肉强食 第六章 稳定性模型