1实验昆虫繁殖问题试验目的:建立昆虫繁殖数学模型,掌握简单程序设计,运用线性代数知识解决实际问题2一、问题描述有一种昆虫,最长寿命为六周,将其分为三组:第一组0~2周龄;第二组2~4周龄;第三组4~6周龄。第一组为幼虫(不产卵);第二组每只两周内平均产卵100个;第三组每只两周内平均产卵150个。每个卵的孵化成活率为9%,第一组能成长为第二组昆虫的存活率为10%,第二组能成长为第三组昆虫的存活率20%。320%10%存活率9%9%无孵化率150100不产卵产卵数4~6周龄2~4周龄0~2周龄4现有三个组昆虫各100只,计算第2周、第4周、第6周后各组昆虫数目,并考虑下面问题:(1)以两周为一时间段,分析各组昆虫数目变化规律。昆虫数目是无限增长还是趋于灭亡?(2)如果使用一种除虫剂,可以使得各组昆虫的成活率以及孵化成活率减半,问这种除虫剂是否有效?5二、问题分析与与数学模型()()()()123()0123246kkkkkXxxxkkX以两周为一个时间段,用表示第个时间段各组昆虫数目。当、、、时,分别为初始时刻、第周后、第周后、第周后的各组昆虫数目。初始时刻两周以后X(0)=[100,100,100]TX(1)=?6x1(0)x1(0)x2(1)=0.1x1(0)x2(1)=0.1x1(0)x2(0)x2(0)x3(1)=0.2x2(0)x3(1)=0.2x2(0)x3(0)x3(0)x1(1)=0.09(100x2(0)+150x3(0))x1(1)=0.09(100x2(0)+150x3(0))7一般,有(1)(0)(0)123(1)(0)21(1)(0)320.091001500.10.2xxxxxxx(1)(0)11(1)(0)22(1)(0)330913.50.10000.20xxxxxx(1)(0)XAX(1)()2(1)(1)(0)==(0,1,2,3)kkkkXAXAXAXk8根据数学模型计算数据结果如下45220100x33622510100x220503602250100x1k=36周后k=24周后k=12周后现在k9根据模型X(k)=LkX(0),当k时X(k)的变化是逐渐增大的趋势MATLAB程序n=input(‘inputn(>3)’)A=[0,9,13.5;0.1,0,0;0,0.2,0],X0=[100;100;100],X1=A*X0,X2=A*X1,X3=A*X2,Xn=A^n*X0%若要取整:round(A^n*X0)10使用除虫剂,可以使得各组昆虫的成活率以及孵化成活率减半,此时数学模型修改为X(k+1)=LX(k)=L(k+1)X(0)(k=0,1,2,3,······)(1)()()123(1)21(1)()320.045(100150)0.05,0.1,(0,1,2,3)kkkkkkkxxxxxxxk04.56.750.050000.10L11MATLAB程序如下n=input(‘inputn(>3)’)L=[0,4.5,6.75;0.05,0,0;0,0.1,0],X0=[100;100;100],X1=L*X0,X2=L*X1,X3=L*X2,Xn=L^...