《大学数学实验》由会员分享，可在线阅读，更多相关《大学数学实验（18页珍藏版）》请在金锄头文库上搜索。
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）鹿离开草无法独立生存，因此鹿独立生存时的模型为 ,但是草的存在会使得鹿的死亡率得到补偿，则满足如下差分方程: ( ） (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)*
