我国以前一直以来都是世界上大豆生产的第一大国。
但由于各国的日益强大，导致我国豆种植面积和产量持续缩减。
因此，预测我国的大豆产量对中国未来的经济发展有着极其重要的作用。
我们为一位客户进行了短暂的咨询工作，他正在构建一个主要基于ARIMA的大豆产量预测应用程序
运用SAS与eviews软件对全国1957年到2009年大豆产量的历史数据建立时间序列ARIMA模型，通过判断其稳定性与是否通过白噪声检验，建立AR（2）模型。最后，因通过残差白噪声检验，对数据进行分析并预测未来五年的我国的大豆产量。
一、问题分析
预测未来5年我国的大豆产量及其发展趋势，对国民经济与人民生活都是百利而无一害的。同时，还可以提前让国家了解未来的情况，及时作出应对措施。
我国的大豆出产量的数据，在此假设所有数据真实可靠（详见下表1），并假设在预测期内不发生任何影响我国的大豆出产量的突发事件。
difx=dif(x);
year=intnx('year', **'1jan1957'd**,_n_-**1**);
二、模型识别
首先，根据原始数据绘制时间序列图形，观察序列特征。图形如下：
plot x*year difx*year;
视频
时间序列分析模型 ARIMA-ARCH GARCH模型分析股票价格数据
视频
在Python和R语言中建立EWMA，ARIMA模型预测时间序列
图1 我国的大豆产量时序图
由原始数据的时间序列图可以看出，资料数据呈现明显的上升趋势，其平均数不稳定，是非平稳序列。
图2 原始序列的自相关图
上图是自相关函数的结果，自相关函数衰减到0的速度缓慢，由此可以再次确定序列是非平稳的。
因为原序列呈现出上升的趋势，故选择1阶差分。1阶差分后的时序图如下所示。
图3 1阶差分后的时序图
由该图可以看出差分后的时间序列在均值附近比较稳定地波动。为了进一步确定平稳性，考察差分后序列的自相关图（如下所示）。
接下来，我们对差分后的时间序列进行ARMA模型的建立。季节差分后数据的自相关函数如下：
图4 1阶差分后的自相关系数图
随时关注您喜欢的主题
从上面的分析结果可以看到自相关图显示很强的短期相关性，所以可以初步认为1阶差分后序列平稳。
随后，对1阶差分后序列进行白噪声检验，结果如下图所示。
在检验的显著水平取为0.05时，由于上述所有延迟阶数的P值都小于0.05，所以该差分后的序列不能视为白噪声序列，即差分后序列还蕴藏着不容忽视的相关信息可供提取。
1阶差分后序列的自相关图已经显示该序列自相关系数具有拖尾的性质。再考虑其偏自相关系数的性质（见下图）。
图6 偏自相关系数图
根据自相关图和偏自相关图的特点，进行模型的定阶。偏自相关图显示，延迟1阶和2阶的偏自相关系数显著大于2倍标准差范围内波动，其他阶数的偏自相关系数都比较小。通过多方面的考虑，最后认为AR（2）模型为最优模型。接着，综合考虑前面的差分运算，实际上是对原序列拟合模型ARIMA(2,1,0)。
三、参数估计
在此，本文采用最小二乘法来估计参数，得到未知参数的估计值为：
图7 参数估计图
四、模型确定
由上面的输出结果可知拟合的方程如下：
图8 模型拟合结果图
该输出形式等价于
五、模型检验
本文需要检验残差是否有自相关性，由SAS的分析结果得知，不存在自相关性，即残差序列通过白噪声检验。
图9残差白噪声检验
六、模型的预测
本文给出了后面5年的人口自然增长率预测值以及置信区间。
forecast lead=**5** id=year out=out;
plot x*year=**1** forecast*year=**2** l95*year=**3** u95*year=**3**/overlay;
图10 预测值以及置信区间
通过图示可以直观地看出该模型对序列的拟合效果良好。
七、 总结与建议
从对我国大豆出产量的预测值可以看出，大豆的产量会相对提高，不过提高的速率跟之前的相比也不会相差太大。因此，为了提高我国的大豆产量，我提出了以下几点：
1、为大豆生产提供技术支撑。目前我国水稻、玉米和小麦的单产基本位于国际领先水平，大豆单产却与国际水平有较大差距，这也说明提高我国大豆单产还是很有潜力的。所以，我们应该选育和改良大豆品种，为大豆生产提供优良的种质资源。
2、充分利用我国的自然资源，进一步扩大大豆的种植面积，为扩大我国大豆生产规模和优化区域布局提供保障。
3、提高我国大豆在国际市场的竞争优势，以食品的安全性为主打，向外国销售。
八、 参考文献
[1]应用时间序列分析（第三版），王燕 编著 中国人民大学出版社
可下载资源
关于作者
Kaizong Ye是拓端研究室（TRL）的研究员。在此对他对本文所作的贡献表示诚挚感谢，他在上海财经大学完成了统计学专业的硕士学位，专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。
本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。
非常感谢您阅读本文，如需帮助请联系我们！
MATLAB深度学习Transformer神经网络量化金融时间序列预测交易策略回测
MATLAB图注意力网络GAT多标签图分类预测可视化
Python用CEEMDAN-LSTM-VMD金融股价数据预测及SVR、AR、HAR对比可视化
Python注意力机制Attention下CNN-LSTM-ARIMA混合模型预测中国银行股票价格|附数据代码
