野兔生长问题摘要本文根据已知的野兔连续十年的统计情况，探讨野兔的合理的存活率并推测当前的发展趋势，针对不同情况给出方法推算出野兔数量的走向的目的。
首先，充分利用给出的前两年来野兔的数量变化，分析近两年来的野兔群落的情况，建立一个线性方程组的数学模型，通过求解方程组得出不同年份野兔的数量的数学关系，并且求出了平均增长率为：1.718%；所以通过一些比例之间的关系得到这个野兔群落的T＝10的数量（见表1）。
然后，建立一个种群增长的差分方程模型，求出的野兔生长规律。
求解当前野兔对应的Leslie矩阵的特征根，发现该特征根大于1，根据Leslie矩阵的稳定性理论知道：如果不进行避孕注射该野兔种群将无限增长（如果环境允许）；据此，利用Leslie矩阵稳定的充要条件求出应该保持多大的繁殖率才能使种群保持稳定，求解的主要思路是：特征根取为1、把繁殖率当成未知数，将此时的各年龄段的存活率代入方程⑥即可。
最后，只需将野兔的存活率代入那个以繁殖率为未知数的方程（方程⑥），求出在哪些年内野兔的增长有异常现象，。
考虑到求解的数据比较多，采取计算机模拟的方法来确定移走野兔后所需要进行避孕的母兔头数为了检验计算机模拟的正确性，用理论去验证。
问题重述位于某国的国家公园中栖息着近10000头野兔。
管理者要求有一个健康自由的环境以便观察这个10000头野兔的数量变化情况。
管理者逐年统计了野兔的数量，发现在过去的10年中，野兔的生长变化并不稳定，呈现波浪式起伏，根据这些信息我们需要解决以下问题：1. 探讨年龄在1岁到10岁之间的野兔的合理的存活率的模型，推测这个野兔群落的当前的年龄结构。
2. 知道哪些环境和内部因素对野兔生长数量的影响，并测算出各个影响的程度如何。
3. 探求偶然突发事件对野兔生长数量的巨大影响和它的规律性。
4. 根据野兔的生长变化，对野兔的生长特点进行分析。
问题假设1、假设野兔的性别比近似认为1：1，并且采用措施维持这个性别比；2、假设母兔可以怀孕的年龄为1岁—6岁、最高年龄为10岁，10岁的死亡率为100%，并且6—10岁的野兔的只数呈线性递减；3、假设野兔在各年龄段中的分布率不变，即年龄结构不变，并采用各种措施维持这一结构；4、假设兔子的内部因素对其生存率的影响不大5、假设0岁野兔能够活到1岁的比例为75%；6、假设各个环境因素对野兔生长的影响是互不影响的。
符号说明X ： 表示一年中野兔的头数（i=0表示0岁野兔的头数，i=1表示1--10岁大象头数，i=2表示1—10岁野兔的头数）；p ： 表示存活率（0p 表示0岁野兔的存活率，1p 表示1—60岁野兔的存活率，p表示61岁—70岁野兔的存活率）；2x k：表示时段k第i年龄组的野兔数量；()b：第i年龄组每个（母兔）个体在1个时段内平均繁殖的数量；s：第i年龄组的存活率；L： Leslie矩阵；：L矩阵的那个唯一正特征根；1n：表示移出野兔的头数；问题分析对于问题一，利用给出的近两年来运出的野兔的数量统计表，可以分析近十年来的野兔群落的变化情况，比如各个年龄段的野兔占总的野兔的头数的比例是多少，还可以根据十年野兔总数的不规则变化来建立方程，用于求解野兔生存规律。
对于问题二，因为考虑的是公园在未来很长一段时间的野兔生长变化问题，所以可以建立一个按年龄分组的种群增长的差分方程模型，根据差分方程的Leslie矩阵的特征根，结合 Leslie矩阵的稳定性理论对当前野兔种群的情况进行分析。
为了保持野兔种群的稳定，必须使得Leslie矩阵的最大特征根为1，而这样，特征根取为1、把繁殖率当成未知数，将此时的各年龄段的存活率代入方程特征方程，求解这个以繁殖率为未知数的方程可以得到野兔在某一年的繁殖率的取值；根据一条件建立方程来求解野兔在哪些年数量变化有异常情况。
对于问题三，可以建立一个按年龄分组的种群增长的差分方程模型结合Leslie矩阵的稳定性理论对当前野兔种群的情况进行分析，然后求解这个以繁殖率为未知数的方程可以得到野兔在T=10年的繁殖率的取值，接着把取值代入线性方程求出野兔在这一年的数量。
探讨野兔的存活率和当前野兔的年龄结构下面将根据给出的近两年来运出的野兔的数量与性别统计表，分析近两年来的野兔群落的情况，建立一个线性方程组数学模型，通过求解方程组得到年龄在1岁到10岁之间的野兔的存活率，并给出野兔各年龄所占的比例，进而得到这个野兔群落的当前的年龄结构。
1、线性方程组模型的建立（1）首先，计算一年中野兔的只数。
野兔群是由0岁，1—0岁，10岁—12岁组成 ,且稳定在10000头。
设0岁的头数为X 0，1—10岁野兔头数为X 1，10岁—12岁野兔头数为X 2。
所以得到第一个方程：X 0+X 1+X 2=10000 ①（2）其次，考虑到前一年野兔的总数等于前两年存活下来的野兔加上新生的幼儿再减去运出的野兔数。
设0岁野兔的存活率为0p ，1—10岁野兔的存活率为1p ，10岁—12岁野兔的存活率为2p 。
则经过一年后，新生的野兔存活下来的头数为X 0⨯0p ；1到10岁的野兔存活下来的头数为X 11p ⨯；10岁——12岁的野兔能存活下来的头数为X 22p ⨯，因此得到第二个方程：（X 0⨯0p + X 11p ⨯+ X 22p ⨯）+ X 0-622=10000 ②联立①②得到方程组：0120011220X +X +X =11000X + X + X + X -622=11000 p p p ⎧⎨⨯⨯⨯⎩⑴ 2、模型的求解根据近两年来运出的野兔的数量与性别统计表，得到如下分析结果： （1）计算0岁的野兔头数由表中统计，1岁—3岁的野兔占1岁—10岁的野兔比例为：（67/620+169/876）/2=15.05%所以得到：3岁—10岁能生小野兔的母野兔占1岁—10岁的野兔比例为：（1-15.032%）⨯0.5 =42.48%因为能生小野兔的母野兔每3个月生一头小野兔，且十胞胎的机会为1.35%，相当于每年生126头 ，所以0岁的野兔占1岁—10岁的野兔比例为：0.4248⨯0.2896=0.12303这样0岁的野兔共有： 0X=0.12303⨯1 X ③（2）计算10岁—12岁的野兔头数从表中计算运出的9岁的野兔占运出的总野兔比率为：（14/622+22/876）/2=0.0238由于运出的野兔都是1岁—10岁的，所以0.0238也可看为9岁的野兔占1—10岁的野兔的头数比例，得到10岁的野兔占的比例为0.02381p ⨯，由假设可以知道：10岁—12岁的野兔头数为：2X=1/2⨯10⨯0.02381p ⨯⨯X 1 ④10岁——12岁的野兔经过一年能存活下来的头数为：2 211X =(1/2)90.0238X p p ⨯⨯⨯⨯⨯ ⑤ （3）、将③⑤和④两个式子代入上面方程组⑴得：111110111110.12303X +X +(1/2 )100.0238 X =110000.12303X + X + (1/2)90.0238X +0.12303X -622=11000 p p p p ⨯⨯⨯⨯⨯⎧⎨⨯⨯⨯⨯⨯⨯⨯⨯⎩又由假设知道，0岁野兔的存活率为0p =75%代入上述方程组，然后用Mathematica 解之得：110.9897198864.85p X =⎧⎨=⎩ 再依次将1X 、s 代入③⑤和④求得：201044.07X 1090.66X =⎧⎨=⎩所以， 0岁野兔的总头数为1091（头）；1—10岁的野兔的存活率为98.9719%，总头数为8865（头）；10岁—12岁的野兔头数为1091（头）。
把0—10岁的野兔分为八个年龄段，由假设知道，各个年龄段占总数可以用各个年龄段移出的头数除以移出的总头数来衡量。
下面以1—10年龄段的野兔头数计算为例：前一年总共移出622头，其中1—10岁移出为67头；前两年总共移出876头，其中1—10岁移出169头。
故1—10年龄段的野兔头数可以这样计算：11X =671698865 [()/2]622 876⨯+=1332（头）3、结果分析（1）由结果可以知道，2—10岁野兔的存活率为98.9718%，这与题目给出的大于95%是相一致的，所以可以认为结果是合理的；（2）从图1可以看出，各个年龄段的野兔所占的比例基本上是一样的，2—3岁和4—7岁的野兔比例相对比较大，因为这段大象正处于年龄的黄金时期。
由此，可以认为求出的野兔年龄也是合理的。
估计注射避孕药后野兔数量的变化首先建立一个按年龄分组的种群增长的差分方程模型；然后用Leslie 矩阵稳定的充要条件分析如果不进行避孕注射种群的增长情况；最后仍然利用Leslie 矩阵稳定的充要条件求出应该保持多大的繁殖率才能使种群保持稳定，进而利用一个方程求出每年注射避孕药的母野兔头数。
1、按年龄分组的种群增长的差分方程模型的建立记时段k 第i 年龄组的野兔数量为()i x k ，k=0，1，2，…，i=1，2，…n,第i 年龄组的繁殖率为b ，即第i 年龄组每个（母野兔）个体在1个时段内平均繁殖的数量，第i 年龄的存活率为s ，我们这里假设b 和s 不随时段k 变化，在稳定的环境下这个假设是合理的。
b 和s 可由统计资料获得。
()i x k 的变化规律由以下的基本事实得到：时段k+1第1年龄组种群数量是时段k 各年龄组繁殖数量之和，即11(1)()ni i i x k bx k =+=∑ 时段k+1第i+1年龄组的种群数量是时段k 第i 年龄组存活下来的数量，即 1(1)(),1,2,,1i ii x k s x k i n ++==⋯- 记时段k 种群按年龄组的分布向量为12()[(),(),()]Tnx k x k x k x k =⋯ 由繁殖率b 和存活率s 构成的矩阵为010116061697001910110...0.........0...000......00...000...000......00...0000...000......00...0000...00......00...0000...00......00...00.........00......00...00.........000......00...000b b b b b b b s s s s L s =606159...000......00...0000...000......0...0000...000......0...0000...000......00 (000)............0s b b ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦根据Leslie 矩阵的性质可以得到如下的定理：定理1：L 矩阵有唯一的正特征根1λ，且它是单重根的，1λ对应正特征向量*12111221111...1,,,...,Tn n s s s s s s x λλλ--⎡⎤=⎢⎥⎣⎦L 矩阵的其他n-1个特征根k λ都是满足 1,2,3,...,k k n λλ≤= 该定理表明L 矩阵的特征方程121211122121(.........)0n n n n n n n b b s b s s sb s s s λλλλ-----+++++= 只有一个正根，并且易知，**1Lxx λ=2、如果不进行避孕注射种群的增长情况（1）建立Leslie 矩阵首先，由第一问的求解知道，0岁的野兔的存活率为0.75；1—10岁野兔的存活率为0.989718；根据假设10—12岁野兔头数是线性递减的，而且到70岁所有的野兔都死完了，所以很容易求出存活率为（1-0.1）1s ⨯=0.9⨯0.989718=0.8907；9—10岁野兔的繁殖率为0.1448。