昆虫繁殖问题.pptVIP

可修改 欢迎下载 可修改 欢迎下载 可修改 欢迎下载 可修改 欢迎下载 可修改 欢迎下载 * 欢迎下载 可修改 实验 昆虫繁殖问题 试验目的： 建立昆虫繁殖数学模型，掌握简单程序设计，运用线性代数知识解决实际问题 一、问题描述 有一种昆虫，最长寿命为六周，将其分为三组：第一组0~2周龄；第二组2~4周龄；第三组4~6周龄。 第一组为幼虫(不产卵)；第二组每只两周内平均产卵100个；第三组每只两周内平均产卵150个。每个卵的孵化成活率为9%，第一组能成长为第二组昆虫的存活率为10%，第二组能成长为第三组昆虫的存活率20%。 20% 10% 存活率 9% 9% 无 孵化率 150 100 不产卵 产卵数 4~6周龄 2~4周龄 0~2周龄 现有三个组昆虫各100只，计算第2周、第4周、第6周后各组昆虫数目，并考虑下面问题： (1)以两周为一时间段，分析各组昆虫数目变化规律。 昆虫数目是无限增长还是趋于灭亡？ (2)如果使用一种除虫剂，可以使得各组昆虫的成活率以及孵化成活率减半，问这种除虫剂是否有效？ 二、问题分析与与数学模型 初始时刻 两周以后 X(0) = [100, 100, 100]T X(1) =? x1(0) x2(1)=0.1x1(0) x2(0) x3(1)=0.2x2(0) x3(0) x1(1)=0.09(100x2(0)+150x3(0)) 一般, 有 根据数学模型计算数据结果如下 45 2 20 100 x3 36 225 10 100 x2 2050 360 2250 100 x1 k=3 6周后 k=2 4周后 k=1 2周后 现在 k 根据模型X(k) = Lk X(0) ,当k ? ? 时X(k)的变化是逐渐增大的趋势 MATLAB程序 n=input(‘input n(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) 使用除虫剂，可以使得各组昆虫的成活率以及孵化成活率减半，此时数学模型修改为 X(k+1) = L X(k) = L(k+1) X(0) (k = 0,1,2,3,······) MATLAB程序如下 n=input(‘input n(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^n*X0 %round(L^n*X0) 对不同的n进行计算实验，得数据结果如下 0.0105 0.0596 0.6138 k=12 24周后 5.6250 4.5000 256.50 k=3 6周后 0.50 56.25 90 k=2 4周后 10 5 1125 k=1 2周后 100 100 100 现在 x3 x2 x1 k 对问题的进一步思考 利用eig(A)和eig(L)命令可得特征值分别为: d1= 1.0731 d2= -0.7271 d3= -0.3460 d1= 0.5366 d2= -0.3635 d3= -0.1730 由此可知，An趋一无穷大，而Ln趋于0 数学原理:设A的三个特征值为?1, ?2, ?3,对应特征向量为?1, ?2, ?3 , 则有 令 P=[?1, ?2, ?3 ], D=diag([?1, ?2, ?3]),则 AP=PD = A= PDP-1 = An = (PDP-1 )······(PDP-1 )= PDnP-1 三、试验内容 1、建立三个周龄组的昆虫的种群动态分析的数学模型 建立数学模型并写出X(k)和X(k+1)的递推关系式，以及莱斯利矩阵L 2、根据数学模型计算并将数周后各组昆虫数填写到下表空白单元处 100 x3 100 x2 100 x1 k=25 50周后 k=2 4周后 k=1 2周后 现在 k 3、写出计算所用程序 4、在除虫剂作用下，情况将有所变化，重建数学模型 建立新的数学模型并写出X(k)和X(k+1)的递推关系式，以及新的莱斯利矩阵L 5、记录结果，分析除虫剂是否有效 k=10 10周后 100 x3 100 x2 100 x1 k=11 20周后 k=2 4周后 k=1 2周后 现在 k 结论： 1、最困难的事就是认识自己。20.6.286.28.202020:1120:11:15