数学模型讲座艾世猛种群增长模型在自然界和人类社会中.docVIP

数学模型讲座 艾世猛 第三讲 种群增长模型 在自然界和人类社会中，细菌、病毒的增殖；动植物的生长、繁殖；人口的增长等等是常见的现象。这些现象用生物学的术语来讲属于生态系统中种群数量的变化现象。这一讲就来研究一下这些数量变化的内在规律。对这些规律的了解和掌握可以使我们对不利的种群增长带来的风险加以规避；对有利的种群增长进行有效的保护，开发和利用。首先我们考虑人口的增长问题。 一、人口增长模型 世界人口在不断的增长，现在全世界每秒钟大约出生4.3人、每分钟大约出生259人、每小时大约出生15540人、每天大约出生37万人、每年增长约8296万人 年份 1625 1804 1927 1960 1974 1987 1999 2015 人口/亿 5 10 20 30 40 50 60 73 中国 人口 年份 1949 1953 1965 1982 1990 2000 2010 人口/亿 5.42 5.88 7.25 10.17 11.43 12.67 13.40 世界人口数据来源于《数学模型》（第四版），姜启源等编著，高等教育出版社 中国人口数据来源于《中国统计年鉴》，不含港澳台。 世界人口和中国人口的变化情况见下图1。从世界人口和中国人口的变化情况（表1和图1）可以看出，二者都在增长，而且很多人口是现代和当代以来迅速增加的。控制人口过快增长，提高人们的生活水平已经成为所有人的共识。下面就来研究一下人口的变化规律。 模型1：指数增长模型 最简单的人口增长模型类似于经济学中的产品增长模型。记今年的人口数量为，人口的年增长率为常数，则年以后的人口数为： （1） 模型（1）称为离散形式的人口增长模型，这是常用的一种简单预报人口的模型。从这个模型可以看出，随着时间趋向于无穷大，人口数量将趋近于无穷大。 图1、世界和中国人口增长图 模型（1）还有另一个连续型的版本，即著名的马尔萨斯[Malthus（1766-1834）]指数增长模型。这个模型是200多年前由英国人口学家Malthus提出来的。记时刻的人口为，并将函数视为连续可微的函数，记初始时刻（=0）的人口为。假设人口增长率为常数，即单位时间内人口的增量为乘以人口数。考虑时间段内人口的变化情况。这个时间段内的人口增长量即函数的增量为: 根据假设，这个增量等于增长率乘以人口数再乘以所用时间，从而有： = （2） 上式两边同时除以，再令取极限，得： = 根据导数的定义，上式左端就是函数关于变量的导数，而上式右端被取极限的函数与无关，取极限等于它本身，所以可得： 这是一个简单的微分方程，利用初始条件可以解得： （3） 上式（3）就是著名的马尔萨斯人口模型！是一个指数增长模型。 （思考题：模型（1）和模型（3）的关系是什么？） 指数增长模型的变化规律也是随着时间趋近于无穷，人口数量将增长为无穷大。这显然是不符合实际的，因为地球资源是有限的，养活不了无穷多的人口。指数增长模型只能用来预测短期内人口快速增长的情况。对于长期人口的预测就不合理了。 例如，指数增长模型可以较好的吻合美国1790-1900期间的人口变化情况（见表2和图2）。 表2、美国人口统计 年份 1790 1800 1810 1820 1830 1840 1850 1860 人口/百万 3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 年份 1870 1880 1890 1900 1910 1920 1930 1940 人口/百万 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 年份 1950 1960 1970 1980 1990 2000 人口/百万 150.7 179.3 204.0 226.5 251.4 281.4 数据来源于《数学模型》（第四版），姜启源等编著，高等教育出版社。 但是对于全部数据的拟合就存在较大的误差，如图3。指数增长模型不符合长期人口增长，所以要寻找新的模型来预测长期人口增长。 图2、指数增长模型拟合美国人口（1790年-1900年），参数取值x0 = 4.1884, r = 0.2743。（如何得到参数的取值？参见“最小二乘估计”） 图3、指数增长模型拟合美国人口（1790年-2000年），参数取值x0 = 6.045, r = 0.2022，存在较大误差。 模型2：逻辑增长模型（logistic 模型） 指数增长模型一个不符合实际的原因是假设增长率是常数，实际上增长率是变化的，由于资源，环境等的限制，人口增长到一定数量以后，增长率会下降（这一现象称为受到密度的制约）。所以增长率应该是人口数