尹曉燕 王旭陽 史澳 何曉豐 王雪 劉云

摘要：為了進(jìn)一步提高新疆棉花產(chǎn)量的預(yù)測精度，運(yùn)用GM模型與ARIMA模型對新疆棉花產(chǎn)量進(jìn)行了預(yù)測，進(jìn)而賦予合理權(quán)重建立了基于灰色理論和時(shí)間序列的組合預(yù)測模型對新疆棉花產(chǎn)量進(jìn)行了預(yù)測。實(shí)證分析結(jié)果表明：組合模型的預(yù)測結(jié)果更加逼近于真實(shí)數(shù)據(jù)，預(yù)測精確度較單一預(yù)測模型的結(jié)果準(zhǔn)確性更高、誤差更小。因此， GM模型與ARIMA模型組合預(yù)測新疆棉花產(chǎn)量是可行的。

關(guān)鍵詞：灰色理論;時(shí)間序列;組合模型;棉花產(chǎn)量;預(yù)測

中圖分類號(hào)：S562

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：2095-3143（2021）01-0015-07

DOI：10.3969/j.issn.2095-3143.2021.01.003

Abstract：In order to improve the forecast precision of Xinjiang cotton yield， the GM model and ARIMA model are used to predict Xinjiang cotton yield， the combination forecasting model based on grey theory and time series was established to predict the cotton yield in Xinjiang. The results of empirical analysis showed that the forecasting accuracy of the combined model is higher and the error is smaller than that of the single model. Therefore， it is feasible to combin the GM model and the ARIMA model to forecast cotton yield in Xinjiang.

Key words：Grey theory; Time series; Combination model; Cotton yield; Forest

新疆維吾爾自治區(qū)是我國主要的植棉區(qū)之一，同時(shí)也是我國最大的優(yōu)質(zhì)棉和唯一的長絨棉生產(chǎn)基地。在新疆農(nóng)業(yè)中，棉花產(chǎn)業(yè)占居主導(dǎo)地位，其產(chǎn)量的變動(dòng)對棉花行業(yè)相關(guān)從業(yè)者的經(jīng)營和發(fā)展有重大的影響。棉花產(chǎn)業(yè)的健康發(fā)展對新疆農(nóng)業(yè)及紡織工業(yè)的發(fā)展具有極其重要的意義[1]。相對于一熟旱地作物，棉花生長周期長，受自然環(huán)境的影響大。長期以來，棉花產(chǎn)量面臨較大的風(fēng)險(xiǎn)，波動(dòng)頻繁，需要進(jìn)行產(chǎn)量預(yù)測，以改變管理措施和調(diào)控市場。

近年來，基于不同的研究對象眾多學(xué)者運(yùn)用不同的方法對新疆棉花產(chǎn)量進(jìn)行了預(yù)測。目前學(xué)者研究的熱點(diǎn)是系統(tǒng)分析法，主要包括灰色理論預(yù)測法、人工神經(jīng)網(wǎng)絡(luò)預(yù)測法、時(shí)間序列預(yù)測法以及系統(tǒng)動(dòng)力學(xué)法等。這些方法對平穩(wěn)的時(shí)間序列取得較好的預(yù)測效果，得出了一系列重要的結(jié)論[2]。為了能夠合理的利用各種預(yù)測模型的優(yōu)點(diǎn)，Bates，等[3]將不同的預(yù)測模型進(jìn)行適當(dāng)?shù)慕M合，于1969年提出了組合預(yù)測的方法，該方法能夠綜合利用各種預(yù)測模型提供的信息，盡可能地提高預(yù)測的精度。近年來已有許多對組合預(yù)測模型得到了推廣。邢棉[4]提出了灰色神經(jīng)網(wǎng)絡(luò)組合模型，為較復(fù)雜的季節(jié)性預(yù)測問題，提供了一種新的有效的方法，并對電網(wǎng)供電量進(jìn)行了預(yù)測。朱輝[5]運(yùn)用了灰色預(yù)測模型，對新疆生產(chǎn)建設(shè)兵團(tuán)棉花產(chǎn)量進(jìn)行了預(yù)測，得到了精確度更高的預(yù)測值。陳黎明，等[6]運(yùn)用灰色預(yù)測和回歸模型建立組合預(yù)測模型對GDP統(tǒng)計(jì)數(shù)據(jù)質(zhì)量進(jìn)行了評(píng)估研究。單銳，等[7]為了更加精確的預(yù)測出人均GDP，采用了ARIMA模型與GM模型的組合模型，結(jié)果顯示組合預(yù)測結(jié)果更加接近于真實(shí)值。佟長福，等[2]為了減少隨機(jī)性，提出了基于小波分析理論的農(nóng)業(yè)需水量組合預(yù)測模型。李曉彤，等[8]通過建立灰色模型和ARIMA-GARCH模型的組合模型對我國的股票市場進(jìn)行了分析。張利平，等[3]通過建立六種時(shí)間序列組合模型并結(jié)合它們的優(yōu)點(diǎn)，使其擬合及預(yù)測效果達(dá)到最優(yōu)。為了對棉花的產(chǎn)量進(jìn)行動(dòng)態(tài)、科學(xué)、精準(zhǔn)地預(yù)測，作者運(yùn)用灰色理論與ARIMA模型對新疆棉花產(chǎn)量進(jìn)行預(yù)測，運(yùn)用相關(guān)統(tǒng)計(jì)軟件分析得出單一灰色GM（1， 1）模型與ARIMA模型的預(yù)測結(jié)果，之后結(jié)合二者各自的優(yōu)點(diǎn)賦予合理權(quán)重建立基于灰色理論和時(shí)間序列的組合預(yù)測模型，以期提高棉花產(chǎn)量預(yù)測的精度，為政府及棉花相關(guān)從業(yè)者制定合理的生產(chǎn)、銷售、儲(chǔ)存和加工等產(chǎn)銷策略提供參考與借鑒。

