差分模型(虫子)080.docVIP

按年龄分组的种群增长 摘要：为了研究种群的结构增长模型，如果不考虑种群的年龄结构，种群的数量主要由总量的固有增长率决定。从一个实际问题着手，建立差分方程模型，讨论稳定状况下种群的增长规律，即其总数量与各年龄结构模型。 【关键字】 种群增长规律 差分方程 矩阵 繁殖率 死亡率 问题重述 一种昆虫每两周产卵一次，6周以后死亡，孵化后幼虫2周后成熟，平均产卵100个；4周龄的成虫平均产卵150个。假设只有0.09的卵可以孵化。孵化出来的幼虫有0.09可以发育成2周龄的成虫，2周龄的成虫发育成4周龄的成虫的概率是0.2。 建立差分方程模型，讨论稳定状况下（时间充分长）种群的增长规律（总数量与各年龄结构），然后分析各年龄组的繁殖率与死亡率对种群总量增减的影响。 模型假设 （1）种群按年龄大小等分为 n 个年龄组，记 i 1,2,… , n （2）时间离散为时段, 长度与年龄组区间相等, 记k 1,2,… （3）第i 年龄组1个体在1时段内的繁殖率为 bi （4）第i 年龄组在1时段内的死亡率为 di , 存活率为 si 1－ di 问题分析 对于本文所述的问题，利用昆虫不同时期的产卵率和死亡率，可以分析出该昆虫的一些情况，并可根据所得情况进行分析，算出该昆虫多年以后的生长情况和年龄组成，可以建立一个按年龄分组的种群增长的差分方程模型，根据差分方程的 Leslie矩阵的特征根，结合 Leslie矩阵的稳定性理论对当前昆虫种群的情况进行分析。为了保持该昆虫种群的稳定，必须使得Leslie矩阵的最大特征根为1，这样就可求出所求。 模型的建立和求解 假设时段第年龄组的种群数量，则 可得矩阵 按年龄分组可得分布向量~ 预测任意时段种群按年龄组的分布 经过对其稳定状态的分析可得L矩阵存在正单特征根 1，其中实部最大就，记为。 其特征向量为 ， 经验证可得若L矩阵存在, 则，且，其中是由决定的常数。 其中种群按年龄组的分布趋向稳定，x*称稳定分布, 与初始分布无关，记为 且各年龄组种群数量按同一倍数增减， 称固有增长率， ，与基本模型相比较，可以推出。从而我们可以确定当时，，即各年龄组种群数量不变，即。 因为当时， ，， 此时， 将其当为1个个体在整个存活期内的繁殖数量为1，可得。 在将年龄组分为3组的模型内记为。 当时，同理可得，其中 反映了种群增长或是衰减。 由可得存活率，即，其中是同一时段的与 之比 现针对本文所述问题建立模型： 由于种群是通过雌性个体的繁殖而增长的，所以忽略雌雄昆虫的区别，把所有雌性个体分成3个年龄组： 根据题设可得： 若记，则 在Matlab命令窗口输入矩阵L，便可求出特征向量U,V，如下： L [0 9 13.5; 0.09 0 0; 0 0.2 0]; [U,V] eig L U 0.9960 -0.9903 0.9603 0.0876 0.1334 -0.2432 0.0171 -0.0399 0.1368 V 1.0234 0 0 0 -0.6680 0 0 0 -0.3554 程序编写： function out worm n out zeros 3,n ； L [0 9 13.5; 0.09 0 0; 0 0.2 0]; out :,1 [100,100,100]; for i 2:n out :,i L*out :,i-1 end 运行结果： （1）利用程序便可知n年后，小虫种群的数量情况，将结果描绘称图像如下所示（以下讨论100年时的结果）： 在命令窗口输入： x worm 100 ; plot x 1,: ,or- ; hold on plot x 2,: ,^- ; hold on plot x 3,: ,dk- ; 可得图像如下： （红线为0~2周的虫，蓝线为2~4周的虫，黑线为4~6周的虫） 从图中我们可以发现，虫未灭绝。当时间足够长时，可观察到不同年龄阶段虫的数量变化如上图所示：0~2周的虫的数量急剧增多，2~4周虫的数量也增多，而4~6周虫的数量相对较少，但三者间无太多比例关系。最总导致整个种群数量增多。 （2）利用程序便可知n年后，小虫种群的数量情况，将结果描绘称图像如下所示（以下讨论1000年时的结果）： 在命令窗口输入： plot x 1,[1:100:end] ,or- ; hold on plot x 2,[1:100:end] ,^- ; hold on plot x 3,[1:100:end] ,dk- ; 可得图像如下： （红线为0~2周的虫，蓝线为2~4周的虫，黑线为4~6周的虫） 从图中我们可以发现，虫未灭绝。当时间足够长时，可观察到不同年龄阶段虫的数量变化： 从图中我们可