大学学建模.docx

大学学建模.docx

《大学学建模.docx》由会员分享，可在线阅读，更多相关《大学学建模.docx（44页珍藏版）》请在冰豆网上搜索。

大学学建模大学学建模动物群落的稳定发展摘要：

本文通过对某公园近两年内被运出的某种动物的年龄和性别的数据进行统计分析，并针对题目的四个问题分别建立了符合实际的数学模型，在模型的求解过程中，应用C语言进行编程调试，通过统计学软件SAS，数学软件MATLAB等计算工具，编写相应的程序，对建立的模型进行求解，得出了符合实际的结果。

问题一：

我们假设新生幼仔的数量为，然后通过对各年龄阶段的存活率、被运走的动物数量以及该动物的总体数量的分析来建立该群落的动态变化模型，利用该群落近两年内被运走的各年龄阶段的个体数量分布，用C语言编程计算，推测出当前该动物的年龄结构（具体结果见7页表一）。

并利用MATLAB软件对得出的数据用图形表示，利用对比分析法，得到该动物群落的基本分布轨迹，最后用统计软件SAS对模型进行相关性的分析检验，求得相关系数R与P的值，验正了模型的稳定性。

问题二：

由于现在采用注射避孕药的方法来维持该种群的稳定，而且已经没有个体被运走或被偷猎的情况，为此我们把该种群的稳定性转化为求目标函数（该种群每年的新生幼仔的数量减去该年死亡个体的数量的差值）；另外从（即年头的数量与该年年底的数量的差值）当趋于0时，即认为该群落的个体数量是稳定的，从而把问题的稳定性问题转化为求单目标的最优化问题建立模型；利用MATLAB对模型进行求得，得出当不考虑不确定性因素影响时要注射药物的雌性动物数量为276头，而当考虑了双胞胎和被重复注射这两个不确定性因素影响后，得到要注射药物的雌性动物数量为352头，其中有110头是被重复注射的。

问题三：

其大致模型与问题二相近，不同之处在于要考虑到被运走的动物的数量（b），即目标函数应考虑上被运走的数量，即只是对问题二的模型进行扩充建立新的目标模型；b和b；利用MATLAB对不同b值进行求解，从而得出相应的避孕措施。

（具体结果见19页表二）问题四：

我们引进了增量加速度的概念，利用c语言进行编程求解，然后用MATLAB软件对得到的数据进行线性回归分析，得到该群落在减少至M时重新壮大该动物群落能力的模型：

M3.9010+0.0047D。

最后应用统计软件SAS对模型进行稳定性分析。

关键字：

存活率年龄结构新生幼仔数稳定性最优目标增量加速度一问题重述与提出位于非洲某国的国家公园中栖息着近11000头某种野生动物。

管理员要求有一个健康稳定的环境以便维持这个11000头该动物的稳定群落。

过去的20年中，整个该动物群是通过一些偷猎枪杀以及转移到外地而稳定下来的。

但是近年来，偷猎被禁止，而且每年要转移这些动物也比较困难，因此，要控制现在的数量就使用了一种避孕注射法。

用这种方法注射一次可以使得一头成熟雌性动物在两年内不会受孕。

要探讨这种避孕注射法的实用性，我们需要完成以下问题：

1探讨该动物年龄在2岁到60岁之间的合理的存活率的模型，推测这个动物群落的当前的年龄结构。

2估计每年在该群落中有多少雌性动物要注射避孕药，可以式群落固定在11000头左右。

这里不免有些不确定性，也要估计这种不确定性的影响。

3假如每年转移50至300头此动物到别处，那么上面的避孕措施将可以有怎样的改变？

4如果由于某种原因，突然使得注射避孕的方法不得不停止（例如由于一场灾难导致大量该动物的死亡），那时重新壮大该动物群的能力如何？

二基本假设与符号说明

（一）模型假设1该公园是非开放式的，它与外界不发生关系，从而构成独立的生物群落，该动物群落不存在与其它动物种群的竞争，或虽有竞争，但其影响只局限于该动物群落的死亡率内。

