- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学建模兔生长问题完整论文
一、问题重述和分析
(二)野兔生长问题
在某地区野兔的数量在连续十年的统计数量(单位十万)如下:
T=0
T=1
T=2
T=3
T=4
T=5
T=6
T=7
T=8
T=9
1
2.31969
4.50853
6.90568
6.00512
5.56495
5.32807
7.56101
8.9392
9.5817
分析该数据,得出野兔的生长规律。 并指出在哪些年内野兔的增长有异常现象,预测T=10 时野兔的数量。
根据数据中野兔生长数量增长规律, 对于生物增长模型, 我们可以考虑到logistic模型,因为此种模型曲线是单调递增的,但是表格中明显不是单调的,于是可以
分三段讨论,由统计数据可以客观得到如下结果:
T=0、1、2、3时种群数量单调上升,对于生物增长模型可考虑到logistic模型
T=3、4、5、6时种群数量单调递减,是一种反常现象,仍可考虑logistic模型
T=6、7、8、9时种群数量单调上升,对于生物增长模型可考虑到logistic模型
野兔在自然条件不变下,野兔的种群应该保持不变。然而通过读数据的观察发现。野兔的数量并不是单一地增长,T=3,6.90568;T=4,6.00512;T=5,5.56495;T=6,5.32807。第三年到第六年野兔的增长有异常现象,这四年野兔的数量不增反降,说明其间有影响野兔生长的因素存在。
我们探讨了其中的因素:
1、兔子的内部矛盾,兔子之间因为食物的减少而引发争斗
2、天敌大量地捕食使野兔生存受到威胁
3、疾病的侵扰,在野兔种群中蔓延并流行疾病
4、人类的捕杀与破坏。
二、模型的假设
上述野兔生长问题,我们作出以下假设:
1.假设各个环境因素对野兔生长的影响是互不影响的
2.假设兔子没有受到传染性疾病的影响
3.假设它使处于自然的情况(没有人的作用),人类活动对其生存不产生影响
4.假设野兔性别比接近1:1,且采用措施维持这个比列
三、符号说明
连续三年中第一年兔子的数量
连续三年中第二年兔子的数量
连续三年中第二年兔子的数量
: 表示兔子的数量
: 表示兔子的出生率
: 表示兔子的死亡率
: 表示年份
模型的分析与建立
对于生物模型,首先考虑的是logistic模型,考虑到logistic模型的增长曲线是单调的,而题目所给的数据中有一段是下降的,这是反常的情况,而正常情况应当是单调上升的。考虑到可能在这段时间内有使野兔减少的因素,因此不能在整个时间段内进行拟合,所以我们应当在每个单调区间上进行拟合。
第一单调增区间
T=0
T=1
T=2
T=3
1
2.31969
4.50853
6.90568
第二单调减区间
T=3
T=4
T=5
T=6
6.90568
6.00512
5.56495
5.32807
第三单调增区间
T=6
T=7
T=8
T=9
5.32807
7.56101
8.9392
9.5817
我们把野兔生长情况分成了上表三个区间,建立野兔生长的logistic模型。
模型的求解
= 1 \* GB3 ①.首先求解
对于logistic连续模型,设微分方程为
,,
其中参数a,b 需要通过拟合得到。
将整理求得:
. (2)
= 2 \* GB3 ②.再求出变量a和变量b
设已知连续三年的数据,其中
则由(2)得方程组
(3)
这三个方程中有三个未知量,可以解出a,b如下: 将(3)中第1式代入第2、3式消去x0, 得
(4)
消去a 后得b 满足的方程
(5)
解得:
. (6)
代入(4) 的第1式得a 满足的方程
(3)
求参数a,b的MATLAB程序:
Function [a,b, c]=hare(c,t)%
输入单调的连续三年数量c和时间间隔T(本题T=1), 输出参数 a, b和下一年的数量d
a=log(c(3)*(c(2)-c(1))/(c(1)*(c(3)-c(2))));
b=(c(2)^2-c(3)*c(1))/(c(3)*c(2)+c(1)*c(2)-2*c(1)*c(3))/c(2)
d=1/(b+(1/c(3)-b))*exp(-a*T));
在第一个上升阶段, 对于连续三年时间分别为(0,1,2)和(1,2,3)分别计算得到两组 a,b值:
0.9999957823 0.0999975822
1.0000035509
文档评论(0)