大学数学实验

1、4.在某种环境下猫头鹰的主要食物是田鼠,设田鼠的年平均增长率为 ,猫头鹰的存在引起的田鼠增长率的减少与猫头鹰的数量成正比，比例系数为 ；猫头鹰的年平均减少率为 ;田鼠的存在引起的猫头鹰减少率的增加与田鼠的数量成正比,比例系数为 。建立差分方程模型描述田鼠和猫头鹰共处时的数量变化规律，对以下情况作图给出50年的变化过程。(）设开始时有100只田鼠和50只猫头鹰。（2) 同上，开始时有10只田鼠 和20只猫头鹰。（)适当改变参数 （初始值同上）（)求差分方程的平衡点,它们稳定吗？解: 记第k代田鼠数量为k，第k代猫头鹰数量为yk,则可列出下列方程: xk+1=xk+(r-k*1)xkk+1=yk+

2、(r2+ka）*y 运用matlb计算，程序如下:funcion z=disiti(x,0,a,a2，r1,r2)=x0;y=y；fo =:4 x(k+1)=x(k）+(r-y(k)*a)*x(k); y（k+1)=（)+(-r2+x(k)*a2)*y(k);endz=x,；（)z=dsti(100,5,000,0.002,0.,0.3)plot(1:5,z(:,1);ol on; plot（1:0,z(:，2),r)得到如下结果:z = 100.00 50.00 115.000 5.0000 32.5 41.8500 153.831 0.412 178.30 4.722 2067930 43

3、036 239.255 4.9216 26376 5.4831518 70.6668 35.969 94.0149 393.6908 132.74 42.196 9383 41.259 304.12 373850 49.057 8285 6.408 142757684.819 .3682 41.8092 17.832672.36 9.0 44443 6.5828 35.3402 5.5602 53.414 52318.296 5.3674 12.036 5.36 91.0051 .387 64.7508 .437 4.8 8.3581 32.949 .7541 23.6336 11.474 1

