完整word版8动物繁殖问题.docx

1、完整word版8动物繁殖问题一、实验目的 了解矩阵的方幂和矩阵的特征值的应用 二、实验内容 1. 动物繁殖的规律问题 2. 商品的市场占有率问题 3. 常染色体遗传问题 三、实验仪器和设备 1. 计算机若干台（装有matlab6.5及以上版本软件） 2. 打印机 四、实验要求 1. 独立完成各个实验任务； 2. 实验的过程保存成 .m 文件，以备检查； 3. 实验结果保存成 .mat 文件 五、实验原理 当矩阵的列数与某一个列向量元素个数一致时，用矩阵乘以向量将得到另一向量，这就是向量的线性变换。当矩阵是方阵时，线性变换可持续进行。即，用矩阵乘以一个向量得一个新的向量，用同一矩阵再乘以新的向量

2、又获得另一新的向量， ，这种运算的本质是用矩阵的方幂乘以最早的哪一个向量。在线性代数应用中称为矩阵的方幂问题，它和矩阵的特征值问题有密切关系。对它的研究导致了矩阵对角化方法，这类方法在生物学研究等方面应用广泛。本章介绍的几个实际问题包括动物繁殖的规律问题、商品的市场占有率问题、常染色体遗传特征问题。 （一）动物繁殖的规律问题 某农场饲养的某种动物所能达到的最大年龄为15岁，将其分为三个年龄组：第一组05岁；第二组610岁；第三组1115岁。动物从第二个年龄组开始繁殖后代，第二个年龄组的动物在其年龄段平均繁殖4个后代，第三年龄组的动物在其年龄段平均繁殖3个后代。第一年龄组和第二年龄组的动物能顺利

3、进入下一个年龄组的存活率分别为0.5和0.25。假设农场现有三个年龄段的动物各1000头，计算5年后、10年后、15年后各年龄段动物数量。20年后农场三个年龄段的动物的情况会怎样？ 根据有关生物学研究结果，对于足够大的时间值k，有 （ 是莱斯利矩阵L的唯一正特征值）。请检验这一结果是否正确，如果正确给出适当的k的值。 如果每五年平均向市场供应动物数c=s s sT，在20年后农场动物不至灭绝的前提下，c应取多少为好？ 1 .问题分析和数学模型 由题设，在初始时刻05岁、610岁、1115岁的三个年龄段动物数量分别为： =1000， =1000， =1000以五年为一个年龄段，则某一时刻三个年龄

4、段的动物数量可以用一个向量 X=x1 x2 x3T 表示。以五年为一个时间段，记 X(k) = x1(k) x2(k) x3(k)T 为第k个时间段动物数分布向量。当k= 0，1，2，3时，X(k) 分别表示现在、五年后、十年后、十五年后的动物数分布向量。根据第二年龄组和第三年龄组动物的繁殖能力，在第k个时间段，第二年龄组动物在其年龄段平均繁殖4个后代，第三年龄组动物在其年龄段平均繁殖3个后代。由此得第一个年龄组在第k+1个时间段的数量如下 同理，根据第一年龄组和第二年龄组的存活率，可得等式 建立数学模型如下 (k = 0，1，2，3) （1） 或写成矩阵形式 (k = 0，1，2，3) （2

5、） 由此得向量X(k)和X(k+1)的递推关系式 X(k+1) = LX(k) (3) 其中，矩阵 称为莱斯利矩阵。由式(3)可得 X(k+1) = L k+1X (0) 2 . 程序和计算结果 为了计算5年后、十年后、十五年后农场中动物的数量，输入初始数据和莱斯利矩阵在MATLAB中键入下面命令 x0=1000;1000;1000;A=0 4 3;1/2 0 0;0 1/4 0;x1=A*x0x2=A*x1x3=A*x2x4=A*x3可得数据结果 x1 = 7000 500 250x2 = 2750 3500 125x3 = 14375 1375 875x4 = 1.0e+003 *8.12

