实验昆虫繁殖问题.pptVIP

*/13实验昆虫繁殖问题第1页，共13页，星期日，2025年，2月5日莱斯利于1945年提出用于预测单种群生物数量增长的矩阵模型。将一个生物种群按年龄分为m个年龄组。设xk(t)表示t时刻第k个年龄组的生物数量,xk(0)是初始时刻数量。生物数量向量随时间t=0,t1,t2,t3,······变化规律用矩阵描述。即P.H.Leslie1900-1974第2页，共13页，星期日，2025年，2月5日一种昆虫最长寿命六周.分为三个年龄组,第一组幼虫(不产卵),第二组每只成虫平均两周产卵100，第三组每只成虫平均两周产卵150。假设每个卵的成活率为9%，第一组和第二组的昆虫两周后成为下一年龄组昆虫的存活率分别为10%和20%。?X(k+1)=LX(k)以两周为一时间段,设t0=0,t1=2,t2=4,t3=6,······.各年龄组昆虫数量为：x1(k)=x1(tk),x2(k)=x2(tk),x3(k)=x3(tk)第3页，共13页，星期日，2025年，2月5日x1(0)x2(0)x3(0)x2(1)=0.1x1(0)x3(1)=0.2x2(0)x1(1)=0.09(100x2(0)+150x3(0))?第4页，共13页，星期日，2025年，2月5日实验任务:以两周为一时间段，分析昆虫各年龄组数量向量变化规律.昆虫数量向量变化趋势是无限增长还是趋于灭亡？一种刹虫剂使用后可以控制各年龄组昆虫数量，使得各年龄组昆虫成活率减半，这种除虫剂是否有效？现有三个组的昆虫各100只,计算第2周、第4周、第6周后各个周龄的昆虫数目开始时刻X(0)=[100,100,100]T第5页，共13页，星期日，2025年，2月5日实验任务一:昆虫数量变化规律计算FunctionX=insect(n)X=[100;100;100];L=[0913.5;0.100;00.20];P=X;fork=1:nX=L*X;P=[P,X];endfigure(1),bar(P(1,:))figure(2),bar(P(2,:))figure(3),bar(P(3,:))调用函数X=insect(27)?X=7368.05686.52127.97第6页，共13页，星期日，2025年，2月5日L=[0913.5;0.100;00.20];[P,lamda]=eig(L)A的主特征值实验任务二:主特征值的特征向量试验p=-P(:,1);D=sum(p);X=[p(1)/D,p(2)/D,p(3)/D]*300P=-1.000.99-0.95-0.09-0.140.27-0.020.04-0.16Lamda=1.07000-0.73000-0.35X=270.1425.174.69昆虫数量向量按年龄显示出倒金字塔结构270.1425.174.69第7页，共13页，星期日，2025年，2月5日取n=27取n=3?270292310336252729315555functionP=insect(n)X=[270;25;5];P=X;L=[0913.5;0.100;00.20];fork=1:nX=L*X;P=[P,X];endK=0:n;figure(1),bar(K,P(1,:))figure(2),bar(K,P(2,:))figure(3),bar(K,P(3,:))第8页，共13页，星期日，2025年，2月5日主特征值:三个线性无关特征向量:??取初始时刻:通项:?···········第9页，共13页，星期日，2025年，2月5日实验任务三：使用除虫剂,各组昆虫的成活率以及孵化成活率减半,数学模型X(k+1)=L1X(k)(k=0,1,2,3,······)L1=？X(1)=?X(2)=?X(6)=?L1的