4、70 3.5742 12.2935 16.221 8.92 19.2024 6.557 22.172 4.8333 739 64 32761 2.79 9238 2.82 46.9988 1.621 56.322 1.875 7.13 .46380.946 0873 9.051 0.75 11.051 063 19077 .6280 67.441 0.650 20269 636 24045 0.711 288.85 08286 346.4029 1.58 415.3167 1.475497777 2.2570红线为猫头鹰数量曲线，蓝线为田鼠曲线(2) z=dt(00，20，0.01,0.02,

5、02，.）po(1:50,z(:,1）)；hod on； plt(1:5,(：,2),r） = 0.0000 20.000 00.0000 18.000 02.000012.00 15870 146.4811.559 13.243 189610 123.55 .0906 .6037 1.9010 1.5381 51.373 0884 15.2367 1037 180.3771 111.6737 93091 1.84 12.3165 19.4298 22797 145.610 29.37 68.0647 247.3278 19806 2771 23.7886 23.661 230909 18.8

6、20 333.2864 18.6595 379.163155.693 409.238 123.052 13.8872 6.7745 391.624 78.2302 349.36 66.5007 99.7060 598702 249655 56.8973 204.57 56.6326 1662 58.272 15.445 2.304 10.663 67.8714 92566 75.19 76261 84.5631 6.84 5.9905 56.76 1.77 50.5 126.211 .212 145613 44.1 6887 43.655319.461 45.256 224765 42965

7、58.39 56.66 2571 68.807 33.547 9.02 3113 21.9469 40.5 175.988 101726 264.167 853 01.23.4586 589.4677 18.1029 77.222 79.7204 81.31 红线为猫头鹰数量曲线，蓝线为田鼠曲线(3)当a,a分别取0002,02时,得到如下图像：可见,当a1,a参数在一定范围内改变时，猫头鹰与田鼠数量在一定范围内震荡，且不灭绝。()令xk=k1x； yk=yk1=y解方程得到如下结果:x=150y=0经matlab验证如下：isiti(50，0,0.0,0.002,.2,0.3）lo(1:50

8、,(：,);hold on； plt(1:50,z(:，2),r)z = 150 00 10 200 150 200 1 20 150 200 0 200 150 200 150 200 10 0 15 00 15 00 150 20 15 150 00 150 20 150 200 15 20 10 00 0 20 10 00 150 20 150 200 50 200 150 20 10 0 150 200 5 00 150 200 150 20 150 200 50 200 0 0 50 200 15 200 50 150 15 00 150 20 10 0 150 00 10 200

9、15 2 150 20 150 2 15 20 10 200 5 00 150 00 0 10 00 由此可知:平衡点为：x10 y200研究将鹿群放入草场后草和鹿两种群的相互作用。草的生长遵从logistic规律，年固有增长率08，最大密度为000(密度单位），在草最茂盛时每只鹿每年可吃掉16（密度单位)的草。若没有草，鹿群的年死亡率高达0.9，而草的存在可使鹿的死亡得以补偿,在草最茂盛时补偿率为15。作一些简化假设,用差分方程模型描述草和鹿两种群数量的变化过程，就以下情况进行讨论:(1)比较将100只鹿放入密度为1000和密度为3000的草场两种情况。（)适当改变参数，观察变化趋势。解:设

10、1.草独立生存,独立生存规律遵从logisc规律;草场上除了鹿以外，没有其他以草为食的生物;.鹿无法独立生存。没有草的情况下,鹿的年死亡率一定；4.假定草对鹿的补偿率是草场密度的线性函数；每只鹿每年的食草能力是草场密度的线性函数。记草的固有增长率为r，草的最大密度为n,鹿独立生存时的年死亡率为d,草最茂盛时鹿的食草能力为a，草对鹿的年补偿作用为b;第k年草的密度为 ，鹿的数量为 ,第k年草的密度为 ，鹿的数量为 。草独立生存时，按照lgisti规律增长，则此时草的增长差分模型为，但是由于鹿对草的捕食作用，草的数量会减少,则满足如下方程：（） (1）鹿离开草无法独立生存，因此鹿独立生存时的模型为

11、 ,但是草的存在会使得鹿的死亡率得到补偿，则满足如下差分方程: ( ） (2）另外，记初始状态鹿的数量为 ,草场密度初值为 。各个参数值为:， , , , 利用tlab编程序分析计算该差分方程模型,源程序如下：%定义函数wui,实现wuiogisti综合模型的计算，计算结果返回种群量unctio b =iw(x0,y0,r,n，a，d，n） 描述dui-ost综合模型的函数 x(1) = x0; 草场密度赋初值 （) = y0; 鹿群数量赋初值for k ： n; x(k+1) = x(） r*（1-(k)/n)*(k) - *x(k)*y(k)n; y（k+) y(k)+ (-d + b*(

12、k）)*y(k); end b = x;y;%lerallc1 =diwuti(100,10,300,1.5,1.6，.9,50）；c2 = diti(3000,0,08,000,1，.6，0.9,50); = 0 ： 50;plot（,c1（1,:),b,，c1(2,:),b,k，c2(1,:),r,c2（2,:),r,),.axs（0 50 300）;xabl（时间/年)ylel（种群量草场:单位密度,鹿：头)ttle(图1草和鹿两种群数量变化对比曲线)gtext（=1000)text(x0=00)gtet(草场密度)gte(鹿群数量）比较将1只鹿放入密度为100和密度为300的草场两种情

13、况（绘制曲如图1所示):由图中可以看到,蓝色曲线代表草场密度的初始值为1000时,两种群变化情况；而红色曲线则代表草场密度的初始值为300时,两种群的变化情况。观察两种情况下曲线的演变情况,可以发现大约0年左右时间后，两种群的数量将达到稳定。使用atb计算可以得到，当 ，即两种群数量的平衡点为(00,600)。为进一步验证此结论,下面通过改变相关参数，研究两种群变化情况，找到影响平衡点的因素：（1)改变草场密度初始值；从图2中可以看到,改变草场的初始密度不会对两种群数量的平衡点造成影响。（2）改变鹿的数量初值由图2可以看到,鹿初始的数量的改变在理论上也不会改变最终种群数量的平衡值。但是,我们可

14、以看到,0200的那条曲线(紫色曲线)，在5-区间内降低到了非常小的值，这显然是不符合鹿的现实繁殖规律的，因为鹿的种群可持续繁殖的最小数量是存在域值的。当种群数量低于这个值时，在实际情况下，鹿的种群就要灭绝。同样道理，草场的密度也存在一个最小量的域值,低于这个阈值，草也将灭绝。综合上面分析，可以在此得出一个结论:最大密度一定的草场所能承载的鹿的数量存在上限。(3)改变草场的最大密度n，画图比较结果；如图4所示,如果草场密度的最大值n发生变化,则最终两种群数量的平衡点也会发生相应的变化。结论:值越大,平衡点两种群的数量就越大；n越小，平衡点两种群的数量就越小。（4)改变鹿群独立生存时的死亡率实验

