(完整版)Logistic人口阻滞增长模型.doc

Logistic 人口阻滞增长模型 一、模型的准备 阻滞增长模型的原理：阻滞增长模型是考虑到自然资源、环境条件等因素对人口增 长的阻滞作用，对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口 增长率 r 的影响上，使得 r 随着人口数量 x 的增加而下降。若将 r 表示为 x 的函数 r (x) 。 则它应是减函数。于是有： dx dt r (x)x , x (0) x 0 （1） 对r (x) 的一个最简单的假定是，设 r(x) 为x 的线性函数，即 r (x) r sx (r 0 , s 0) （2） 设自然资源和环境条件所能容纳的最大人口数量 x ，当x xm 时人口不再增长， 即 m 增长率 r ( ) 0，代入（ 2）式得 x m r s ，于是（ 2）式为 x m x r( x) r (1 ) （3） x m 将（3）代入方程（ 1）得： dx dt rx (1 x ) m （4） x x(0) x 0 解方程（4）可得： x(t ) 1 ( x m x m x 0 1)e rt （5） 二、模型的建立 我国从 1954 年到 2005 年全国总人口的数据如表 1 表 1 各年份全国总人口数（单位：千万） 年份 1954 1955 1956 1957 1958 1959 1960 1961 1962 总人口 60.2 61.5 62.8 64.6 66.0 67.2 66.2 65.9 67.3 年份 1963 1964 1965 1966 1967 1968 1969 1970 1971 总人口 69.1 70.4 72.5 74.5 76.3 78.5 80.7 83.0 85.2 年份 1972 1973 1974 1975 1976 1977 1978 1979 1980 总人口 87.1 89.2 90.9 92.4 93.7 95.0 96.259 97.5 98.705 年份 1981 1982 1983 1984 1985 1986 1987 1988 1989 1 总人口 100.1 101.654 103.008 104.357 105.851 107.5 109.3 111.026 112.704 年份 1990 1991 1992 1993 1994 1995 1996 1997 1998 总人口 114.333 115.823 117.171 118.517 119.850 121.121 122.389 123.626 124.761 年份 1999 2000 2001 2002 2003 2004 2005 总人口 125.786 126.743 127.627 128.453 129.227 129.988 130.756 1、将1954 年看成初始时刻即 t 0 ，则 1955为t 1，以次类推， 以 2005 年为t 51 作为终时刻。用函数（ 5）对表 1 中的数据进行非线性拟合，运用 Matlab 编程得到相关 的参数 xm 180.9871 , r -0.0336 ，可以算出可决系数 （可决系数是判别曲线拟合效果的 一个指标）： 5 ( y i y? i ) 2 2 R 1 i 1 5 ( y y) i 2 0. 9959 i 1 由可决系数来看拟合的效果比较理想。 所以得到中国各年份人口变化趋势的拟合曲 线： 180.9871 x(t) （6） 180.9871 0.0.0336t 1 ( 1)e 60.2 根据曲线（6）我们可以对 2010 年（t 56）、2020年（t 66）、及2033 年（t 79） 进行预测得（单位：千万） ： x (56) 138.6161, x (66) 148. 5400, x( 79) 158.6028 结果分析：从所给信息可知从 1951年至 1958年为我国第一次出生人口高峰，形成 了中国人口规模“由缓到快”的增长基础；因此这段时期人口波动较大，可能影响模型 结果的准确性。 1959、1960、1961年为三年自然灾害时期，这段时期人口的增长受到很 大影响， 1962 年处于这种影响的滞后期，人口的增长也受到很大影响。总的来说 1951-1962 年的人口增长的随机误差不是服从正态分布， 程序： 2 结果： 2、 将 1963年看成初始时刻即 t 0，以 2005年为 t 32作为终时刻。运用 Matlab 2 编程得到相关的参数 xm 151.4513 ,r 0.0484 ，可以算出可决系数 0. 9994 R 得到中国 各年份人口变化趋势的另一拟合曲线： 151.4513 x(t ) （7） 151.4513 0.0484t 1 ( 1)e 69.1 根据曲线（7）我们可以对 2010 年（t 47）、202