食饵—捕食者模型.docVIP

精心整理 楚雄师范学院数学系《数学模型》课 程 食饵—捕食者模型 议论拥有自己阻滞作用的两种群食饵-捕食者模型，第一依据该两种群的相互关系建立模型，讲解参数的意义，尔后进行牢固性解析，讲解平衡点牢固的实质意义，对模型进行相轨线解析来验 证理论解析的正确性，并用matlab软件画出图形。 自然界中不同样种群之间还存在着一种特别幽默的既有相互依存、又有相互限制的生活方式：种群甲靠丰富的天然资源生长，而种群乙靠捕食甲为生，形成鱼和鲨鱼，美洲兔和山猫，落叶松和蚜虫等等都是这类生计方式的典型，生态学称种群甲为食饵，种群乙为捕食者。两者共同组成食饵—捕食者系统。 一食饵—捕食者 采用食饵(食用鱼)和捕食者(鲨鱼)为研究对象,设x(t)/x1(t)为食饵(食用鱼)在时刻t的数目， y(t)/x2(t)为捕食者(鲨鱼)在时刻t的数目，r1为食饵(食用鱼)的相对增加率，r2为捕食者(鲨鱼) 的相对增加率；N1为大海中能容纳的食饵(食用鱼)的最大容量，N2为大海中能容纳的捕食者(鲨 鱼)的最大容量，1为单位数目捕食者（相对于N2）供给的供养食饵的实物量为单位数目捕食者(相 对于N1)耗资的供养甲实物量的1倍；2为单位数目食饵（相对于N1）供给的供养捕食者的实物 量为单位数目捕食者(相对于N2)耗资的供养食饵实物量的2倍；d为捕食者走开食饵独立生计时 的死亡率 二模型假设 假设捕食者（鲨鱼）走开食饵没法生计； 假设大海中资源丰富，食饵独立生计时以指数规律增加； 三模型建立 精心整理 食饵(食用鱼)独立生计时以指数规律增加，且食饵(食用鱼)的相对增加率为r1，即xrx，而 捕食者的存在使食饵的增加率减小，设减小的程度与捕食者数目成正比，于是x(t)知足方程 x(t)x(ray)rxaxy(1) 比率系数a反响捕食者打劫食饵的能力。 因为捕食者走开食饵没法生计，且它独立生计时死亡率为d，即ydy，而食饵的存在为捕 食者供给了食品，相当于使捕食者的死亡率降低，且促使其增加。设这类作用与食饵数目成正比， 于是y(t)知足 y(t)y(dbx)dybxy(2) 比率系数b反响食饵对捕食者的供养能力。 方程(1)、(2)是在自然环境中食饵和捕食者之间依存和限制的关系，这里没有考虑种群自己的 阻滞作用，是Volterra提出的最简单的模型。结果以下。 不考虑自己阻滞作用 令x(0)=x0,y(0)=0, 求解以下 1）先建立M文件  :数值解 设r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2  使用  Matlab  求解 functionxdot=shier(t,x) r=1;d=0.5;a=0.1;b=0.02; xdot=[(r-a*x(2)).*x(1);(-d+b*x(1)).*x(2)]; 在命令窗口输入以下命令： ts=0:0.1:15; x0=[25,2]; [t,x]=ode45(shier,ts,x0);[t,x], ts=0:0.1:15; x0=[25,2]; [t,x]=ode45(shier,ts,x0);[t,x], ans= 精心整理 省略 plot(t,x),grid,gtext(x(t)),gtext(y(t)), pause plot(x(:,1),x(:,2)),grid, （可以猜想，x(t),y(t)是周期函数,与此相应地相轨线y(x)封闭曲线，从数值解近似定出周期为10.7，x的最大最小值分别为99.3，2.0，y的最大，最小值分别为28.4和2.0，简单算出x(t),y(t)再一个周期的平均值为25，10.） 考虑阻滞作用 前面我们没有考虑种群自己的阻滞作用，接下来我们考虑种群自己的阻滞作用，在上边（1），（2）两式中加入Logistic项，即建立以下数学模型： x(t) rx1 x1 x2 (3) 1 11 N 1N 2 1 x2 (t) r2x2 1 x1 x2(4) 2N N 2 1 四平衡点进行理论解析 下边对（3）（4）进行平衡点牢固性解析： 由微分方程(3)、(4) 令f(x1,x2)=0,g(x1,x2)=0 获得以下平衡点: P1(N1,0),P2( N1( 1 1) , N2( 2 1) ),P3(0,0) 1 1 2 1 1 2 因为仅当平衡点位于平面坐标系的第一象限时 (x1,x20)才有意义，因此，对P2而言要求20。 依据判断平衡点牢固性的方法计算： 依据p等于主对角线元素之和的相反数，而 q为其行列式的值，我们获得下表： 平衡点 牢固条件 21 1 精心整理 不牢固 五模型解析与查验 平衡点牢固性的解析及其实质意义： 1)对P1(N1,0)而言，有p=r1r2(21)，q=r1r2(21)，故当21时，平衡点P1(N1,0)是 牢固的。 意义：若是P(N,0