基于MATLAB的人口预测模型

基于MATLAB的人口预测模型

摘要

本文以1980-2014年中国年终总人口数据资料为依据，分别使用了一次拟合、灰色预测模型和时间序列模型进行拟合，最终得出时间序列模型的效果最优，得到了中国人口数量逐年增长，但同时增长速度逐渐放缓的结论，为政府制定人口、经济政策提供了一定的依据。

关键词：人口数量；一次拟合；灰色预测；时间序列

前言

世界人口的迅猛增长引起了许多问题。特别是一些经济不发达国家的人口过度增长，影响了整个国家的经济发展、社会安定和人民生活水平的提高，给人类生活带来许多问题。为了解决人口增长过快的问题，人类必须控制自己，做到有计划地生育，使人口的增长与社会、经济的发展相适应，与环境、资源相协调。我国是世界上人口最多的发展中国家。人口数量多、增长快、可耕地少、国家底子薄，这是我国的基本国情。人口增长过快，严重制约着我国经济和社会发展的进程，影响着人民生活的改善和民族素质的提高。从而造成社会再生产投入不足，严重影响国民经济的可持续发展。认真分析我国目前的人口现状和特点，采取切实可行的措施控制人口的高速增长，提高人口的整体素质，已成为我国目前经济发展中需要解决的首要问题。

本文以中国近35年的人口数据尝试建立模型，分别建立了一次模型、灰色预测GM（1,1）模型和时间序列AR模型，最终选取了拟合效果最好的时间序列模型，用于说明我国人口问题以及预测短期内人口数量变化，以及为我国即将面临的人口问题提供一些建议。

概念与引理

定义1]1[：人口问题，是由于人口在数量、结构、分布等方面快速变化，造成人口与经济、社会以及资源、环境之间的矛盾冲突。

人口数量问题，主要由非均衡生育（多子化和少子化）以及人口迁移造成，只有通过均衡生育（发达国家2.17胎，发展中国家2.3胎）和调控迁移来解决。

人口结构问题，主要包括年龄、性别、收入、人种、民族、宗教、教育程度、职业、家庭人数等人口结构问题；其中最为突出的是年龄（多子化、少子高龄化）、性别（男女比例失调）和收入（基尼系数高、中产塌陷）结构问题。人口年龄结构问题，只有通过均衡生育（发达国家2.17胎，发展中国家2.3胎）来解决。人口性别结构问题，只有通过限制堕胎来解决。人口收入结构问题，原因较复杂，但最终都只有通过壮大中产阶层，使中产阶层成为社会主体才能真正解决。

人口分布问题，主要包括大城市病、高密度连绵城市群的环境污染问题、大片乡村缺少就近（200公里内）特大城市辐射带动的发展难题、生态气候等自然条件恶劣地区人口的生存困境、以及高密度大流量的人口迁移等问题。人口分布问题，主要是通过城镇化的合理布局，构建合理的城镇体系来解决。

定义2]2[：最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。

原理：在我们研究两个变量 x,y)之间的相互关系时，通常可以得到一系列成对的数据 x1,y1,x2,y2…x m,y m)；将这些数据描绘在x−y直角坐标系中，若发现这些点在一条直线附近，可以令这条直线方程如（式1-1）。

y1x（式1-1）其中：、1是任意实数

为建立这直线方程就要确定和1，应用《最小二乘法原理》，将实测值与利用计算值1x)（式1-1）的离差− )的平方和∑− )2最小为“优化判据”。

令：∑− )2（式1-2) 把（式1-1）代入（式1-2）中得：

∑−−1)2（式1-3) 当∑− )2最小时，可用函数对、1求偏导数，令这两个偏导数等于零。

∑1−)（式1-4)

∑1−)（式1-5)

亦即：

∑)1（式1-6)

∑)(∑2) 1 ∑)（式1-7)

得到的两个关于、1为未知数的两个方程组，解这两个方程组得出：

∑)−

1 ∑)

（式1-8)

1 ∑)∑∑

∑∑∑

（式1-9)

定义3]3[：灰色预测是一种对含有不确定因素的系统进行预测的方法。灰色预测通过鉴别系统因素之间发展趋势的相异程度，即进行关联分析，并对原始数据进行生成处理来寻找

系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程模型，从而预测事物未来发展趋势的状况。其用等时距观测到的反应预测对象特征的一系列数量值构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。

GM （1，1）模型

令x ) x ) 1),x ) ),…,x ) )),x(0)=(x ⑴，x ⑵，…，x(n)) 作一次累加, 生成消除数据的随机性和波动性,即

x 1) k ) ∑x

)

),k =1k 1, ,…,n 有x 1) (x 1) 1),x 1) ),…,x 1) ))

x ) 1),x ) 1) x ) ),…,x ) −1) x ) )) x 可建立白化方程：dx

dt ax u 即gm 1,1) 该方程的解为： x 1) k ) x 1) 1) u

a

)e ak u

a

x(k+1)=(x ⑴-u/a)exp( )+u/a

其中：α称为发展灰数；μ称为内生控制灰数

定义4]4[顺序排列而成的数列。时间序列分析的主要目的是根据已有的历史数据对未来进行预测。

AR 模型是一种线性预测，即已知N 个数据，可由模型推出第N 点前面或后面的数据（设推出P 点），所以其本质类似于插值，其目的都是为了增加有效数据，只是AR 模型是由N 点递推，而插值是由两点（或少数几点）去推导多点，所以AR 模型要比插值方法效果更好。

主要结果

数据：