1、.淮 阴 工 学 院数学建模课程设计班 级: 计科1091 姓 名: 刘红斌 学 号: 1094101109 选 题: A 组 第 09 题教 师: 王小才胡 平 姜红燕 教 师王小才胡 平姜红燕答辩分总 分得 分数 数 理 院 2011年12月一年生植物的繁殖摘要本文研究生植物的繁殖问题,根据生植物的繁殖规律建立了一个多年后该植物繁殖数量变化情况的三阶线性常系数差分方程模型。实验利用数学软件采用一维搜索的方法,最终确定了1岁种子,2岁种子和3岁种子的比例的取值范围,得到了当时该植物就能一直繁殖下去,而当时就不能繁殖的结果。此模型能够很好地解决类似此类预计某项事物发展规律的问题,具有较强的规律
2、性。关键词:MATLAB,三阶差分方程,一维搜索一 、问题重述1.1背景资料与条件一年生植物春季发芽,夏天开花,秋季产种,不考虑腐烂,被人为掠取。这些种子如果可以活过冬天,其中一部分能在第二年春季发芽,然后开花,产种,其中的另一部分虽未能发芽,但如又能活过一个冬天,则其中一部分可在第三年春季发芽,然后开花,产种,如此继续,一年生植物只能活1年,而近似的认为,种子最多可以活过三个冬天。现在在一片空地上种上500颗某种该植物。记一棵植物春季产种的平均数为,种子能活过一个冬天的(1岁种子)比例为b,活过一个冬天没有发芽又活过一个冬天的(2岁种子)比例仍为b, 活过两个冬天没有发芽又活过一个冬天的(3
3、岁种子)比例仍也为b,1岁种子发芽率,2岁种子发芽率,3岁种子发芽率,为固定值,是变量。1.2需要解决的问题试建立数学模型研究这种植物数量变化的规律,及它能一直繁殖下去的条件。二、问题的假设1. 不考虑恶劣的气候环境影响种子春季的种量;2. 不考虑食物链对该植物的影响;3. 不存在自然灾害的破坏;4. 不考虑外部作用使该植物发生突变或变异。三、符号的说明:第年植物数量:一棵植物春季产种的平均数:空地上初始的植物数量:1岁种子发芽率:2岁种子发芽率:3岁种子发芽率:1岁种子,2岁种子和3岁种子的比例四、问题分析根据所给条件,能够将年之后植物的数量表示出来,但是1岁种子,2岁种子和3岁种子的比例不
4、能确定,所以本题是在其他条件都确定的情况下,比较在的不同取值下,植物数量的变化规律。记第年植物数量为,显然与有关,由决定的部分是,由决定的部分是;由决定的部分是。 记: 实际上,就是 ,我们需要知道, 考察不同时,种子繁殖的情况,而,下面就可以利用软件求解了。五、模型建立本题目提供了该植物初始数量,活过1个冬天的种子的比例,一棵植物春季产种的平均数以及1岁种子的发芽率,因此不难求出1年后该植物的数量,2年后该植物的数量,以此类推。但是种子最多可以活3个冬天,所以年后的植物总数,所以建立模型为:六、模型的求解针对这个问题,我们利用数学软件求解该植物在未来20年的数量变化,并从中找出其变化规律。为
5、了获得比例的合适范围,我利用数学软件采用一维搜索的方法,以步长为0.001确定了。为了获得理想的答案数据及可观美观的关系图像,发现当选择时,所得到的图像更能清晰地反映出该植物的数量变化规律及繁殖条件,更显示了在不同的的作用下,数量变化趋势的不同,这样易于理解与分析(程序见附录一)。植物数量变化规律图(分别取0.13、0.14和0.15):图(1)植物数量变化规律图由运行的结果(见附录二)及植物数量变化规律图得:当时,20年里该植物数量从500下降到135;当时,20年里该植物数量从500增加到591;当时,20年里该植物数量从500增加到2329,并且,运行结果中,因此,该植物的变化规律可以总
6、结为:当时,该植物的总数量会不断地增加;当时,该植物的总数量会保持不变;当时,该植物的总数量会不断减少。 七、模型的评价与推广7.1模型的评价1 本文所建立的模型成功地找到了该一年生植物在未来20里的数量变化规律以及找到了能使之繁殖下去的条件;2 在模型的建立与求解过程中,文中利用图表结合,使文章所要表达的思想简洁明了,更形象、直观;3 在模型中,可以很明确地发现与的关系,通过模型将复杂的高阶线性常系数差分方程问题简单化,更易于理解;4 在模型的求解过程中,利用数学软件来求解,所得结果清晰明了,大大减少了计算量;5 该模型不足之处就是单纯地找出了该植物20年内数量变化规律及繁殖条件,20年以后
7、是否会发生变化没能进行检验。7.2模型的推广本模型的建立为解决高阶线性常系数差分方程提供了一个合理的方案,可以应用于其他类似的问题。由此可见,该模型的应用非常广泛,对于实际问题,具有很强的适应性,例如可以推广到蜘蛛网模型上,汽车租赁公司的运营问题上,动物按年龄分组的种群增长问题上等。参考文献1 胡运权,运筹学习题集(修订版),北京:清华大学出版社,1994年.2 胡运权,运筹学教程(第三版),北京:清华大学出版社,2006年.3 白峰杉,数值计算引论(第二版),北京:高等教育出版社,2010年.4 王沫然,与科学计算,北京:电子工业出版社,2003年.附录附录一:程序代码:fun函数:func
8、tion x=zwfz(x0,n,b)c=10;a1=0.7;a2=0.4;a3=0.2;k1=a1*b*c;k2=a2*b*(1-a1)*b*c;k3=a3*b*(1-a1)*(1-a2)*b*c;x(1)=x0;x(2)=k1*x(1);x(3)=k1*x(2)+k2*x(1);for i=4:n x(i)=k1*x(i-1)+k2*x(i-2)+k3*x(i-3);end主函数:clear allclcS=500; %初始野生植物数量N=800; %模拟次数b=0; flag=1;while flag x=zwfz(S,N,b); if x(end)1 flag=0; end b=b+0
9、.001;endbk=(0:20);y1=zwfz(500,21,0.13);y2=zwfz(500,21,0.14);y3=zwfz(500,21,0.15);round(k,y1,y2,y3)plot(k,y1,k,y2,r-,k,y3,g-),gtext(b=0.13),gtext(b=0.14),gtext(b=0.15)xlabel(时间);ylabel(一年生植物数量);title(一年生植物的繁殖);附录二:程序结果:b = 0.1390ans = 0 500 500 500 1 455 490 525 2 424 492 565 3 398 497 611 4 374 502 661 5 351 507 715 6 329 513 774 7 309 518 837 8 290 523 906 9 272 528 980 10 255 534 1060 11 240 539 1147 12 225 545 1241 13 211 550 1343 14 198 556 1453 15 186 562 1572 16 175 567 1700 17 164 573 1839 18 154 579 1990 19 144 585 2153 20 135 591 2329