微分方程模型3.ppt

食饵-捕食者模型 Ancona 无法解释这个现象，于是求助于著名的意大利数学家V.Volterra，希望建立一个食饵—捕食系统的数学模型，定量地回答这个问题. 模型建立 我们得到 针对一组具体的数据用Matlab软件进行计算. 设食饵和捕食者的初始数量分别为25和2，r=1,d=0.5,a=0.1,b=0.02 MATLAB 程序 r=1;d=0.5;a=0.1;b=0.02; volte=@(t,x)[r*x(1)-a*x(1)*x(2), -*x(2)+b*x(1)*x(2)]; ts=[0 50];x0=[25 2]; [t,x]=ode45(volte,ts,x0); plot(t,x(:,1),r,t,x(:,2),b); figure; plot(x(:,1),x(:,2)); 图形结果 结构稳定的生态系统 MATLAB程序 r1=1;n1=30;sigma1=0.1;w=0.2;r2=0.5;sigma2=0.3; fun=@(t,x)[r1*x(1)*(1-x(1)/n1-sigma1*x(2)/(1+w*x(1))); r2*x(2)*(-1+sigma2*x(1)/(1+w*x(1)))]; [t,x]=ode45(fun,[0,100],[15,10]); plot(t,x(:,1),r);hold on; plot(t,x(:,2),b); figure; plot(x(:,1),x(:,2)); 稳定性问题 意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究，他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中，发现鲨鱼等的比例有明显增加（见下表），而供其捕食的食用鱼的百分比却明显下降.显然战争使捕鱼量下降，食用鱼增加，鲨鱼等也随之增加，但为何鲨鱼的比例大幅增加呢？ Volterra将鱼划分为两类。一类为食用鱼（食饵），数量记为x(t)，另一类为食肉鱼（捕食者），数量记为y(t)，建立双房室系统模型。 对于食饵（Prey）系统 ： 大海中有食用鱼生存的足够资源，可假设食用鱼独立生存将按增长率为r的指数律增长 由于捕食者的存在，食用鱼数量因而减少，设减少的速率与两者数量的乘积成正比 a反映捕食者的能力 对于捕食者（Predator）系统 ： 捕食者离开食饵独立存在时的死亡率为d 但食饵提供了食物，使捕食者增长，设增长作用与食饵数量 成正比 b反映食饵对捕食者的供养能力 方程组无解析解 猜测 x(t), y(t)是周期函数，相图(x,y)是封闭曲线 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项 有稳定平衡点 相轨线是封闭曲线，结构不稳定——一旦离开某一条闭轨线，就进入另一条闭轨线，不恢复原状。 自然界存在的周期性平衡生态系统是结构稳定的，即偏离周