2种群是通过雌性个体的繁殖而增长的，所以用雌性个体数量的变化为主要研究对象。

3为了讨论的必要，我们把新生的幼儿的存活率定为75，而其后的存活率为95，直到60岁为止。

各年龄组的该动物经过一年后即进入高一级的年龄组，而龄超过60即认为全部死亡，退出该系统。

4由于该公园加强了对该动物群落的保护，我们认为该动物没有再被偷猎射杀。

而该动物群落个体数量的减少只是因为自然死亡以及被运走。

5假设同一年龄组的动物个体之间是同质的，我们只考虑其平均水平，不讨论个别差异。

6题设该动物在1012岁开始怀孕，我们这里设定为11岁开始，经过22个月（约两年）的怀孕期后生幼仔，即可认为该雌性动物在1360岁的时间内可以生幼仔。

7该群落的自然死亡是在生完幼仔后才发生的，产幼仔只发生在每年的年初时段，而被运走只发生在年底时段。

（二）符号说明：

新生幼儿的存活率，其值为0.75；：

160岁个体的存活率，其值为0.95；：

双胞胎出生的几率，其值为0.0135；：

该动物第k年时刻的数量；：

该动物第k年初i龄动物的数量；：

该动物第k年初底i龄动物的数量：

第j年被运走的动物的数量；：

表示该动物第k年初时的总数量；:

表示每年没有注射避孕药的雌性动物生幼仔的几率，其值为；：

表示被注射过避孕药但在两年内不再被注射的雌性动物生幼仔的几率，其值为；：

表示被注射过避孕药但在两年内被重复注射的雌性动物生幼仔的几率，其值为；：

表示从1360岁该动物的雌性个体的总数；：

表示从159岁该动物的个数总和；：

表示60岁该动物的个体总和；：

表示1360岁雌性动物没有被注射避孕药部分的数量；：

表示1360岁雌性动物被注射过避孕药但在两年内不再被注射部分的数量；：

表示1360岁雌性动物被注射过避孕药但在两年内被重复注射部分的数量；：

表示每年出生幼仔的数量与该年个体死亡的数量的差值；：

表示该种群每年的新生幼仔的数量减去该年死亡个体的数量与运走个体数量的和的差值；：

表示该动物群落在年底时的总数量与年初的数量加上被运走的个体数量b的差值。

三问题分析与模型建立问题一：

1我们要研究该动物群落的稳定性问题，首先要根据存活率确定其当前的年龄结构。

该动物的新生幼仔存活率较低，题设是70到80之间，为了讨论的需要，我们这里设定为75。

在1岁后的存活率比较高，在这里设为95，直到60岁，而超过60岁则认为退出该系统。

因此，我们先建立出该动物群落中年龄在2岁到60岁之间的合理的存活率的模型。

模型一：

式

（1）表示该动物第k年增长的数量；式

（2）表示该动物第k年初时的总数量，可由已有的数据计算出来；式（3）表示该动物被运走的数量；式（4）和（5）表示该动物第i龄到了年底全部转化为（i1）龄；式（6）和（7）表示该动物各年龄段的变化；式（8）表示该动物新生的幼仔数量。

2通过对该公园近两年内从这个地区运出的该动物的年龄和性别的数据进行统计分析，并利用编程工具TurboC2.0对该模型进行编程计算（源程序及计算过程见附录1），可得到当前该动物群落的年龄结构，如下表所示：

