斑马的种群分析 食饵与捕食者.docx

B 题：斑马的种群分析摘要本文针对斑马与草场种群的相互影响这一自然界最常见的食饵—捕食者问题，在 Volterra 模型的基础上，提出了完善方案，使之更切合实际本文共建立三个模型：模 型一假设草的生长遵从 Logistic 规律，斑马的增长不受限制，建立了单阻滞模型；模 型二假设斑马的增长也遵从Logistic规律，建立了双阻滞模型同时，运用MATLAB软 件画出不同情况下草场密度和斑马数量随着时间变化的图像以及草场密度和斑马数量 变化的相轨线模型三则是把线性功能反应函数替换为Holling II功能反应函数下的 Leslie 捕食者-食饵模型，定性的讨论了其长久性及稳定性x(t) = xr (1■上1 x y mm• ◎・x y (t) = y gd+t )xm针对问题一， 模型一中考虑到食饵自身的阻滞作用， 建立了单阻滞模型 )；模型二进一步考虑捕食者的自身阻滞作用，建立了双阻滞模型X (t) = f (x, y) = rx(l—丄一Q —)1 x i y然后，分别求出平衡点，进行稳定性分析，< mmy(t) = g(x, y) = r y(—1—丄 +q —) 、 2 y 2 xmm再运用 MATLAB 软件分析验证，并且通过图像分析，对比模型一和模型二的不同点。

针对问题二，在模型一、二中适当改变参数r , d , x，观察各参数变化对两个种群 m数量的影响，得出结论但通过作图只能直观地描述其变化趋势，缺乏说服力因此建立模型三x(t) = x(1-x)-xy/(ay + x)，以进行定性分析当邮+1卄， y (t) = 5 y (卩-y / x)(«P + 2)p < (8p+ 1)(ap +1)2时，捕食者和食饵最终能达到平衡状态，种族能持续生存 当邮+1 >卩，(ap+ 2) p> (8p+ 1)(ap +1)2时，平衡点不稳定，围绕平衡点出现极限环, 且极限环唯一，捕食者和食饵振荡共存最后，本文提出了考虑建立自然环境受时间影响而产生周期性变化的条件下的食饵 —捕食者模型等改进建议关键词：Vol terra模型 单阻滞 双阻滞Holling II功能反应函数 极限环一、问题重述自然界中不同种群之间存在着一种非常有趣的既依存又有制约的生存方式：种群甲 靠丰富的自然资源生长，而种群乙靠捕食种群甲为生生态学上称种群甲为食饵，种群 乙为捕食者，两者共组成食饵一捕食者关系（简称P-P系统），本文的问题就是围绕食饵 —草，捕食者—斑马，而提出的现假设将一定数量的斑马放入草场，研究草和斑马两种群的相互作用，草的生长遵从Logistic规律，年固有增长率r =0.8，最大密度为x二3000 （密度单位），有草时每 m只斑马每年可吃掉b 1 = 1.6 （密度单位）的草；若没有草，斑马的年死亡率高达d = 0.9, 而草的存在可使斑马的死亡得以补偿，有草时补偿率为b 2 = 1.5,草场中最多容许4000 只斑马生存。

要求根据以上资料，作出一些简化假设，用微分方程模型描述草和斑马两 种群数量的变化过程，回答以下两小问：1、 将y0 = 100只斑马放入密度x为1000和密度为3000的草场两种情况，编程计算， 画出相应的草场密度和斑马数量随着时间变化的图形，草场密度和斑马数量变化的相轨 线2、 适当改变参数r，d，x，观察两种群数量新的变化趋势m二、问题分析本文斑马与草场种群问题是典型的捕食者 -食饵模型，由于其广泛的存在性和重要 性，长久以来一直是生态学和生物数学的研究论点原有的 Volterra 模型有其局限性， 许多生态学家指出，多数食饵—捕食者系统都观察不到 Volterra 模型显示的那种周期震 荡，而是趋向某种平衡状态，即系统存在稳定平衡点实际上，只要在Vol terra模型中 加入考虑自身的阻滞作用的 Logistic 项就可以描述现实中的现象我们在Vol terra食饵一捕食者模型的基础上，进一步改进和完善首先考虑草的生 长遵从Logistic规律，斑马的增长不受限制的单阻滞食饵一捕食者（P-P）模型，并运 用matlab软件画出不同情况下草场密度和斑马数量随着时间变化的图形，草场密度和斑 马数量变化的相轨线；接着考虑斑马的增长也遵从Logistic规律的双阻滞食P-P模型， 并画出相应图形。

在模型一、模型二下适当改变参数r , d , x，观察两种群数量变化趋势但通过作m 只能直观地描述其变化趋势，缺乏说服力因此建立模型三，把线性功能反应函数替换 为Hollingll功能反应函数，进行定性分析求出平衡点，讨论其长久性及稳定性与 现实生活对照，应该有捕食者和食饵最终能达到平衡状态，种族能持续生存或者捕食 者和食饵振荡共存最后我们就本文所有模型的优缺点及用途进行评价，给出进一步优 化方向，使模型更符合实际三、符号说明草遵从Logistic规律增长这一条件下的年固有增长率，已知r =0.8r2 斑马遵从Logistic规律增长这一条件下的年固有增长率有草时，每只斑马每年的食草量，已知b = 1.6 (密度单位) 1有草时，草对斑马死亡率的补偿，使得斑马数量有所增长，已知b = 1.52没有草时，斑马的年死亡率，已知d二0.9x 自然环境下草的最大密度，已知x二3000 (密度单位)mmy 斑马的最大环境容纳量，已知y二4000 (只)mm四、模型假设1. 假设自然环境中只有草和斑马两个物种2. 假设自然环境不随时间的变化而变化3．草独立生存，独立生存规律遵从 Logistic 规律4. 在模型二中，斑马种群的生存规律也遵从 Logistic 规律 5．鹿无法独立生存。

