代数方程与差分方程模型（三）：按年龄分组的人口模型

考虑的问题

• 不同年龄组的繁殖率和死亡率不同
• 以雌性个体数量为对象

假设与建模

• 种群按年龄大小等分为n个年龄组，记i=1,2,…,n
• 时间离散为时段，长度和年龄组区间相等，记k=1,2,…
• 第i个年龄组1雌性个体在1时间段内发繁殖率为b_i
• 第i年龄组在1时段内的死亡率为d_i，存活率为s_i=1-d_i

$x_i(k):时间段k第i个年龄组的种群数量$
$x_i(k+1)={\sum }_{i=1}^n{b}_i{x}_i(k)(设至少1个b_i>0)$
$x_{i+1}(k+1)=s_i{x}_i(k),i=1,2,...,n-1$
把n个年龄组k时刻的向量记成：$x(k)=[x_1(k),x_2(k),...,x_n(k){]}^T$

可以表示为$x(k+1)=Lx(k)$
$$x(k)=L^{k}x(0)$$

稳定状态分析的数学知识

两个定理：

• L矩阵存在正单特征根λ₁，|λ_k|≤λ₁，k=2,3,…n(比其他所有特征根的模都不会小)
  特征向量$x^{\ast }=[1,\frac{s_1}{{\lambda }_1},\frac{s_1{s}_2}{{\lambda }_1^2},...,\frac{s_1{s}_2...s_{n-1}}{{\lambda }_1^{n-1}}]$
• 若L矩阵存在$b_i,b_{i+1}>0,则|{\lambda }_k|<{\lambda }_1,k=2,3,...,n$
  且$li{m}_{k\to \infty }\frac{x(k)}{{\lambda }_1^k}=c{x}^{\ast }，c是由b_i,s_i,x(0)决定的常数$

当k充分大，种群按年龄组的分布

1. $x(k)\approx c{\lambda }^k{x}^{\ast }$:种群按年龄组的分布趋向稳定，x^*称稳定分布，与初始分布无关。
2. $x(k+1)\approx \lambda x(k)$:各年龄组种群数量按同一倍数增减，λ称固有增长率
   推导出$x_i(k+1)\approx \lambda x_i(k)$
   与基本模型$x(k+1)=Lx(k)$比较发现k充分大，L矩阵可以用特征根λ代替
3. $\lambda =1时x(k+1)\approx x(k)\approx c{x}^{\ast }$各个年龄组的分布都是稳定的
   $x^{\ast }=[1,s_1,s_1{s}_2,...,s_1{s}_2...s_{n-1}{]}^T$,后一个年龄组的数量与前一个年龄组的数量之比是他们的存活率，各年龄组种群数量不变
   将L矩阵与x^*代入可以得到
   $$b_1+b_2{s}_1+...+b_n{s}_1{s}_2...s_{n-1}=1$$
   此式子表明1个个体在整个存活期内的繁殖数量为1
4. $x_{i+1}(k)\approx s_i{x}_i(k),i=1,2,...n-1$:表示存活率s_i是同一时段的x_i+1与x_i之比