表一该动物的年龄结构统计表年龄（岁）前一年数量（头）前一年运走数量（头）前一年剩下数量（头）前兩年数量（头）前兩年运走数量（头）前两年剩下数量（头）假设无运走数量（头）今年数量（头）079278680680674680015940594604060455960025640564573205535315693535053554421523504540450835055161350347851254824478490124784544866457745046513452431461743420414441224194094378412940341814404388415939115376397403573683941037193623771436334937411352223303582633233135512334333134013327314337133172329432214308298320143015296305272782833031528513272289328626828716270212492741425025427217256025626012248241258182432222124620226228245192301421623325208216232202185213221172042052202120713194209141951942082219610186198101881841972318601861880188174187241761316317821761651772516730137169316615616826158141441604156148159271501213815141471401512814201421433140132143291342011413521331251353012761211283125118128311203117121131081121213211351081141698106114331078991081395100108341011289102109294102359510859610868996369038791127984913785778861670798638801466811269758139751065761066717640711655721260677241672146681949636842631350641351596443591049602436566044561244561739535645536475316375053465034750252547504747641471235444748449354445-1414449411328412318384150381028383443638513633336132334365234628341618323453322111321022303254301515301713283055284242813152628562613132613132426572410142412122224582232-1022319202259201462022-2182060180181820-21718生幼仔的雌性数量7847927347891360岁雌性数量27352322503277230224702569总数量117146221109211876876110001100611808注1：

0岁表示新生幼仔。

注2：

由于每个年龄段的数据均为推测值，而实际上运走的各年龄段的数量不一定全部与预测值相符，故表中“剩下数量”两组数据中出现负数可认为是独异点，不影响模型整体的准确性。

每年新生幼仔的数量（）减去生幼仔的雌性的数量（，由于雄性与雌性的数量比接近1：

1，我们可近似地认为1360岁个体的雌雄数量相等），其差值即为双胞胎的数量，这个差值与生幼仔的雌性数量之比即为双胞胎的几率（）1.35。

由表中数据可得，这些比例都基本上接近题设的双胞胎的几率0.0135，说明以上推测得出的数据是准确的。

利用Matlab软件对以上四组数据用图形表示，并进行比较，得到该动物群落的基本分布情况图（源程序见附录2），如下图所示图1分析该图，可以看出，这四组曲线的轨迹、分布情况基本相同。

由于“预测当前的年龄结构情况（无运走）”一组数据没有减去被运走的个体数量，故其每个年龄层的数量都略多于前三组的数量，因此其曲线比前三组的曲线略高一点，利用SAS软件对模型进行相关性的分析检验（源程序见附录3），得到如下结果：

图2程序的分析及统计结论：

程序中的x1是前一年的该动物群落的年龄结构，x2是前两年该动物群落的年龄结构，x3是该动物群落没有被运出是的年龄结构，x4是预测的当前的该动物群落的年龄结构。

过程中的PROCCORR是分析变量中两两变量之间的PEAROS简单相关的。

输出结果中的结果1是一些基本的描述统计量，结果2是两两变量之间的相关矩阵，其中包括相关系数和显著性检验的概率。

由结果可知前一年的该动物群落的年龄结构（x1）与前两年的该动物群落的年龄结构（x2）的相关系数R=0.99744,P=0.00010.01，所以前一年的该动物群落的年龄结构（x1）与前两年的该动物群落的年龄结构（x2）之间存在着极显著的正相关；前一年的该动物群落的年龄结构（x1）与没有运走是的该动物群落的年龄结构（x3）的相关系数R=0.99896,P=0.00010.01，所以前一年的该动物群落的年龄结构（x1）与没有运走时的该动物群落的年龄结构（x3）之间存在着极显著的正相关；同理可知x1与x4的相关系数R=0.99896,P=0.00010.01;x2与x3的相关系数R=0.99846,P=0.00010.01；x2与x4的相关系数R=0.99854,P=0.00010.01；x3与x4的相关系数R=0.99999,P=0.00010.01。

由以上的分析可知，x1,x2,x3,x4之间的相关系数接近1，可见模型一的稳定行很强，而且由公式推出的前一两年的数据与该公园已有的数据基本相符合，可见模型是很优的。

问题二：

由于目前该动物已经很少被移出或移入，而且偷猎枪杀的情况微乎其微，所以暂时不列入考虑范围内。

因此对该动物群落若不采用人工手段控制，则其在一定时间范围内会大幅度增加，从而破坏该种群的动态平衡。

