大象群落的稳定发展

3、周晓阳 数学实验与Matlab [M].武汉：华中科技大学出版社 2002 4、郑 谏 当代数学的若干理论与方法 [M].上海：华东理工大学出版社 2002 5、李尚志 数学建模竞赛教程 [M].江苏：江苏教育出版社 1996 6、北峰数模网 http://mcm.zjnu.net.cn/ 2006-8-13

附录：

1、计算机模拟Matlab程序代码 function [bys,b0]=moni(n)

b=[12 16 10 12 17 14]; c=b/sum(b); b=cumsum(c)'; a=rand(n,1);

k1=0;k2=0;k3=0;k4=0;k5=0;k6=0; i=1;

for i=1:n

if(a(i)<=b(1)) k1=k1+1;

elseif(a(i)>b(1)&a(i)<=b(2)) k2=k2+1;

elseif(a(i)>b(2)&a(i)<=b(3)) k3=k3+1;

elseif(a(i)>b(3)&a(i)<=b(4)) k4=k4+1;

elseif(a(i)>b(4)&a(i)<=b(5)) k5=k5+1;

else(a(i)>b(5)&a(i)<=b(6)) k6=k6+1; end end

s10=1-(11000*0.12*(1-0.9897)+k1)/(11000*0.12); s20=1-(11000*0.16*(1-0.9897)+k2)/(11000*0.16); s30=1-(11000*0.1*(1-0.9897)+k3)/(11000*0.1); s40=1-(11000*0.12*(1-0.9897)+k4)/(11000*0.12); s50=1-(11000*0.17*(1-0.9897)+k5)/(11000*0.17); s60=1-(11000*0.14*(1-0.9897)+k6)/(11000*0.14);

b0=1/(s10^10*(1+s20+s20^2+s20^3+s20^4+s20^5+s20^6+s20^7+s20^8+s20^9)+s10^10*s20^10*(1+s30+s30^2+s30^3+s30^4+s30^5+s30^6+s30^7+s30^8+s30^9)+...

13

s10^10*s20^10*s30^10*(1+s40+s40^2+s40^3+s40^4+s40^5+s40^6+s40^7+s40^8+s40^9)+s10^10*s20^10*s30^10*s40^10*(1+s50+s50^2+s50^3+s50^4+s50^5+s50^6+s50^7+s50^8+s50^9)+...

s10^10*s20^10*s30^10*s40^10*s50^10*(1+s60+s60^2+s60^3+s60^4+s60^5+s60^6+s60^7+s60^8+s60^9))*4/3;

bys=0.5*0.85*8864*(0.1448-b0)/(2*0.1448);

2、求大象年龄各个年龄头数及饼图代码

a1=[ 0 0 0 3 4 7 20 9 15 9 22 3 23 5 13 21 0 22 14 5 13 10...

0 13 30 14 12 0 20 6 3 5 8 12 10 3 7 14 10 16 21 13 10 12...

6 3 6 9 13 10 3 6 21 15 4 13 10 32 14 0];

a2=[0 20 21 13 12 13 22 14 40 14 26 13 14 27 3 14 12 20 25 17 14 10...

0 2 3 4 4 3 2 3 13 16 13 10 10 12 16 12 10 12 19 13 24 17...

16 25 12 45 23 34 13 16 10 17 13 13 12 3 22 20];

a11=sum(a1([1:10]))/sum(a1,2); a12=sum(a1([11:20]))/sum(a1,2); a13=sum(a1([21:30]))/sum(a1,2); a14=sum(a1([31:40]))/sum(a1,2); a15=sum(a1([41:50]))/sum(a1,2); a16=sum(a1([51:60]))/sum(a1,2); a21=sum(a2([1:10]))/sum(a2,2); a22=sum(a2([11:20]))/sum(a2,2); a23=sum(a2([21:30]))/sum(a2,2); a24=sum(a2([31:40]))/sum(a2,2); a25=sum(a2([41:50]))/sum(a2,2); a26=sum(a2([51:60]))/sum(a2,2); c1=(a11+a21)/2*8865 c2=(a12+a22)/2*8865 c3=(a13+a23)/2*8865 c4=(a14+a24)/2*8865 c5=(a15+a25)/2*8865 c6=(a16+a26)/2*8865

d=[1091 c1 c2 c3 c4 c5 c6 1044]; pie(d);

14

附录：有关数据

表1：该动物的年龄和性别统计表

前一年的情况 该动年龄 物头数 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

前两年的情况 该动雌性年龄 物头头数 数 0 2 5 3 4 1 3 2 0 12 11 4 6 4 4 2 0 3 2 4 1 4 11 6 1 4 5 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 15

该动雌性年龄 物头头数 数 0 0 0 2 1 2 10 2 7 3 12 1 13 2 6 10 0 12 6 4 7 5 0 5 12 6 5 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 3 5 8 12 10 3 7 14 10 16 21 13 10 12 6 3 6 9 13 10 3 6 21 15 4 13 10 该动雌性年龄 物头头数 数 0 10 12 3 5 4 4 7 21 8 10 10 4 12 1 3 8 10 11 14 10 7 0 0 0 2 2 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 13 16 13 10 10 12 16 12 10 12 19 13 24 17 16 25 12 45 23 34 13 16 10 17 13 13 12 雌性头数 7 12 6 6 0 8 5 2 4 6 10 7 10 10 4 12 3 23 12 10 9 4 4 7 3 6 3 0 0 0 3 4 7 20 9 15 9 22 3 23 5 13 21 0 22 14 5 13 10 0 13 30 14 12 0 20 21 13 12 13 22 14 40 14 26 13 14 27 3 14 12 20 25 17 14 10 0 2 3 4 4

28 29 30

0 20 6 0 10 5 58 59 60 32 14 0 12 8 0 28 29 30 3 2 3 1 1 0 58 59 60 3 22 20 2 11 10 16