2022年8动物繁殖问题

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

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

3、 和 0.25。 假设农场现有三个年龄段的动物各1000头，计算 5 年后、10 年后、15 年后各年龄段动物数量。 20 年后农场三个年龄段的动物的情况会怎样？根据有关生物学研究结果，对于足够大的时间值k，有（是精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 12 页 - - - - - - - - -莱斯利矩阵 l 的唯一正特征值） 。请检验这一结果是否正确， 如果正确给出适当的 k 的值

4、。如果每五年平均向市场供应动物数c=s s st，在 20 年后农场动物不至灭绝的前提下， c 应取多少为好？1 . 问题分析和数学模型由题设，在初始时刻05 岁、610岁、1115岁的三个年龄段动物数量分别为：=1000，=1000，=1000 以五年为一个年龄段，则某一时刻三个年龄段的动物数量可以用一个向量x= x1 x2 x3 t 表示。以五年为一个时间段，记x( k) = x1( k) x2(k) x3(k) t为第 k 个时间段动物数分布向量。当k= 0，1，2，3 时，x(k) 分别表示现在、五年后、十年后、十五年后的动物数分布向量。根据第二年龄组和第三年龄组动物的繁殖能力，在第

5、k 个时间段，第二年龄组动物在其年龄段平均繁殖4 个后代，第三年龄组动物在其年龄段平均繁殖3 个后代。由此得第一个年龄组在第k+1个时间段的数量如下同理，根据第一年龄组和第二年龄组的存活率，可得等式建立数学模型如下(k = 0，1，2，3) （1）或写成矩阵形式 ( k = 0，1，2，3) （2）由此得向量 x(k) 和 x( k+1)的递推关系式x( k+1) = lx(k) (3) 其中，矩阵精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - -

6、 - - - - - - - - - 第 2 页，共 12 页 - - - - - - - - -称为莱斯利矩阵。由式 (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*x0 x2=a*x1 x3=a*x2 x4=a*x3 可得数据结果x1 = 7000 500 250 x2 = 2750 3500 125 x3 = 14375 1375 875 x4 = 1.

7、0e+003 * 8.1250 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).1 x=y; 精品学习资料 可选择p d f - - - - - - - - - - - -

8、 - - 第 3 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 12 页 - - - - - - - - -y=a*x; y1=d1*x; k=k+1; end 可知，当 k=291时，有结论成立。3 . 问题解答与进一步思考根据数学模型计算将数据结果填写如下表 1 k ( 现在) k=1 ( 五年后 ) k=2(十年后 ) k=3(十五年后 ) k=4(20 年后) x1 1000 7000 2750 14375 8125 x2 1000 500 3500 1375 7187.5 x

9、3 1000 250 125 875 343.8 从表中数据变化，如果没有其它的原因，估计农场的动物总数量会逐步增加。在验证生物学研究的结论时，当 k=291可以得到如下结论x = 1.0e+054 * 3.3121 1.1040 0.1840 这说明多年以后，动物数量是大得非常惊人的。如果每个五年平均向市场供应动物c=s s st，分析动物数分布向量变化规律可知x(1) = ax (0) cx(2) = ax (1) c x(3) = ax (2) c x(4) = ax (3) c 所以有x(4) = a4x(0) ( a3 + a2 + a + i ) c 考虑二十年后动物不灭绝，应有x

10、(4) 0 即( a3 + a2 + a + i ) c a 4x(0) 由于 c 是常数向量，故可简单求解不等式组，可取c= 152 152 152 t这说明当五年平均向市场供应三个年龄段的动物各152 头可以使 20 年后有各年龄段的动物生存。 如果将这一限制作为约束条件，而求 c 的各分量之和最大， 这将是一个线性规划问题，可用单纯形法求解（或利用数学软件见第七章）。精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 4 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - -

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

12、初始分配数据可得如下向量所以一年后，市场分配为两年后，市场分配为以向量 xn记第 n 年后市场分配的分额，则（n = 1，2，）设有数据 a 和 b 作为 r公司和 s公司的初始市场分额，则有a + b = 1 为了使以后每年的市场分配不变，根据顾客数量转移的规律，有即这是一个齐次方程组问题。 如果方程组有解， 则应该在非零解的集合中选取正数精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 12

13、 页 - - - - - - - - -解作为市场稳定的初始分额。2、程序和计算结果为了知道两年、五年、十年后市场分配的情况，在matlab 中键入下面命令a=1/4 1/3;3/4 2/3 x0=3/5;2/5 x2=a2*x0 x5=a5*x0 x10=a10*x0 可得数据结果x2 = 0.3097 0.6903 x5 = 0.3077 0.6923 x10 = 0.3077 0.6923 由此得下表表 2 r公司的市场分额s公司的市场分额两年后31% 69% 五年后31% 69% 十年后31% 69% 为了求 a 和 b 作为 r公司和 s公司稳定的初始市场分额，需要求解齐次方程组。键

