目录收集数据的简单分析结果
代码
1、matlab代码
我参考了BP神经网络的一个实例，这个例子的数据我都给了
%% I .环境变量clear allclc%% II .生成培训集/测试集%1.导入数据x=xlsread ('人口. xlsx '，1，' A2:A72 ' )； y2=xlsread ('人口. xlsx '，1，' C2:C72 ' )； y3=xlsread ('人口. xlsx '，1，' D2:D72 ' )； y4=xlsread ('人口. xlsx '，1，' E2:E72 ' )； y5=xlsread ('人口. xlsx '，1，' F2:F72 ' )； y6=xlsread (人口. xlsx )，1，G2:G72 )； y7=xlsread (人口. xlsx )，1，H2:H72 )； y8=xlsread ('人口. xlsx '，1，' I2:I72 ' )； y9=xlsread ('人口. xlsx '，1，' J2:J72 ' )； y=[y2 y3 y4 y5 y6 y7 y8 y9]； y1=xlsread (人口. xlsx )，1，(B2:B72 ) )由8个因素测定； %输出的结果是人口数%% %2.随机生成训练集和测试集，预测2020年出生人口数。 此时，从第2行至第73行%temp=randperm(size(y，1 )-1 )读取。 %temp(72 )=72； % %训练集--62个样本，如果没有足够的代表性，训练模型就不好。 这里采用随机生成的方法。 %p_train=y(temp(1:62 )、 )； %P为输入，t为输出%t_train=y1(temp(1:62 )、 ) '； % %测试集--10个样品%year=x(temp(6:end )； %p_test=y(temp(6:end )， )’； %t_test=y1(temp(6:end )， )； %n=size(p_test，2 )； %固定训练% p _ train=y (1:61， )； % t _ train=y1 (1:61， )； %year=x(623360end； %p_test=y(62:end， )； %t_test=y1(62:end， )； %n=size(p_test，2 )； %随机数训练temp=randperm(size(y，1 )； p_train=y(temp(1:61 )， )； t_train=y1(temp(1:61 )， )； year=x (时间) 623360end )； p_test=y(temp(62:end )， )； t_test=y1(temp(62:end )， ) '； n=size(p_test，2 )； %% III .数据规范化[p_train，PS_input]=mapminmax(p_train，0，1 )； P_test=mapminmax(apply )、p_test、ps_input )； [t_train，PS_output]=mapminmax(t_train，0，1 )； % % iv.BP (后退传播)神经网络构建、培训和模拟测试(1.创建网络net=newff(p_train、t_train、9 ) %信号前向传播、误差反向传播、随机生成权重和阈值%单抑制层神经元%%2.设置训练参数net.trainParam.epochs=1000； net.trainParam.goal=1e-3； %培训要求精度net.trainParam.lr=0.001； 学习速度%%%3.训练网络net=train(net，p_train，t_train )； %%%4.模拟测试t_sim=sim(net，p_test )； %%%5.数据的反规范化t_sim=mapminmax('reverse '，t_sim，ps_output )； %% V .性能评估%%% 1.相对误差errorerror=ABS(t_sim-t_test )./T_test； %%% 2.决策系数r^2r2=(n*sum(t_sim.*t_test )-sum (t _ sim ) sum (t _ test ) ) (2/) (n*sum ) ) %% VI .绘图figureplot(1:N，T_test，' b:* '，1:n，T_sim，' r-o ' )； 20年真实值为用matlab拟合函数工具箱的拟合值，2020年预测值为神经网络模型训练预测值legend ('真实值)、'预测值)； %右上有xlabel (预测样本) (x轴坐标描述ylabel ) (人数(万) ) (y轴坐标描述set ) GCA、(XTickLabels )、Year )； string={ '测试人口预测结果比较'； ['r^2='num2str(R2 ) ]； title (字符串)； 数据收集2. 数据来源
C-J列我认为是影响人口每年出生率的因素1949-2019的数据，是我花了一天从国家统计局、年鉴、查阅论文或者是推算得到的数据，我就不给了，收集数据是要下功夫的，比如一年国赛中的出租车问题就需要去爬取数据，可以考虑爬虫~
简要分析
3. 分析
这个例子里面，由于我建模时间比较匆忙，所以其中输入神经元的八个特征值，我理想化认为都直接影响中国的每年出生人口的量纲。正常来说，应该对这八个因素进行逻辑回归分析和显著性水平检验，筛选掉一些显著性水平不达标的因素/特征值(例如就业人口等)。可以用多种模型和指标来进行探究。例如说AIC、BIC指标，引入激励函数，损失函数等方法。总而言之，还是很有意思的。
首先对这70年的人口数据进行分析，结合历史背景，我们可以清晰地看到在5年-10年的小时间段内，影响人口的主要因素。例如:
1)建国初期，国内外环境相对稳定，中国迎来了第一波婴儿潮；
2)1959年,1960年,1961年我国遭受三年自然灾害，图表1中可以看到这三年出生人口骤减。1961年-1963年是中国建国后70年人口增幅最大的两年，中国迎来第二波婴儿潮。
3)然而在2000-2015年的十五年之间，中国每年出生人口比较平缓。直到2016年之后，人口开始呈现明显的下降趋势。
图表1.中国每年出生人口的折线图
通过中国国家统计局和快易理财网的数据，共收集到九组数据，分别是中国每年出生人口、总人口统计(万/人)、人均GPA(按照美元算)、性别比例(按照女生等于100人计算)、人口自然增长率、城镇人口和乡村人口(两者加起来等于100)、美元兑换人民币汇率、总就业人口(万/人)。数据的范围从1949年不间断统计直到2019年为止。
我采用了BP 神经网络模型来预测中国每年的出生人口数据，表格中(详见“人口.xlsx”)每一行，分别利用总人口、人均GPA、性别比例、自然增长率、城镇人口、乡村人口、美元兑换人民币汇率、就业人口共八列数据作为特征值，来影响中国每年出生人口的数量。
为获得2020年特征值中的拟合数据
首先对八个特征值近10年-20年的数据进行分析，得到两种拟合的方向:
1、初步分析呈线性的数据:
包括“人口”表中的总人口、人均GPA、性别比例、城镇人口、乡村人口共5列的数据。我通过灰色模型GM(1,1)来对这五项数据进行预测，目的是得到未知的2020年的5项人口特征值。图表在GM(1,1)模型求解中展示
通过以上的分析，可以得到近似地拟合得到2020年的总人口数为140908万人，人均GDP为11158美元，男女比例为104.58:100，城乡村人口比例为61.65:38.35。之所以采用近十年来预测，是因为近十年我国在科技、军事、文化等方面飞速发展，综合国力和文化软实力显著提升，采用近十年的数据能较好地反应出中国在这一段时间内数据的总体趋势。
图表2.GPA
2、初步分析呈非线性的数据:
通过观察自然增长率、美元兑换人民币汇率和中国就业人口数据，可以发现近10年-20年，这些数据都是上下波动的。为此，我利用matlab数据拟合的工具箱Curve Fitting进行的拟合。
图表3.自然增长率
以上图表我就不一一给出!!!
4.得出结论
matalb代码段的第三部分，我分别采用两种方式训练模型:
第一种是固定训练样本，采用了1949年-2010年共61/62年的训练样本，然后训练得到的结果非常不理想。我这里的评价指标采用了决定系数R^2以及相对误差矩阵。
图表4.评价指标
图表5.R^2太小了
图表6.2020年预测样本
当我们要预测2020年出生人口数据的时候，要加多2020年这一列，因为要输入八个特征值以及拟合预测结果:
导入数据的代码块全部要改为:
由于写的比较快，不足之处希望可以指正，希望大家国赛、美赛取得好成绩，考研保研出国都能一切顺利！
—————————————————————————————更新—————————————————————————————
之前有朋友在评论或者私信我说如何预测2060年或者2070年的数据，我在这里做一下解释。
在文章中我们要预测2020年的出生人口数据，首先我们要假定1950-2019年收集的数据是真实的，事实上也是我们从网站上爬取或者从文献中找到的数据，八九不离十。我们将这69年的数据喂入网络，用matlab的黑盒训练后，这个bp神经网络就有了一定的能力。这是什么能力呢？就是你输入一行八个特征值（上文中有提到，这也是我们定义的网络的输入神经元个数），网络就能给你输出一个数，这个数就是bp网络预测的出生人口数。
有不少人在看到这篇博文有点疑惑，我要预测2050年的数据，是不是要将1950年到2049年的数据都进行输入呢，也就是将 100 × 8 100times8 100×8的矩阵输入到网络中，让网络学习。答案是否定的。这是为什么呢？因为你要获得2021-2049年的数据（现在是2021年底，今年的数据应该还没统计出来），只能通过拟合值来估计。就好像上文中我将每一列的数据通过matlab工具箱预测，得到2020的八个特征值输入数据（之前写这篇博文的时候是2020年，还没有2020年的数据），但这不一定准确，很有可能是fake的，存在比较大的出入。我用2000-2019预测2020的八个特征值都不一定准确了，更加何况预测今后30年的 30 × 8 30times8 30×8的数据呢？那可能是灾难性的错误，如果将这些特征值输入网络，会导致网络的理解能力有很大偏差，甚至是不稳定。
但毫无疑问，你要让网络输出2050年的出生人口数，你得有2050年的八个特征值数据呀（人均GPA、性别比例等）。最简单的方法就是通过matlab工具箱预测，选择一个合适的拟合曲线去拟合数值。虽然这个数据不一定准确，但是我们的bp神经网络的理解能力是基本上准确的，因为输入的数据都是没有问题的。到这里，你应该明白了吧。最后一点我还想提一下，就是如果你要预测未来n年的数据，最好不要选择1978之后的年份数据输入网络中训练，因为在那一年的春天之后，我们的生活发生了很大的变化，陈旧的数据（生产力的差异太大）对预测未来还是有一定的影响。