1灰色時(shí)間序列組合預(yù)測模型

1.1 GM（1， 1）模型

灰色預(yù)測是一種對具有不確定因素的系統(tǒng)進(jìn)行預(yù)測的方法，能有效解決數(shù)據(jù)少、序列的完整性及可靠性低的問題。GM（1， 1）模型是一種較為常用的灰色模型，GM（1， 1）預(yù)測模型的建立實(shí)質(zhì)上就是對原始數(shù)據(jù)序列作一次累加生成，使生成數(shù)據(jù)序列呈顯出一定規(guī)律，然后通過建立微分方程模型，求得擬合曲線，進(jìn)而對系統(tǒng)進(jìn)行預(yù)測。該模型構(gòu)建思路如下。

1.2 ARMIA模型

ARIMA模型（自回歸積分滑動(dòng)平均模型）是一種根據(jù)系統(tǒng)觀察得到的初始時(shí)間序列數(shù)據(jù)，通過曲線擬合和參數(shù)估計(jì)建立數(shù)學(xué)模型的理論和方法。ARIMA模型認(rèn)為依照時(shí)間順序進(jìn)行排列的所有觀測值之間都具有自相關(guān)性，這種自相關(guān)性延續(xù)了變量的發(fā)展趨勢。若將這種自相關(guān)性用定量的方法描述，就可以依據(jù)時(shí)間序列的過去值預(yù)測其將來值。該模型構(gòu)建思路如下。

（1）序列預(yù)處理：原始時(shí)間序列經(jīng)過預(yù)處理，可以被識(shí)別為平穩(wěn)非白噪聲序列，就能夠利用模型對該數(shù)據(jù)序列進(jìn)行建模。

（2）模型識(shí)別：由給出的序列樣本，根據(jù)樣本自相關(guān)系數(shù)和偏自相關(guān)系數(shù)的性質(zhì)，選擇階數(shù)適當(dāng)?shù)腁RMA（p， q）模型進(jìn)行擬合。

（3）參數(shù)估計(jì)：根據(jù)識(shí)別的模型及階數(shù)，對模型進(jìn)行參數(shù)估計(jì)和顯著性檢驗(yàn)。

（4）模型的構(gòu)建：記原始時(shí)間序列數(shù)據(jù)為at（t=1，2，3，...，n），記at的秩為Rt=R（at），考慮變量對（t，Rt）（t=1，2，3，...，n）的Spearman的相關(guān)系數(shù)qs，則有（14）數(shù)據(jù)。

（5）平穩(wěn)性檢驗(yàn)：假設(shè)H0為平穩(wěn)序列，H1為非平穩(wěn)序列。對于顯著水平a，由時(shí)間序列at計(jì)算（t，Rt）的Spearman秩相關(guān)系數(shù)qs，若T>ta2（n-2），則拒絕H0，認(rèn)為序列非平穩(wěn)。且當(dāng)qs>0時(shí)，認(rèn)為序列有下降趨勢。又當(dāng)T>ta2（n-2）時(shí)，接受H0，可以認(rèn)為序列是平穩(wěn)序列。

1.3組合預(yù)測模型

