Logistic人口阻滞增长模型.docVIP

(完满版)Logistic人口阻滞增加模型 (完满版)Logistic人口阻滞增加模型 (完满版)Logistic人口阻滞增加模型 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 总人口 年份 1963 1964 1965 1966 1967 1968 1969 1970 1971 总人口 年份 1972 1973 1974 1975 1976 1977 1978 1979 1980 总人口 年份 1981 1982 1983 1984 1985 1986 1987 1988 1989 1 总人口 年份 1990 1991 1992 1993 1994 1995 1996 1997 1998 总人口 年份 1999 2000 2001 2002 2003 2004 2005 总人口 1、将1954 年看作初始时辰即 t 0 ，那么 1955为t 1，以次类推， 以 2005 年为t 51 作为终时辰。用函数〔 5〕对表 1 中的数据进行非线性拟合，运用 Matlab 编程获取相关 的参数 xm , r ，可以算出可决系数 〔可决系数是鉴识曲线拟合收效的 一个指标〕： 5 ( y i y? i ) 2 2 R 1 i 1 5 ( y y) i 2 0. 9959 i 1 由可决系数来看拟合的收效比较理想。 因此获取中国各年份人口变化趋势的拟合曲 线： x(t) 〔6〕 1 ( 1)e 依照曲线〔6〕我们可以对 2021 年〔t 56〕、2021年〔t 66〕、及2033 年〔t 79〕 进行展望得〔单位：千万〕 ： x (56) 138.6161, x (66) 148. 5400, x( 79) 结果解析：从所给信息可知从 1951年至 1958年为我国第一次出生人口顶峰，形成 了中国人口规模“由缓到快〞的增加基础；因此这段时期人口颠簸较大，可能影响模型 结果的正确性。 1959、1960、1961年为三年自然灾害时期，这段时期人口的增加碰到很 大影响， 1962 年处于这种影响的滞后期，人口的增加也碰到很大影响。总的来说 1951-1962 年的人口增加的随机误差不是遵从正态分布， 程序： 2 结果： 2、 将 1963年看作初始时辰即 t 0，以 2005年为 t 32作为终时辰。运用 Matlab 2 编程获取相关的参数 xm ,r ，可以算出可决系数 0. 9994 R 获取中国 各年份人口变化趋势的另一拟合曲线： x(t ) 〔7〕 1 ( 1)e 依照曲线〔7〕我们可以对 2021 年〔t 47〕、2021年〔t 57〕、及2033 年〔t 70〕 进行展望得〔单位：千万〕 ： x( 47) , x( 57) , x (70) 结果解析：1963年-1979 年此间，人口的增加根本上是依照自然的规律增加，特别 是在农村是这样，城市碰到收入的影响，生育率较低，但都有规律可寻。总的来说，人 口增加的外界大的搅乱因素根本上没有，可以认为这一阶段随机误差遵从正态分布； 1980-2005 年这一时间段，诚然人口的增加碰到国家方案生育政策的控制，但方案生育 的政策是根本牢固的，这一阶段随机误差也应遵从正态分布，因此用最小二乘法拟合所 获取的结果应有较大的可信度。 程序： 3 结果： 3、从 1980-2005 年，国家