




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数学实验设计课题:两种群相互竞争模型如下:x(t)r1x(1xs1y)n1n2y(t)r2y(1s2xy)n1n2其中x(t),y(t)分别是甲乙两种群`的数量,r1,r2为它们的固有增长率,n1,n2为它们的最大容量。s1的含义是,对于供养甲的资源而言,单位数量乙(相对n2)的消耗量为单位数量甲(相对n1)消耗的s1倍,对于s2也可做相应的解释。分析:这里用x(t)表示甲种群在时刻 t的数量,即一定区域内的数量。y(t)表示乙种群在时刻 t的数量。假设甲种群独立生活时的增长率(固有增长率)为 r1,则x(t)/x=r1,而种群乙的存在会使甲的增长率减小,且甲种群数量的增长也会抑制本身数量的增长,即存在种间竞争。这里,我们设增长率的一部分减少量和种群乙的数量与最大容纳量的比值成正比,与s1(s1表示最大容纳量乙消耗的供养甲的资源是最大容纳量甲消耗该资源的 s1倍)成正比。另一部分的减少量和种群甲的数量与甲的最大容纳量的比值成正比。则我们可以得到如下模型:x(t)=r1*x*(1-x/n1-s1*y/n2)精品文档交流同样,我们可以得到乙种群在 t时刻的数量表达式:y(t)=r2*y*(1-s2*x/n1-y/n2)如果给定甲、乙种群的初始值,我们就可以知道甲、乙种群数量随时间的演变过程。对于上述的模型,我们先设定好参数以后,就可以用所学的龙格库塔方法及 MATLAB 软件求其数值解;问题一:设r1=r2=1,n1=n1=100,s1=0.5,s2=2,初值x0=y0=10,计算x(t),y(t),画出它们的图形及相图 (x,y),说明时间 t充分大以后x(t),y(t)的变化趋势(人民今天看到的已经是自然界长期演变的结局)。编写如下 M文件:function xdot=jingzhong(t,x)r1=1;r2=1;n1=100;n2=100;s1=0.5;s2=2;xdot=diag([r1*(1-x(1)/n1-s1*x(2)/n2),r2*(1-s2*x(1)/n1-x(2)/n2)])*x;然后运行以下程序:ts=0:0.1:10;x0=[10,10];[t,x]=ode45(@jingzhong,ts,x0);[t,x]plot(t,x),grid,精品文档交流gtext( '\fontsize{12}x(t)' ),gtext( '\fontsize{12}y(t)' ),pause,plot(x(:,1),x(:,2)),grid,xlabel( 'x' ),ylabel( 'y' )得到10年间甲、乙两种群数量变化的图象为:100908070x(t)605040302010y(t)0123456789100相图为:精品文档交流252015y1050203040506070809010010x结论:当t充分大时,x和y的数量悬殊变大,最终是一方灭绝,一方繁荣。如上述模型中,甲种群繁荣下去,乙种群很快灭绝。问题二:改变r1,r2,n1,n2,x0,y0,但s1,s2不变,(或保持s1<1,s2>1),计算并分析所得结果;若 s1=1.5(>1),s2=0.7(<1)再分析结果,由此你的得到什么结论,请用各参数生态学上的含义作出解释。分析:当s1,s2不变(或保持s1<1,s2>1)时精品文档交流5045 x(t)4035302520151050
y(t)0 1 2 3 4 5 6 7 8 9 10当s1=1.5(>1),s2=0.7(<1)时100908070605040302010
y(t)x(t)00 1 2 3 4 5 6 7 8 9 10当s1,s2不变(或保持s1<1,s2>1)时总有甲种精品文档交流群繁荣,乙种群灭绝。而当 s1=1.5(>1),s2=0.7(<1)时,有乙种群繁荣,甲种群灭绝。因此我们得到:在两个种群的相互竞争中 s1,s2是两个关键指标.从上面对它们的解释可知, s1>l,s2<1表示在消耗供养甲的资源中,乙的消耗多于甲,因而对甲增长的阻滞作用乙大于甲,即乙的竞争力强于甲.问题三:实验当s1=0.8(<1),s2=0.7(<1)时会有什么样的结果:当s1=1.5(>1),s2=1.7(>1)时又会有什么样的结果。能解释这些结果吗?分析:当s1=0.8(<1),s2=0.7(<1)时有如图:70y(t)60x(t)50403020100 1 2 3 4 5 6 7 8 9 10即甲、乙竞争激烈程度加剧,没有一方有明显优势;精品文档交流当s1=1.5(>1),s2=1.7(>1)时又会有如图:8070x(t)60504030y(t)20100 1 2 3 4 5 6 7 8
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 专利使用保证协议书
- 假货和解协议书范本
- 物业合作买菜协议书
- 物品撞坏赔偿协议书
- 奖励股权分配协议书
- 商事和解协议书范本
- 试制模具技术协议书
- 物品收集转运协议书
- 土地确权协议书范本
- 养鸡回收协议书范本
- “条令条例学习月”主题授课课件
- 海洋生态环境监测技术-全面剖析
- 2024年湖北省中学教师招聘考试真题
- 卫星科普知识
- 北京市朝阳区2025届高三一模质量检测一 语文试题(含答案)
- 新教材高中生物选择性必修2课件:1 2 种群数量的变化(人教版)
- 车辆租赁服务保障计划
- 《装配式混凝土建筑》全套教学课件
- 《敖丙》幼儿园小学少儿美术教育绘画课件创意教程教案
- 高职单招职业技能测试题库
- 2025年正德职业技术学院单招职业技能测试题库一套
评论
0/150
提交评论