灰色模型是通過對原始數(shù)據(jù)加工處理來弱化隨機(jī)性的，若數(shù)據(jù)存在較大的波動(dòng)性，預(yù)測出來的結(jié)果可能會(huì)存在較大誤差。ARIMA模型對于預(yù)測的模型比較理想，要求時(shí)序數(shù)據(jù)是穩(wěn)定的，或者是進(jìn)行差分后是穩(wěn)定的，往往需要對數(shù)據(jù)進(jìn)行預(yù)處理，導(dǎo)致預(yù)測的精度不夠準(zhǔn)確。單一的預(yù)測模型均存在一定的局限性，很難把握預(yù)測結(jié)果的準(zhǔn)確性，從而使得誤差較大。組合模型可以通過對各個(gè)預(yù)測模型加權(quán)平均，把單一預(yù)測模型的優(yōu)點(diǎn)集合起來，形成預(yù)測精度更高的預(yù)測模型，可以大大減小誤差，使得結(jié)果更加準(zhǔn)確。

設(shè)原始序列的為X（0），預(yù)測序列為X^（0），觀測樣本數(shù)列為N，預(yù)測樣本數(shù)為k，組合模型為（15），其中ω1，ω2的公式為（16）[4]。

式中：X^（N+k）（0）表示組合模型的k個(gè)預(yù)測值，X^（N+k）（1）表示灰色模型的k個(gè)預(yù)測值，X^（n+k）（2）表示ARIMA模型的k個(gè)預(yù)測值。ω1，ω2分別表示灰色模型、ARIMA預(yù)測模型的權(quán)重，σ1，σ2分別表示灰色模型、ARIMA預(yù)測模型的標(biāo)準(zhǔn)差。

檢測結(jié)果表明，該組合預(yù)測方法優(yōu)于單一預(yù)測方法，由此可得最優(yōu)的組合預(yù)測模型。

2實(shí)證分析

選取2005～2019年的《新疆統(tǒng)計(jì)年鑒》中新疆維吾爾自治區(qū)棉花產(chǎn)量（單位為萬噸）年度數(shù)據(jù)為研究對象，對棉花產(chǎn)量進(jìn)行短期預(yù)測。

2.1 GM（1， 1）預(yù)測模型

按照GM（1， 1）的模型求解思路，建立傳統(tǒng)GM（1， 1）模型，得到時(shí)間響應(yīng)公式（17）。

利用Matlab軟件求解GM（1， 1）模型，將新疆歷年棉花產(chǎn)量作為訓(xùn)練數(shù)據(jù)，預(yù)測值與真實(shí)值的曲線擬合結(jié)果如圖1所示。

一般情況，如果絕對誤差百分比MAPE≤10%，同時(shí)遠(yuǎn)點(diǎn)誤差小于2%時(shí)可以認(rèn)定研究結(jié)果滿足精度的要求。由模型擬合結(jié)果可得：MAPE=0.0764SymbolcB@

0.1認(rèn)為模型精度誤差滿足要求。后驗(yàn)差比值C=0.2828<0.35，小誤差概率p=1>0.95，認(rèn)為精度檢驗(yàn)合格。綜上可知，本次GM（1， 1）灰色預(yù)測模型精度檢驗(yàn)合格。模型精度等級(jí)為1級(jí)，模型擬合結(jié)果良好，建立的模型所得的預(yù)測值誤差較小，其模型方程為（18）。將數(shù)據(jù)代入公式（18）中，可預(yù)測出棉花近幾年的總產(chǎn)量，具體預(yù)測數(shù)據(jù)見表2。

2.2時(shí)間序列預(yù)測模型

ARIMA模型要求序列是平穩(wěn)序列，因此要對數(shù)據(jù)進(jìn)行平穩(wěn)性分析。對樣本數(shù)據(jù)進(jìn)行ADF檢驗(yàn)，發(fā)現(xiàn)原始序列是非平穩(wěn)序列，經(jīng)二階差分后的P值<5%，在5%的顯著性水平下是平穩(wěn)序列。棉花產(chǎn)量序列經(jīng)二階差分后的自相關(guān)圖ACF和偏自相關(guān)圖PACF，如圖1所示。

由圖2也可以看出，二階差分后，數(shù)據(jù)逐步趨于平穩(wěn)，且滿足置信區(qū)間，故差分次數(shù)d=2。差分后的序列可以通過自相關(guān)圖和偏自相關(guān)圖初步確定p和q的值，由圖2可知，延遲5階后，ACF和PACF均落在二倍標(biāo)準(zhǔn)差范圍內(nèi)，因此該數(shù)據(jù)序列滿足ARIMA（5，2，1）模型。

2.3組合預(yù)測模型

