大象群落的稳定发展

大象群落的稳定发展摘要本文根据非洲某国的国家公园近两年内从公园运出的大象的大致年龄和性别的统计情况，探讨大象的合理的存活率并推测当前的年龄结构，针对不同情况给出如何进行避孕注射以达到控制大象数量的目的。首先，充分利用给出的近两年来运出的大象的数量与性别统计表，分析近两年来的大象群落的情况，建立一个线性方程组的数学模型，通过求解方程组得到年龄在2岁到60岁之间的大象的总数，并且求出了存活率为：98.9718%；因为假设公园内2岁到60岁之间的大象占总大象的比例等于运出的2岁到60岁之间的大象占总移出大象的比例，所以通过一些比例之间的关系得到这个大象群落的当前的年龄结构（见表1）。最后，假设被转移的大象只考虑处于1—60岁之间，这样可以认为转移后的大象看成每年多死了这么多头大象，即意味着死亡率将增加，存活率将减少；仍然按照解决第二问的模型，只需将此时不同的各年龄段大象的存活率代入那个以繁殖率为未知数的方程（方程VI），求出应该保持多大的繁殖率才能使种群保持稳定。考虑到求解的数据比较多，采取计算机模拟的方法来确定移出大象后所需要进行避孕的母象头数（见表2），为了检验计算机模拟的正确性，用理论去验证。模拟的思路方法见计算机模拟流程图—图2。关键词：关键字：线性方程组、差分方程模型、Lelie矩阵、计算机模拟问题重述位于非洲某国的国家公园中栖息着近11000头大象。管理者要求有一个健康稳定的环境以便维持这个11000头大象的稳定群落。管理者逐年统计了大象的数量，发现在过去的20年中，整个大象群经过一些偷猎枪杀以及转移到外地还能保持在11000头的数量，而其中每年大约有近600头到800头是被转移的。由于近年来，偷猎被禁止，而且每年要转移这些大象也比较困难，现决定采取避孕注射法以维持大象数量的平衡。我们已知此公园近两年内从这个地区运出的大象的大致年龄和性别的统计。根据这些信息我们需要解决以下问题：探讨年龄在2岁到60岁之间的象的合理的存活率的模型，推测这个大象群落的当前的年龄结构。估计每年有多少母象要注射避孕药，可以使象群固定在11000头左右。这里不免有些不确定性，是否能估计这种不确定性的影响。假如每年转移50至300头象到别处，那么上面的避孕措施将可以有怎样的改变？问题假设1、假设大象的性别比近似认为1：1，并且采用措施维持这个性别比2、假设母象可以怀孕的年龄为11岁—60岁、最高年龄为70岁，70岁的死亡率为100%，并且61—70岁的大象的头数呈线性递减；3、假设大象在各年龄段中的分布率不变，即年龄结构不变，并采用各种措施维持这一结构；4、假设被转移的大象只考虑处于1—60岁之间，转移后的大象看成每年多死了这么多头大象；5、假设0岁大象能够活到1岁的比例为75%；符号说明某i：表示一年中大象的头数（i=0表示0岁大象的头数，i=l表示1--60岁大象头数，i=2表示61—70岁大象的头数）；pi：表示存活率（p0表示0岁大象的存活率，pl表示1—60岁大象的存活率，p2表示61岁—70岁大象的存活率）；某i（k）:表示时段k第i年龄组的大象数量；bi：第i年龄组每个（母象）个体在1个时段内平均繁殖的数量；i：第i年龄组的存活率；L：Lelie矩阵；21：L矩阵的那个唯一正特征根；n表示移出大象的头数；问题分析对于问题一，利用给出的近两年来运出的大象的数量与性别统计表，可以分析近两年来的大象群落的情况，比如移出的各个年龄段的大象占移出的总的大象的头数的比例是多少，还可以根据两年移出大象后大象总数都是11000来建立方程，用于求解存活率。对于问题二，因为考虑的是公园在未来很长一段时间的大象种群控制问题，所以可以建立一个按年龄分组的种群增长的差分方程模型，根据差分方程的Lelie矩阵的特征根，结合Lelie矩阵的稳定性理论对当前大象种群的情况进行分析。为了保持大象种群的稳定，必须使得Lelie矩阵的最大特征根为1，而这样，特征根取为1、把繁殖率当成未知数，将此时的各年龄段的存活率代入方程特征方程，求解这个以繁殖率为未知数的方程可以得到要使种群保持稳定繁殖率的取值；根据需要避孕掉母象所生的幼象的数目等于注射避孕药使得母象没有繁殖幼象的数目这一条件建立方程来求解应该对多少头母象进行避孕。对于问题三，由于假设被转移的大象只考虑处于1—60岁之间，故可以认为转移后的大象看成每年多死了这么多头大象，即意味着死亡率将增加，存活率将减少。按照解决第二问的模型，只需将此时不同的各年龄段大象的存活率代替原来的存活率，就可以求出此时应该保持多大的繁殖率才能使种群保持稳定。为了方便，可用采用计算机模拟的方法来确定移出的大象在哪个年龄段，考虑到计算机模拟的不确定性，可以对模拟结果进行检验。探讨大象的存活率和当前大象的年龄结构下面将根据给出的近两年来运出的大象的数量与性别统计表，分析近两年来的大象群落的情况，建立一个线性方程组数学模型，通过求解方程组得到年龄在2岁到60岁之间的大象的存活率，并给出大象各年龄所占的比例，进而得到这个大象群落的当前的年龄结构。1、线性方程组模型的建立（1）首先，计算一年中大象的头数。大象群是由0岁，1—60岁，61岁—70岁组成,且稳定在11000头。设0岁的头数为某0，1—60岁大象头数为某1，61岁—70岁大象头数为某2。所以得到第一个方程：某0+某1+某2=11000（I）（2）其次，考虑到前一年大象的总数等于前两年存活下来的大象加上新生的幼儿再减去运出的大象数。设0岁大象的存活率为p0,l—60岁大象的存活率为pl,61岁一70岁大象的存活率为p2。则经过一年后，新生的大象存活下来的头数为某0P0;1到60岁的大象存活3下来的头数为某1p1；61岁一一70岁的大象能存活下来的头数为某2p2，因此得到第二个方程：（某0P0+某1p1+某2p2）+某0-622=11000（II）联立（I）、（II）得到方程组：某0+某1+某2=11000（某）某p+某p+某p+某-622=1100011220002、模型的求解根据近两年来运出的大象的数量与性别统计表，得到如下分析结果：(1)计算0岁的大象头数由表中统计，1岁—10岁的大象占1岁—60岁的大象比例为：(67/622+169/876)/2=15.032%所以得到：11岁—60岁能生小象的母象占1岁—60岁的大象比例为(1-15.032%)0.5=42.48%因为能生小象的母象每3.5年生一头小象，且双胞胎的机会为1.35%相当于每年生0.2896头，所以0岁的大象占1岁—60岁的大象比例为：0.42480.2896=0.12303这样0岁的大象共有：某0=0.12303某1(111)(2)计算61岁一70岁的大象头数从表中计算运出的59岁的大象占运出的总大象比率为：(14/622+22/876)/2=0.0238由于运出的大象都是1岁—60岁的，所以0.0238也可看为59岁的大象占1—60岁的大象的头数比例，得到60岁的大象占的比例为0.0238p1，由假设可以知道：61岁—70岁的大象头数为：某2=1/2100.0238p1某1(IV)61岁——70岁的大象经过一年能存活下来的头数为：某2p2=(1/2)90.0238p1某1(V)（3）、将（III）、（V）和（^）两个式子代入上面方程组（某）得：0.12303某1+某l+（l/2）100.0238pl某1=110000.12303某1p0+某lpl+（l/2）90.0238pl某1+0.12303某1-622=11000又由假设知道，0岁大象的存活率为p0=75%代入上述方程组，然后用Mathematica解之得：p10.989719某18864.854再依次将某1、1代入（III）、（V）和（旦）求得：某21044.07某1090.660所以，0岁大象的总头数为1091（头）；1—60岁的大象的存活率为98.9719%，总头数为8865（头）；61岁—70岁的大象头数为1091（头）。把0—70岁的大象分为八个年龄段，由假设知道，各个年龄段占总数可以用各个年龄段移出的头数除以移出的总头数来衡量。下面以1—10年龄段的大象头数计算为例：前一年总共移出622头，其中1—10岁移出为67头；前两年总共移出876头，其中1—10岁移出169头。故1—10年龄段的大象头数可以这样计算：67169）/2]=1333（头）某11=8865[（622876其他的年龄段用同样的方法计算，得到如下表（附饼形图）：表1（大象年龄结构）年龄头数010911-10133311-20177721-30106931-40125541-5018881751-6015441461-7010449比例1012161011（%）图1（大象年龄结构饼图）3、结果分析(1)由结果可以知道，2—60岁大象的存活率为98.9718%，这与题目给出的大于95%是相一致的，所以可以认为结果是合理的；(2)从图1可以看出，各个年龄段的大象所占的比例基本上是一样的，21—30岁和41—50岁的大象比例相对比较大，因为这段大象正处于年龄的黄金时期。由此，可以认为求出的大象年龄也是合理的。估计每年注射避孕药的母象头数为了估计每年注射避孕药的母象头数，首先建立一个按年龄分组的种群增长的差分5方程模型；然后用Lelie矩阵稳定的充要条件分析如果不进行避孕注射种群的增长情况；最后仍然利用Lelie矩阵稳定的充要条件求出应该保持多大的繁殖率才能使种群保持稳定，进而利用一个方程求出每年注射避孕药的母象头数。1、按年龄分组的种群增长的差分方程模型的建立,l,2,nl某(k)[某1(k),某2(k),某n(k)]T由繁殖率bi和存活率i构成的矩阵为b0b10001000000L 0000000000...b9...0...0...0...9...0...0...0...0...0...0...0...0b100000100000000b110000011000000.. b60...0...0...0...0...0...0...0...0...60...0...0...0b610000000006100...b69...0...0...0...0...0...0...0...0...0...0...0...69b70000000000000根据Lelie矩阵的性质可以得到如下的定理：定理1：L矩阵有唯一的正特征根1,且它是单重根的，1对应正特征向量6...某某1,1,122,...,12n1n1111L矩阵的其他n-1个特征根k都是满足k1,k2,3,...,nT该定理表明L矩阵的特征方程n（b1n1b21n2...bn112...n2bn12...n1）0只有一个正根，并且易知，L某某1某某2、如果不进行避孕注射种群的增长情况（1）建立Lelie矩阵首先，由第一问的求解知道，0岁的大象的存活率为0.75；1—60岁大象的存活率为0.989718；根据假设61—70岁大象头数是线性递减的，而且到70岁所有的大象都死完了，所以很容易求出存活率为（1-0.1）p1=0.90.989718=0.8907；11—60岁大象的繁殖率为0.1448。然后根据上面的矩阵L建立起如果不进行避孕注射种群的增长的Lelie矩阵如下所示：0...00.14480.1448 0.14480...0000.750...000 00...0000.9897...000 00...0000...000 00...0000...0.989700 ..00...0000...00.98970 00...00 000.9897 00...00这是一L7171= 000 00...00...00...000 00...0000...000. 0.98970...0000...000 00.8907...0000...000 00...0000...000 00...0.89070个7171的矩阵。(2)讨论L7171的特征根，分析种群增长规律用Matlab软件求得特征根为R=1.0481，根据定理1知道，如果不进行避孕注射，该大象种群将无限增长下去(如果环境允许)，所以要进行避孕注射。3、求出每年注射避孕药的母象头数根据Lelie矩阵的性质知道，要保持种群稳定，必须使得特征根r=1即使得下面式子成立：b0b10b201...0bn112...n2bn012...n11(VI)具体到本题来就是使得如下成立：2b010(1++7+4.9..+)=1解这个方程求出要保持大象种群的稳定，繁殖率应该为b0=0.0377保持大象种群数量不变的繁殖率b0与没采取避孕时的繁殖率b有一定的差距，所以需要避孕掉具有(b-b0)繁殖率母象所生的幼象。假设每年要避孕n0头大象，由于一次注射可以使得一头成熟的母象在两年内不会受孕，所以每年实际上共有2n0头大象处于避孕期。这样根据需要避孕掉具有(b-b0)繁殖率母象所生的幼象的数目等于注射避孕药使得母象没有繁殖幼象的数目这个条件得到一个方程：(1/2)85%某(0.1448-0.0377)=2n0.14484、分析不确定因素的影响（1）最初一两年避孕母象发情期增多，与未避孕母象产生竞争求偶的公象，使部分能怀孕的母象不能怀孕。而避孕的母象每月发情一次，会扰乱了正常求偶的母象，这样会造成未避孕母象的繁殖率出现下降，避孕的母象数量应该减少。（2）随着时间的增长，如果持续使用避孕药，会使象的年龄结构发生变化，象的结构呈老龄化，所以随着时间的增长，要保证象群的稳定，避孕药的使用量必定会逐年减少直至禁用。考虑转移大象时母象的避孕策略被转移的大象只考虑处于1—60岁之间，转移后的大象看成每年多死了这么多头大象，即意味着死亡率将增加，存活率将减少。下面首先通过计算机模拟来确定移出大象后所需要进行避孕的母象头数；然后用理论去验证计算机模拟的正确性。1、通过计算机模拟确定需要进行避孕的母象头数产生了n（可自己指定）个0~1的随机数。具体算法如下页图所示。图2（计算机模拟流程图）初始化：按照各个年龄段大象所占比重，给每个年龄段分配一个区间，每个年龄运产生n个0~1随机数记为向量a（i）No0.148令n=100,进行10次计算机模拟，得到当运出大象头数为100头时，要使大象头数稳定在11000头需避孕的大象头数如下：第一次模拟：1164（头）第二次模拟：1195（头）第三次模拟：1192（头）第四次模拟：1167（头）第五次模拟：1206（头）第六次模拟：1163（头）第七次模拟：1207（头）第八次模拟：1190（头）第九次模拟：1173（头）第十次模拟：1178（头）次模拟得到需避孕的大象头数的平均数为1184头。因此可以认为当运出大象头数为100头时，要使大象头数稳定在11000头需避孕的大象头数为1184头。同理，可以得到题目中要求的当运出大象50—300头要使大象头数稳定在11000头需避孕的大象头数分别为：表2（移出大象与对应得避孕母象头数）移出头数避孕头数移出头数避孕头数18019020021022023024025026027028029030012971277125312311212118411641150111210961036103010105060708090100110120130140150160170970932917860838822764738714.654653591570图3（运出大象数目与应该避孕母象数目的关系）10图4（原始数据关系与拟合关系对比图）2、用理论去验证计算机模拟的正确性11假设所转移的大象是有目的的挑选的，即挑选大象时是按照第i（i=l,2,3,4,5,6）年龄组所占比例进行的，这也是符合情理。设每年被转移的大象共有M头，1――60岁的大