《种群相互竞争模型》.docxVIP

数学实验设计 课题： 两种群相互竞争模型如下： x(t) r1x(1 — slX) n1 n2 y(t) 5 吒 * 其中x (t), y(t)分别是甲乙两种群、的数虽，r1, r2为它们 的固有增长率，n1 , n2为它们的最大容虽。s1的含义是， 对于供养甲的资源而言，单位数虽乙(相对 n2)的消耗虽 为单位数虽甲(相对 n1)消耗的s1倍，对于s2也可做相 应的解释 分析: 这里用x (t)表示甲种群在时刻t的数虽，即一定区域 内的数H。y(t)表示乙种群在时刻t的数H。假设甲种群独 立生活时的增长率(固有增长率)为 「1,则x (t)/ x=r1 ,而 种群乙的存在会使甲的增长率减小，且甲种群数虽的增长 也会抑制本身数H:的增长，即存在种间竞争。这里，我们 设增长率的一部分减少虽和种群乙的数虽与最大容纳虽的 比值成正比，与s1 (s1表示最大容纳虽乙消耗的供养甲的 资源是最大容纳虽甲消耗该资源的 s1倍)成正比。另一部 分的减少虽和种群甲的数虽与甲的最大容纳虽的比值成正 比。则我们可以得到如下模型： x(t)=r1*x*(1-x/n1-s1*y/n2) 同样，我们可以得到乙种群在 t时刻的数虽表达式: y(t)=r2*y*(1-s2*x/n1-y/n2) 如果给定甲、乙种群的初始值，我们就可以知道甲、乙种 群数H随时间的演变过程。 对于上述的模型，我们先设定好参数以后，就可以用 所学的龙格库塔方法及 MATLAB 软件求其数值解； 问题一： 设 r1=r2=1,n1=n1=100,s1=0.5,s2=2,初值 x0=y0=10 ,计算 x(t),y(t),画出它们的图形及相图(x,y)，说明时「d 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),r 2*(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 gtext( \fontsize(12}x(t) ),gtext( \fontsize (12}y(t) ), pause,plot(x(:,1),x(:,2)),grid, xlabel( x ),ylabel( y) 得到10年间甲、乙两种群数虽变化的图象为: 100 90 80 70 60 50 40 30 20 10 0 x(t) / / y (t) 0110 相图为: 25 20 ■ A \ 15 y Ti nn^\ / \ I 10 \ 5 0 10 20 30 40 50 60 70 80 90 100 x 结论：当t充分大时，x和y的数量悬殊变大，最终 是一方灭绝，一方繁荣。如上述模型中，甲种群繁 荣下去，乙种群很快灭绝。 问题二： 改变 r1,r2,n1,n2,x0,y0,但 s1,s2 不变，(或保持 s11,s21),计算并分析所得结果；若s1=1.5(1), s2=0.7(1)B分析结果，由此你的得到什么结论，请 用各参数生态学上的含义作出解释。 分析：当s1,s2不变(或保持s11,s21)时 当 s1=1.5(1), s2=0.7(1)时 当s1,s2不变（或保持s11,s21）时总有甲种 群繁荣，乙种群灭绝。而当 s1=1.5(1), s2=0.7(1) 时，有乙种群繁荣，甲种群灭绝。因此我们得到： 在两个种群的相互竞争中s1, s2是两个关键指标.从 上面对它们的解释可知，s1 l, s21表示在消耗供 养甲的资源中，乙的消耗多于甲，因而对甲增长的 阻滞作用乙大于甲，即乙的竞争力强于甲. 问题三： 实验当s1=0.8(1),s2=0.7(1)0寸会有什么样的结果： 当s1=1.5(1),s2=1.7(1)0寸乂会有什么样的结果。能 解释这些结果吗？ 分析：当 s1=0.8(1),s2=0.7(1)3寸有如图： 70 60 50 40 30 y( ) _一 “ x (t) 20 / / 0110 即甲、乙竞争激烈程度加剧，没有一方有明显优势; 当 s1=1.5(1),s2=1.7(1时乂会有如图: 80 100 x(t y(t) / 70 60 50