2018年动物的繁殖问题.docVIP

动物的繁殖问题 　　实验六 动物的繁殖问题 　　【实验目的】 　　1．加深对矩阵线性变换、方幂、特征根、对角化等概念的理解。 2．掌握用矩阵变换在实际问题中的应用。 　　3．学习使用MATLAB软件中与矩阵有关的命令。 　　【实验内容】 　　某农场饲养的某种动物所能达到的最大年龄为15岁，将其分为三个年龄组：第一组0～5岁；第二组6～10岁；第三组11～15岁。动物从第二个年龄组开始繁殖后代，第二个年龄组的动物在其年龄段平均繁殖4个后代，第三个年龄组的动物在其年龄段平均繁殖3个后代。第一年龄组和第二年龄组的动物能顺利进入下一个年龄组的存活率分别为0.5和0.25。假设农场现有三个年龄段的动物各有1000头，计算5年后、10年后、15年后各年龄段动物数量。20年后农场三个年龄段的动物的情况会怎样？ 　　根据有关生物学研究结果，对于足够大的时间值k，有X如果每五年平均向市场供应动物数c＝?s提下，c应取多少为好？ 　　(k?1) 　　≈?1X 　　(k) 　　（?1是莱斯利 　　矩阵L的惟一正特征值）。请检验这一结果是否正确，如果正确给出适当的k的值。 　　T 　　ss?，在20年后农场动物不至灭绝的前 　　【实验准备】 　　1．矩阵知识回顾 　　当矩阵的列数与某一个列向量元素个数一致时，用矩阵乘以向量将得到另一个向量，这就是向量的线性变换。当矩阵是方阵时，线性变换可以持续进行。即用矩阵乘以一个向量得到一个新的向量，用同一矩阵再乘以新的向量又获得另一个的向量??，这种运算的本质就是用矩阵的方幂乘以最初的向量。在线性代数应用中称为矩阵的方幂问题，它和矩阵的特征根问题有密切关系，对它的研究导致了矩阵对角化方法，这类方法在生物学研究等领域有着广泛应用。 　　设A＝(aij)是数域F上的一个n阶矩阵，行列式 　　x?a11 　　?a21 　　fA(x)＝xI?A＝ 　　... 　　?an1 　　叫作矩阵A的特征多项式，若有 fA(?)＝0 　　?a12x?a22 　　...?an2 　　......... 　　?a1n 　　?a2n 　　（1） ... 　　...x?ann 　　那么?就是矩阵A的特征多项式的特征根，那么方程 　　AX＝?X （2） 的一个非零解叫做矩阵A的属于特征根?的特征向量。 　　定理1 令A是数域F上的一个n阶矩阵，如果A的特征多项式fA(x)在F内有n个 根，那么存在一个n阶可逆矩阵T，使得 　　??10? 　　?0?2?1 　　TAT＝? 　　......? 　　?00? 　　5．求矩阵特征根的MATLAB命令 　　... 　　0??...0? 　　（3） ?......?...?n?? 　　【实验方法与步骤】 　　1．引例问题的分析 　　由题设，在初始时刻0～5岁、6～10岁、11～15岁的三个年龄段动物数量分别为： 　　(0)(0)(0) 　　x1＝1000， x2＝1000， x3＝1000 　　以五年为一个年龄段，则某一时刻三个年龄段的动物数量可以用一个向量 X＝[x1表示。以五年为一个时间段，记 X 　　(k) 　　x2 　　x3]T 　　(k) 　　x2 　　(k)Tx3] 　　(k) 　　(k) 　　＝[x1 　　为第k个时段动物数分布向量。当k＝0，1，2，3时，X分别表示现在、五年后、十年 后、十五年后的动物数分布向量。根据第二年龄组和第三年龄组动物的繁殖能力，在第k个时间段，第二年龄组动物在其年龄段平均繁殖4个，第三年龄组动物在其年龄段平均繁殖3个后代。由此得第一年龄组在第k＋1个时间段的数量如下 x1 　　(k?1) 　　(k) 　　＝4x2＋3x3 　　(k) 　　同理，根据第一年龄组和第二年龄组的存活率，可得等式 x2建立数学模型如下 x1 　　(k?1)(k?1)(k?1) 　　(k?1)＝0.5x1， x3＝0.25x2 　　(k)(k) 　　(k)＝4x2＋3x3 　　(k) 　　x2改写成矩阵形式 　　＝0.5x1 （k＝0，1，2，3） 　　(k) 　　(k) 　　(k?1) 　　x3＝0.25x2 　　?x1(k?1)??043??x1(k)? 　　?(k?1)????x(k)? （k＝0，1，2，3） 　　00 ?x2??0.5???2? 　　(k?1)?(k)??x3?x3??00.250?????? 　　由此得向量X 　　(k) 　　和X 　　(k?1) 　　的递推关系式 　　(k?1) 　　X其中矩阵 　　＝LX 　　(k) 　　　　43??0 　　?? 　　00