14、入下面命令format rat rref(a-eye(2) 得数据结果ans = 1 -4/9 0 0 由此得化简后的方程a 4/9 b = 0 结合约束条件a + b = 1 得a= 4/1331% b= 9/1369% 这是使市场稳定的两家公司的初始分额，也正好与表4-1 中的数据相吻合。3、问题的解答和进一步思考精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 6 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 6 页，共 12 页 - - - - -

15、 - - - -在 r公司和 s公司的市场初始分额分别为60% 和 40% 的情况下，根据计算结果，两年后情况变化较大： s公司大约占 31% ，r公司大约占 69% 。而五年以后与两年以后比较变化不大： s公司大约占 30.8%，r公司大约占 69.2%。十年后的情况与五年后的情况比较大约不变。市场已经趋于稳定。是否所有市场初始分配分额，在经过若干年后均会趋于稳定状态。（三）常染色体遗传问题假定所考虑的遗传特性由两个基因a和 a 来支配，人类的眼睛染色体是通过常染色体遗传来控制，例如aa及 aa型产生棕色眼睛， aa 型的是兰色眼睛。在常染色体遗传中， 一个个体从它的亲本的每一基因对中遗传一

16、个基因，以形成它自己特殊的基因对： aa ，aa ，aa。亲本的两个基因中的哪一个传给后代纯属机会问题，如果一个亲本是 aa型，后代从这个亲本遗传获得a基因或 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

17、 1/2 0 aa 0 0 0 1/4 1/2 1 例 假定一个农民有一大片作物，它由三种可能基因型aa ，aa及 aa 的某种分布所组成。农民要采用的育种方案是：作物总体中的每种作物都总是用基因型aa的作物来授粉，我们要导出在任何一个后代总体中三种可能基因型的分布表达式。解 记 a n (n = 0，1，2，)为在第 n 代中 aa基因型作物所占的分数， bn为在第 n 代中 aa基因型作物所占的分数，cn 为在第 n 代中 aa 基因型作物所占的分数。 a0，b0，c0 表示基因型的原始分布，且a0 + b0 + c0 = 1 由于用基因型 aa的作物来授粉，分析基因表（前三列数据）可知，

18、从上一代的基因型分布产生的下一代的基因型分布可用下列递推公式求出：精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 7 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 7 页，共 12 页 - - - - - - - - -其中，第一式表明，基因型aa的所有后代都是 aa型基因，基因型 aa的后代，有一半是 aa型。这一递推公式的矩阵表示为x( n) = mx ( n - 1) ( n = 1，2，) 其中，由递推公式可得x( n) = mx ( n - 1

19、) = m 2x( n - 2) = = m n x(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 985/1393 -881/1079 0 0 881/215

20、8 d = 1 0 0 0 1/2 0 0 0 0 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 8 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 8 页，共 12 页 - - - - - - - - -这表明， m的三个特征值为d1 = 1 ，d2 = 1/2 ，d3 = 0 因为特征向量乘一非零数仍是特征向量，所以可取三个特征值对应的特征向量分别为，于是可逆矩阵为了求逆矩阵，使用命令p=1 1 1 ；0 -1 -2；0 0 1 ；inv(p) 可得

21、数据结果ans = 1 1 1 0 -1 -2 0 0 1 所以由前面递推公式，得x( n) = pd np1x(0) 而精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 9 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 9 页，共 12 页 - - - - - - - - -故所以这是原始基因型分数表示第n 代作物总体中三种基因型分数。 显然，当 n时，有an(a0 + b0 + c0)=1，bn0，cn0这说明在极限情况下，总体中所有作物都将是基因型a

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

23、实验性生产线每年一月份进行熟练工人的人数统计；然后将其1/6 的熟练精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 10 页，共 12 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 10 页，共 12 页 - - - - - - - - -工人支援其他生产部门， 缺额由招收非熟练工人补齐。 新、老非熟练工人经过培训及实践至年终考核有2/5 成为熟练工人。设第 n 年一月份统计的熟练工人和非熟练工人所占的百分比分别为xn 和 yn 记为向量 xn ，yn t(1) 试推

24、导向量 xn ，yn t和xn+1，yn+1t 的关系，并写成矩阵形式；(2) 当x1 ，y1 t=0.5 ，0.5 t时，求第 10 年一月份统计的熟练工人和非熟练工人所占的百分比。5 在一城市的某商业区内， 有两家有名的快餐店“肯德基”分店和“麦当劳”分店。据统计每年“肯德基”保有其上一年老顾客的1/3 ，而另外的 2/3 顾客转移到“麦当劳”；每年“麦当劳”保有其上一年的老顾客的1/2 ，而另外的 1/2 顾客转移到“肯德基”。用二维向量 x k = xk yk t 表示两个快餐店市场分配的情况，初始的市场分配为x0 = 1/3 2/3t如果有矩阵l 存在，使得xk+1 = lxk，则称

25、l为状态转移矩阵。(1) 写出xk=xk ykt和xk+1=xk+1 yk+1t的递推关系式，以及状态转移矩阵l。(2) 根据递推关系计算近几年的市场分配情况；(3) 求可逆矩阵 p和对角矩阵 使得 l=pp 1。对于足够大的时间值k，有（是状态转移矩阵 l 的唯一正特征值）。请检验这一结果是否正确，并给出适当的k 的值6 某厂生产 a， b两种品牌的味精， 顾客的喜好决定了这两种味精的市场占有率。在生产中可根据占有率调整比例，获得最佳收益。 该厂做市场调查后发现， 一般情况下，顾客若购买a牌，下次有 80% 的可能性购买 a牌；若购买了 b 牌，下次有 60% 的可能性购买 b牌。开始时，两种品牌的市场占有率分别为50% ，顾客每一次的购买必将改变二者市场占有率。（1）预测某一个顾客经过前四次购买之后，他可能第五次购买哪一个品牌的味精。（2）预测 100 个顾客经过前四次购买之后，两种品牌的可能市场占有率各为多少？7足球比赛排名问题：下表给出了我国 8 支足球队在 19881989 年全国足球甲级队联赛中的成绩，请依据现有的数据，对各足球队的实力给予科学客观的评价；给出 8 支球队按实力排序的名次，并将你的排名算法推广到任意 n 个队的情形。t1 t2 t3 t4 t5 t6 t7 t8 t1 0:1