种群年龄结构的估算.ppt

种群年龄结构的估算 上海交大数学科学学院 数学实验 ■ 建立种群年龄结构的数学模型－Leslie矩阵 ■ 矩阵的特征值和特征向量知识的回顾 ■ Leslie矩阵有关性质的介绍和讨论 ■ 利用Matlab 进行矩阵特征值和特征向量 的运算 实验目的 ■ 进一步了解矩阵特征值在实际中的应用 一种昆虫每两周产卵一次，六周后死亡，孵化 的幼虫2周后成熟平均产卵100个，4周龄的成虫平均 平均产卵150个，设虫卵发育为2周龄成虫的概率为 0.09，2周龄成虫发育为4周龄成虫的概率为0.2 （1）假定开始时，0～2，2～4，4～6周龄昆虫数 相同，计算2周，4周，6周，8周后各周龄的昆虫数 （2）讨论各周龄昆虫数的变化趋势 实际问题 （3）若使用一种除虫剂控制昆虫数，已知药剂 将使各周龄该昆虫的成活率减半，问这种除虫剂是 否有效？ ? 这种昆虫在一个地区繁衍，其雌、雄数目之 比基本恒定，讨论时可将虫数设为雌虫数 ? 假设开始时三种周龄：C0:[0,2),C1:[2,4) C2:[4,6) 的昆虫数均为一个单位 情况假定 设 x0(t),x1(t),x2(t)分别为在t 时刻周龄 C0 , C1,C2 的昆虫数，以2周为时间单位，那么 引进矩阵 数学模型 导出关系式 投影矩阵 显然在时刻 2k 周( k 个时间单位) 2k 周后昆虫数 键入 A=[0 9 13.5 ; 0.09 0 0; 0 0.2 0]; X0=[1 1 1]; A*X0, A^2*X0, A^3*X0 结果情况 2周后 22.5 0.09 0.2 C0 C1 C2 利用 Matlab 看一看 4周后 3.51 2.025 0.018 6周后 18.468 0.3159 0.405 ? 可以用数值计算观察结果(作为练习) ? 理论分析 昆虫数的趋势 k越来越大时，X(k)=AkX(0)变化趋势如何？ 若A有特征值λ0 ,λ1,λ2 ,对应特征向量α0,α1,α2 线性无关，其中λ0的模严格大于其它特征值的模 那么可表示 X(0)=c0α0+c1α1+c2α2 ? AkX(0)=c0λ0kα0+c1λ1k α1+c2k λ2α2 ? AkX(0) 这意味着k充分大时，昆虫数 记 (s0, s1, s2 是 α0的分量) 各周龄昆虫数占昆虫总数之比趋于定值 [v,d]=eig(A) Matlab 实现 求A的特征值和特征向量 有最大模的特征值为 ?0 ＝1.0234 对应的特征向量?0 ＝ 0.9960 0.0976 0.0171 0.9049 0.0796 0.0156 归一化 ? 2) 最终各周龄的昆虫数成稳定的比例 随着时间增长 1) 昆虫数成几何级数增长 ( ?0 1 ) (恰好是对应主特征值的特征向量各分量比例数) 结论 分别所占比例： 90.49%，7.96%，1.56% ? 经多少时间才能达到稳定的增长和结构比例？ 依然可以采取观察数值计算的结果的方法 t X0 X1 X2 ∑ q1 q2 q3 0 1 1 1 3 0.3333 0.3333 0.3333 1 22.5000 0.0900 0.2000 22.7900 0.9873 0.0039 0.0088 2 3.5100 2.0250 0.0180 5.5530 0.6321 0.3647 0.00352 3 18.4680 0.3159 0.4050 19.1889 0.9624 0.0165 0.0221 4 8.3106 1.6621 0.0632 10.0359 0.8281 0.1656 0.0063 5 15.8120 0.7480 0.3324 16.8924 0.9360 0.0443 0.0197 10 14.0763 1.3263 0.2327 15.6353 0.9003 0.0848 0.0149 15 16.3044 1.4220 0.2814 18.0078 0.9054 0.0790 0.0156 20 18.2407 1.6056 0.3133 20.1597 0.9148 0.0796 0.0155 28 21.9640 1.9315 0.3774 24.2729 0.9049 0.0796 0.0155 29 22.4793 1.9768 0.3863 24.8424 0.9049 0.0796 0.0156 30 23.0060 2.0231 9.3954 25.4245 0.9049 0.0796 0.0156 在Matlab程序设定误差，计算到某步停止 由于昆虫成活率减半引起投影矩阵变化