《大学数学实验》由会员分享，可在线阅读，更多相关《大学数学实验（24页珍藏版）》请在金锄头文库上搜索。
1、4.在某种环境下猫头鹰的主要食物是田鼠，设田鼠的年平均增长率为 ，猫头鹰的存在引起的田鼠增长率的减少与猫头鹰的数量成正比，比例系数为 ；猫头鹰的年平均减少率为 ；田鼠的存在引起的猫头鹰减少率的增加与田鼠的数量成正比，比例系数为 。建立差分方程模型描述田鼠和猫头鹰共处时的数量变化规律，对以下情况作图给出50年的变化过程。（1）设开始时有100只田鼠和50只猫头鹰。（2） 同上，开始时有100只田鼠 和200只猫头鹰。（3）适当改变参数 （初始值同上）（4）求差分方程的平衡点，它们稳定吗？解： 记第k代田鼠数量为xk ，第k代猫头鹰数量为yk，则可列出下列方程： xK+1=xK+(r1-yK*a1
2、)*xK K+1=yK+(-r2+xK*a2)*yK 运用matlab计算，程序如下：function z=disiti(x0,y0,a1,a2,r1,r2)x=x0;y=y0;for k=1:49 x(k+1)=x(k)+(r1-y(k)*a1)*x(k); y(k+1)=y(k)+(-r2+x(k)*a2)*y(k);endz=x,y;（1）z=disiti(100,50,0.001,0.002,0.2,0.3)plot(1:50,z(:,1);hold on; plot(1:50,z(:,2),r)得到如下结果：z = 100.0000 50.0000 115.0000 45.0000
3、132.8250 41.8500 153.8313 40.4125 178.3808 40.7221 206.7930 43.0336 239.2525 47.9216 275.6376 56.4758 315.1983 70.6668 355.9639 94.0149 393.6908 132.7422 420.1696 197.4383 421.2459 304.1220 377.3850 469.1057 275.8285 682.4408 142.7576 854.1819 49.3682 841.8092 17.6832 672.3836 9.3300 494.4483 6.5828
4、 355.3402 5.5602 253.4164 5.2632 180.2096 5.3674 128.0436 5.7536 91.0051 6.3807 64.7508 7.2437 46.1518 8.3581 32.9749 9.7541 23.6336 11.4744 17.0046 13.5742 12.2935 16.1221 8.9392 19.2024 6.5457 22.9172 4.8333 27.3899 3.6049 32.7691 2.7209 39.2338 2.0829 46.9988 1.6215 56.3224 1.2875 67.5143 1.0463
5、80.9465 0.8737 97.0651 0.7530 116.4051 0.6733 139.6077 0.6280 167.4416 0.6150 200.8269 0.6364 240.8645 0.7011 288.8685 0.8286 346.4029 1.0587 415.3167 1.4745 497.7677 2.2570 红线为猫头鹰数量曲线,蓝线为田鼠曲线(2) z=disiti(100,200,0.001,0.002,0.2,0.3)plot(1:50,z(:,1);hold on; plot(1:50,z(:,2),r)z = 100.0000 200.0000
6、100.0000 180.0000 102.0000 162.0000 105.8760 146.4480 111.5459 133.5243 118.9610 123.2551 128.0906 115.6037 138.9010 110.5381 151.3273 108.0844 165.2367 108.3713 180.3771 111.6737 196.3091 118.4584 212.3165 129.4298 227.2997 145.5610 239.6737 168.0647 247.3278 198.2066 247.7713 236.7886 238.6561 283
7、.0909 218.8260 333.2864 189.6595 379.1639 155.6793 409.2388 123.1052 413.8872 96.7745 391.6244 78.2302 349.9356 66.5007 299.7060 59.8702 249.6555 56.8973 204.6527 56.6326 166.5452 58.5272 135.4454 62.3054 110.6663 67.8714 91.2566 75.2519 76.2671 84.5631 64.8654 95.9905 56.3762 109.7770 50.2865 126.2
8、121 46.2412 145.6183 44.0412 168.3287 43.6553 194.6461 45.2556 224.7665 49.2965 258.6395 56.6680 295.7109 68.9807 334.4547 89.0832 371.5513 121.9469 400.5521 175.9818 410.1726 264.1671 383.8530 401.6251 306.4586 589.4677 187.1029 773.9222 79.7204 831.3517 红线为猫头鹰数量曲线,蓝线为田鼠曲线(3)当a1,a2分别取0.002,0.002时，得
9、到如下图像：可见，当a1,a2参数在一定范围内改变时，猫头鹰与田鼠数量在一定范围内震荡，且不灭绝。(4)令xK=xK+1=x; yK=yK+1=y解方程得到如下结果：x=150y=200经matlab验证如下：z=disiti(150,200,0.001,0.002,0.2,0.3)plot(1:50,z(:,1);hold on; plot(1:50,z(:,2),r)z = 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 2
10、00 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 150 200 由此可知：平衡点为
11、：x=150 y=2005.研究将鹿群放入草场后草和鹿两种群的相互作用。草的生长遵从Logistic规律，年固有增长率0.8，最大密度为3000（密度单位），在草最茂盛时每只鹿每年可吃掉1.6（密度单位）的草。若没有草，鹿群的年死亡率高达0.9，而草的存在可使鹿的死亡得以补偿，在草最茂盛时补偿率为1.5。作一些简化假设，用差分方程模型描述草和鹿两种群数量的变化过程，就以下情况进行讨论：（1）比较将100只鹿放入密度为1000和密度为3000的草场两种情况。（2）适当改变参数，观察变化趋势。解：设1草独立生存，独立生存规律遵从Logistic规律；2草场上除了鹿以外，没有其他以草为食的生物；3鹿
12、无法独立生存。没有草的情况下，鹿的年死亡率一定；4假定草对鹿的补偿率是草场密度的线性函数；5每只鹿每年的食草能力是草场密度的线性函数。记草的固有增长率为r，草的最大密度为N，鹿独立生存时的年死亡率为d，草最茂盛时鹿的食草能力为a，草对鹿的年补偿作用为b；第k1年草的密度为 ，鹿的数量为 ，第k年草的密度为 ，鹿的数量为 。草独立生存时，按照Logistic规律增长，则此时草的增长差分模型为 ，但是由于鹿对草的捕食作用，草的数量会减少，则满足如下方程： （ ） （1）鹿离开草无法独立生存，因此鹿独立生存时的模型为 ，但是草的存在会使得鹿的死亡率得到补偿，则满足如下差分方程： （ ） （2）另外，记初始状态鹿的数量为 ，草场密度初值为 。各个参数值为：， ， ， ， 利用MATLAB编程序分析计算该差分方程模型，源程序如下：%定义函数diwuti，实现diwuti-Logistic综合模型的计算，计算结果返回种群量function B =diwuti(x0,y0,r,N,b,a,d,n) % 描述diwuti-Logistic综合模型的函数 x(1) = x0; % 草场密度赋初值 y(1) = y0; % 鹿群数量赋初值 for k = 1 : n;
