动物繁殖问题数学建模实验报告.doc
《动物繁殖问题数学建模实验报告.doc》由会员分享,可在线阅读,更多相关《动物繁殖问题数学建模实验报告.doc(8页珍藏版)》请在冰点文库上搜索。
科学技术学院
上机报告
课程名称数学建模
上机项目动物繁殖的规律问题
专业班级姓名学号
一、问题提出
谋农场饲养的某种动物所能达到的最大年龄为15岁,将其分为三个年龄组:
第一组0~5岁;第二组6~10岁;第三组11~15岁。
动物从第二个年龄组开始繁殖后代,第二个年龄组的动物在其年龄段平均繁殖4个后代,第三个年龄组的动物在其年龄段平均繁殖3个后代。
第一年龄组和第二年龄组的动物能顺利进入下一个年龄组的存活率分别为0.5和0.25。
假设农场现有三个年龄段的动物各有1000头。
(1)编程,计算5年后、10年后、15年后、20年后各年龄段动物数量,50年后农场三个年龄段的动物的情况会怎样?
(2)根据有关生物学研究结果,对于足够大的时间值,有≈(是莱斯利矩阵的唯一正特征值)。
请检验这一结果是否正确,如果正确给出适当的的值。
(3)如果每五年平均向市场供应动物数=,在20年后农场动物不至灭绝的前提下,计算应取多少为好?
二、问题分析
在初始时刻0~5岁、6~10岁、11~15岁的三个年龄段动物数量分别为:
=1000, =1000, =1000
以五年为一个年龄段,则某一时刻三个年龄段的动物数量可以用一个向量
=
表示。
以五年为一个时间段,记
=
为第个时段动物数分布向量。
当=0,1,2,3时,分别表示现在、五年后、十年
后、十五年后的动物数分布向量。
根据第二年龄组和第三年龄组动物的繁殖能力,在第个时间段,第二年龄组动物在其年龄段平均繁殖4个,第三年龄组动物在其年龄段平均繁殖3个后代。
由此得第一年龄组在第+1个时间段的数量如下:
=4+3
同理,根据第一年龄组和第二年龄组的存活率,可得等式
=0.5, =0.25
三、模型假设
四、模型建立
建立数学模型如下:
=4+3
=0.5(=0,1,2,3)
=0.25
改写成矩阵形式
(=0,1,2,3)
由此得向量和的递推关系式
=
其中矩阵
=
称为莱斯利矩阵,进一步有
=
五、模型求解
(1)编程,计算5年后、10年后、15年后、20年后各年龄段动物数量,50年后农场三个年龄段的动物的情况会怎样?
程序为:
x0=[1000;1000;1000];
L=[043;1/200;01/40];
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*
8.125000000000000
7.187500000000000
0.343750000000000
x10=
1.0e+005*
1.323769531250000
0.666669921875000
0.068061523437500
(2)根据有关生物学研究结果,对于足够大的时间值,有≈(是莱斯利矩阵的唯一正特征值)。
请检验这一结果是否正确,如果正确给出适当的的值。
>>eig(L)%计算Leslie矩阵的特征值
ans=
1.5000
-1.3090
-0.1910
即矩阵L的唯一正特征值
程序为:
%exam01_17.m
x=[1000;1000;1000];d1=1.5;
L=[043;1/200;01/40];
y=L*x;
y1=d1*x;
k=1;
whilemax(abs(y-y1))>0.1
x=y;
y=L*x;
y1=d1*x;
k=k+1;
end
k
%Ö´ÐÐexam01_17.m
:
k=
285
即当k=285时,有结论成立。
(3)如果每五年平均向市场供应动物数=,在20年后农场动物不至灭绝的前提下,计算应取多少为好?
解:
如果每个五年平均向市场供应动物c=[sss]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由于c是常数向量,故可简单求解不等式组,可取c=[152152152]T
这说明当五年平均向市场供应三个年龄段的动物各152头可以使20年后有各年龄段的动物生存。