线性代数方法建模2种群年龄结构的估算--数学建模案例分析_百重点

§2 种群年龄结构的估算

英国的生物学家莱斯利（P.H.Leslie）在1945年给了种群的年龄结构的离散数学模型，它是利用矩阵的特征值和特征向量来预测种群年龄结构及种群的数量变化趋势`，这种方法在种群动态分析及人口统计学中得到了广泛的应用，并得以发展。

在某动物群体中，仅考虑雌性动物的年龄和数量。设雌性动物的最大生存年龄为（年或其它单位）。把等分成个年龄组，每一个年龄组的长度为：

设第个年龄组的生育率为，存活率为。应注意表示第个年龄组的每一雌性动物平均生育的雌性幼体个数；表示第个年龄组中可存活到第

+1年龄组的雌性数与该年龄组总数之比。在不发生意外事件（灾害等）的条件下，，均为常数，且，。同时，假设至少有一个年龄组的雌性动物具有生育能力，即至少有一个。

利用统计资料可获得基年该种群在各年龄组的雌性动物数量。记为

时第个年龄组雌性动物的数量，就得到初始时刻年龄分布向量

如果以年龄组的间隔作为时间单位，记

并统计时各年龄组雌性动物的数量，可得到时的年龄分布向量

随着时间的变化，由于出生、死亡以及年龄的增长，该种群中每一年龄组的雌性动

物数量都将发生变化。实际上，在时，种群中第一年龄组的雌性个数应等于在

和之间出生的所有雌性幼体的总和，即

（1）

同时，在时，第+1年龄组中雌性动物的数量应等于在时第个年龄组中雌性动物数量乘以存活率，即

（2）

综合上述分析，由（1）和（2）可得到与时各年龄组中雌性动物数量间的关系：

（3）

记矩阵

则（3）可写成

（4）

其中矩阵称为莱斯利矩阵。显然

（5）

要研究种群年龄结构的变化趋势，需讨论随的演变情况。设矩阵的特征值为

不妨设将它们按模由大到小排列，依次为

称为的主特征值。若，称为的严格主特征值。从前面的分析可见，如果有个线性无关的特征向量，且有严格主特征值，则

记为常向量，则有

（6）

当充分大时，，这说明各年龄组的个体个数可以近似地看成以比率

按几何级数“增长”。时增多；时稳定；时减少。若记

则随着时间的延续，各年龄组的个体数目在总数中所占比率由的各分量确定，与初始年龄分布无关。由此可见，不管总数是增加还是减少，不管初始时刻年龄分布如何，各年龄组的个体总数占总数的比例趋向稳定。这是种群的一种内禀性质。用莱斯利矩阵模型描述的种群年龄结构，其年龄分布是否稳定，稳定值是多少，达到稳定分布的时间以及种群总数的变化与下列问题有关：1、矩阵是否有一个严格

主特征值；2、主特征值的大小及对应的特征向量；3、按模第二大的特征值

与之比的模。即需要对矩阵的特征值与特征向量进行分析。矩阵

的每一个元素都是非负的，称这种实矩阵为非负矩阵。对莱斯利矩阵有许多理论上的研究成果，下面仅列出与解决上述问题有关的两个结论：

命题1矩阵有唯一一个正的单特征值，且为主特征值。

命题2 如果是矩阵的一个非零特征值，则

为与对应的一个特征向量（注意：的各分量为正数）。

例：已知一种昆虫每2周产卵一次，6周以后即死亡，孵化后的幼虫2周后成熟，平均产卵100个，4周龄的成虫平均产卵150个。假设每个卵发育为2周龄成虫的概率为0.09（成活率），2 周龄的成虫发育成4周龄成虫的概率为0.2。又假设在开始时，0~2，2~4，4~6周龄的昆虫数目相同。则有

（这里不妨设）

这样，只要用矩阵乘法就能得到各周龄的昆虫数目分别为：

，，

计算出的特征值为，，

一组对应的特征向量为，

，

由此可算出，第58周以后，各周龄组的昆虫比例趋于稳定值。