捕食者-被捕食者模型稳定性分析

捕食者—捕食者模型稳定性分析
【摘要】自然界中不同种之间还存在着一种非常有趣的既有相互依存、又有相互制约的生活方式:种甲靠丰富的天然资源生存,种乙靠捕食甲为生,形成食饵-捕食者系统,如食用鱼和鲨鱼,美洲兔和山猫,害虫和益虫等。本文是基于食饵—捕食者之间的有关规律,建立具有自身阻滞作用的两种食饵—捕食者模型,分析平衡点的稳定性,进行相轨线分析,并用数值模拟方法验证理论分析的正确性。
【关键词】食饵—捕食者模型  相轨线  平衡点  稳定性 
一、问题重述
在自然界中,存在这种食饵—捕食者关系模型的物种很多。下面讨论具有自身阻滞作用的两种食饵-捕食者模型,首先根据该两种的相互关系建立模型,解释参数的意义,然后进行稳定性分析,解释平衡点稳定的实际意义,对模型进行相轨线分析来验证理论分析的正确性。
二、问题分析
  本文选择渔场中的食饵(食用鱼)和捕食者(鲨鱼)为研究对象,建立微分方程,并利用数学软件MATLAB求出微分方程的数值解,通过对数值结果和图形的观察,猜测出它的解析解构造。然后,从理论上研究其平衡点及相轨线的形状,验证前面的猜测。
三、模型假设
1.假设捕食者(鲨鱼)离开食饵无法生存;
2.假设大海中资源丰富,食饵独立生存时以指数规律增长;
四、符号说明
/——食饵(食用鱼)在时刻的数量;
/——捕食者(鲨鱼)在时刻的数量;
——食饵(食用鱼)的相对增长率;
——捕食者(鲨鱼)的相对增长率;
——大海中能容纳的食饵(食用鱼)的最大容量;
——大海中能容纳的捕食者(鲨鱼)的罪的容量;
——单位数量捕食者(相对于)提供的供养食饵的实物量为单位数量捕食者(相对于)消耗的供养甲实物量的倍;
——单位数量食饵(相对于舵角指示器)提供的供养捕食者的实物量为单位数量捕食者(相对于)消耗的供养食饵实物量的倍;
——捕食者离开食饵独立生存时的死亡率。
五、模型建立
食饵独立生存时以指数规律增长,且食饵(食用鱼)的相对增长率为,即,而捕食者的存在使食饵的增长率减小,设减小的程度与捕食者数量成正比,于是满足方程
        (1)
比例系数反映捕食者掠取食饵的能力。
由于捕食者离开食饵无法生存,且它独立生存时死亡率为,即,而食饵的存在为捕食者提供了食物,相当于使捕食者的死亡率降低,且促使其增长。设这种作用与食饵数量成正比,于是满足
            (2)
比例系数反映食饵对捕食者的供养能力。
方程(1)、(2)是在自然环境中食饵和捕食者之间依存和制约的关系,这里没有考虑种自身的阻滞作用,是Volterra提出的最简单的模型。
下面,我们加入种自身的阻滞作用,在上两式中加入Logistic项,即建立以下数学模型:
    (3)
    (4)
六、模型求解
在此,我们采用MATLAB软件求解此微分方程组中的金刚石绳锯的图形及相轨线图形。
,,使用MATLAB软件求解,程序代码如下:
1)建立M文件
function y=fun(t,x)
y=[x(1).*(1-x(1)./3500-1.5*x(2)./500),0.4.*x(2).*(-1+4.*x(1)./3500-x(2)./500)]';
2)在命令窗口输入如下命令:
[t,x]=ode45('fun1',[0,40],[2000,35])
得到数值解如下:
t
(x(1),x(2)) 1.0e+003 * (单位:千克)
        0
    0.1033
    0.2066
    0.3099
    0.4132
    0.8079
    1.2026
    1.5973
    1.9919
    2.3806
    2.7693
    3.1579
    3.5466
    3.9353
    4.3239
    4.7126
    5.1012
    5.4167
    5.7322
    6.0477
    6.3631
    6.7446
    7.1261
    7.5076
    7.8891
    8.2967
    8.7042
    9.1117
    9.5193
  10.0173
  10.5153
  11.0133
  11.5113
  12.0453
  12.5793
  13.1133
  13.6473
  14.2231
  14.7989
  15.3747
  15.9505
  16.5523
  17.1541
  17.7559
  18.3577
  19.0403
  19.7229
  20.4054
  21.0880
  21.8574
  22.6268
  23.3962
  24.1656
  25.0656
  25.9657
  26.8657
  27.7657
  28.7657
  29.7657
  30.7657
  31.7657
  32.7657
  33.7657
  34.7657
皮画  35.7657
  36.8243
8700g  37.8828
  38.9414
  40.0000
2.0000    0.0350
    2.0654    0.0369
    2.1276    0.0389
    2.1863    0.0412
    2.2412    0.0438
    2.4113    0.0560
    2.5111    0.0732
    2.5358    0.0961
    2.4870    0.1248
    2.3741    0.1577
    2.2127    0.1922
    2.0228    0.2246
    1.8247    0.2514
    1.6360    0.2697
    1.4727    0.2793
    1.3431    0.2813
    1.2427    0.2775
    1.1780    0.2717
    1.1294    0.2642
    1.0951    0.2561
    1.0728    0.2476
    1.0599    0.2377
    1.0593    0.2286
    1.0684    0.2206
    1.0851    0.2139
    1.1094    0.2081
    1.1376    0.2038
    1.1676    0.2009
    1.1975    0.1993
    1.2313    0.1990
    1.2599    0.2000
    1.2815    0.2020
    1.2954    0.2047
    1.3022    0.2079
    1.3017    0.2111
    1.2958    0.2138
    1.2864    0.2159
    1.2745    0.2174
    1.2627    0.2181
    1.2529    0.2181
    1.2455    0.2177
    1.2402    0.2169
    1.2377    0.2160
    1.2375    0.2151
    1.2391    0.2144
    1.2420    0.2138
    1.2454    0.2134
    1.2484    0.2134
    1.2506    0.2135
    1.2524    0.2137
    1.2530    0.2140
    1.2527    0.2142
    1.2520    0.2144
    1.2509    0.2145
工业盐水    1.2499    0.2145
    1.2495    0.2144
    1.2495    0.2143
    1.2496    0.2143
    1.2498    0.2142
    1.2500    0.2142
    1.2501    0.2143
    1.2501    0.2143
    1.2501    0.2143
    1.2500    0.2143
    1.2500    0.2143
    1.2500    0.2143
    1.2500    0.2143
    1.2500    0.2143
    1.2500    0.2143
语音系统

本文发布于:2024-09-22 17:40:06,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/3/243644.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:捕食者   食饵   分析   生存   模型   种群   作用
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议