大象群落的稳定发展

大象群落的稳定发展摘要本文根据非洲某国的国家公园近两年从公园运出的大象的大致年龄和性别的统计情况，探讨大象的合理的存活率并推测当前的年龄结构，针对不同情况给出如何进行避孕注射以达到控制大象数量的目的。首先，充分利用给出的近两年来运出的大象的数量与性别统计表，分析近两年来的大象群落的情况，建立一个线性方程组的数学模型，通过求解方程组得到年龄在2 岁到 60 岁之间的大象的总数，并且求出了存活率为：98. 9718%;因为假设公园2岁到60岁 之间的大象占总大象的比例等于运出的2 岁到 60 岁之间的大象占总移出大象的比例，所以通过一些比例之间的关系得到这个大象群落的当前的年龄结构（见表1） 。然后，建立一个按年龄分组的种群增长的差分方程模型，运用第一问求出的各年龄段大象的存活率以及繁殖率，求解当前大象群落对应的Leslie 矩阵的特征根，发现该 特征根大于1,根据Leslie矩阵的稳定性理论知道：如果不进行避孕注射该大象种群将无限增长（如果环境允许）；据此，利用 Leslie 矩阵稳定的充要条件求出应该保持多大的繁殖率才能使种群保持稳定，求解的主要思路是：特征根取为1、把繁殖率当成未知数，将此时的各年龄段的存活率代入方程VI, 求解这个以繁殖率为未知数的方程可以得到要使种群保持稳定繁殖率的取值；根据需要避孕掉母象所生的幼象的数目等于注射避孕药使得母象没有繁殖幼象的数目这一条件建立一个方程，最后求得每年注射避孕药的母象头数为: 1393 （头）。最后，假设被转移的大象只考虑处于1 一 60 岁之间，这样可以认为转移后的大象看成每年多死了这么多头大象，即意味着死亡率将增加，存活率将减少；仍然按照解决第二问的模型，只需将此时不同的各年龄段大象的存活率代入那个以繁殖率为未知数的方程（方程 VI ） ,求出应该保持多大的繁殖率才能使种群保持稳定。考虑到求解的数据比较多，采取计算机模拟的方法来确定移出大象后所需要进行避孕的母象头数（见表 2）,为了检验计算机模拟的正确性，用理论去验证。模拟的思路方法见计算机模拟流程图一图2O关键词：关键字：线性方程组、差分方程模型、Leslie 矩阵、计算机模拟问题重述位于非洲某国的国家公园中栖息着近11000 头大象。管理者要求有一个健康稳定的环境以便维持这个11000 头大象的稳定群落。管理者逐年统计了大象的数量，发现在过去的 20年中，整个大象群经过一些偷猎枪杀以及转移到外地还能保持在11000 头的数 量，而其中每年大约有近 600 头到 800 头是被转移的。由于近年来，偷猎被禁止，而且每年要转移这些大象也比较困难，现决定采取避孕注射法以维持大象数量的平衡。我们已知此公园近两年从这个地区运出的大象的大致年龄和性别的统计。根据这些信息我们需要解决以下问题：1. 探讨年龄在2 岁到 60 岁之间的象的合理的存活率的模型，推测这个大象群落的当前的年龄结构。2. 估计每年有多少母象要注射避孕药，可以使象群固定在11000头左右。这里不免有些不确定性，是否能估计这种不确定性的影响。3. 假如每年转移50至 300头象到别处，那么上面的避孕措施将可以有怎样的改变？问题假设1 、 假设大象的性别比近似认为1: 1,并且采用措施维持这个性别比；2、 假设母象可以怀孕的年龄为11 岁一 60岁、最高年龄为70岁， 70岁的死亡率为100%, 并且 6170 岁的大象的头数呈线性递减；3、 假设大象在各年龄段中的分布率不变，即年龄结构不变，并采用各种措施维持这一结构；4、 假设被转移的大象只考虑处于1 60岁之间，转移后的大象看成每年多死了这么多头大象；5、 假设0岁大象能够活到1 岁的比例为75%；符号说明Xj ：表示一年象的头数（匚 。表示。岁大象的头数，i=l表示160岁大象头数，i=2表示61 70岁大象的头数）；Pi： 表示存活率（几表示0岁大象的存活率，门表示1 60岁大象的存活率，必表示 61 岁一 70 岁大象的存活率）；讹）：表示时段k 第 i 年龄组的大象数量；$ : 第 i 年龄组每个（母象）个体在1 个时段平均繁殖的数量；? ：第 i 年龄组的存活率；L : Leslie 矩阵；A : L 矩阵的那个唯一正特征根；n : 表示移出大象的头数；问题分析对于问题一，利用绐出的近两年来运出的大象的数量与性别统计表，可以分析近两年来的大象群落的情况，比如移出的各个年龄段的大象占移出的总的大象的头数的比例是多少，还可以根据两年移出大象后大象总数都是11000 来建立方程，用于求解存活率。对于问题二，因为考虑的是公园在未来很长一段时间的大象种群控制问题，所以可以建立一个按年龄分组的种群增长的差分方程模型，根据差分方程的Leslie 矩阵的特征根，结合Leslie 矩阵的稳定性理论对当前大象种群的情况进行分析。为了保持大象种群的稳定，必须使得Leslie 矩阵的最大特征根为1,而这样， 特征根取为1、 把繁殖 率当成未知数，将此时的各年龄段的存活率代入方程特征方程，求解这个以繁殖率为未知数的方程可以得到要使种群保持稳定繁殖率的取值；根据需要避孕掉母象所生的幼象的数目等于注射避孕药使得母象没有繁殖幼象的数目这一条件建立方程来求解应该对多少头母象进行避孕。对于问题三，由于假设被转移的大象只考虑处于1 60 岁之间，故可以认为转移后的大象看成每年多死了这么多头大象，即意味着死亡率将增加，存活率将减少。按照解决第二问的模型，只需将此时不同的各年龄段大象的存活率代替原来的存活率，就可以求出此时应该保持多大的繁殖率才能使种群保持稳定。为了方便，可用采用计算机模拟的方法来确定移出的大象在哪个年龄段，考虑到计算机模拟的不确定性，可以对模拟结果进行检验。探讨大象的存活率和当前大象的年龄结构下面将根据给出的近两年来运出的大象的数量与性别统计表，分析近两年来的大象群落的情况， 建立一个线性方程组数学模型，通过求解方程组得到年龄在2 岁到 60 岁 之间的大象的存活率，并给出大象各年龄所占的比例，进而得到这个大象群落的当前的年龄结构。1 、线性方程组模型的建立(1) 首先，计算一年象的头数。大象群是由0岁，1 60岁， 61 岁一 70岁组成，且稳定在11000头。设 0岁的头数为 Xo, 1-60岁大象头数为Xi, 61岁一 70岁大象头数为X?。所以得到第一个方程：Xo+X I+X 2= 11000(I)(2) 其次，考虑到前一年大象的总数等于前两年存活下来的大象加上新生的幼儿再减去运出的大象数。设。岁大象的存活率为几，1-60岁大象的存活率为广，61岁 70岁大象的存活 率为“2则经过一年后，新生的大象存活下来的头数为Xoxpo； 1 到 60岁的大象存活下来的头数为X,xp 1； 61 岁一一 70岁的大象能存活下来的头数为X2X/A,因此得到第二 个方程：(Xox 几 +X 2xp,) + Xo-622= 11000(II)联立(I)、 (II) 得到方程组：Xo+Xj+X,=11000* / (*)Xo x p0+ X)x Pi +X2x p2+ Xo-622=110002、模型的求解根据近两年来运出的大象的数量与性别统计表，得到如下分析结果：(1) 计算 0 岁的大象头数由表中统计，1 岁一 10岁的大象占1 岁一 60岁的大象比例为：(67/622+169/876)/2=15. 032%所以得到：11 岁一 60岁能生小象的母象占1 岁一 60岁的大象比例为：(1-15.032%) x0. 5 =42.48%因为能生小象的母象每3.5年生一头小象，且双胞胎的机会为1.35%,相当于每年生 0. 2896 头，所以0岁的大象占1 岁一 60岁的大象比例为：0. 4248x0. 2896=0. 12303这样 0 岁的大象共有：Xo=0. 12303x X,(III)(2) 计算 61 岁一 70 岁的大象头数从表中计算运出的59 岁的大象占运出的总大象比率为：(14/622+22/876) /2=0. 0238由于运出的大象都是1 岁一 60 岁的，所以0. 0238 也可看为59 岁的大象占1 60岁的大象的头数比例，得到60岁的大象占的比例为0. 0238X/7,由假设可以知道：61 岁一 70 岁的大象头数为：X2=l/2Xl0X0. 0238XpXX( IV)61 岁一一 70 岁的大象经过一年能存活下来的头数为：X2 X p2 =( 1/2) X 9 X 0.0238 X pt X X,(V)、将(III)、(V)和(IV)两个式子代入上面方程组(*)得：0.12303xX.+Xj+C 1/2 )x10x0.0238 x p.xXNlOOO0.12303XX X p0+ Xl Xp、 + (1/2)X9X0.0238Xp, XX,+0.12303X X-622= 11000又由假设知道，。岁大象的存活率为P0=75%代入上述方程组，然后用 Mathematica解之 得：p, =0.989719 w,?ww, nu Li就3、结果分析（1） 由结果可以知道，2-60岁大象的存活率为9& 9718%,这与题目给出的大于 95%是相一致的，所以可以认为结果是合理的；（2） 从图1可以看出，各个年龄段的大象所占的比例基本上是一样的，21-30岁 和41-50岁的大象比例相对比较大，因为这段大象正处于年龄的黄金时期。由此，可以认为求出的大象年龄也是合理的。估计每年注射避孕药的母象头数为了估计每年注射避孕药的母象头数，首先建立一个按年龄分组的种群增长的差分方程模型；然后用Leslie矩阵稳定的充要条件分析如果不进行避孕注射种群的增长情况；最后仍然利用Leslie矩阵稳定的充要条件求出应该保持多大的繁殖率才能使种群保持稳定，进而利用一个方程求出每年注射避孕药的母象头数。1、按年龄分组的种群增长的差分方程模型的建立记时段k第i年龄组的大象数M为 兀为，k二0,1,2, ；- i=l, 2,门,第i年龄 组的繁殖率为也，即第i年龄组每个（母象）个体在 1个时段平均繁殖的数第i年 龄的存活率为，我们这里假设切和耳不随时段k变化，在稳定的环境下这个假设是合理的。勺和耳可由统计资料获得。舌伙）的变化规律由以下的基本事实得到：时段 k+1第 1年龄组种群数M是时段 k各年龄组繁殖数M之和，即X伙+1)=工g伙)/-1时段k+1第i+1年龄组的种群数M是时段k第 兀年龄组存活下来的数星，即?+（k +1） = s （ k） J = 1,2,? ?,舁一 1 记时段k种群按年龄组的分布向M为x伙）=片伙）*2伙），兀伙）由繁殖率勺和存活率构成的矩阵为5,b910 %500.0000图.00000.00000$9 0000.00? ? ? ? ?.00几? ? ? ? ?.000A600000000%” 69”700.000.000.000.000.000.000.0000.000.000.000.000 00056000 000 00.000.00阳.000.0000.00000$690根据Leslie矩阵的性质可以得到如下的定理定理1: L矩阵有唯一的正特征根入,且它是单重根的，入对应正特征向MA 1x71- 00.00.750.000.9897.00.00.0.98970.0? ? ? .0? ? ? . 000.000.000.000.000.00.1448 1448 0.144800 000 000 000 00.98970.000.9897.000 000 000 0.989700 000 000 00.8907.0.0.0.0.0.0.0.0.0.0.0.0. 0.89070_这是一个71x71的矩阵。如果不进行避孕注射，该大象种群将无限增长下去(如果环(2)讨论A呦的特征根境允许)，所以要进行避孕注射。r2?L矩阵的其他n-l个特征根人都是满足闲5人*=2,3,中 该定理表明L矩阵的特征方程入”(片九+ b-AsAA91 _+?+ 仇- 比为一翻+) = 0只有一个正根，并且易知，Lx = 2,/2、如果不进行避孕注射种群的增长情况(1)建立Leslie矩阵首先，由第一问的求解知道，0岁的大象的存活率为0.75; 1-60岁大象的存活率 为0.989718;根据假设61-70岁大象头数是线性递减的，而且到70岁所有的大象都死 完了，所以很容易求出存活率为(1-0? 1) X7A=0. 9x0. 989718-0. 8907 ; 1160 岁大象的繁殖率为0. 1448。然后根据上面的矩阵L建立起如果不进行避孕注射种群的增长的Leslie矩阵如下所(2)讨论A ”的特征根，分析种群增长规律示:用Matlab软件求得特征根为 R=1.0481，根据定理1知道,3、求出每年注射避孕药的母象头数根据Leslie矩阵的性质知道，要保持种群稳定，必须使得特征根*1,即使得下面 式子成立:% +b*0 + A2 八0八1 +?+ SO?*2? %_2= 1具体到本题来就是使得如下成立：(VI)bo ? s - sIG ( 1 +s+s)+.+s ) = 1解这个方程求出要保持大象种群的稳定，繁殖率应该为ao=o. 0377保持大象种群数M不变的繁殖率 5与没采取避孕时的繁殖率 b有一定的差距，所以 需要避孕 掉具有(b- b。)繁殖率母象所生的幼象。 假设每年要避孕头大象，由于一 次注射可以使得一头成熟的母象在两年不会受孕，所以每年实际上共有2xn。头大象处 于避孕期。这样根据需要避孕掉具有（ b-bQ 繁殖率母象所生的幼象的数目等于注射避孕药使得母象没有繁殖幼象的数目这个条件得到一个方程：（ 1/2） X85%XX1 x（ 0.1448-0.0377） =2nQ x 0.1448解之得n0=1393所以每年注射避孕药的母象头数为：1393 （头）4、分析不确定因素的影响（ 1） 最初一两年避孕母象发情期增多，与未避孕母象产生竞争求偶的公象，使部分能怀孕的母象不能怀孕。而避孕的母象每月发情一次，会扰乱了正常求偶的母象，这样会造成未避孕母象的繁殖率出现下降，避孕的母象数量应该减少。（ 2） 随着时间的增长，如果持续使用避孕药，会使象的年龄结构发生变化，象的结构呈老龄化，所以随着时间的增长，要保证象群的稳定，避孕药的使用量必定会逐年减少直至禁用。考虑转移大象时母象的避孕策略被转移的大象只考虑处于1 60 岁之间，转移后的大象看成每年多死了这么多头大象，即意后所需要进行避孕的母象头数；然后用理论去验证计算机模拟的正确性。1 、通过计算机模拟确定需要进行避孕的母象头数产生了 n （可自己指定）个01 的随机数。具体算法如下页图所示。图 2 （计算机模拟流程图）0J48a(i)0.346a(i)0.469a( i)0.617a( i)0,827a( i)Ixsliei = i+1k6=k6+k!=kI+k4=k4+下面议n=100为例进行计算机模拟。令n=100，进行10次计算机模拟，得到当运出大象头数为100头时，要使大象头数 稳定在11000头需避孕的大象头数如下：第一次模拟：1164 （头） 第二次模拟：1195 （头） 第三次模拟：1192 （头） 第四次模拟：1167 （头） 第五次模拟：1206 （头） 第六次模拟：1163 （头） 第七次模拟：1207 （头）第八次模拟：1190 （头）第九次模拟：1173 （头）第十次模拟：1178 （头）10次模拟得到需避孕的大象头数的平均数为1184头。因此可以认为当运出大象头数为100头时，要使大象头数稳定在 11000头需避孕的大象头数为1184头。同理，可以得到题目中要求的当运出大象50-300头要使大象头数稳定在11000头 需避孕的大象头数分别为：表2 （移出大象与对应得避孕母象头数）移 出 头 数5060708090100110120130140150160170避 孕 头 数1297127712531231121211841164115011121096103610301010移 出 头 数180190200210220230240250260270280290300遂 孕头 数970932917860838822764738714.654653591570图3 （运出大象数目与应该避孕母象数目的关系）从该图的变化趋势可以看出为了使大象总数稳定在11000头而应该避孕的母象数目与运出大象的头数近似成二次函数关系。通过最小二乘进行拟合二次函数得到移出大象头数X与避孕母象头数y关系：y = -0.00344X 2 - 1.7558x +1396.2画出图像为：图4 （原始数据关系与拟合关系对比图）原始关系与拟合得到关系对比图原始美系与拟吉得到关系对比图13001200110010009008D070060050050ICC1502325030Q移出大象头数2 .用理论去验证计算机模拟的正确性假设所转移的大象是有目的的挑选的, 即挑选大象时是按照第i （ i=l, 2, 3, 4, 5,6） 年龄组所占比例进行的，这也是符合情理。设每年被转移的大象共有M 头， 1 60 岁的大象共有 W 头。设第 i 年龄组占1 一一 60 岁的比例为叱，第i 年龄组的存活率九，则有几 =（ X, XSlXWi -M XMA） /（ Xj XMA） = （ X X - M） /X则要保持大象群落稳定，如第二问的做法有：Z? X X X （ 1 + S + . + 护 ） 1当移出大象头数帖50 时，解得b=0. 0455o又由第三问，把这个方程（1/2） X85%XXi x （0.1448-0.0377） =2n0 xO.1448 的亿换为 n、0. 0377换为0. 0455,求得n=129U所以，现在只需避孕1291 头，由计算机模拟得到的是1297 头，非常地接近。同样的道理可以验证当H二60, 70, 80的时候也是与计算机模拟很接近的。由此，可以说明我们用计算机模拟的方法是有效的。模型的评价和改进方法1、 模型的优点（ 1） 本文解决问题的模型都是比较简单的，但是这并不影响得到的结果的准确性，因为这些简单的模型都有很强的理论依据；（ 2） 在求解第二问的时候，充分利用Leslie 矩阵稳定性理论来求解应该让多少母象进行避孕注射，这些理论在差分方程中都是经典的理论，经得起许多事实的考验;（ 3） 第三问的求解中运用了计算机模拟方法来模拟移出大象属于哪个年龄段，这样不 仅求解方便、简洁（只需要把算法程序写好就可以得到结果），得到的结果与实际也更接近；（ 4） 第三问用计算机模拟得到数据后，又用理论去验证，这样使得结果更具有说服力；2、 模型需要改进的地方（ 1） 因为假设了大象性别是严格地1: 1 关系，而实际中不一定那么地严格是这样，所以如果能够把各个年龄段大象的性别比例分别计算，那么模型的结果可能更接近实际；（ 2） 在进行计算机模拟时，最开始的随机数的产生个数只有几十个，这几十个随机数不能很好的反映各个年龄段的大象所占的比重，这样势必会对结果造成一定的误差；参考文献：1、 启源数学模型（第三版）M. ： 高等教育20032、 静数学建模与数学实验（第2 版） M. ： 高等教育20033、 周晓阳数学实验与Mat lab M.: 华中科技大学20024、 谏 当代数学的若干理论与方法M. ： 华东理工大学20025、 尚志数学建模竞赛教程M. ： 教育 19966、 北峰数模网 mem, zjnu. net, cn/2006-8-13附录：1、计算机模拟Mat lab程序代码function bys,b0=moni(n)b=121610121714;c=b/sum(b);b=cumsum(c)* ;a=rand(n. 1);kl 二 0; k2=0; k3=0; k4=0; k5=0; k6=0;i=l;for i=l : n if (a(i)b(l)&a(i)b(2)&a(i)b(3)&a(i)b(4)&a(i)b(5)&a(i) 1-0置=数 -1_| |_|数性 数雌 头叫V 土 头数11OO31301A0031372OO3 25222011 03211 6123OO3 38532112333643234123413334C IY654?13 51A 045125350 IXC06723 63116134361287201A 03 7737224371658923 8M28473812291573,911 009,0 0213911 04C IX9310 1A 612?一 ?08 aI 0126O IX1122121121?126C 1Y1A 111 911 0123?1423412130 IX4237323313 1A 0613M0 1443240145214124142712447IA 01513615641531A4511 641 1 621 IX 016321643，4 6251217O04I 76011 712817123182212489318200 zI 845 :3 31 1 91461932I 925119 92312205450 1A 04，201 74503402137513IX211411 0513922 11 O55i 2642211 075211 6423O05i 32123005311 04243554562! 42054177253012 5i 54IX25305533261 16561312612561362712557105274257123280058321228315832292010591482921592211306560003030602010