食饵—捕食者模型.pdfVIP

《数学模型》课程 食饵—捕食者模型 3. 讨论具有自身阻滞作用的两种群食饵-捕食者模型，首先根据 该两种群的相互关系建立模型，解释参数的意义，然后进行稳定性分 析，解释平衡点稳定的实际意义，对模型进行相轨线分析来验证理论 分析的正确性，并用 matlab 软件画出图形。 自然界中不同种群之间还存在着一种非常有趣的既有相互依存、 又有相互制约的生活方式：种群甲靠丰富的天然资源生长，而种群乙 靠捕食甲为生，形成鱼和鲨鱼，美洲兔和山猫，落叶松和蚜虫等等都 是这种生存方式的典型，生态学称种群甲为食饵，种群乙为捕食者。二 者共同组成食饵—捕食者系统。 一食饵—捕食者 用食饵(食用鱼)和捕食者(鲨鱼)为研究对象,设x (t) / x (t) 为食饵 1 (食用鱼)在时刻t 的数量，y (t) / x 2 (t) 为捕食者(鲨鱼)在时刻t 的数量，r1 为食饵(食用鱼)的相对增长率，r 为捕食者(鲨鱼)的相对增长率；N 为 2 1 大海中能容纳的食饵(食用鱼)的最大容量，N 2 为大海中能容纳的捕食 者(鲨鱼)的最大容量， 1 为单位数量捕食者（相对于N 2 ）提供的供养食 饵的实物量为单位数量捕食者(相对于N 1 )消耗的供养甲实物量的 1 倍； 2 为单位数量食饵（相对于N 1 ）提供的供养捕食者的实物量为单位 数量捕食者(相对于N 2 )消耗的供养食饵实物量的 2 倍；d 为捕食者离 开食饵独立生存时的死亡率 二模型假设 1.假设捕食者（鲨鱼）离开食饵无法生存； 2.假设大海中资源丰富，食饵独立生存时以指数规律增长； 三模型建立 食饵(食用鱼)独立生存时以指数规律增长，且食饵(食用鱼)的相对 增长率为r ，即x   rx ，而捕食者的存在使食饵的增长率减小，设减小 1 的程度与捕食者数量成正比，于是x (t) 满足方程  x (t)  x (r  ay )  rx  axy (1) 比例系数a 反映捕食者掠取食饵的能力。 由于捕食者离开食饵无法生存，且它独立生存时死亡率为 ，即 d y   dy ，而食饵的存在为捕食者提供了食物，相当于使捕食者的死亡 率降低，且促使其增长。设这种作用与食饵数量成正比，于是y (t) 满足  y (t)  y (d  bx)  dy  bxy (2) 比例系数b 反映食饵对捕食者的供养能力。 方程(1)、(2)是在自然环境中食饵和捕食者之间依存和制约的关 系，这里没有考虑种群自身的阻滞作用，是 Volterra 提出的最简单的 模型。结果如下。 不考虑自身阻滞作用:数值解 令 x(0)=x0,y(0)=0,设r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2 使用 Matlab 求解 求解如下 1）先建立 M 文件 function xdot=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)]; 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]