對于上述兩種單一預(yù)測模型的預(yù)測精度不高、數(shù)據(jù)波動(dòng)較大等問題，若結(jié)合以上兩種單一預(yù)測模型的特征，構(gòu)建組合預(yù)測模型，經(jīng)計(jì)算可得，GM（1，1）模型與ARIMA模型的組合權(quán)重系數(shù)分別為0.46497和0.53503。預(yù)測值與真實(shí)值的具體數(shù)據(jù)見表2。

基于灰色預(yù)測模型與ARIMA模型得組合預(yù)測模型，較好地結(jié)合了兩種單一模型的優(yōu)點(diǎn)，能夠有效地?cái)M合非線性波動(dòng)及非線性增長趨勢。從表2可以看出，組合預(yù)測模型的平均相對誤差（MAPE）小于GM（1，1）模型與ARIMA模型這兩種單一預(yù)測模型。說明組合預(yù)測模型有較高的預(yù)測精度，取得很好的預(yù)測效果，預(yù)測精度較兩種單一預(yù)測模型有了很大的提高，實(shí)證分析結(jié)果表明該方法是有效和實(shí)用的，灰色預(yù)測模型與ARIMA模型的組合預(yù)測模型可以為相關(guān)從業(yè)者科學(xué)準(zhǔn)確地預(yù)測新疆棉花產(chǎn)量提供依據(jù)和實(shí)踐參考。

3結(jié)論

單獨(dú)采用灰色理論預(yù)測模型和時(shí)間序列模型進(jìn)行預(yù)測時(shí)，由于各個(gè)單一預(yù)測模型自身?xiàng)l件的限制，不能全面準(zhǔn)確地掌握原始數(shù)據(jù)的信息，影響預(yù)測結(jié)果的準(zhǔn)確性。而通過加權(quán)組合方式對兩個(gè)單一預(yù)測模型進(jìn)行組合，獲到的組合預(yù)測模型彌補(bǔ)了兩種單一預(yù)測模型自身的不足。本文基于灰色GM（1， 1）模型與時(shí)間序列ARIMA模型，通過加權(quán)組合方式將兩種單一模型很好地融合在一起，并選取2005～2019年的新疆維吾爾自治區(qū)棉花產(chǎn)量年度數(shù)據(jù)為研究對象，對棉花產(chǎn)量進(jìn)行了短期預(yù)測。實(shí)證分析結(jié)果表明組合模型預(yù)測的精度和穩(wěn)定性高于單項(xiàng)預(yù)測模型。因此，這種組合預(yù)測模型用于棉花產(chǎn)量預(yù)測是有效可行的。這為復(fù)雜性時(shí)間序列的建模和預(yù)測提供了切實(shí)可行的方法。

參考文獻(xiàn)

[1]李娟麗， 許英. 基于ARIMA模型的新疆棉花產(chǎn)量的預(yù)測分析[J]. 教育現(xiàn)代化，2018，5（17）：149-150.

[2]佟長福， 史海濱， 包小慶， 等. 基于小波分析理論組合模型的農(nóng)業(yè)需水量預(yù)測[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2011，27（05）：93-98.

[3]張利平， 于貞杰， 張建華， 等. 六種時(shí)間序列組合建模方案在衛(wèi)生費(fèi)用趨勢擬合中的應(yīng)用研究[J]. 中國衛(wèi)生經(jīng)濟(jì)， 2015，34（08）：56-58.

[4]邢棉. 季節(jié)性預(yù)測的組合灰色神經(jīng)網(wǎng)絡(luò)模型研究[J]. 系統(tǒng)工程理論與實(shí)踐，2001， 21（01）：31-35.

[5]朱輝.GM（1， 1）模型在新疆生產(chǎn)建設(shè)兵團(tuán)棉花總產(chǎn)量預(yù)測中的應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué)，2008（20）：8399-8400.

[6]陳黎明， 傅珊. 基于組合預(yù)測模型的GDP統(tǒng)計(jì)數(shù)據(jù)質(zhì)量評(píng)估研究[J]. 統(tǒng)計(jì)與決策，2013（08）：8-11.

[7]單銳， 王淑花， 高東蓮， 等. 基于時(shí)間序列模型與灰色模型的組合預(yù)測模型的研究[J]. 燕山大學(xué)學(xué)報(bào)，2012，36（01）：79-83.

[8]李曉彤. 基于小波分析的灰色模型與ARMA-GARCH模型的組合預(yù)測[D]. 濟(jì)南： 山東大學(xué)，2018.

[9]楊中旭. 棉花產(chǎn)量的Fuzzy預(yù)測法[J]. 江西棉花，1995，17（01）：18-21.

[10]鄧聚龍. 灰色系統(tǒng)理論教程[M]. 武漢： 華中理工大學(xué)出版社， 1990.