1、假如出现疾病或者失控的偷猎,使大象总数突然大幅度下降,即使停止避孕,总数恢复到期望值也需要很长时间,并且会对大象群的种群结构产生很大影响,对于恢复存在不良影响。最后,对所设立的方案模型通过蒙特卡罗随机模拟法进行计算机模拟,验证以上计算的理论结果,模拟结果表明,结果是合理的。问题背景:一家大型自然公园散养了大约11000头大象,管理部门希望为大象创造一个健康的生存环境,将大象的总数控制在11000头左右。每年,公园的管理人员都要统计当年大象的总数。过去20年里,公园每年都要移去一些大象,以便保持大象总数维持在11000头左右,通常都是采用捕杀或者迁移的方法来实现的。统计表明,每年约处理600-8
2、00头大象。近年来,公众强烈反对捕杀大象行为,而且即使是迁移少量的大象也是不允许的。但是一种新的给大象打避孕针的方法也被研制成功。一只成年母象打了避孕针后,两年内不再怀孕。1 建立并利用模型推算2-60岁大象可能的存活率,以及目前的大象年龄结构.2估计每年需要避孕多少大象,才能保证大象总数控制在11000头左右,说明数据不确定性对你的结论的影响。3假设每年可以移出50-300头大象,避孕大象数可以减少多少。4有一些反对观点认为,假如出现疾病或者失控的偷猎,使大象总数突然大幅度下降,即使停止避孕,也会对大象群的恢复存在不良影响,研究并回答这个问题.模型假设:1几乎没有大象迁入或迁出;2性别比接近
3、1:1,采取控制后,也维持这个比例;3初生象的性别比也是大约1:1;4母象初次怀孕大约在10-12岁,一直到60岁大约每3.5年怀胎一次,60岁后不再受孕,怀孕期为22个月,可以假设母象均在11岁怀孕,且从13岁开始生出小象。5取按年循坏的方案;6避孕针对母象没有副作用,打了避孕针的母象2年内不再受孕;7假设初生象存活到1岁的比例为75%,此后,直至60岁前,存活率都比较均匀,大约在95%以上, 大象一般只活到70岁,设其在6070之间的存活率线性递减,而70岁往后的死亡率为100%。8公园里不存在捕杀行为,偷猎可以不考虑。符号说明:1、 :第i年龄组母象个体在1个时段内平均繁殖的数量。2、
4、第i年龄组母象个体在1个时段内的存活率。3、 L : leslie 矩阵。4、 n: 移出大象的头数。5、 r : 特征值。6、 q1: 母象的总数。7、 为岁数为t的大象在第i年时的个数问题一建立并利用模型推算2-60岁大象可能的存活率,以及目前的大象年龄结构.模型1:过去两年迁出的大象时从随机抽样中来的,所以它的结构可以反应向群总体的年龄结构。将过去两年迁移出的总的大象的数目两个向量表示如下X1=103 77 71 0 2X2=98 74 690 0令x为x1与x2的和(或平均值,效果一样)X=X1+X2则x的结构即可以表示目前的大象年龄结构。将x中各值的范围控制在合理的范围内利于输出观察
5、,令y0=x/norm(x,1);利用matlab显示其年龄结构即为:a=0:70;bar(a,y0,stacked);则年龄结构如图所示:其260岁大象可能的存活率可以根据结构向量的后项与前项比得到,本题中,具体做法是,将260岁年龄的大象分为前项为29岁,1019岁,2029岁,3039岁,4049岁,5060岁,求出大致的存活率,再求出平均值,可以得到:求得ans平均值=(0.9672+0.9851+0.9962+0.9789+0.9749+0.9859)/6=0.9814所以说,260岁的存活率为98.14%,与题目中所给的95%一致。所以说,260岁的存活率为98.14%,与题目中所
6、给的95%接近。(程序代码见附录)问题2:在问题二中我们分为不考虑重复打针与考虑重复打针两种情况在这一种情况中我们建立了三个模型序号为2,3,模型2讨论了打针时区分有效年龄,模型3讨论了不区分年龄段的情况;第二种情况我们考虑了重复打针的情况模型2:1不考虑两年内被重复注射的雌性数量。(重复的稍复杂,下面再分析)2 假设打避孕针的时候能够区分有效年龄段 13611,1360岁母象生幼年母象率=1/3.5*(1+0.0135)/2=0.1448可得,leslie模型中的leslie矩阵为:=通过matlab求其特征值: l=zeros(71,71); l(2,1)=0.75;for i=14:61
7、l(1,i)=0.1448;endfor j=3:l(j,j-1)=0.9814;end; l;for k=62:71l(k,k-1)=0.9814-0.9814*(k-61)/10;eig(l)ans = 0 1.0333 如上,求得特征根为1.0333,大于1,如果不进行避孕注射,该大象种群将无限增长下去,所以要进行避孕注射。2、求避孕繁殖率根据Leslie矩阵的性质知道,要保持种群稳定,必须使得特征根为1,即使得下面式子成立: 而此题中 , 带入数据:解得:b=0.0523所以打完避孕针的繁殖率为0.05233、验证b的正确性:l=zeros(71,71);l(1,i)=0.0523;a
8、ns = 0 1.0000 0.9501 + 0.1205i 0.9501 - 0.1205i4、求生育母象的比例解得特征向量为: n1=zeros(1,71);n1(1)=1;n1(2)=0.75;for i=3:n1(i)=n1(i-1)*0.9814;for i=62:n1(i)=n1(61)*(1-(i-61)/10); n1 a=norm(n1)a = 3.8515 a=norm(n1,1) 29.3663 b=n1(:,14:61) norm(b,1) 19.1174 c=ans/ac =0.6510知:稳定后,可生孕的母象的比例为:65.10 %5、求每年需要避孕的母象数量(不考
9、虑重复打针,在有效年龄打针)由以上知道:打避孕药后leslie矩阵中第一行的所有0.1448应该替换为0.0523,而这样的调整需要对母象大避孕药后实现,设每年被打避孕针的1360岁的母象数为n。一次注射可以使得一头成熟的母象在两年内不会受孕,所以实际上每年共有2n头大象处于避孕期方案1:设此系统中1361岁避孕的母象的比例为k,1361母象总数为5500*65.10%(N),则,因为一次注射可以使得一头成熟的母象在两年内不会受孕,所以实际上每年共有2kN头大象处于避孕期。所以新的出生率应该为:=0.0523即为b可以求得:k=31.941%,每年要避孕的大象总数n为: 65.10%*5500
10、*k=1143头。方案2:具有(b- b0)繁殖率母象所生的幼象的数目应该等于注射避孕药使得母象没有繁殖幼象的数目:解之得=1143头。即每年大约需要给1143头母象大避孕针。分析:用两种不同的方案得到的结果是一致的故该数据是正确的模型3两点假设:2抽取打避孕药的大象完全是随机的,可以是任意年龄的大象 如果不区分小象大象老象,直接抽取所有母象中百分比为k的象打避孕药,母象总数q1,则,其他分析均同模型1,只有求出的n值不一样:可以解得,k=0.319,所以n=q1*0.319=1755头。2、假设我们在打避孕针时是随机的,即大象的年龄=12岁时也有打到避孕针的可能性时用利用蒙特卡罗随机模拟法来
11、模拟十次 避孕大象为1755头的结果如下第一次 b=0.0522第二次 b=0.0522第三次b=0.0525第四次 b=0.0526第五次 b=0.0523第六次 b=0.0524第七次 b=0.0524第八次 b=0.0524第九次 b=0.0526第十次 b=0.0523取十次的平均值得 b=0.0524 特征值近似为1这里的b为避孕有效年龄段的避孕率,与理论值中0.0523非常接近,说明当避孕1755头母象时能让b的值为0.0523,即可以让leslie矩阵的特征值为1,能实现让大象的总数控制在11000头左右的目标。模拟程序见附录6、(1)实行避孕政策前的年龄结构:m1=zeros(
12、1,71); m1(1)=1; m1(2)=0.75/1.0333; for i=3:m1(i)=m1(i-1)*0.98/1.0333; for i=62:m1(i)=m1(61)*(1-(i-61)/10); a=0:71; plot(a,m1,r-(2)施行避孕政策后的年龄结构: n1(2)=0.75; plot(a,n1,b-分析结果可以知道打避孕药后的大象的年龄结构中6070的数目有所增加模型4:考虑发生有雌性个体被重复注射的情况时建立模型:每个被注射的个体被重复注射的几率设为k,根据其随机性可以知道,两次注射比率为k2,则实际情况下每年的避孕比率应该为 2*k-k2 分析过程可有如
13、下饼图清晰显示:如图,假设这两幅饼图表示连续两年被打避孕药的情况,设每年被打避孕药的母象的比例为k,所以如图中右上角两块加起来为k,而红的那部分表示前一年大避孕药后,第二年又被打避孕药了,其比例为,所以黑的部分和紫的部分均占(k-)所以说,事实上每一年处于避孕状态的母象的比例为上一年的黑色部分加上这一年的紫色部分和红色部分,即为:2*(k-)+从而可以得到如下计算式0.1448*(1-2K+K2)*N/N=0.0523得到k=0.399N=11000/2*0.399=2195此时算出来n应该为:2195头。即每年大约需要给2195头母象大避孕针。有些幼象一看就能看出来,可以不给它打避孕药,但是
14、得记数目,打针总计数到2195时即可。数据不确定性对结论的影响:1、避孕针可能引起大象每个月都发情,但不受孕,而公象的数目有限,可以推断,打了避孕针后,其他母象与公象配偶受孕的几率会减小,就是说,b可能还不足0.0523。 这样以后,可以假设,b减为0.05则,将0.05带入leslie矩阵后可算得,l(1,i)=0.05;ans = 0 0.9987 0.9496 + 0.1204i其特征值为:0.99871可能会导致大象种群数最终减小。2、随着时间的增长,如果持续使用避孕药,会使象的年龄结构发生变化,幼象变少,老象变多,呈现老年化倾向,如图:幼年不能生育能生育老年年数0岁112岁13-60
15、岁61-70岁1(当前)460330871775492439371083964713420354480224504401338776644305369311470503966383323673234117366309269973938.N(稳定)37131367095398问题3问题三我们分为考虑重复打针跟不考虑重复打针两种情况求解不考虑重复打针时我们分别建立两种不同的模型,而这两种模型对问题的理解是不同的,结果也不一样。A 不考虑重复打针模型5:1、问题理解我们可以这样理解:每年移出50300头大象意思是,稳定后每年的增长率为:50/11000300/11000,在年终移出多余的300头大象
16、,刚好使得大象的总数控制在11000左右。即为:0.0045450.02727,而事实上,趋于稳定状态时,设 leslie矩阵的特征值是r,则各年龄的大象数会近似地按照r-1的比例增长,所以说,此时可得,r的范围即为:1.0045451.027272、求出繁殖率b的范围:对于此式,根据leslie矩阵的性质可知:所以有:将r的值的范围带入,可以解得,b的范围为:0.06090.1224。3、经matlab验算:l(1,i)=0.0609;ans = 1.0045验证,则有b值计算正确l(1,i)=0.1224; 1.0273 知,b值计算正确。4,移出大象后的年龄结构此时r对应的特征向量为:,
17、可以求得,当r=1.004545时,1360岁母象占总象数的63.32%,r=1.02727时,1360岁母象占大象总数的51.02%。方案1:r=1.004545时=0.0609,k=0.2897 所以n=5500*0.6332*0.2897=1009r=1.02727时=0.1224,k=0.07735n=5500*0.5102*0.07735=217n的范围为2171009得:n=1009得,n=217减少的n值可以为:134926头两种方案算得结果完全一样模型6:1 对问题的理解考虑即使移走大象后,leslie模型的特征值是始终仍为1,但其存活率因为大象的移走而不断变化,对此分析考虑。
18、2 理论计算求解避孕大象头数1当移走50头大象时,有a=norm(n1,1)a = 29.3663,1)b = 1 c=b/ac = 0.0341 则,0岁象的等效存活率降低,变为161岁大象的存活率变为:0.9769与问题2的分析相同,有:b=0.0630从1361岁的母象中大避孕针,有(1360岁母象共N头)k=0.28245,所以n=N*k=5500*0.651*0.28245=1012头 移走50头,有1012头需要避孕2、当移走300头时:同以上的方法,有,0岁大象存活率变为:0.7227,161的存活率变为:0.9541,计算得b=0.127,所以k=0.06146,最终有,n=5
19、500*0.651*0.06146=220移走300头 需要避孕220头所以应该避孕的母象头数为:2201012头。3用蒙特卡罗随机模拟法进行计算机模拟在命令方式下输入如下指令:for i=1:26z(i)=smm(50+(i-1)*10)t=50:10:300;plot(t,z)得到z =1.0e+003 *Columns 1 through 12 1.0146 0.9920 0.9670 0.9436 0.9170 0.8918 0.8651 0.8372 0.8110 0.7827 0.7533 0.7247Columns 13 through 24 0.6922 0.6607 0.62
20、96 0.5956 0.5613 0.5262 0.4947 0.4553 0.4239 0.3791 0.3439 0.3022Columns 25 through 26 0.2192 分别对应的n为2191014头与理论计算的2201012几乎一致可见两者的数据是非常精确的,模型也是正确的问题三A的总结: 之所以与模型一有较大差别是因为两者建立的基础不同,模型一种的特征值是改变的,出生率是不变的, 模型二中的特征值始终为1 而出生率则是改变的故这两种不同的方案每年需要避孕的头数是不一样的,但结果都是大象的数目不变。问题三B考虑重复打针的情况如下:模型7:1、搬运走50头大象时,r=1.00
21、4545,b=0.0609同问题二,有:0.1448*(1-2K+K2)*N/N=0.0609得到k=0.3515N=11000/2*0.3515=1933避孕大象减少数目s=2195-1933=2621933头。即每年大约需要给1933头母象大避孕针。有些幼象一看就能看出来,可以不给它打避孕药,但是得记数目,打针总计数到1933时即可。2、搬运走300头时,r=1.02727,b=0.12240.1448*(1-2K+K2)*N/N=0.1224得到k=0.0806N=11000/2*0.0806=443避孕大象减少数目s=2195-443=1752443头。即每年大约需要给443头母象大避
22、孕针。有些幼象一看就能看出来,可以不给它打避孕药,但是得记数目,打针总计数到443时即可。综上,若每年可以移出50-300头大象,避孕大象数可以减少2621752头问题四:一、求解一下当前种群分布,采取避孕措施至到稳定之间的的种群分布,注意观察种群结构的变化。1、先求经过相当长时间稳定后的种群分布,以便于分析对照n1=1n1(2)=0.75*n1(1);61;n1(i)=n1(i-1)*0.98;n1(1)=0;n1(1)=n1(1)+0.5023*n1(i);由此得到稳定时结构向量:n=1.0000 0.7500 0.7361 0.7224 0.7089 0.6957 0.6828 0.67
23、01 0.6576 0.6454 0.6334 0.6216 0.6101 0.5987 0.5876 0.5766 0.5659 0.5554 0.5451 0.5349 0.5250 0.5152 0.5056 0.4962 0.4870 0.4779 0.4690 0.4603 0.4518 0.4434 0.4351 0.4270 0.4191 0.4113 0.4036 0.3961 0.3888 0.3815 0.3744 0.3675 0.3606 0.3539 0.3473 0.3409 0.3345 0.3283 0.3222 0.3162 0.3103 0.3046 0.2989 0.2933 0.2879 0.2825 0.2773 0.2721
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1