数学建模野兔生长问题完整论文

PAGEPAGE13数学建模一周论文论文题目：野兔生长问题2011

摘要假设野兔生长的条件是在无外界干扰的完美条件下（即不考虑外界因素对野兔繁殖的影响），该种群的成长曲线应该为对数型增长。但依题意可知，野兔增长先是成对数增长后来趋于平缓，变化幅度不断降低，这说明野兔生长并不是处于理想的情况下的，考虑到自然的各种原因，诸如，环境条件因为兔群激增而变得恶劣，天气的变化，天敌的增多等等。对于这种种群生态学问题，我们可以用Logistic（逻辑斯蒂方程）模型来模拟。Logistic模型是种群生态学的核心理论之一，它可以很好的表示生物种群的生长规律，动态的表示生物种群的增减情况，例如兔子。由于野兔生长问题相对简单，其涉及的内容和有求也相对较少，并且该问题概过了种群在生态中生长问题。根据逻辑斯蒂方程，以及建立一只双曲线右支可以预测出在T=10时，野兔数量为10.8156十万只。在此，我们结合过去九年野兔数量的历史数据，建立了逻辑斯谛增长模型，得到野兔的生长规律如下：野兔初始于该地方生存时，野兔的生长繁殖有充分的保障，数量增多。随着野兔的不断繁殖，其有限生存空间日趋减小，其数量趋向于某一极值。而当野兔数量超过环境容纳量时,野兔种群的增长受到抑制，数量下降。当野兔种群数量降低到环境容纳量以下时,野兔种群的出生率上升，死亡率下降，自然资源与食物资源较为充裕，种内与种间竞争有所缓解，从而野兔种群增长加快。通过建立Logistic模型，我们小组得出当T=10时，野兔数量为10.8156（十万）只左右。该结果比较符合客观规律。利用Logistic模型可以表征种群的数量动态；如昆虫类种群的增长，收获与时间关系的确定。描述某一研究对象的增长过程如生态旅游区环境容量的确定，森林资源的管理以及耐用消费品社会拥有量的预测、国民生产总值的预测等；也可作为其它复杂模型的理论基础如Lotka-Volterra两种群竞争模型；以上的大多数的工作都是拿逻辑斯蒂模型来用，但也由此可看出逻辑斯蒂方程不管在自然科学领域还是在社会科学中都具有非常广泛的用途。关键字:分析数据异常现象预测数量逻辑斯蒂方程模型问题重述野兔生长问题：在某地区野兔的数量在连续十年的统计数量(单位十万)如下：T=0T=1T=2T=3T=4T=5T=6T=7T=8T=912.319694.508536.905686.005125.564955.328077.561018.93929.5817分析该数据，得出野兔的生长规律。并指出在哪些年内野兔的增长有异常现象，预测T=10时野兔的数量。首先，野兔是生长在自然环境中的。自然很复杂，存在着许多影响种群发展的因素。我们知道，假如给野兔一个理想的环境，野兔数量是呈对数增长的。现实情况中，种群一般是呈S型增长的，从题中表格看出，野兔的数量并不是单一地增长，T=1,2.31969；T=3，6.90568；T=4，6.00512；T=5，5.56495，呈类J型增长，说明兔子数量不多受内外因素的因数影响不明显。第四年到第七年，这三年野兔的数量不增反降呈类S型增长，说明其间有影响野兔生长的因素存在。我们探讨了其中的因素：（1），兔子的内部矛盾，兔子之间因为食物的减少而引发争斗等。（2），自然环境的恶化，比如说兔子的激增使粪便数量大大增加是环境变得恶劣，不在适合兔子的生存；再如气候反常，使野兔的产卵，交配受影响。（3），天敌的捕食，狼，狐狸等天敌大量地捕食使野兔生存受到威胁。（4），疾病的侵扰，野兔种群中，蔓延并流行疾病，必然使野兔存活率下降。。（5），人类的影响，城市扩建，使其栖息地面积减少；捕杀。考虑到上述因素，野兔的生长就不能完全用一个Logistic模型来模拟模型假设因为所学知识有限，所以我们做出以下假设以方便猜想：（1）,假设它使处于自然的情况（没有人的作用）,人类活动对其生存不产生影响。（2），假设各个环境因素对野兔生长的影响是互不影响的。兔子种族内部生存空间足够多，不存在对生存空间需求问题。（3），假设兔子的内部因素对其生存率的影响不大。（4），假设野兔在各年龄段中的分布率不变，即年龄结构不变，并采用各种措施维持这一结构；（5）假设野兔性别比接近1：1，且采用措施维持这个比列。（6）假设母兔可以怀孕的年龄为1到6岁，最高年龄为10岁，10岁的死亡率为100%，并且6到10岁的野兔个数成线性递减趋势。在以上条件成立的前提下，用Logistic模型来模拟野兔的增长情况。建立模型对于生物模型，首先考虑的是logistic模型，考虑到logistic模型的增长曲线是单调的，而题目所给的数据中有一段是下降的，这是反常的情况，而正常情况应当是单调上升的。考虑到可能在这段时间内有使野兔减少的因素。不能在整个时间段进行拟合，我们应当在每个单调区间上进行拟合。第一单调增区间T=0T=1T=2T=312.319694.508536.90568第一单调减区间T=3T=4T=5T=66.905686.005125.564955.32807第二单调增区间T=6T=7T=8T=95.328077.561018.93929.5817我们把野兔生长情况分成了上表三个区间，建立野兔生长的logistic模型。我们先假设在一个小的单位时间间隔内新出生的兔子百分比为b，类似的兔子死亡率的百分比为c。换句话，新的兔子数P(t+Δt)是原有兔子数P(t)加上在Δt时间内新增兔子数减去死亡兔子数，即P(t+Δt)=P(t)+bP(t)Δt-cP(t)Δt这样我们把问题化归到如何确定k。一旦k被确定，通过已知数据，我们解这个微分方程，就可以得到一个野兔数量随时间变化的函数了rM(ttrM(tt)MpperM(tt)pMeP(t)这个模型就成为logistic模型。模型求解对于logistic连续模型，设微分方程为，,（1）其中参数a，b需要通过拟合得到。的解为.（2）设已知连续三年的数据，其中，则由（2）得方程组(3)这三个方程中有三个未知量可以解出a,b如下:将（3）中第一式代入第二、三式消去x0，得(4)消去a后得b满足的方程(5)解得.(6)代入(4)的第一式得a满足的方程(3)求参数a,b的MATLAB程序function[a,b,q]=hare(p,T)%输入单调的连续三年数量p和时间间隔T(本题T=1),输出参数a,b和下一年的数量qa=log(p(3)*(p(2)-p(1))/(p(1)*(p(3)-p(2))));b=(p(2)^2-p(3)*p(1))/(p(3)*p(2)+p(1)*p(2)-2*p(1)*p(3))/p(2);q=1/(b+(1/p(3)-b))*exp(-a*T)）;在第一个上升阶段,对于连续三年（0，1，2）和（1，2，3）分别计算得到二组a，b值0.999996295432800.099998990654181.000001896730560.10000006995945在下降阶段，对于连续三年（3，4，5）和（4，5，6）分别计算得到的二组a，b值0.499999514703010.200000053216010.499983964746560.20000085565547在第二个上升阶段，对于连续三年（6，7，8）和（7，8，9）分别计算得到的二组a，b值1.000005087174110.100000057968451.000009756401800.10000014562299当取a,b为最后一组数据时，T=10时由(2)得到预测数为10.1（十万）,当取a=1，b=0.1时，预测数为10.1（十万）.（M为兔子饱和值）结论是：在T=0到T=3之间增长规律为logistic模型：.在T=3到T=6之间增长规律有异常情况,但仍为logistic模型；：.在T=6到T=9之间增长规律恢复为logistic模型：.在T=10时,在正常情况下,野兔数量为10.8156（十万）只.模型的误差分析这里我们需要讨论几个造成模型误差的问题，它们都很有实际意义。一个是模型对初值的依赖性问题，另一个是模型对参数选取的依赖性问题。模型的优缺点模型的评价与推广模型评价一开始我们根据生物学的知识，对种群数量的饱和值取了一个估计值，而事实上，我们在估计了M=10左右的几个值后，又分析了不同的M值对模型解的影响，从而得出我们对M值的估计具有一定的现实性和真实性。在考虑初值误差对模型解影响过程中，我们采用常微分知识，考虑解对初值的依赖性和解对参数的依赖性，这种思路对于实际研究具有极其重大的现实意义。该模型的使用范围是：在一定的空间区域内，不考虑环境的重大变化对生物种群的影响，而是按照该区域的一般自然规律。同样，这也是本模型的不足之处，因为现实中，我们研究的生物种群可能对环境的依赖性很强，那么环境的影响就是一个重要的因素，但是我们也意识到要用一个很具体的模型，该模型考虑现实中的每一个因素，并且表示这些因素的表达式又是简单的数学式，要找这样的模型是很困难的，而我们采用的模型在绝大部分同类问题中所体现的优越性是易见的，并且我们在以上的分析中也给出了相对较好的模型。模型的推广这个模型是针对一定空间区域内，在一般规律下，对生物种群数量随时间的变化范围作出的估计。这种估计对于研究某一区域内的种群有直接的好处，也利于人类研究人类自身。对于一定区域内的生物种群，如果考虑种群生存环境的一般现象，我们可以从事实上看到我们应用的模型对于该类问题的解决有重大的预见效果。现在，在关于生物种群数量的研究中，有很多种方法，而logistic模型对于这类与现实搭配的事实更趋于成功。由以上的分析，我们看到了该模型对于相对较小误差的测量，是可以得到相当准确的预测值的。于是，我们认为，在环境不发生重大改变的情况下，研究生物种群数量与时间的关系就可以采用我们在解决本问题时的解法参考资料：[1]FrankR.Giordano,MauriceD.Weir,WilliamP.Fox,AFirstCourseinMathematicalModeling,(ThirdEdition),北京：机械工业出版社，2005.1[2]郑阿奇等，MATLAB实用教程，北京：电子工业出版社，2004.5[3]洪伟，吴承祯等，对种群增长模型的改进，应用与环境生物学报，2004.10目录TOC\o"1-2"\h\z\u第一章总论 11.1项目名称与承办单位 11.2研究工作的依据、内容及范围 11.3编制原则 31.4项目概况 31.5技术经济指标 51.6结论 6第二章项目背景及建设必要性 82.1项目背景 82.2建设的必要性 9第三章建设条件 113.1项目区概况 113.2建设地点选择 错误！未定义书签。3.3项目建设条件优劣势分析 错误！未定义书签。第四章市场分析与销售方案 134.1市场分析 134.2营销策略、方案、模式 14第五章建设方案 155.1建设规模和产品方案 155.2建设规划和布局 155.3运输 185.4建设标准 185.5公用工程 205.6工艺技术方案 215.7设备方案 215.8节能减排措施 24第六章环境影响评价 256.1环境影响 256.2环境保护与治理措施 266.3评价与审批 28第七章项目组织与管理 297.1组织机构与职能划分 297.2劳动定员 297.3经营管理措施 307.4技术培训 30第八章劳动、安全、卫生与消防 318.1编制依据及采用的标准 318.2安全卫生防护原则 318.3自然灾害危害因素分析及防范措施 328.4生产过程中产生的危害因素分析及防范措施 328.5消防编制依据及采用的标准 348.6消防设计原则 358.7火灾隐患分析 358.8总平面消防设计 358.9消防给水设计 368.10建筑防火 368.11火灾检测报警系统 378.12预期效果 37第九章项目实施进度 389.1实施进度计划 389.2项目实施建议 38第十章项目招投标方案 4010.1招标原则 4010.2项目招标范围 4010.3投标、开标、评标和中标程序 4010.4评标委员会的人员组成和资格要求 42第十一章投资估算和资金筹措 4311.1投资估算