动物群体的微分方程模型分解

1动物群体的微分方程模型§1引言§2单种群模型与人口问题§3进行开发的单种群模型§4弱肉强食模型§5竞争排斥模型§6竞争排斥原理的数学分析§7无管理的捞鱼模型2ACM-85试题A的标题是“动物群体的管理”，题文曰：“一种资源有限（即有限的食物、空间、水等）的环境里发现天然存在的动物群体，试选择一种鱼类或哺乳动物（例如北美矮种马、鹿、兔、鲑鱼、带条纹的欧洲鲈鱼等）以及一个你能获得适当数据的环境，并建立一个对该动物群体捕获量的最佳方案。与这一试题有相同或相似数学模型问§0引言3题非常之多，例如人口问题，生态与动植物保护的问题，种群之间的竞争排斥问题，等等，这些涉及人口与社会发展、生态与社会发展的重要问题，理应成为数学建模当中急需考虑的内容。本讲用常微分方程这一数学模型定量地或定性地讨论此类问题的建模思想与方法。4动植物种群本身是离散变量，谈不上可微，但由于突然增加或减少的只是单一个体或少数几个个体，与全体数目相比，这种增量是很微小的，所以，可以近似地假设大规模种群随时间是连续地甚至可微地在变化，进而可以引用微分方程这一数学工具来研究。§1单种群模型与人口问题5英国人马尔萨斯（Malthus，1766-1834）认为人口的净增长率为常数，即单位时间内人口增量与人口总量成正比，设t时刻人口数位p(t)，则有Malthus人口模型00,.0()dpapaconstdtptp6用此模型估算1700—1961年间的人口数目，计算结果与人口实际情况竟然惊人地相似。但是，当t+∞，计算结果p(t)+∞，具体地说，此模型可以求得2510年的人口总数为2000亿左右，可见，这一模型必须进行修正。问题出这个Cauchy问题的解为0()0()attptpe7在Malthus只看到繁衍增长的一面，未看到种群内竞争（如人类战争）对种族发展的抑制作用。1837年荷兰生物数学家Verhulst考虑单了种群成员间冲突乃至残害现象，得出容易理解的下述单种族数学模型：00()()dpabppdtptp8其解为：（2）美国和法国都曾用这个公式预报过人口变化，结果相当符合实际。显然（3）0()000attappabpape(t)limtapb9a=0.029，b可以如下求得：1980年5月1日，我国公布的人口总数1979年底为97092万人，当时人口增长率为1.45%，于是a-b×9.7092×108=0.0145，从而求得：b，及=19.42(亿），即，我国的人口极限约为19.42亿人。ba+∞10养鱼场从鱼池中捞鱼出售，每次捕捞得太少不合算，一方面销售收入少，而且池中鱼过多也不利于鱼群生长繁衍，但每次捞得过多，“竭泽而渔”，显然也不可取，应怎样控制捕捞率，使得总经济效益最优？设单位时间内捕捞h条鱼，t时刻池中鱼数为N(t)，则N(t)满足下列数学模型：§3进行开发的单种群模型11hNKKrNdtdN)(（4）02hNKrrN其中K是鱼池中鱼数的最大值（受池子条件限制，此最大值是存在的。h称为收获率。考虑dN/dt=0时，即，124rKh当4rKh得到时，dN/dt0，此时，池中鱼数单调递减，长此下去将无鱼可捞，所以，是最大可承受的产量。4rKh13当时，有两个正的平衡点4rKhrhKKKN42121rhKKKN42122（5）这样，模型（4）可以写成))((21NNNNKrdtdN14可见，当t增加时，N=N1附近的N=N(t)远离N=N1这一水平线（在Nt平面，t为横轴），而在N=N2附近N=N(t)趋近于N=N2这一水平线，N=N1，N=N2是平凡解，即，解N=N1是不稳定的，N=N2是稳定的。;0dtdN当NN1（N2）时，;0dtdN当N2NN1时，当NN2时，;0dtdN15初始时刻，池中鱼数N(t0)N1，则单调下降趋于零，池中鱼会捞净灭绝；而N(t0)N1时，则池中鱼数量将自动调节随时间之增加趋于N2条鱼，又由dtdNrhKKKN42121可见h越小，N1越小所以，),40(rKh16一般要用小收获率h来开发低密度的种群，而用大收获率去开发高密度的种群。反之由可以解得)(01tNNKtNKtrNKKtNKrh))()((4]))(2([00020217即应控制收获率h不要超过KtNKtrN))()((00否则，将无鱼可捕。从上面讨论知，收获率h与种群密度是相关的，密度小时收获率亦应小。令收获率h=kN，k称为捕捞率。由（5）知，是（4）的平凡解，此时221KNN18收获率是最大可承受的单位时间内的产量。可见，欲使池中鱼不至于随时间之增加而趋于灭绝，又使产量最大，仅当池中鱼是最大可能鱼数之半时才可能。这时，从4rKh0)(kNNKKrN得平衡点为0KrkrN19KrkrK2即得r=2k，即鱼的增长率是捕捞率的2倍时，才达到最大收获量.4rK（r≤k则是“败家式”捕捞，不可行），于是下面分析在多大捕捞量时净利润最大。假设价为p元，又开支与捕捞率k成正比，则净利润为：20ckpkNz(6)在池鱼数稳定的条件下，即0)(kNNKKrN)1(KNrk时的利润可写为（上式代入（6））：21)1()()(KNrcpNNz（7）求函数（7）的最大值得知当pcKN22时（7）取最大值。这时捕捞量为：)1(4)22)](22(11[)22(222KpcrKpcKpcKKrpcKk22这时的捕捞量比最大捕捞量小，要少捞一些，少捕4rKKprc224)1(4222KpcrK捕捞开支c越大，越应该少捞一些，鱼价越高，越应该多捞一些，总之，欲使净收入最大，单位时间捞鱼量为23生活在同一环境中的各类生物之间，进行残酷的生存竞争，一类动物靠捕食另一类动物为生，被捕食者只能靠又多又快地繁殖后代和逃跑等方式求生存发展，如此等等。设想一海岛，居住着狐狸和野兔，狐吃兔，兔吃草，青草如此之茂盛，兔子们无无食之忧，于是大量繁殖。兔子一多，狐易得食，狐量亦增。而由于狐狸数目增§4弱肉强食模型24多吃掉大量的兔子，狐群又进入饥饿状态而使其总数下降，这时兔子相对安全些，于是兔子总数回升。这样，狐兔数量交替增减，无休止地循环，遂形成生态的动态平衡。意大利著名生物数学家沃特拉（Volterra)对上述现象建立了下述模型)()(dxcydtdybyaxdtdx（8）25其中x(t)表示t时刻兔子的数目，y(t)是狐狸数，ax项表示兔子繁殖速度与兔子现存总数比例，-bxy项表示狐兔相遇兔子被吃的速度，-cy项表示狐狸因为同类竞争食物造成的死亡速度与狐狸数成正比，+dxy项表示狐兔相遇对狐狸有好处而使狐狸繁衍增加的速度。看来这一模型表达了达尔文主义思想，而且数学分析之后还会充实和精确表达上述直观思想。26方程组等价于xdxdxcydybya)()(积分得KdxxcbyyalnlnlndxcbyaeKxey（9）27从（9）解不出y=f(x)这种显式解，沃特拉发明了一种巧妙的办法：在xOy平面上画出x(t)与y(t)变化相关性的相图。令dxcbyaeKxweyz,)00(00bydxaceyxK其中K由初始值x0，y0定出为于是绘出图5-128图5-129在L4上，随t的增加，动点（x(t)，y(t))依逆时针而动，事实上，点s是使L1：z=w;L2：z=yae-by；L3：w=Kx-cedx；L4：狐兔曲线。0yx30的平衡点（或称奇点），考虑点P2，P2的横坐标大于,故在P2点，，y增加，在P2处向上运动，可见是逆时针运动。dc0y现在考虑对两个物种同时进行捕捉，既抓兔子也捉狐狸，于是，模型（8）变成修正模型：dxyycydxycydtdybxyxaxbxyaxdtdx)()(（10）31从图5-1中已经看到，x(t)，y(t)是周期为T的周期函数，同理（10）的解x(t)、y(t)也是周期函数。对于（8），x(t)，y(t)的平均值为：yx,TTdttyTydttxTx00)(1)(132byaxxTTdttbyaTdtxxT00))((11TxTxdtxx00)0(ln)(ln又得：而33TTaadtTdttbyT001)(1故于是baydcx同理可得34对于（10）则得由（11）可知，当捕捉率不超过兔子的繁殖率a时，兔子反而会增加，狐狸要减少，反过来，捕捉率降低，平均而言，会增加狐狸的数目，而减少兔子的数目。)0(,abaydcx（11）35意大利生物学家棣安奇纳（D.Ancona）发现，第一次世界大战那些年代，地中海各港口捕鱼量百分比表明，掠肉鱼（例如鲨鱼）的百分比急剧增加，从上述数学分析中，对这种现象已经有了理论上的解释。事实上，那时战火连天，渔民大量停业，使捕捉率下降，所以相当于狐狸的掠肉鱼明显增加。这种结论在农业防治病虫害上有很大意义，例如，有两个物种（可能是两36种昆虫或害虫与青蛙等），一者是作物的害虫，一者是害虫的天敌，若施农药不当，虽然可以杀灭一些害虫，但同时也杀死了害虫的天敌，这一“捕捉行为”的实施，由上述结论知，可能造成天敌的减少，害虫的增多，事与愿违，与其施用少量农药治虫，不如采用生物治虫的办法。37§5竞争排斥模型在自然界中不难发现这种现象，两种生物为了争夺有限的同一食物、生活空间或配偶，进行着激烈的斗争。达尔文在《物种起源》一书中明确指出：“最剧烈的斗争，差不多总是发生在同种的个体，因为它们居住在同一地域，需要同样食物，遭受同样威胁。在同种的变种之间，其斗争之剧烈，大体如此，且有时在短期内即见胜负。”这里用数学模型及其解的定性分析来论证达尔文的上述思想。两种相似的38生物之间为争夺生存条件而斗争，直至其中一种生物物种完全灭绝才会中止的现象称为“竞争排斥原理”。这一原理的生物学解释是：已知生物群体在群落中有何种习性、食物和生活繁衍方式等，叫这一种群体“生态龛”。两种同类群体，难以占有同一生态龛。事实上，如果两个群体力图持有同一个生态龛，那么他们之间的生存竞争将是异常之激烈，且以弱者灭亡而告终。生态龛也可称为“小环境”。39在单种群模型中2bNaNdtdN且当t时，batN)(.Kba记这个极限可以认为是这个环境中可以承受的生物体最大数量。又40)()1()1(2KNKaNKNaNNabaNbNaNdtdN（12）（12）可以解释如下：当N很小时，N(t)按照马尔萨斯定律aNdtdN41增长，aN叫“生物势”，它是理想条件下，物种的可能增长率。只要对食物、配偶和空间不加限制，又无各个成员因排泄等造成的对环境毒化引起流行病害，这种增长率是可以实现的。但是，随着总数的增加，随的减少而减少。dtdN)(KNK今设N1(t)，N2(t)分别为物种A和物种B在时刻t的数量，K1和K2分别是A与B在小天地中最大可能的个数，那么，N1(t)，N2(t)满足下面的数学模型（设K1≠K2):42)()(21222221211111KmNKNadtdNKmNKNadtdN(13)其中m2为第二物种B占据A的位置的数量，m1为A占据B的位置的数量。m2=αN2，m1=βN1，如果A和B占有不同的生态龛，利害不冲突。当α=β=1，这时（13）变成：43)()(22122221211111KNNKNadtdNKNNKNadtdN（14）44§6竞争排斥原理的数学分析为了从数学上分析（14）中N1(t)，N2(t)的渐近性态，先介绍一些常微分方程定性理论的概念和结论。称方程组),(),(yxQdtdyyxPdtdx（15）为平面自治系统