实验22按年龄分组的种群增长模型..docVIP

实验　个年龄组，记，时段记作，且年龄组区间与时段长度相等（若5岁为一个年龄组，则5年为一个时段）。以雌性个体为研究对象比较方便，以下种群数量均指其中的雌性。 记第年龄组在时段的数量为；第年龄组的繁殖率为，表示每个（雌性）个体在一个时段内繁殖的数量；第年龄组的死亡率为，表示一个时段内死亡数与总数的比。是存活率。 为建立的变化规律，我们注意到:第1年龄组在时候的数量为各年龄组在第时段繁殖的数量之和，即 (22.1) 而第年龄组在时段的数量是第年龄组在时段存活的数量，即 (22.2) 记在时段种群各年龄组的数量为 。 (22.3) 这样，有 (22.4) 将归一化后的向量记做,称种群按年龄的分布向量。给定在0时段，各年龄组的初始数量，就可以预测任意时段各年龄组的数量。 设一种群分成5个年龄组, 已知繁殖率存活率 。各年龄组现有数量都是100只，下面我们用MATLAB计算。 % 按年龄分组的种群增长 clear all b=[0,0.2,1.8,0.8,0.2]; s=diag([0.5,0.8,0.8,0.1]); % 对角阵，对角元素为0.5,0.8,0.8,0.1 L=[b;s,zeros(4,1)]; % 构造矩阵L x(:,1)=100*ones(5,1); % 赋初值 K=45; for k=1:K x(:,k+1)=L*x(:,k); %迭代计算 end round (x), y=diag([1./sum(x)]); % 为向量x归一化做的计算 z=x*y, % z是向量x的归一化 k=0:K; subplot(1,2,1), plot(k,x),grid % 在一个图形窗内画两张图 subplot(1,2,2), plot(k,z),grid 将归一化后的向量记做，称为种群按年龄分组的分布向量，即各年龄组在时段在数量上占总数的百分比。 y=diag(1./sum(x)); %sum(x) 对列求和 Z=x*y subplot(1,2,2),plot(k,z),grid subplot(1,2,2),plot(k,z),grid 得到的结果见表22-1、表22-2和图22.1、图22.2。 表22-1 的部分计算结果 1 2 3 4 5 6 …… 40 41 42 43 44 45 100 300 220 155 265 251 …… 544 557 572 586 601 616 100 50 150 110 77 132 …… 265 272 279 286 293 301 100 80 40 120 88 62 …… 207 212 217 223 229 234 100 80 64 32 96 70 …… 161 165 170 174 178 183 100 10 8 6 3 10 …… 16 16 17 17 17 18 表22-2 的部分计算结果 1 2 3 …… 43 44 45 0.2000 0.5769 0.4564 …… 0.4559 0.4559 0.4558 0.2000 0.0962 0.3112 …… 0.2222 0.2223 0.2223 0.2000 0.1538 0.0830 …… 0.1734 0.1734 0.1734 0.2000 0.1538 0.1328 …… 0.1353 0.1353 0.1353 0.2000 0.0192 0.0166 …… 0.0132 0.0132 0.0132 结果分析 从上述图表可以看出，时间充分长以后种群按年龄分组的分布向量趋于稳定，这种状况与Leslie矩阵的如下性质有关（设矩阵第一行有两个顺序的大于零）： 矩阵有单特征根，对应特征向量为 （22.5） 对于的其他特征根有，且由（22.4）式确定的满足 ， （22.6） 其中是与，，有关的常数（请读者在矩阵可对角化的条件下证明（22.6）式）。 图22.1 的图形 图22.2从上到下依次为到的图形 由上述性质可以对时间充分长以后的，做出如下分析（以下记作）； 记归一化的特征向量为，则 （22.7） 与无关，