没有草的情况下，鹿的年死亡率一定 6．假定草对鹿的补偿率是草场密度的线性函数 7．每只鹿每年的食草能力是草场密度的线性函数五、模型的建立与数据处理模型一 单阻滞模型・5.1.1微分方程的建立现将一定数量的斑马放入草场，基于模型假设，参照 Volterra 的“食饵—捕食者 模型”在其中加入考虑草自身阻滞作用的Logistic项加以优化斑马和草在时刻t的数量分别记作x(t)，y(t)，此时假设草以Logistic规律增长, 固定增长率为r，即有X(t) = x[r (1-上)]，而斑马的存在使草的增长率减小，设减小的 1xm程度与斑马数量成正比，于是x(t)满足方程：・ /、 - x b . yxx(t)二 xr (1— (1)1 x y mm 比例系数b反映每只斑马每年的食草量1斑马没有食物来源也无法生存，没有草时它们的死亡率为d,即有y (t)二-dy,而 草的存在为斑马提供了食物，相当于使斑马的死亡率降低，且促使其增长设这种作用 与草的数量成正比，于是 y(t) 满足：2)y (t )= y — )xm比列系数b反映草对斑马的供养能力2・5.1.2模型分析 方程(1)，(2)没有解析解，我们分两步对这个模型所描述的现象进行分析。

首先从理论上研究不同初始值下其平衡点的值并进行稳定性分析然后，利用数学软件 MATLAB 编程，画出相应的草场密度和斑马数量随着时间变化的图形及相轨线・5.1.2.1平衡点的求解及稳定性分析根据以上模型二阶方程(1)，(2)可用两个一阶方程表示为x(t) = f (x, y)3)Vy (t) = g (x. y)右端不显含t，是自治方程代数方程组］/(x，y)二0的实根x二x0 , y二y0称为(3)I g (x. y)二 0的平衡点，记作 P (x0, y0) 0如果存在某个领域，使方程( 6)的解 x(t) ， y(t) 从这个领域内的某个 (x(0), y(0)) 出 发，满足4)limx(t)二 xo, limy(t)二 y0t fg t fg则称平衡点 P 是稳定的；0否则,称 P 是不稳定的0求P0(x0,y0),解方程如下:0.8x(1—x30001.6y4000)— 0y (—0.9 +1.5 x3000解得三组解：1 x o = o 1 x o = 3ooo x o = 18ooV 1 ， V 2 ， A 3y o = o y o = o 1 y o = looo1 2 3即有平衡点 P (o,o) ，1P (3ooo,o) ， P (18oo,1ooo)23对于非线性方程，我们用近似线性方法判断其平衡点的稳定性。

在P点将f (x,y)和 0g(x.y)作Taylor展开，只取一次项，得(3)的近似线方程x(t) = f (xo, yo)(x-xo) + f (xo, yo)(y - yo)xyy(t) = g (xo, y0)(x-xo) + g (xo, y0)(y - y0)y6)系数矩阵记作xg'x7)带入f (x)和g (x)得到:ry r - 2 亠 x - r a - 1 x 1 1 y mm-ra112x将平衡点带入A中计算相应的p和q，其中p = 一(f' + g')1、q二det(AI ) (i二1,2,3)，当 x y Pi Pip > o,q >o时稳定经计算得到在各个平衡点稳定性如表1平衡点pq稳定条件p =(o,o)—r + r-rr不稳定11 21 2P = (x ,o)r - r (a -1)-rr (a -1)a < 12 m1 2 21 2 22(d .x y (1-d / a )、P — ( m , ― 2—)r-1-rr (a -1)a > 13 a a2 1a21 2 22表1 具有自身阻滞作用的食饵—捕食者模型的平衡点及稳点性根据表1可知，p = (0,0)不稳定；当& < 1时，因食饵不能为捕食者提供足够食物,12P稳定，即捕食者将灭绝而食饵趋向最大容量；而当b > 1时，P稳定，草能够为捕食2 2 3 者提供足够的食物，二者共存，分别趋向非零的有限值，这也是食饵—捕食者保持共存 的最大数量。

・5.1.2.2利用MATLAB编程绘图记草和斑马的初始数量分别为：x(0)二 x , y (0)二 y00为求微分方程(1), (2)及初始条件(3)的数值解x(t) , y (t)及相轨线y(x), 代入相关系数的值，以及|Xo二100和|Xo二3000分别绘图如下：I y 二 100 I y 二 10000t值46.31447.81449.31450.81452.31453.81455.31456.48557.657x值1800.118001799.9180018001800180018001800y值999.97999.99100010001000999.9910001000100058.828180010006018001000图像分析：假设斑马和草的初始种群数量分别为 1000。