动物繁殖问题数学建模实验m精编b程序修订稿

动物繁殖问题数学建模实验m精编b程序

内部编号：（YUUT-TBBY-MMUT-URRUY-UOOY-DBUYI-0128）

问题：

谋农场饲养的某种动物所能达到的最大年龄为15岁，将其分为三个年龄组：第一组0～5岁；第二组6～10岁；第三组11～15岁。动物从第二个年龄组开始繁殖后代，第二个年龄组的动物在其年龄段平均繁殖4个后代，第三个年龄组的动物在其年龄段平均繁殖3个后代。第一年龄组和第二年龄组的动物能顺利进入下一个年龄组的存活率分别为0.5和0.25。假设农场现有三个年龄段的动物各有1000头。

（1）编程，计算5年后、10年后、15年后、20年后各年龄段动物数量，50年后农场三个年龄段的动物的情况会怎样？

x0=[1000; 1000; 1000];

L=[0 4 3; 1/2 0 0; 0 1/4 0];

x1=L*x0 %计算5年后农场中三个年龄段动物的数量

x2=L*x1 %计算10年后农场中三个年龄段动物的数量

x3=L*x2 %计算15年后农场中三个年龄段动物的数量

x4=L*x3 %计算20年后农场中三个年龄段动物的数量

x1 =

7000

500

250

x2 =

2750

3500

125

x3 =

14375

1375

875

x4 =

1.0e+003 *

（2）根据有关生物学研究结果，对于足够大的时间值k ，有)1(+k X ≈1λ)(k X （1λ是莱斯利矩阵L 的唯一正特征值）。请检验这一结果是否正确，如果正确给出适当的k 的值。

>> eig(L) %计算Leslie 矩阵的特征值

ans =

1.5000

-1.3090

-0.1910

即矩阵L 的唯一正特征值 1.5λ=

%exam01_17.m

x=[1000; 1000; 1000]; d1=1.5;

L=[0 4 3; 1/2 0 0; 0 1/4 0];

y=L*x;

y1=d1*x;

k=1;

while max(abs(y-y1))>0.1

x=y;

y=L*x;

y1=d1*x;

k=k+1;

end

k

%′DDexam01_17.m

k =

285

（3）如果每五年平均向市场供应动物数c＝[]T s

s，在20年后农场动物

s

不至灭绝的前提下，计算c应取多少为好？

如果每个五年平均向市场供应动物c=[s s s]T，分析动物数分布向量变化规律可知

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年后有各年龄段的动物生存。

提示：现在给大家作出如下问题分析：

在初始时刻0～5岁、6～10岁、11～15岁的三个年龄段动物数量分别为：

)0(1x ＝1000， )0(2x ＝1000， )0(3x ＝1000

以五年为一个年龄段，则某一时刻三个年龄段的动物数量可以用一个向量 X ＝T x x x ][321

表示。以五年为一个时间段，记

)(k X ＝T k k k x x x ][)(3)(2)(1 为第k 个时段动物数分布向量。当k ＝0，1，2，3时，)(k X 分别表示现在、五年后、十年

后、十五年后的动物数分布向量。根据第二年龄组和第三年龄组动物的繁殖能力，在第k 个时间段，第二年龄组动物在其年龄段平均繁殖4个，第三年龄组动物在其年龄段平均繁殖3个后代。由此得第一年龄组在第k ＋1个时间段的数量如下：

)1(1+k x ＝4)(2k x ＋3)(3k x

同理，根据第一年龄组和第二年龄组的存活率，可得等式

)1(2+k x ＝0.5)(1k x ， )1(3+k x ＝0.25)(2k x

建立数学模型如下：

)1(1+k x ＝4)(2k x ＋3)(3k x

)1(2+k x ＝0.5)(1k x （k ＝0，1，2，3）

)1(3+k x ＝0.25)(2k x

改写成矩阵形式

⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣

⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++)(3)(2)(1)1(3)1(2)1(1025.00005.0340k k k k k k x x x x x x （k ＝0，1，2，3） 由此得向量)(k X 和)1(+k X 的递推关系式 )1(+k X ＝L )(k X 其中矩阵

L ＝⎥⎥⎥⎦

⎤⎢⎢⎢⎣⎡025.00005.0340 称为莱斯利矩阵，进一步有 )1(+k X ＝1+k L (0)X