实验四动物鳘殖.pptVIP

*/13 实验四:动物养殖问题 莱斯利矩阵模型 实验任务与操作 思考题与练习题 ? ? ? 某种动物最大年龄为15岁,将其分为三个年龄组：第一组0—5岁；第二组6—10岁；第三组11—15岁。第二组在其年龄段平均繁殖4个后代,第三组平均繁殖3个后代。第一和第二组五年的存活率分别为0.5和0.25。现有三个年龄组动物各1000，计算5年后、10年后、15年后各年龄组动物数量。 ? X(k+1)=L X(k) 设 t0 = 0, t1 = 5, t2 = 10, t3 = 15. 各年龄组动物数量 x1(k)=x1(tk), x2(k)= x2(tk), x3(k)= x3(tk) x1(0) x2(0) x3(0) x2(1)=0.5x1(0) x3(1)=0.25x2(0) x1(1)=4x2(0)+3x3(0)) ? 其中L 称为莱斯利矩阵。 莱斯利于1945年提出用于预测单种群生物数量增长的矩阵模型。 将一个生物种群按年龄分为 m 个年龄组。设 xk( t ) 表示 t 时刻第 k 个年龄组的生物数量, xk(0)是初始时刻数量。生物数量向量 随时间 t = 0, t1, t2, t3, ······ 变化规律用矩阵 描述。即 P.H.Leslie 1900-1974 莱斯利矩阵一般形式： 实验任务: 以五年为一时间段，分析动物各年龄组数量变化规律. 动物总数量变化趋势是无限增长还是趋于灭亡？ 3*.如果每五年向其它养殖场输送动物C=[s1 s2 s3]T 要求20年后本养殖场动物不灭绝，C 取多少为好？ 4*.设 为L的主特征值，则k取多大时，有近似式 （要求误差不大于0.0001）。 现有三个组的动物各1000,计算第5年、第10年、第15年后各个周龄的动物数量 开始时刻 X(0) = [1000, 1000, 1000]T 实验任务一:动物数量变化规律计算 function X=animal(n) L=[0 4 3;0.5 0 0;0 0.25 0]; X=[1000;1000;1000]; P=X; for k=1:n X=L*X;P=[P,X]; end figure(1),bar(P(1,:)) figure(2),bar(P(2,:)) figure(3),bar(P(3,:)) 调用函数 X=animal(12) ? X =314754.15 143543.21 16547.12 主特征值: 三个线性无关特征向量: ? ? 取 初始时刻: 通项: ? ··········· 为了验证 　　　　　　　　　　 d=1.5; x=[1000;1000;1000]; L=[0 4 3;1/2 0 0;0 1/4 0]; y=L*x; y1=d*x; k=1; while max(abs(y-y1))0.01 x=y; y=L*x; y1=d*x; k=k+1; end k 编辑M文件animal.m： 实验任务三：每五年平均向市场供应动物 C = [s1 s2 s3]T 修正数学模型 X(k+1) = L X(k) – C , (k = 0, 1,2,3) X(1) = L X(0) – C, X(2) = L X(1) – C X(3) = L X(2) – C, X(4) = L X(3) – C X(2) = L2 X(0) – LC – C X(3) = L3 X(0) – L2C – LC – C X(4) = L4 X(0) – L3C – L2C – LC – C X(4) = L4 X(0) – (L3 + L2 + L + I ) C L=[0 4 3;0.5 0 0;0 0.25 0]; [P,lamda]=eig(L) L 的主特征值 主特征值特征向量试验与注记 p1=P(:,1);d=sum(p1); p=p1/d X0=p*3000 P= 0.95 0.93 0.23 0.32 -0.36 -0.59 0.05 0.07 0.77 Lamda= 1.50 0 0 0 -1.31 0 0 0 -0.19 X0= 2160.00 720.00 120.00 动物数量按年龄显示出倒金字塔结构 2160 720 120 取 n=3 ? 2160. 3240. 4860. 7290 720. 1080. 1620.