15、中，改变了鹿单独生存的死亡率得到如图5.1和5.2两幅图,可以得出结论：鹿单独生存的死亡率越大，则两种群数量达到平衡点的时间越短;相反，鹿单独生存的死亡率越小，则两种群数量达到平衡点的时间越长（甚至有可能会出现分叉、混沌)。（5)草场密度对鹿数量的补偿作用变化(变化）从图中可以看到,如果b增大,则达到稳定点的时间会加长，但如果b减小则会有一个域值,当b低于域值时,草鹿种群数量的平衡时将不收敛于同一个平衡点，出现多值性。6.leli种群年龄结构的差分方程模型 已知一种昆虫每两周产卵一次,六周以后死亡(给出了变化过程的基本规律)。孵化后的幼虫周后成熟,平均产卵10个,四周龄的成虫平均产卵10个。假

16、设每个卵发育成2周龄成虫的概率为0.09,（称为成活率）,2周龄成虫发育成4周龄成虫的概率为02。假设开始时，2,4,周龄的昆虫数目相同,计算2周、4周、6周后各种周龄的昆虫数目; 讨论这种昆虫各种周龄的昆虫数目的演变趋势：各周龄的昆虫比例是否有一个稳定值?昆虫是无限地增长还是趋于灭亡？ 假设使用了除虫剂，已知使用了除虫剂后各周龄的成活率减半,问这种除虫剂是否有效？解:将两周分成一个时段，设k时段周后幼虫数量为:x1(）, 2到周虫的数量为:x2(k）, 4到6周虫数量为：3()。据题意可列出下列差分方程:1(k+)=2(k)*100+x3()*150 x2（k+1)=x1()*0.09 x(

17、1)=x（k)*0.2运用matlab编写的程序如下:funcin z=dili(a,r1，r2,n)x(1)=;(1）=a；w(）=a；for =:n x(k1)=y(k)*100+(k)*15; y(+1)=(k)r1; w(k+1)=y(k)*r2;ndz=x,，w;fr k=1:n+ =（k）y（k)+（）edplot（1：,x）;old olot(1:n+1,r)；hold npot(:n+,w,)，grid计算前三年的结果为： =diliuti（10，0.0,2,）m = 30m= 2.5021e+00m = .312e+00z = .0e+00 * 00100 0.0100 00

18、10 2.50 00 0.0200.300 .025 0.0000(蓝线为02周的虫,红线为24周的虫，黑线为46周的虫）其中，m表示三个不同生长周期的虫的总数，可见虫并未灭绝。当年份足够长时,可观察到各年龄段虫的数量变化： =iliut(10,.09，0.,20)m = 30m 2.502e+004m = 3.352e003m = 2260e+004m = 739e00 = 2.1235e+04m= 1.4867+0m .14e004m = 1914e004m= 2.38e+04m =33e+00m = 2.384e+00m 2713e+04m = 2.9975e+04 .153e+004m 3433+04 3.42e04m = 3.389+004= 4.095e04m = 4.5301+00 = 4.521e04z = 10e+04 * 0.0100 010 0.0100 250 0001 0.0 .3090 .022 0.000 2.52 0.0 005 . 0203 0.6 2.1 0.006 0041 .460 .190 .1 215 0132 0.038 1.893 .14 0.002 2.3372 00170 0.0039 .28 0.210 .00 2.136 0.005 0042 685