食饵捕食模型.docVIP

楚雄师范学院数学系《数学建模》课程

教学论文

题目：具有自身阻滞作用的两种群食饵—捕食模型

专业：信息与计算科学

班级：08级3班

学号：20081022152

学生姓名：罗文枢

完成日期：2011年6月

具有自身阻滞作用的两种群食饵—捕食模型

摘要：在自然界中，更多的生物是杂居在一起的，各种生物根据其生理特点、食物来源分成了不同的层次，各层次之间及同一层次的生物种群之间有着各样的联系，尤其是相互之间影响非常大的生物种群，需要放在一起讨论，在这里，我们一两种群为例进行建模和讨论，具有自身阻滞作用的两种群食饵—捕食者模型。捕食—食饵模型是数学生态学研究的重要内容,影响种群波动的因素很多,自身阻滞作用就是其中重要的一种因素。因为资源环境是有限的，相互竞争是不可防止的，所以自身阻滞也是影响平衡位置的不稳定性和周期波动现象的主要因素。时滞可以对生态系统的性质产生相当大的影响,理论生态学家们普遍认为在种群的相互作用中,自身阻滞作用是不可防止的。本文主要通过对两类具有自身阻滞作用的典型的捕食-食饵模型的研究,通过分析发现时滞对模型的稳定性有非常重要的作用。事实上只要在Volterra模型参加考虑自身阻滞作用的Logsitic项就可以得到这种现象了。

关键字：自身阻滞，稳定性分析，相轨线分析，平衡点分析，Logistic模型；

一.问题重述：

讨论具有自身阻滞作用的两种群食饵—捕食者模型，首先根据两种群的相互关系建立模型，解释参数的意义，然后进行稳定性分析，解释平衡点稳定的实际意义，对模型进行相轨线分析来验证理论分析的正确性。

二．问题分析：

本论文主要是讨论具有自身阻滞作用的食饵—捕食者模型。我们用Logistic模型来描述这个种群数量的演变过程，即食饵会受到自然界中的资源所限制，它不仅会无限的增大，而且捕食者也会受到食饵的数量的影响。此种情况下会出现以下的3种现象：

1.当捕食者灭绝时，食饵也不会无限的增长，即指数函数型增长，因为有自身的阻滞作用，它到达某个数量就不在会增长而趋于稳定了；

2.当食饵受到自然资源的影响的灭绝时，捕食者也会因食物而灭绝；

3.当两种群都不灭绝时，它们会趋于某个非零的有限值，从而到达稳定状态。

三．模型假设：

1.假设在某特定环境中只存在食饵和捕食者两种群；

2.假设食饵和捕食者均能正常生长，没有疾病等原因促使死亡；

3.假设两种群的增长率不变；

4.食饵由于捕食者的存在使增长率降低，假设降低的程度与捕食者数量成正比；

5.捕食者由于食饵为它提供食物的作用使其死亡率降低或使之增长，假设增长的程度与食饵数量成正比。

四．符号说明：

：食饵在时刻的数量；

：捕食者在时刻的数量；

：食饵独立生存时以指数规律增长，相对增长率；

：捕食者独立生存时以指数规律增长，相对增长率；

：食饵生存的最大容量；

：捕食者生存的最大容量；

?1：食饵自身的竞争能力；

?2：捕食者自身的竞争能力；

五．模型的建立及求解：

1.模型建立：

将食饵和捕食者在时刻的数量分别记为和，我们根据自身的阻滞增长作用，知道存在一定的规律，在Volterra模型中参加考虑自身阻滞作用的Logistic项，那么根据前面的假设建立以下模型：

2.模型求解：

根据数学建模的相关知识，利用数学软件MATLAB求解。下面对这两个个方程进行平衡点稳定性分析，即设

;

通过解这个方程组得到3个平衡点,分别是：

(0,0),(N,0),;

当时，才有意义

按照平衡点稳定性的方法计算可知：

|,

将3个平衡点的p,q德结果及稳定条件情况如下表所示：

表一：平衡点p,q的结果及稳定条件

平衡点

p

q

稳定条件

(0,0)

不稳定

(N,0)

3.下面根据表一对平衡点稳定性进行评价：

1〕对点：当时，因为食饵不能为捕食者提供足够的食物，而使捕食者灭绝，从而食饵到达最大容量，那么使到达稳定；

2〕对点：当时，因为二者共存，分别趋向非零的有限值，从而到达稳定。

六．模型的检验：

首先先来看相轨线的图像，从图像中进行分析；

并设=,=m;那么得到相应的MATLAB代码如下：

1.先写M文本：

functionxdot=shier(t,x)

=2.5;=1.8;=0.5;=1.6;=1.6;=1;

xdot=[r1*(1-x(1)/-m1*x(2)/).*x(1);R2*(-1+m2*x(1)/N1-x(2)/)];

2.MATLAB程序为：

ts=0:0.1:15;

x0=[25,2];

[t,x]=ode45(shier,ts,x0);[t,x]

plot(t,