摘 要：本文首先简要介绍Malthus和Logistic两种单种增长模型，然后详细介绍双种竞争的Volterra模型，最后介绍了多种的Gause-Lotka-Volterra 和三种的RPS博弈模型，对其做了比较和分析，得出了一些有益的启示。
为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。本文首先简要介绍Malthus和Logistic两种单种增长模型，然后详细介绍双种竞争的Volterra模型，最后介绍了三种的Gause-Lotka-Volterra 和RPS博弈模型。
一般生态系统的分析可以通过一些简单模型的复合来研究，根据生态系统的特征建立相应的模型。种的数量本应取离散值，但由于种数量一般较大，为建立微分方程模型，可将种数量看作连续变量，甚至允许它为可微变量，由此引起的误差将是十分微小的。
1.1 马尔萨斯（Malthus）模型
马尔萨斯在分析人口出生与死亡情况的资料后发现，人口净增长率r基本上是一常数，（r=b-d, b为出生率，d为死亡率），既：
或 (1)
其解为
(2)
其中N0=N(t0)为初始时刻t0时的种数。
马尔萨斯模型的一个显著特点：种数量翻一番所需的时间是固定的。
令种数量翻一番所需的时间为T，则有：
(3)
(4)
人口统计数据与Malthus模型计算数据对比：
年 | 1625 | 1830 | 1930 | 1960 | 1974 | 1987 | 1999 |
人口(亿) | 5 | 10 | 20 | 30 | 40 | 50 | 60 |
| | | | | | | |
表1 世界人口数量统计数据
年 | 1908 | 1933 | 1953 | 1964 | 1982 | 1990 | 2000 |
人口(亿) | 3.0 | 4.7 | 6.0 | 7.2 | 10.3 | 11.3 | 12.95 |
| | | | | | | |
表2 中国人口数量统计数据
比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如，1961年世界人口数为30.6亿（即3.06×1010），人口增长率约为2%，人口数大约每35年增加一倍。查1700年至1961年共260年的人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每34.6年增加一倍，两者也几乎相同。
卢允忠假如人口数真能保持每34.6年增加一倍，那么人口数将以几何级数的方式增长（如图1）。例如，到湿法炼锌2515年，人口约达2×1014人，即使海洋全部变成陆地，每人也只有9.3平方英尺的活动范围，而到2665年，人口约达4×1015人，只好一个人站在另一人的肩上排成二层了。故马尔萨斯模型是不完善的。
图1 马尔萨斯模型人口预测呈几何级数增长
Malthus模型实际上只有在体总数不太大时才合理，到总数增大时，生物体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象，后续将展开讨论。所以Malthus模型假设的人口净增长率不可能始终保持常数，它应当与人口数量有关。
1.2 Logistic模型
人口净增长率应当与人口数量有关，即：r=r(N)，
从而有：
(3)
对马尔萨斯模型引入一次项（竞争项），令 r(N)=r-aN
此时得到微分方程：
或 (4)
水培鱼缸 (4)式被称为Logistic模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（Verhulst）首先提出的。一次项系数是负的，因为当种数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。
(4)式可改写为
(5)
(5)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种个体，当种数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种数量的上界为K（近似地将K看成常数），N表示当前的种数量，K-N恰为环境还能供养的种数量，(5)式指出，种增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是(5)式也被称为统计筹算律的原因。
对(5)式分离变量：
(6)
两边积分并整理得：
(7)
令N(0)=N0，求得：
(8)
故(5)的满足初始条件N(0)=N0的解为：醋酸乙酯生产
(9)
易见：
N(0)=N芳纶头盔0， (10)
N(t)的Logistic模型曲线如图2
图2 Logistic模型曲线图
历史上已有多个种增长的实验验证了Logistic模型所描述的规律。1945年克朗皮克（Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（E·F·Gauss）也做了一个原生物草履虫实验，实验结果都和Logistic曲线十分吻合。
大量实验资料表明用Logistic模型来描述种的增长，效果还是相当不错的。例如，高斯把5只草履虫放进一个盛有0.5cm3营养液的小试管，他发现，开始时草履虫以每天230.9%的速率增长，此后增长速度不断减慢，到第五天达到最大量375个。
实验数据与r=2.309，a=0.006157，N(0)=5的Logistic曲线：，几乎完全吻合，见图3。
图3 r=2.309，a=0.006157，N(0)=5的Logistic曲线图
1.3 Lotka-Volterra模型
20世纪40年代，Lotka和Volterra分别奠定了种间竞争关系的理论基础，他们提出的种间竞争方程对现代生态学理论的发展有着重大影响，该模型因此用两人的名字命名为Lotka-Volterra模型，该模型是对Logistic模型的延伸。
意大利生物学家D’Ancona曾致力于鱼类种相互制约关系的研究，在研究过程中他无意中发现了一些第一次世界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很理想的鱼类占总渔获量的百分比。在 1914-1923年期间，意大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加，如表3：
年份 | 1914 | 1915 | 1916 | 1917 | 1918 | 1919 | 1920 | 1921 | 1922 | 1923 |
比例 | 11.9 | 21.4 | 22.1 | 21.2 | 36.4 | 27.3 | 16.0 | 15.9 | 14.8 | 透水混凝土施工方案10.7 |
| | | | | | | | | | |
表3 意大利阜姆港食肉鱼占总渔货量的比例
捕获的各种鱼的比例可以近似地反映地中海里各种鱼类的比例。战争期间捕鱼量大幅下降，但捕获量的下降为什么会导致鲨鱼、鳐鱼等食肉鱼比例的上升，即对捕食者有利而不是对食饵有利呢？D’Ancona百思不得其解，无法解释这一现象，就去求教当时著名的意大利数学家V.Volterra，希望他能建立一个数学模型研究这一问题。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议