- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
三个种群关系的数学模型
摘要:自然界中不同种群之间存在着一种非常有趣的既有依存、乂有制约的牛 存方式:种群甲靠丰富的自然资源生长.而种群乙靠捕食种群甲为生.种群丙又 靠捕食种群为生.类似的现象还存在很多。假设这种情况:在草原上栖息着羚羊 和狮子,羚羊以草为食,狮子又以羚羊为食。由此建立描述三个种群数量变化规 律的微分方程模型,分析平衡点稳定性。
关键词:种群数解值平衡点
一、问题重述
自然界中不同种群之间存在着一种非常有趣的既有依存、又有制约的生存方 式:种群甲靠丰富的自然资源生长。而种群乙靠捕食种群甲为生。种群丙又靠捕 食种群为生.类似的现象还存在很多。假设这种情况:在草原上栖息着羚羊和狮 子,羚羊以草为食,狮子乂以羚羊为食。由此建立描述三个种群数量变化规律的 微分方程模型,分析平衡点稳定性。
二、模型假设
1、 假设在草原上只存在了草、羚羊和狮子三个物种。并且羚羊离开草,狮子离 开羚羊无法生存。
2、 不考虑自然资源对草的生长的限制。羚羊和狮子都可以正常生长,没有其他 原因可以使狮子和羚羊死亡。
3、 草的种子可以正常的生长,小羚羊也可以正常的生长。
三 符号说明
西(/):草在吋刻t的数量;
x2 (0 :羚羊在时刻t的数量;
(r):狮子在时刻t的数量;
中 草独立生存时以指数规律增长,相对增长率;
r2:羚羊离开植物无法生存,它独自存在时的死亡率;
厂3:狮子离开羚羊无法生存,它独自存在吋的死亡率;
羚羊食草的能力;
入:草对羚羊的供养能力;
希:羚羊对狮子的供养能力;
A:狮子掠取羚羊的能力。
环境容许的草的最人数量;
n2:羚羊的最大容量;
n3:狮子的最大容量;
5:单位数量的羚羊(相对于N2而言)掠取5倍的单位数量的草量(相对于M 而言);
t2:单位数量的草(相对于R而言)供养6倍的单位数量的羚羊量(相对于M 而言);
O-3:单位数量的狮子(相对于他而言)掠取6倍的单位数量的羚羊量(相对于
“2而言);
6:单位数量的羚羊(相对于N?而言)供养6倍的单位数量的狮子量(相对于
弘而言);
!1!
把草、羚羊、狮子的数量分别记作坷(/),兀卫),心⑴若不考虑自然资源对 草的限制,草独立生存时以指数规律增长。相对增长率为,即羚羊的存在使草 的增长率减小.设减小的程度与羚羊数量成正比。于是草的模型为:
石(/)=兀](斤一人兀2) ( 1)
比例系数入反映羚羊食草的能力。
羚羊离开草无法生存.设它独自存在吋死亡率为厂2,即:
x2(t) = -r2x2
而草的存在乂为羚羊提供了食物。草的存在相当于使羚羊的死亡率降低。且 促使羚羊增长。设这种作用与草的数量成正比。则有:
兀2(=无2(一G +入兀J
比例系数乙反映草对羚羊的供养能力;
羚羊又为狮子提供了食物。狮子的存在使羚羊的增长率减小。设减小的程度 与羚羊的数量成正比,于是上式右端应减去狮子对羚羊增长的阻滞作用。于是羚
羊的模型应为:
(2)X2 (Z) = X2 (-?2 +兄2石一严3)
(2)
比例系数“反映狮子掠取羚羊的能力。
狮子离开羚羊无法生存。设它独自存在时死亡率为心,即
x3(t) = -r3x3
而羚羊的存在又为狮子提供了食物。相当于使狮子的死亡率降低。且促使狮 子的增长。于是狮子的模型为:
七,(0 =心(一E + 几3*2 ) (3)
比例系数希反映羚羊对狮子的供养能力;
方程(1)、(2)、(3)构成草、羚羊、狮子三者既依存又制约的数学模型:
兀](0 =兀1 (-斤+人兀2)
X2 (/) = X2 (-,2 + 久2兀1 一 %尢3) (4)
x3 (/) = x3 (-r3 + A3x3)
五、模型求解
微分方程组(4)没有解析解,可Mat lab求微分方程组(4)的数值解,通过 对数值结果和图形的观察,猜测它的解析解的构造;
兀](0)=花,兀2(°)=他,兀3(°)= X1 (5)
求方程组⑷及⑸的数值解x}(“,尢2(0,七⑴,并作图。
设厂]=1 , r2 =0. 5 , r3 =0. 6 , \ =0. 1 , A2 =0. 02 , 入二°? 06. “二0. 1 ,
兀5 = 1OQ 兀6 = 40, Xj — 6 o
Mat lab程序:
M程序:
function f=funl(t, x);
rl=l;r2=0. 5;r3=0. 6;lambdal二0. 1;lambda2=0. 02;lambda3=0. 06;mu二0. 1;
f二[x (1) * (rITambdal*x (2)) ; x (2) * (-r2+lambda2*x (1) 一mu*x (3)) ; x (3) * (-r3+l ambda3*x(2))];
在Mat lab窗口输入:随时间的变化,三个种群的数量的变化曲线。
? [t,x]=
原创力文档


文档评论(0)