数学实验2优选.docx
《数学实验2优选.docx》由会员分享,可在线阅读,更多相关《数学实验2优选.docx(26页珍藏版)》请在冰豆网上搜索。
数学实验2优选
数学实验.实验2
4.某日凌晨一住所发生一件凶杀案,警方于6时到达现场后测得尸温26℃,室温17℃,2小时后尸温下降了3℃,试根据冷却定理建立差分方程,估计凶杀案发生的时间(可设正常体温为37℃)
解:
设t0为正常体温,y为死亡后人的体温,x为时间。
差分方程由冷却定理得。
冷却定律:
T=C+(T0-C)*e^-kt
C是室温,T0是人体常温37℃,T为死亡后第t小时的体温,k为可求常数。
程序:
函数
functiony=swsj1(x)
%带参数求得k=0.2027,令a=k
t0=37;
a=0.2027
y=17+(t0-17)*exp(-a*x)
end
n=6;
x=0:
n;
y=swsj1(x)
plot(x,y','r*')
gtext('死亡时间'),
运行结果:
a=
0.2027
y=
37.000033.330530.334227.887725.890024.258922.9271
y=
37.000033.330530.334227.887725.890024.258922.9271
由结果知此人的死亡时间为2 :
00.
5.据报道,某种山猫在较好、中等及较差的自然环境下,年平均增长率分别是1.68%,0.55%和-4.50%,假定开始时有100只山猫,按以下情况讨论山猫数量逐年变化过程及趋势:
(1)3种自然环境下25年的变化过程(作图);
(2)如果每年捕获三只,会发生什么情况?
山猫会灭绝吗?
如果每年捕获一只呢?
(3)在较差的自然环境下,如果想使山猫数量稳定在60只左右,每年需要人工繁殖多少只?
(1)3种自然环境下25年的变化过程(作图)
程序:
functionx=cat(x0,n,r)
x=x0;
fork=1:
n
x(k+1)=(1+r)*x(k);
end
k=(0:
25)';
y1=cat(100,25,0.0168);
y2=cat(100,25,0.0055);
y3=cat(100,25,-0.0450);
round([k,y1',y2',y3'])
plot(k,y1,k,y2,'r:
',k,y3,'g--'),grid
gtext('较好环境'),gtext('中等环境'),gtext('较差环境')
结果:
ans=
0100100100
110210196
210310191
310510287
410710283
510910379
611110376
711210472
811410469
911610566
1011810663
1112010660
1212210758
1312410755
1412610852
1512810950
1613110948
1713311046
1813511044
1913711142
2014011240
2114211238
2214411336
2314711335
2414911433
2515211532
(2)如果每年捕获三只,会发生什么情况?
山猫会灭绝吗?
如果每年捕获一只呢?
(i)每年捕获三只
程序:
functionx=cat(x0,n,r)
x=x0;
fork=1:
n
x(k+1)=(1+r)*x(k)-3;
end
command窗口程序同第一问
结果:
ans=
0100100100
1999893
2979585
3969378
4959072
5938866
6928560
7908354
8898049
9877743
10867539
11847234
12837029
13816725
14796421
15786217
16765913
17745610
1873546
1971513
2069480
216746-3
226543-6
236340-9
246137-11
255935-14
(ii)每年捕获一只
程序:
functionx=cat(x0,n,r)
x=x0;
fork=1:
n
x(k+1)=(1+r)*x(k)-1;
end
command窗口程序同第一问
结果:
ans=
0100100100
110110095
21019989
31029984
41039879
51049875
61049770
71059766
81069662
91079659
101079555
111089551
121099448
131109445
141119342
151119339
161129236
171139234
181149231
191159129
201169126
211179024
221189022
231198920
241208818
251218816
答:
若每年捕获三只,无论在什么环境情况下,均会灭绝;若每年捕获一只,则在中等和较差的环境下会灭绝.
(3)在较差的自然环境下,如果想使山猫数量稳定在60只左右,每年需要人工繁殖多少只?
程序:
symsb
solve(60-(1-0.0450)*60-b)
结果:
ans=
27/10
答:
在较差的自然环境下,如果想使山猫数量稳定在60只左右,每年需要人工繁殖3只.
8.在某种环境下猫头鹰的主要食物来源是田鼠,设田鼠的年平均增长率为r1,猫头鹰的存在引起的田鼠增长率的减少与猫头鹰数量成正比,比例系数a1;猫头鹰的年平均减少率为r2;田鼠的存在引起的猫头鹰减少率增加与田鼠的数量成正比,比例系数为a2,建立差分方程模型描述田鼠和猫头鹰共处时数量变化规律,对以下情况作图给出50年的变化过程。
解:
设x为田鼠数,y为猫头鹰数。
程序如下:
函数:
functionz=ms(x0,y0,r1,r2,a1,a2,n)
x=x0;
y=y0;
fork=1:
n
x(k+1)=(1+r1-a1*y(k))*x(k);
y(k+1)=(1-r2+a2*x(k))*y(k);
end
z=[x',y'];
运行结果;
ans=
37.2001
z=
2
y=
0
37.2001
33.5035
30.4834
28.0159
26.0000
(1)设r1=0.2,r2=0.3,a1=0.001,a2=0.002,开始时有100只田鼠和50只猫头鹰。
z=ms(100,50,0.2,0.3,0.001,0.002,50);
k=(0:
50)';
[k,z]
plot(k,z(:
1),k,z(:
2)),grid,
gtext('x(k)'),gtext('y(k)'),
(2)r1,r2,a1,a2同上,开始时有100只田鼠和200只猫头鹰。
z=ms(100,200,0.2,0.3,0.001,0.002,50);
k=(0:
50)';
[k,z]
plot(k,z(:
1),k,z(:
2)),grid,
gtext('x(k)'),gtext('y(k)'),
(3)适当改变参数a1,a2(初始值同上)。
令a1=0.003,a2=0.001,有如下数据图像
z=ms(100,200,0.2,0.3,0.003,0.001,50);
k=(0:
50)';
[k,z]
plot(k,z(:
1),k,z(:
2)),grid,
gtext('x(k)'),gtext('y(k)'),
(4)求差分方程的平衡点,它们稳定吗?
令x(k)=x(k+1)=x,y(k)=y(k+1)=y,得平衡点是(150,200)。
但是由图像可以看出它们并不稳定。
9.研究将鹿群放入草场后草和鹿两种群的相互作用,草的生长遵从Logistic规律,年固有增长率0.8,最大密度为3000(密度单位),在草最茂盛时每只鹿每年可吃掉1.6(密度单位)的草,若没有草,鹿群的年死亡率高达0.9,而草的存在可使鹿的死亡的以补偿,在草最茂盛时补偿率为1.5,作出一些简化假设,用差分方程模型描述草和鹿两种群数量的变化过程,就一下情况进行讨论,
(1)比较将100只鹿放入密度为1000和密度3000的草场两种情况。
(2)适当改变参数,观察变化趋势。
解:
设x为鹿群数,y为草的数量,H为草场的最大密度。
(1)
程序:
functionz=chl(x0,y0,n)
H=3000;
x=zeros(1,n);y=x;
a=x;b=x;
a
(1)=1.6;b
(1)=1.5;
x
(1)=x0;
y
(1)=y0;
fork=1:
n
y(k+1)=y(k)+0.8*(1-y(k)/H)*y(k)-a(k)*x(k);
x(k+1)=x(k)+(b(k)-0.9)*x(k);
a(k+1)=a
(1)*(y(k+1)/H);
b(k+1)=b
(1)*(y(k+1)/H);
end
z=[x',y'];
end
z1=chl(100,3000,60)
z2=chl(100,1000,60)
k=0:
60;
plot(k,z1(:
1),'b',k,z1(:
2),'b',k,z2(:
1),'r',k,z2(:
2),'r'),grid,
gtext('x,“鹿群”'),gtext('y“草”'),
运行结果:
z1=
1.0e+003*
0.10003.0000
0.16002.8400
0.24322.7188
0.35492.5700
0.49162.3782
0.63372.1490
0.74431.9104
0.78541.7071
0.74891.5806
0.66681.5476
0.58261.5966
0.52341.6980
0.49671.8136
0.50011.9070
0.52681.9542
0.56741.9501
0.61001.9059
0.64231.8419
0.65581.7797
0.64911.7364
0.62851.7203
0.60351.7307
0.58261.7595
0.57081.7949
0.56931.8253
0.57651.8429
0.58891.8449
0.60211.8338
0.61221.8152
0.61691.7960
0.61571.7817
0.61001.7755
0.60261.7776
0.59581.7858
0.59161.7965
0.59061.8062
0.59241.8123
0.59611.8137
0.60011.8109
0.60341.8055
0.60511.7996
0.60491.7949
0.60341.7926
0.60121.7929
0.59901.7952
0.59761.7984
0.59711.8015
0.59761.8035
0.59861.8042
0.59991.8035
0.60091.8019
0.60151.8001
0.60151.7986
0.60111.7978
0.60051.7978
0.59981.7984
0.59931.7994
0.59911.8003
0.59921.8010
0.59951.8013
0.59991.8011
z2=
1.0e+003*
0.10001.0000
0.16001.3733
0.12591.8519
0.12912.2945
0.16112.5682
0.22292.6433
0.31692.5805
0.44062.4330
0.58002.2292
0.70451.9978
0.77421.7811
0.76681.6246
0.69961.5560
0.61431.5746
0.54501.6573
0.50611.7689
0.49831.8721
0.51631.9377
0.55181.9531
0.59401.9236
0.63071.8663
0.65161.8027
0.65251.7518
0.63681.7252
0.61301.7258
0.59021.7480
0.57491.7813
0.56951.8141
0.57351.8368
0.58411.8447
0.59711.8384
0.60861.8224
0.61541.8032
0.61641.7868
0.61231.7775
0.60541.7765
0.59831.7825
0.59311.7924
0.59081.8027
0.59161.8102
0.59461.8134
0.59861.8121
0.60221.8076
0.60451.8018
0.60511.7966
0.60401.7934
0.60201.7927
0.59981.7942
0.59811.7972
0.59731.8003
0.59741.8028
0.59821.8040
0.59941.8038
0.60051.8026
0.60131.8008
0.60151.7992
0.60131.7981
0.60071.7978
0.60001.7981
0.59951.7990
0.59921.8000
由图可知,鹿群与草的数量符合寄生模型,最终鹿和草的数量都趋于一个稳定的值,由差分方程(y(k+1)=y(k)+0.8*(1-y(k)/H)*y(k)-a
(1)*(y(k+1)/H)*x(k);
x(k+1)=x(k)+(b
(1)*(y(k+1)/H)-0.9)*x(k);
令x=x(i),y=y(i),i=1,2…n)
可算出稳定值约为(600,1800)。
(2)
functionz=chl(x0,y0,a1,b1,c,d,n)
H=3000;
x=zeros(1,n);y=x;
a=x;b=x;
a
(1)=a1;b
(1)=b1;
x
(1)=x0;
y
(1)=y0;
fork=1:
n
y(k+1)=y(k)+c*(1-y(k)/H)*y(k)-a(k)*x(k);
x(k+1)=x(k)+(b(k)-d)*x(k);
a(k+1)=a
(1)*(y(k+1)/H);
b(k+1)=b
(1)*(y(k+1)/H);
end
z=[x',y'];
end
z1=chl1(100,3000,1.6,1.8,0.8,0.9,60)
z2=chl1(100,1000,1.6,1.8,0.8,0.9,60)
k=0:
60;
plot(k,z1(:
1),'b',k,z1(:
2),'b',k,z2(:
1),'r',k,z2(:
2),'r'),grid,
gtext('x“鹿群”'),gtext('y“草”'),
运行结果:
z1=
1.0e+003*
0.10003.0000
0.19002.8400
0.34282.6734
0.58412.4175
0.90562.0400
1.19901.5769
1.25431.1669
1.00370.9567
0.67650.9659
0.45971.1413
0.36081.4272
0.34501.7512
0.39702.0121
0.51902.1162
0.71082.0292
0.93651.7852
1.09681.4718
1.07831.2106
0.89111.0921
0.67301.1287
0.52311.2868
0.45621.5157
0.46051.7469
0.52871.9016
0.65611.9224
0.82241.8022
0.97151.5874
1.02241.3629
0.93831.2147
0.77771.1851
0.63081.2671
0.54271.4264
0.51871.6121
0.55361.7628
0.64091.8239
0.76541.7725
0.89061.6291
0.95961.4509
0.93131.3077
0.82381.2483
0.69941.2829
0.60831.3918
0.56881.5371
0.58151.6704
0.64101.7446
0.73501.7323
0.83751.6388
0.90721.5017
0.90811.3751
0.84011.3049
0.74181.3101
0.65721.3822
0.61081.4940
0.60861.6073
0.64781.6825
0.71871.6924
0.80171.6338
0.86601.5305
0.88191.4233
0.84131.3523
0.76671.3398
z2=
1.0e+003*
0.10001.0000
0.19001.3733
0.17561.8299
0.21032.2295
0.30242.4375
0.47252.4101
0.73042.1819
1.02931.8079
1.21951.3902
1.13911.0828
0.85400.9786
0.58681.0604
0.43201.2770
0.37421.5695
0.38981.8550
0.47282.0357
0.62482.0458
0.82941.8846
1.02091.6115
1.08911.3308
0.97861.1501
0.77311.1172
0.59561.2175
0.49461.4095
0.46781.6355
0.50581.8226
0.60371.9032
0.74971.8471
0.90581.6764
1.00171.4582
0.97661.2787
0.84691.1996
0.69431.2337
0.58341.3580
0.53361.5301
0.54331.6944
0.60661.7934
0.71341.7902
0.83761.6866
0.93141.5239
0.94471.3667
0.86921.2734
0.75101.2694
0.64711.3468
0.58761.4757
0.57901.6131
0.61831.7116
0.69681.7352
0.79511.6756
0.87891.5568
0.90891.4262
0.86861.3334
0.78181.3083
0.69191.3530
0.63081.4480
0.61121.5601
0.63321.6506
0.69041.6871
0.76791.6566
0.84011.5716
0.87611.4661
z1=chl1(100,3000,1,1.5,0.8,0.85,60)
z2=chl1(100,1000,1,1.5,0.8,0.85,60)
k=0:
60;
plot(k,z1(:
1),'b',k,z1(:
2),'b',k,z2(:
1),'r',k,z2(:
2),'r'),grid,
gtext('x“鹿群”'),gtext('y“草”'),
改变参数可看出a1明显影响着两者数量的平衡性;b1明显影响其震荡的频率及次数。
最新文件----------------仅供参考--------------------已改成-----------word文本---------------------方便更改