6、50 7.1875 0.3438（x4的数据结果中，1.0e+003是科学计数法。表示用103乘后面的每一个数） 为了计算莱斯利矩阵的特征值，键入下面命令 eig(A)得数据 ans =1.5000 -1.3090 -0.1910这说明矩阵A的唯一正特征值为 =1.5为了验证 运行下面程序 x=1000;1000;1000;d1=1.5;A=0 4 3;1/2 0 0;0 1/4 0;y=A*x;y1=d1*x;k=1;while max(abs(y-y1).1x=y;y=A*x;y1=d1*x;k=k+1;end可知，当K=291时，有结论 成立。 3 . 问题解答与进一步思考 根据数学模型

7、计算将数据结果填写如下 表1k (现在) k=1 (五年后) k=2(十年后) k=3(十五年后) k=4(20年后) x1 1000 7000 2750 14375 8125 x2 1000 500 3500 1375 7187.5 x3 1000 250 125 875 343.8 从表中数据变化，如果没有其它的原因，估计农场的动物总数量会逐步增加。 在验证生物学研究的结论 时，当k=291可以得到如下结论 x = 1.0e+054 *3.3121 1.1040 0.1840这说明多年以后，动物数量是大得非常惊人的。 如果每个五年平均向市场供应动物c=s s sT，分析动物数分布向量变化规

8、律可知 X(1) = AX(0) c X(2) = AX(1) c X(3) = AX(2) c X(4) = AX(3) c 所以有 X(4) = A4X(0) ( A3 + A2 + A + I )c 考虑二十年后动物不灭绝，应有 X(4) 0 即 ( A3 + A2 + A + I )c A4X(0) 由于c是常数向量，故可简单求解不等式组，可取c= 152 152 152 T 这说明当五年平均向市场供应三个年龄段的动物各152头可以使20年后有各年龄段的动物生存。如果将这一限制作为约束条件，而求c的各分量之和最大，这将是一个线性规划问题，可用单纯形法求解（或利用数学软件见第七章）。 （

9、二）商品的市场占有率问题 有两家公司R和S经营同类的产品，它们相互竞争。每年R公司保有1/4的顾客，而3/4转移向S公司；每年S公司保有2/3的顾客，而1/3转移向R公司。当产品开始制造时R公司占有3/5的市场分额，而S公司占有2/5的市场分额。问两年后，两家公司所占的市场分额变化怎样，五年以后会怎样？十年以后如何？是否有一组初始市场分额分配数据使以后每年的市场分配成为稳定不变？ 1 、问题分析和数学模型 根据两家公司每年顾客转移的数据资料，形成以下转移矩阵 根据产品制造之初，市场的初始分配数据可得如下向量 所以一年后，市场分配为 两年后，市场分配为 以向量Xn记第n年后市场分配的分额，则 （

10、n = 1，2，） 设有数据a和b作为R公司和S公司的初始市场分额，则有 a + b = 1 为了使以后每年的市场分配不变，根据顾客数量转移的规律，有 即 这是一个齐次方程组问题。如果方程组有解，则应该在非零解的集合中选取正数解作为市场稳定的初始分额。 2、程序和计算结果 为了知道两年、五年、十年后市场分配的情况，在MATLAB中键入下面命令 A=1/4 1/3;3/4 2/3x0=3/5;2/5x2=A2*x0x5=A5*x0x10=A10*x0可得数据结果 x2 = 0.3097 0.6903x5 = 0.3077 0.6923x10 = 0.3077 0.6923由此得下表 表2 R公司

11、的市场分额 S公司的市场分额 两年后 31% 69% 五年后 31% 69% 十年后 31% 69% 为了求a和b作为R公司和S公司稳定的初始市场分额，需要求解齐次方程组。键入下面命令 format ratrref(A-eye(2)得数据结果 ans =1 -4/9 0 0 由此得化简后的方程 a 4/9b = 0 结合约束条件 a + b = 1 得 a= 4/1331% b= 9/1369% 这是使市场稳定的两家公司的初始分额，也正好与表4-1中的数据相吻合。 3、问题的解答和进一步思考 在R公司和S公司的市场初始分额分别为60%和40%的情况下，根据计算结果，两年后情况变化较大：S公司大

12、约占31%，R公司大约占69%。而五年以后与两年以后比较变化不大：S公司大约占30.8%，R公司大约占69.2%。十年后的情况与五年后的情况比较大约不变。市场已经趋于稳定。 是否所有市场初始分配分额，在经过若干年后均会趋于稳定状态。 （三）常染色体遗传问题 假定所考虑的遗传特性由两个基因A和a来支配，人类的眼睛染色体是通过常染色体遗传来控制，例如AA及Aa型产生棕色眼睛，aa型的是兰色眼睛。在常染色体遗传中，一个个体从它的亲本的每一基因对中遗传一个基因，以形成它自己特殊的基因对：AA，Aa，aa。亲本的两个基因中的哪一个传给后代纯属机会问题，如果一个亲本是Aa型，后代从这个亲本遗传获得A基因或

13、a基因的机会是等可能的。例如，一个亲本是aa型，另一个亲本是Aa型，后代总是从aa亲本接受一个a基因，再从Aa亲本以等概率或是接受一个A基因或是接受一个a基因，结果后代为aa型或者Aa型的概率是相同的。对于各种亲本基因型，后代的可能基因型的概率可列表如下 表4.3亲本 后代 AA-AA AA-Aa AA-aa Aa-Aa Aa-aa aa-aa AA 1 1/2 0 1/4 0 0 Aa 0 1/2 1 1/2 1/2 0 aa 0 0 0 1/4 1/2 1 例 假定一个农民有一大片作物，它由三种可能基因型AA，Aa及aa的某种分布所组成。农民要采用的育种方案是：作物总体中的每种作物都总是用

14、基因型AA的作物来授粉，我们要导出在任何一个后代总体中三种可能基因型的分布表达式。 解 记a n (n = 0，1，2，)为在第n代中AA基因型作物所占的分数，bn为在第n代中Aa基因型作物所占的分数，cn为在第n代中aa基因型作物所占的分数。a0，b0，c0表示基因型的原始分布，且 a0 + b0 + c0 = 1 由于用基因型AA的作物来授粉，分析基因表（前三列数据）可知，从上一代的基因型分布产生的下一代的基因型分布可用下列递推公式求出： 其中，第一式表明，基因型AA的所有后代都是AA型基因，基因型Aa的后代，有一半是AA型。这一递推公式的矩阵表示为 X(n) = MX(n - 1) (

15、n = 1，2，) 其中 ， ， 由递推公式可得 X(n) = MX(n - 1) = M 2X(n - 2) = = M nX(0) 计算上式有两种方法，即直接计算和将矩阵对角化的计算方法。对角化方法需要将矩阵M对角化，需要找出一个可逆矩阵P和一个对角阵D，使 M = PDP 1 于是 M n = PD nP 1， (n = 1，2，) 其中 d1，d2，d3，是M的特征值。所以，只需求得M的特征值和对应的特征向量，就可使M对角化。在MATLAB环境中输入命令 M=1 1/2 0；0 1/2 1；0 0 0； p d=eig(M)得数据结果 p =1 -985/1393 881/2158 0

16、 985/1393 -881/1079 0 0 881/2158 d =1 0 0 0 1/2 0 0 0 0 这表明，M的三个特征值为 d1 = 1，d2 = 1/2，d3 = 0 因为特征向量乘一非零数仍是特征向量，所以可取三个特征值对应的特征向量分别为 ， ， 于是 可逆矩阵 为了求逆矩阵，使用命令 P=1 1 1；0 -1 -2；0 0 1； inv(P)可得数据结果 ans =1 1 1 0 -1 -2 0 0 1 所以 由前面递推公式，得 X(n) = PD nP 1X(0) 而 故 所以 这是原始基因型分数表示第n代作物总体中三种基因型分数。显然，当n时，有 an(a0 + b0

17、 + c0)=1，bn0，cn0这说明在极限情况下，总体中所有作物都将是基因型AA的。 假设一片作物是由AA，Aa，及aa基因型的某种分布组成，且作物总体中每种作物不是全部都用基因型AA授粉，而是用每种作物自身的基因型来授粉。求任何一个后代总体中三种可能基因型的分布表达式。 六、实验任务 1某一种甲虫最多可活两年。且其年龄群体分配数的矩阵如下： 如果有600只在第一年龄群体，300只在第二年龄群体，100只在第三年龄群体，则年复一年各年龄群体的甲虫数目是否会改变，从数学上给以解释。 2假设某一个城市的气候不是下雨就是干旱。根据以前所保留下来的记录可知，干旱天之后下雨天为的可能性为1/3，而下雨

18、天之后为下雨天的可能性为1/2。试建立数学模型分析气候变化情况。 3某实验性生产线每年一月份进行熟练工人的人数统计；然后将其1/6的熟练工人支援其他生产部门，缺额由招收非熟练工人补齐。新、老非熟练工人经过培训及实践至年终考核有2/5成为熟练工人。设第n年一月份统计的熟练工人和非熟练工人所占的百分比分别为xn和yn记为向量 xn，ynT(1) 试推导向量xn，ynT和xn+1，yn+1T的关系，并写成矩阵形式； (2) 当x1，y1T=0.5，0.5T时，求第10年一月份统计的熟练工人和非熟练工人所占的百分比。 5在一城市的某商业区内，有两家有名的快餐店“肯德基”分店和“麦当劳”分店。据统计每年

19、“肯德基”保有其上一年老顾客的1/3，而另外的2/3顾客转移到“麦当劳”；每年“麦当劳”保有其上一年的老顾客的1/2，而另外的1/2顾客转移到“肯德基”。 用二维向量X k =xk ykT表示两个快餐店市场分配的情况，初始的市场分配为 X0 = 1/3 2/3T 如果有矩阵L 存在，使得 Xk+1 = LXk，则称 L 为状态转移矩阵。 (1)写出Xk=xk ykT和Xk+1=xk+1 yk+1T的递推关系式，以及状态转移矩阵L。 (2)根据递推关系计算近几年的市场分配情况； (3)求可逆矩阵P和对角矩阵使得L=PP 1。对于足够大的时间值k，有 （ 是状态转移矩阵L的唯一正特征值）。请检验这

20、一结果是否正确，并给出适当的k的值 6某厂生产A，B两种品牌的味精，顾客的喜好决定了这两种味精的市场占有率。在生产中可根据占有率调整比例，获得最佳收益。该厂做市场调查后发现，一般情况下，顾客若购买A牌，下次有80%的可能性购买A牌；若购买了B牌，下次有60%的可能性购买B牌。开始时，两种品牌的市场占有率分别为50%，顾客每一次的购买必将改变二者市场占有率。 （1）预测某一个顾客经过前四次购买之后，他可能第五次购买哪一个品牌的味精。 （2）预测100个顾客经过前四次购买之后，两种品牌的可能市场占有率各为多少？ 7足球比赛排名问题： 下表给出了我国 8 支足球队在 19881989 年全国足球甲级

21、队联赛中的成绩，请依据现有的数据，对各足球队的实力给予科学客观的评价；给出 8 支球队按实力排序的名次，并将你的排名算法推广到任意 n 个队的情形。 T1 T2 T3 T4 T5 T6 T7 T8 T1 0:11:02:2 2:21:00:2 2:03:11:0 0:11:3 0:22:1 1:04:0 1:11:1 T2 1:00:12:2 2:00:11:3 0:02:00:0 1:11:1 0:00:0 2:01:1 0:20:0 T3 2:20:12:0 0:21:03:1 4:21:10:0 1:01:4 0:13:1 1:02:3 0:12:0 T4 0:21:30:1 0:00:20:0 2:41:10:0 0:52:3 2:11:3 0:10:0 0:11:1 T5 1:03:1 1:11:1 0:14:1 5:03:2 1:02:00:0 2:13:01:0 3:13:02:2 T6 2:01:2 0:00:0 1:01:3 1:23:1 0:10:20:0 0:11:22:0 1:11:00:1 T7 0:10:4 0:21:1 0:13:2 1:00:0 1:20:30:1 1:02:10:2 3:01:00:0 T8 1:11:1 2:00:0 1:00:2 1:01:1 1:30:32:2 1:10:11:0 0:30:10:0