第 70A 章 种群生态学¶
前置：非负矩阵与 Perron-Frobenius (Ch17) · 矩阵分析 (Ch14) · 差分方程
本章脉络：生物系统的线性表示 \(\to\) 莱斯利 (Leslie) 矩阵的构造：存活率与生育率 \(\to\) 种群演化的离散模型 \(\mathbf{x}_{t+1} = L \mathbf{x}_t\) \(\to\) 核心判据：Perron 特征值（主增长率） \(\to\) 稳定年龄分布 (Stable Age Distribution) \(\to\) 种群动态的谱分析 \(\to\) 衰减系数与灭绝判定 \(\to\) 应用：濒危物种保护规划、渔业捕捞限额设定、入侵物种预测
延伸：种群生态学是线性代数的“生命脉动”；它将繁衍生息的生物过程抽象为矩阵的迭代幂运算，证明了种群的最终结构是由其生理参数（特征值与特征向量）预先决定的，是数理生物学的代数底座
在自然界中，不同年龄段的生物对种群增长的贡献是不同的。幼体不繁殖但存活率高，成体繁殖力强。种群生态学利用 Leslie 矩阵 将这些生理参数整合进一个统一的线性模型。通过研究矩阵的谱性质，我们不仅可以预测种群未来的总数，还能知道各年龄段比例在长期演化后会如何趋于一个“完美的平衡态”。本章将介绍这一作为生态资源管理核心的代数工具。
70A.1 Leslie 矩阵的构造¶
定义 70A.1 (Leslie 矩阵 \(L\))
对于 \(n\) 个年龄组的种群，\(L\) 具有如下特定结构： $\(L = \begin{pmatrix} f_1 & f_2 & \cdots & f_{n-1} & f_n \\ s_1 & 0 & \cdots & 0 & 0 \\ 0 & s_2 & \cdots & 0 & 0 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \cdots & s_{n-1} & 0 \end{pmatrix}\)$ - \(f_i \ge 0\)：第 \(i\) 年龄组的单位生育率（第一行）。 - \(s_i \in (0, 1]\)：从第 \(i\) 组存活到第 \(i+1\) 组的存活率（次对角线）。
70A.2 主特征值与稳定分布¶
定理 70A.1 (Perron-Frobenius 应用)
由于 Leslie 矩阵是非负且通常不可约的，它存在一个唯一的正主特征值 \(\lambda_1\)。 1. 增长判定：若 \(\lambda_1 > 1\) 增长；若 \(\lambda_1 < 1\) 衰减。 2. 稳定分布：对应的正特征向量 \(\mathbf{v}_1\) 描述了长期演化后的稳定年龄比例。
70A.3 长期演化预测¶
技术：幂法投影
对于任何初始种群 \(\mathbf{x}_0\)，由于 \(\lambda_1\) 的支配性： $\(\mathbf{x}_t \approx c \lambda_1^t \mathbf{v}_1 \quad (t \to \infty)\)$ 这意味着无论初始状态如何，种群最终都会按照 \(\lambda_1\) 的速度扩张，且年龄结构固定为 \(\mathbf{v}_1\)。
练习题¶
1. [基础] 给定一个 \(2 \times 2\) Leslie 矩阵 \(L = \begin{pmatrix} 0 & 2 \\ 0.5 & 0 \end{pmatrix}\)。描述其生物含义。
参考答案
解析： 1. 第一行：幼体不育 (\(f_1=0\))，成体每单位产出 2 个后代 (\(f_2=2\))。 2. 次对角线：幼体到成体的存活率为 0.5 (\(s_1=0.5\))。 3. 成体在下一阶段不再存活（生命周期结束）。
2. [计算] 求上题中种群的增长率 \(\lambda\)。
参考答案
计算步骤： 1. 特征方程：\(\det(L - \lambda I) = \begin{vmatrix} -\lambda & 2 \\ 0.5 & -\lambda \end{vmatrix} = \lambda^2 - 1 = 0\)。 2. 解得 \(\lambda = \pm 1\)。 结论：主特征值为 1。这意味着种群规模长期保持恒定（不增不减）。
3. [稳定分布] 求上题中 \(\lambda=1\) 对应的特征向量（稳定年龄比例）。
参考答案
解方程组 \((L-I)\mathbf{v} = 0\)： 1. \(\begin{pmatrix} -1 & 2 \\ 0.5 & -1 \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \end{pmatrix} = 0 \implies v_1 = 2v_2\)。 2. 取 \(v_2 = 1, v_1 = 2\)。 结论：稳定年龄比例为 \(2:1\)。长期来看，种群中幼体数量是成体的 2 倍。
4. [生存分析] 若由于环境恶化，存活率 \(s_1\) 从 0.5 降至 0.4。种群会灭绝吗？
参考答案
判定： 1. 新矩阵 \(L' = \begin{pmatrix} 0 & 2 \\ 0.4 & 0 \end{pmatrix}\)。 2. 特征值 \(\lambda^2 - 0.8 = 0 \implies \lambda = \sqrt{0.8} \approx 0.89\)。 3. 由于 \(\rho(L') < 1\)。 结论：是的，种群将以每代约 11% 的速度萎缩，最终灭绝。
5. [不变量] 为什么 Leslie 矩阵的第一行不一定全是正数？
参考答案
理由： 生物在幼体期（juvenile）通常没有繁殖能力。因此 Leslie 矩阵的前几个对角元 \(f_1, f_2 \ldots\) 通常为 0。只要最终有某一年龄段能够繁殖，系统仍可能具有正的 Perron 特征值。
6. [计算] 证明：对于 Leslie 矩阵，特征方程可以写为 \(\sum_{i=1}^n f_i \left( \prod_{j=1}^{i-1} s_j \right) \lambda^{-i} = 1\)。
参考答案
证明要点： 利用拉普拉斯展开或直接观察状态转移：每一项代表了一个个体存活到 \(i\) 岁并产生后代的概率回馈。这个方程被称为欧拉-洛特卡方程的矩阵版。
7. [应用] 简述“渔业捕捞限额”如何通过矩阵特征值设定。
参考答案
策略： 1. 设原始增长率为 \(\lambda > 1\)。 2. 捕捞相当于给存活率矩阵乘上一个缩减因子 \(H = \operatorname{diag}(1-h_1, 1-h_2, \ldots)\)。 3. 目标是寻找捕捞率向量 \(\mathbf{h}\)，使得调整后的矩阵 \(L_{new}\) 的主特征值恰好为 1。 这样既能获得最大产量，又保证了种群资源的可持续性（平衡态）。
8. [性质] 判定：若 Leslie 矩阵不可约且是原初的（见 Ch17），种群是否会发生震荡？
参考答案
结论：不会。 理由：根据原初矩阵理论，主特征值 \(\lambda_1\) 是唯一的模最大特征值（即具有绝对支配性）。任何初始波动都会被指数级地抹平，种群平滑地收敛到稳定年龄分布。
9. [计算] 若 \(\mathbf{x}_0 = (100, 0)^T\)，使用题 1 中的矩阵 \(L\) 计算 \(\mathbf{x}_1, \mathbf{x}_2\)。
参考答案
迭代： 1. \(\mathbf{x}_1 = L \mathbf{x}_0 = \begin{pmatrix} 0 & 2 \\ 0.5 & 0 \end{pmatrix} \begin{pmatrix} 100 \\ 0 \end{pmatrix} = \begin{pmatrix} 0 \\ 50 \end{pmatrix}\)。 2. \(\mathbf{x}_2 = L \mathbf{x}_1 = \begin{pmatrix} 0 & 2 \\ 0.5 & 0 \end{pmatrix} \begin{pmatrix} 0 \\ 50 \end{pmatrix} = \begin{pmatrix} 100 \\ 0 \end{pmatrix}\)。 注：此例中 \(\lambda_2 = -1\)，导致种群在两代间发生震荡（非原初矩阵）。
10. [应用] 什么是“左弗科维奇矩阵”（Lefkovitch Matrix）？它与 Leslie 矩阵有何不同？
参考答案
对比： - Leslie 矩阵：基于年龄。个体每年必须且只能移动到下一个年龄组。 - Lefkovitch 矩阵：基于阶段（如幼虫、蛹、成虫）。个体可以在同一阶段停留，或者跳过阶段。 在代数上，Lefkovitch 矩阵的对角元可以是非零的，这为描述发育速度不均的复杂种群提供了更灵活的线性模型。
本章小结¶
线性代数是解析生物繁衍逻辑的“数字基因”：
- 生命的谱分析：通过特征值 \(\lambda_1\)，它将繁杂的生死数据提炼为一个单一的增长指标，确立了物种存续的数学判据。
- 平衡的必然性：特征向量理论证明了稳定年龄结构是环境与生理约束下的必然归宿，展示了动力系统向吸引子收敛的自然美。
- 干预的科学：从渔业配额到森林抚育，矩阵模型提供了对自然资源进行定量管理的唯一解析途径，是人类与生态系统和谐共生的理性指南。