为了保持该种群的平衡，而又不必每年运走一定数量动物，现在使用一种避孕注射法，可使该动物群落的数量固定在一定范围内，用这种方法注射一次可以使得一头成熟雌性动物在两年内不会受孕，但不会引起其它附加的反应。

我们所要做的就是估计出每年在该群落中要注射避孕药的雌性动物的数量，并且要考虑到各种不确定性因素的影响。

为了分析的方便，我们先建立初步模型，该模型暂时不考虑注射避孕药所产生的不确定性因素的影响，即不考虑两年内被重复注射的雌性数量及双胞胎的几率。

在这里我们只认为新生幼仔的数量由两部分组成，一部分为没注射过避孕药的雌性个体所生，另一部分为被注射过避孕药的雌性个体所生。

另外，由于已经没有个体被运走或被偷猎的情况，为此我们把该种群的稳定性转化为求目标函数（该种群每年的新生幼仔的数量减去该年死亡个体的数量的差值），当趋于0时，即认为该群落的个体数量是稳定的，从而把问题的稳定性问题转化为求单目标的最优化问题。

从而建立模型如下:

模型二：

其中：

，表示每年没有注射避孕药的雌性动物生幼仔的几率。

，表示被注射过避孕药但在两年内不再被注射的雌性动物生幼仔的几率；这里设2596，表示从1360岁该动物的雌性个体的总数为2596；设10243，表示从159岁该动物的个数总和为10243；设17，表示60岁该动物的个体总和为17；：

表示该种群每年的新生幼仔的数量与该年死亡个体的数量总和的差值。

把已知的数据代入上述模型，从而得到以下模型：

利用Matlab软件进行编程（源程序见附录4），求得.其中，说明该模型是稳定的，即该动物群落的数量被控制在一定的范围内。

由可知，每年大约有276头雌性动物要注射避孕药，才能使该群落的数量保持在11000头左右。

模型三（模型二的改进）：

由于每年被注射的雌性动物数量一定，所以被注射过后两年有可能又被注射，则将其归入到“被注射过避孕药但在两年内不再被注射的雌性动物”内。

因此，注射避孕药所产生的不确定性因素之一即为“被注射过避孕药但在两年避孕期内被重复注射的雌性动物”。

另外不确定因素之二即双胞胎的几率问题，所以这里我们加入一个参数。

当第一年采用注射避孕药的方法时，是不会发生有雌性个体被重复注射的情况的，故有以下模型：

其中：

，表示每年没有注射避孕药的雌性动物生幼仔的几率；，表示被注射过避孕药但在两年内不再被注射的雌性动物生幼仔的几率；，表示被注射过避孕药但在两年内被重复注射的雌性动物生幼仔的几率；，表示新生幼仔的存活率；0.95,表示160岁动物的存活率；,表示新出生的幼仔中双胞胎的概率。

利用Matlab软件求解（源程序见附录5），求得其中，说明该模型是稳定的，即该动物群落的数量被控制在一定的范围内。

由可知，每年大约有367头雌性动物要注射避孕药，才能使该群落的数量保持在11000头左右。

由于考虑双胞胎的机率，故必须增加注射的数量。

当注射避孕药一年后再次注射时，就会有某些数量的雌性个体被重复注射的情况出现，但这部分一定比前一年注射的雌性个体的数量少，故建立以下模型模型四：

代入已知数据，得到以下目标函数模型：

其中是前一年被注射避孕药的雌性个体的数量。

利用Matlab软件编写程序求解，通过对的不同取值进行调试，（源程序见附录6）求得当242时，模型稳定，此时其中与模型三的解相同，说明新生幼仔的数量比较稳定，而每年注射避孕药的雌性动物则减少到352头，其中有110头是在注射后一年又被重复注射的数量，而则说明了该动物群落的个体数量是稳定的，被控制在一定范围内。

由此可知当每年注射352头时，就能把该动物群落的数量控制在11000头范围内。

问题三:

