[论文解读] Applications of Perron-Frobenius Theory to Population Dynamics

Chi-Kwong Li, Hans Schneider|ArXiv.org|Sep 2, 2001

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 12被引用 4

一句话总结

本文将Perron-Frobenius理论应用于种群动态中的矩阵模型，对人口统计学基本定理提供了简洁的证明，并对Cushing与Yicang关于净繁殖率$ R_0 $的定理进行了改进。研究证明，通过$ q(s) $对生育力矩阵进行缩放，可实现任意期望的增长率$ s $，其中$ R_0 = 1 $对应于$ s = 1 $时的人口稳定状态，并展示了$ R_0 $如何通过下一代矩阵的谱性质决定种群的长期行为。

ABSTRACT

By the use of Perron-Frobenius theory, simple proofs are given of the Fundamental Theorem of Demography and of a theorem of Cushing and Yicang on the net reproductive rate occurring in matrix models of population dynamics. The latter result is further refined with some additional nonnegative matrix theory. When the fertility matrix is scaled by the net reproductive rate, the growth rate of the model is 1. More generally, we show how to achieve a given growth rate for the model by scaling the fertility matrix. Demographic interpretations of the results are given.

研究动机与目标

• 使用Perron-Frobenius理论，为种群动态中的关键结果提供简洁而严谨的证明。
• 利用非负矩阵理论，改进并推广Cushing与Yicang关于净繁殖率$ R_0 $的定理。
• 建立一个函数关系$ q(s) $，通过该关系对生育力矩阵进行缩放，以实现任意指定的增长率$ s $。
• 阐明$ R_0 $、$ q(s) $以及投影矩阵$ P $的谱性质在人口统计学中的解释。
• 证明当$ R_0 = 1 $时，对应于零增长率和人口静止状态，前提是满足适当条件。

提出的方法

• 利用不可约非负矩阵的Perron-Frobenius定理，分析投影矩阵$ P = T + F $的谱性质。
• 定义下一代矩阵$ Q = F(I - T)^{-1} $，其谱半径$ \rho(Q) = R_0 $即为净繁殖率。
• 推导出$ q(s) $作为函数，使得对$ F $按$ q(s) $缩放后，所得模型的增长率为$ s $，其中$ q(s) = \text{矩阵} F(I - T/s)^{-1}/s \text{的首项} $。
• 应用Perron-Frobenius理论的图论推广，以改进不可约矩阵的结果。
• 在$ \rho(T) < 1 $的条件下，利用预解恒等式$ (I - T)^{-1} = \sum_{k=0}^\infty T^k $，将$ Q $表示为一系列生存与繁殖事件的无穷级数。
• 分析特征方程$ q(r) = 1 $，其中$ r $为增长率，表明$ r $是$ s^{-1} $的多项式方程的正根。

实验结果

研究问题

• RQ1 如何利用Perron-Frobenius理论，为本征投影矩阵的种群统计学基本定理提供证明？
• RQ2 净繁殖率$ R_0 $在决定长期种群增长或稳定状态中起何种精确作用？
• RQ3 如何对生育力矩阵进行缩放以实现期望的人口增长率$ s $，以及所需缩放因子$ q(s) $的函数形式是什么？
• RQ4 当$ R_0 = 1 $时，其人口统计学解释是什么，它与种群静止状态有何关联？
• RQ5 当$ P $为不可约但非本原矩阵时，$ Q $与$ P $的谱性质之间有何关系，这对世代稳定性有何影响？

主要发现

• 净繁殖率$ R_0 = \rho(Q) $是下一代矩阵$ Q = F(I - T)^{-1} $的谱半径，表示每个新生个体在其一生中预期产生的后代数量。
• 当$ R_0 = 1 $时，种群增长率$ r = 1 $，在本征投影矩阵下表示种群处于静止状态。
• 对于任意$ s > 0 $，可通过$ q(s) $对生育力矩阵$ F $进行缩放，使得所得模型的增长率为$ s $，其中$ q(s) $为$ F(I - T/s)^{-1}/s $的首项。
• 在Leslie矩阵情形下，$ q(s) = f_1 s^{-1} + f_2 t_1 s^{-2} + \cdots + f_n (t_{n-1} \cdots t_1) s^{-n} $，即$ s^{-1} $的多项式，且增长率$ r $满足$ q(r) = 1 $。
• 在$ T $与$ F $被缩放使得$ \rho(T) = 0 $的示例中，缩放函数为$ q(s) = (1 + 2s^2)/(8s^4) $，对应的$ R_0(s) = 3s^4 / (1 + 2s^2) $。
• 即使$ \rho(P) = 1 $，归一化种群仍可能无限振荡，表明世代稳定性并不意味着种群分布的收敛。