I*者被捕食者模型稳定捕性分析Prepared on 22 November 2020被捕食者—捕食者模型稳定性分析【摘要】自然界中不同种群之间还存在着一种非常有趣的既有相互依存、又 有相互制约的生活方式：种群甲靠丰富的天然资源生存，种群乙靠捕食甲为 生，形成食饵-捕食者系统，如食用鱼和鲨鱼，美洲兔和山猫，害虫和益虫等 本文是基于食饵—捕食者之间的有关规律，建立具有自身阻滞作用的两种群食 饵—捕食者模型，分析平衡点的稳定性，进行相轨线分析，并用数值模拟方法 验证理论分析的正确性关键词】食饵—捕食者模型 相轨线 平衡点 稳定性一、问题重述 在自然界中，存在这种食饵—捕食者关系模型的物种很多下面讨论具有 自身阻滞作用的两种群食饵-捕食者模型，首先根据该两种群的相互关系建立模 型，解释参数的意义，然后进行稳定性分析，解释平衡点稳定的实际意义，对 模型进行相轨线分析来验证理论分析的正确性二、问题分析本文选择渔场中的食饵（食用鱼）和捕食者（鲨鱼）为研究对象，建立微分方程，并利用数学软件MATLAB求出微分方程的数值解，通过对数值结果和图形的 观察，猜测出它的解析解构造然后，从理论上研究其平衡点及相轨线的形 状，验证前面的猜测。
三、模型假设1.假设捕食者（鲨鱼）离开食饵无法生存； 2.假设大海中资源丰富，食饵独立生存时以指数规律增长；四、符号说明x（t） / x （t） 食饵（食用鱼）在时刻t的数量；1y（t） / x （t）——捕食者（鲨鱼）在时刻t的数量；2r ——食饵（食用鱼）的相对增长率；1r ——捕食者（鲨鱼）的相对增长率；2N ——大海中能容纳的食饵（食用鱼）的最大容量 1N ——大海中能容纳的捕食者（鲨鱼）的罪的容量2单位数量捕食者（相对于 N ）提供的供养食饵的实物量为单位数量捕食2者（相对于N ）消耗的供养甲实物量的b倍;11b ——单位数量食饵（相对于 N ）提供的供养捕食者的实物量为单位数量捕 21食者（相对于N ）消耗的供养食饵实物量的b倍；22d ——捕食者离开食饵独立生存时的死亡率五、模型建立食饵独立生存时以指数规律增长，且食饵（食用鱼）的相对增长率为r,即1x'二rx，而捕食者的存在使食饵的增长率减小，设减小的程度与捕食者数量成(1)正比，于是x（t）满足方程x '(t) = x(r - ay) = rx - axy比例系数a反映捕食者掠取食饵的能力由于捕食者离开食饵无法生存，且它独立生存时死亡率为 d ，即y' = -dy，而食饵的存在为捕食者提供了食物，相当于使捕食者的死亡率降(2)低，且促使其增长。
设这种作用与食饵数量成正比，于是y (t)满足y '(t) = y (-d + bx) = -dy + bxy比例系数b反映食饵对捕食者的供养能力方程(1)、(2)是在自然环境中食饵和捕食者之间依存和制约的关系，这里没有考虑种群自身的阻滞作用，是Vol terra提出的最简单的模型下面，我们加入种群自身的阻滞作用，在上两式中加入Logistic项，即建 立以下数学模型：h X1 x2 I N1 1N 2(3)%2(t) = 丫2 x2—1+b(4)六、模型求解在此，我们采用MATLAB软件求解此微分方程组中的x (t)、x (t)的图形及12相轨线图形设 b = 1.5，b = 4，r = 1，r = 0.4，N = 3500, N = 500，使用 MATLAB1 2 1 2 1 2软件求解，程序代码如下：1)建立M文件function y=fun(t，x)y=[x(1).*(1-x(1)./*x(2)./500)，.*x(2).*(-1+4.*x(1)./3500-x(2)./500)]';2)在命令窗口输入如下命令：[t,x]=ode45('fun1',[0,40],[2000,35])得到数值解如下：t(x(l),x(2)) +003 * (单位：千克)0>> plot(t,x),grid,gtext('x(t)'),gtext('y(t)')图1.数值解x (t), x (t)的图形12>> plot(x(:,1),x(:,2)),grid,图 2. 相轨线图形从数值解及x (t) , x (t)的图形可以看出他们的数量变化情况，随着时间的 12推移，都趋于一个稳定的值，从数值解中可以近似的得到稳定值为：(1250,214)。
下面对其平衡点进行稳定性分析：由微分方程(3)、(4)得到如下平衡点:P (N ,0) ,11P(P3(0,0)因为仅当平衡点位于平面坐标系的第一象限时（x ,x > 0）才有意义，所12以，对P而言要求b >022按照判断平衡点稳定性的方法计算：根据p等于主对角线元素之和的相反数，而q为其行列式的值，我们得到下表：平衡点稳定条件b <1b >1不稳定七、模型分析与检验1.平衡点稳定性的分析及其实际意义：q = 一 r r (b -1)，故当 b <11 2 2 2没有种群的共存1)对 P (N ,0)而言，有 p = r 一 r (b 一 1),1 1 1 2 2时，平衡点件（N 1，0）是稳定的意义：如果 P1 （N1 ,0）稳定，则种群乙灭绝2)对P (+ bi)2 1 +b b12N2 °2一-）而言，有p1 + b b12r (1 + b ) + r (b 一1) =-1 + 2 2—1 +b b12r r (1 + b )(b 一1) q = -^2 + 2—1 + b b12故当b >1时，平衡点p（N^,2 ?）是稳定的会互相依意义：如果P ( N 1(1 +b 1) , [ 2(b 2 一 D)稳定，则两物种恒稳发展,2 1 +b b 1 +b b1 2 1 2存生长下去。
3)对 P3(0,0)而言，由于 p = -r + r， q = -rr ，又有题知 r >0, r >0,故3 1 2 1 2 1 2q <0,即P (N ,0)是不稳定的112.平衡点的检验：对于平衡点P (字也,N①—％，把前面给出的初始值带入，在这使 2 1 + Q Q 1 + Q Q1 2 1 2用 MATLAB 软件进行简单的求解，在命令窗口输入如下代码： >> x(1)=(3500.*(1+)./(1+.*4);>> x(2)=(500.*(4-1))./(1+.*4);>> [x(1);x(2)]ans =+003 *把此处求解出的解和前面得出的数值解进行比较可知，平衡点1 + Q Q12p(,工岂)是稳定的1 + Q Q1八、模型的评价与推广1.模型的评价自然界中，任何物种即使是捕食者也有自身的阻滞作用，该模型从原始的 没带自身阻滞作用模型中加入了阻滞项，使得此模型更接近于生态平衡系统 从此模型中，我们知道两物种同时灭绝是不稳定的，也就是不太可能的，但两 种群有一种灭绝一种生存是完全有可能的，两种群共存的可能也是可能的 2.模型的推广本文只考虑两物种模型，我们完全可以把此模型推广到三物种的情形。
自然界里长期存在的呈周期变化的生态平衡系统应该是结构稳定的，即系 统受到不可避免的干扰而偏离原来的周期轨道后，其内部制约作用会使系统自 动回复原状，如恢复原有的周期和振幅，而Vol terra模型描述的周期变化状态却不是结构稳定的要得到能反映周期变化的结构模型，要用到极限环的概念 参考文献[1] 姜启源，谢金星，叶俊．数学模型,高等教育出版社.2003 年[2] 冯杰，黄力伟，王勤.《数学建模原理与案例》科学出版社， 2007年1月。