假如每年转移50至300头此动物到别处，其大致模型与问题二相近，不同之处在于要考虑到被运走的动物的数量（b），即应表示该种群每年的新生幼仔的数量减去该年死亡个体的数量与运走个体数量的和的差值。

相应的避孕措施将改变如下：

模型五：

其中b表示该年被运走的动物的数量（50300）；代入已有数据，并把方程标准化，得到以下模型：

其中b的值在50300之间变化，我们以50头为间距（即b的值分为50，100，150，200，250，300六种情况），利用软件MATLAB6.0计算得到每年转移从50至300头到别处时所应采取的避孕措施（源程序见附录7），当取b50时，得到以下输出结果：

由输出结果的提示可知，在100000次迭代后仍然无法找到最优值；由exitflag=-1可知，以上模型得到的数据是发散的，从而说明上述模型本身并不存在问题，只是计算机无法在有限的迭代次数内找到最优值。

为此解决该问题，我们对以上模型进行修改，把目标函数定义为该动物群落在年底时的总数量与年初的数量加上被运走的数量b的差值，从而得到以下模型：

模型六：

其中式

（1）表示目标函数，其值趋于0；式

（2）表示该动物在第k年初的数量，其值取11000；式（3）和（4）表示该动物第i龄到了年底全部转化为（i1）龄；式（5）和（6）表示该动物各年龄段的变化；式（7）表示该动物新生的幼仔数量；式（8）表示该动物在1360岁的雌性个体数量之和；对b以10为间距，分别从50取到300头，利用MATLAB软件进行求解（源程序见附录8），得到以下表二：

被运走的数量（头）5060708090100110120130140被注射的数量（头）370361351342332322313304294285被运走的数量（头）150160170180190200210220230240被注射的数量（头）275266256247237228218209199190被运走的数量（头）250260270280290300被注射的数量（头）180171161152142133问题四：

用人工手段可以使得该动物群落的数量减少，使其数量控制在一定范围内，但是，我们并不知道用了这种人工手段后，对于一次突发事件后要增加其数量时，该种群的复壮能力如何。

因此，我们建立数学模型，分析研究其复壮能力。

我们假设由于某种原因（例如由于一场大灾难导致该动物的死亡），使得该动物减少至M头，此时不得不停止注射避孕的方法。

在初步模型中，我们先假设灾难发生后该种群虽然数量减少，但其年龄结构没有改变。

建立该动物群落重新壮大的能力的模型如下：

模型七：

对M取不同的值，利用c语言进行编程求解，（源程序见附录9）可以得到以下表：

M值第t年1000090008000700060005000400030002000t=1增量d16651066559495945268526456745638863883225322250425017931791132113当前数量M1t=2增量d270711372638102325629088489794541460823425664267451719133701222235当前数量M2t=3增量d3757恢复699109326029690521846644472463656029286480320535751302365当前数量M3t=4增量d4729恢复64010330557902347277183906419305510821937941402505当前数量M4t=5增量d568911019596961950782254186837325543323540291452650当前数量M5t=6增量d6741恢复6381025753887924457272348578125142801562806当前数量M6t=7增量d76821092958193734727744371615226845481662972当前数量M7t=8增量d8729恢复62399965078251401655328645481783150当前数量M8t=9增量d9665106615438794426697930748351913341当前数量M9t=10增量d10707113685819375456743532751422043545当前数量M10增长加速度：

46.000046.600042.685738.881035.642428.612122.806116.436410.7636该表对M进行不同的取值，分别表示某一年（t0）灾难发生后该动物群落剩下的数量。

我们用第一题的C程序即可算出该群落在下一年（t1）的增量，把这个增量加上上一年群落的总数即为第一年（t1）群落的总数M1,接着我们又用程序算出这一年的增量，加上M1即为第二年（t2）的个体总数M2。

以这样的算法推算下去，一直到群落数量恢复到11000左右或者前十年（t10）为止。

将计算出来的各年增量按时间t（1，210）列在座标上，用MATLAB进行曲线拟合（具体程序见附录