基于BP神經(jīng)網(wǎng)絡(luò)、ARIMA和LS-SVM模型的集成預(yù)測(cè)研究

李鵬飛 王青青 毋建宏 陳華雪

摘要：利用單一預(yù)測(cè)模型進(jìn)行產(chǎn)量預(yù)測(cè)時(shí)，由于各模型的數(shù)學(xué)原理不同，對(duì)同一數(shù)據(jù)的處理只能基于數(shù)據(jù)的部分特征，在預(yù)測(cè)時(shí)無(wú)法深度挖掘數(shù)據(jù)的潛在規(guī)律，易出現(xiàn)較大的預(yù)測(cè)偏差。基于以上問(wèn)題，筆者以1978—2017年陜西省蘋果產(chǎn)量數(shù)據(jù)為研究對(duì)象，將1978—2012年產(chǎn)量作為預(yù)測(cè)模型的訓(xùn)練數(shù)據(jù)，2013—2017年產(chǎn)量作為測(cè)試數(shù)據(jù)，選取BP神經(jīng)網(wǎng)絡(luò)、ARIMA、LS-SVM 3種在數(shù)學(xué)原理上具有明顯差異的預(yù)測(cè)模型，采用集成預(yù)測(cè)策略，依據(jù)3種預(yù)測(cè)模型對(duì)訓(xùn)練數(shù)據(jù)的平均相對(duì)預(yù)測(cè)誤差確定各模型的預(yù)測(cè)權(quán)重，最終對(duì)各模型預(yù)測(cè)的2013—2017年數(shù)據(jù)進(jìn)行加權(quán)獲取集成預(yù)測(cè)值。實(shí)證分析表明，集成預(yù)測(cè)值的平均相對(duì)誤差在2.5%以內(nèi)，其預(yù)測(cè)結(jié)果比單一預(yù)測(cè)模型更加準(zhǔn)確和穩(wěn)健，可有效實(shí)現(xiàn)蘋果產(chǎn)量高精度預(yù)測(cè)。

關(guān)鍵詞：BP神經(jīng)網(wǎng)絡(luò);ARIMA;LS-SVM;陜西蘋果產(chǎn)量;集成預(yù)測(cè)

中圖分類號(hào)： S11+9文獻(xiàn)標(biāo)志碼： A

文章編號(hào)：1002-1302（2020）04-0294-07

收稿日期：2019-07-30

基金項(xiàng)目：國(guó)家社會(huì)科學(xué)基金（編號(hào)：18FGL022）;教育部哲學(xué)社會(huì)科學(xué)研究后期資助項(xiàng)目（編號(hào)：18JHQ082）;陜西省科技廳重大項(xiàng)目（編號(hào)：2018ZDXM-GY-188）;陜西省社會(huì)科學(xué)界2019年重大理論與現(xiàn)實(shí)問(wèn)題研究項(xiàng)目（編號(hào)：2019TJ038）;陜西高校青年創(chuàng)新團(tuán)隊(duì);西安市科技計(jì)劃[編號(hào)：201806117YF05NC13（5）]。

作者簡(jiǎn)介：李鵬飛（1975—），男，陜西戶縣人，教授，碩士生導(dǎo)師，主要從事電商大數(shù)據(jù)分析研究。E-mail：lpf@xupt.edu.cn。

通信作者：毋建宏，副教授，主要從事大數(shù)據(jù)處理研究。E-mail：wujh@xupt.edu.cn。

陜西省因獨(dú)特的地理優(yōu)勢(shì)，已成為中國(guó)蘋果產(chǎn)區(qū)中唯一符合7項(xiàng)氣象指標(biāo)的蘋果生產(chǎn)最適宜區(qū)。蘋果產(chǎn)業(yè)是陜西省農(nóng)業(yè)特色的優(yōu)勢(shì)產(chǎn)業(yè)，不僅是陜西省農(nóng)業(yè)經(jīng)濟(jì)的重要支柱之一，也是農(nóng)民收入的重要來(lái)源[1]。2017年，陜西省蘋果產(chǎn)量達(dá)到1 153.94萬(wàn)t，約占中國(guó)蘋果總產(chǎn)量的1/4，世界蘋果總產(chǎn)量的1/7。蘋果產(chǎn)業(yè)已成為全國(guó)農(nóng)業(yè)結(jié)構(gòu)調(diào)整的典范，因而構(gòu)建符合產(chǎn)量變化趨勢(shì)的預(yù)測(cè)模型，合理準(zhǔn)確預(yù)測(cè)蘋果產(chǎn)量，對(duì)指導(dǎo)陜西省乃至全國(guó)蘋果產(chǎn)業(yè)具有重要的理論價(jià)值和實(shí)際意義。

隨著現(xiàn)代預(yù)測(cè)理論的發(fā)展，傳統(tǒng)的BP神經(jīng)網(wǎng)絡(luò)[2]、ARIMA模型[3]和支持向量機(jī)回歸模型[4]已難以滿足當(dāng)今研究需求。眾多學(xué)者將單一預(yù)測(cè)方法的組合預(yù)測(cè)模型作為研究重點(diǎn)，如ARIMA和BP神經(jīng)網(wǎng)絡(luò)組合[5];經(jīng)驗(yàn)?zāi)B(tài)分解、主成分分析和人工神經(jīng)網(wǎng)絡(luò)組合[6];季節(jié)調(diào)整法和BP神經(jīng)網(wǎng)絡(luò)組合[7];EGA算法和灰色Verhulst計(jì)量組合[8]等。也有部分學(xué)者基于模型預(yù)測(cè)性能的優(yōu)劣，將具有優(yōu)勢(shì)互補(bǔ)的單一模型進(jìn)行集成，如將定量與定性預(yù)測(cè)相結(jié)合的集成預(yù)測(cè)模型（SIF）[9];經(jīng)驗(yàn)?zāi)B(tài)分解和支持向量機(jī)集成[10];灰色GM（1，1）、ARIMA和VAR集成[11]等，集成預(yù)測(cè)模型克服了單一模型預(yù)測(cè)的缺陷，提高了預(yù)測(cè)精度。

為深度挖掘陜西省蘋果產(chǎn)量潛在的變化規(guī)律，同時(shí)多原理、多角度地預(yù)測(cè)蘋果產(chǎn)量，本研究選取BP神經(jīng)網(wǎng)絡(luò)、ARIMA和LS-SVM模型3種在數(shù)學(xué)原理上差異較大的模型作為預(yù)測(cè)模型，由于3種模型對(duì)中短期預(yù)測(cè)效果較好，將1978—2012年數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)，而2013—2017年5年數(shù)據(jù)作為測(cè)試數(shù)據(jù)，采用MATLAB軟件，對(duì)其進(jìn)行擬合和測(cè)試，然后根據(jù)3種預(yù)測(cè)模型訓(xùn)練值的相對(duì)誤差，計(jì)算出各模型的平均相對(duì)誤差，并按照權(quán)重計(jì)算公式確定單一預(yù)測(cè)模型的權(quán)重，最后將加權(quán)的結(jié)果作為最終的集成預(yù)測(cè)值。

1?BP神經(jīng)網(wǎng)絡(luò)、ARIMA、LS-SVM和集成預(yù)測(cè)模型理論

1.1?BP神經(jīng)網(wǎng)絡(luò)

反向傳播（back propagation，簡(jiǎn)稱BP）神經(jīng)網(wǎng)絡(luò)是由輸入層、隱含層、輸出層構(gòu)成的一種前饋型人工神經(jīng)網(wǎng)絡(luò)，通過(guò)約束網(wǎng)絡(luò)神經(jīng)元之間的聯(lián)系，實(shí)現(xiàn)模擬大腦神經(jīng)網(wǎng)絡(luò)智能處理信息的功能[12]，其網(wǎng)絡(luò)的架構(gòu)見圖1。

輸入層和輸出層與外界相連，隱含層承擔(dān)計(jì)算功能。BP神經(jīng)網(wǎng)絡(luò)每組輸入數(shù)據(jù)有對(duì)應(yīng)的目標(biāo)輸出，能夠?qū)崿F(xiàn)信息由輸入層至輸出層的前向傳播，并根據(jù)誤差平方和最小原則實(shí)現(xiàn)誤差的反向傳播，在此過(guò)程中，沿著梯度最速下降的方向[13]，調(diào)整不同層神經(jīng)元之間權(quán)值和閾值，使神經(jīng)網(wǎng)絡(luò)能夠根據(jù)輸入特征實(shí)現(xiàn)對(duì)輸出的逼近。

BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練流程：

1.1.1?正向傳播

設(shè)輸入數(shù)據(jù)為xi，輸入層到隱含層的權(quán)值為vij，隱含層個(gè)數(shù)為m，閾值為αk，激勵(lì)函數(shù)為f1，每個(gè)節(jié)點(diǎn)的輸出值γi，隱含層到輸出層的權(quán)值為vij，輸出層個(gè)數(shù)為n，閾值為αk，激勵(lì)函數(shù)為f2，輸出層每個(gè)節(jié)點(diǎn)的輸出值為yi，則隱含層中第i個(gè)神經(jīng)元的輸出：

γi=f1（∑mj=1xivij+αk）。（1）

輸出層第j個(gè)神經(jīng)元的輸出：

yi=f2（∑nj=1γivij+αk）。（2）

BP神經(jīng)網(wǎng)絡(luò)模型的隱含層和輸出層中神經(jīng)元之間采用的激活函數(shù)為sigmoid函數(shù)：

f（x）=11+e-x。（3）

1.1.2?反向傳播?根據(jù)梯度最速下降法調(diào)整連接系數(shù)：

ΔGij（k+1）=-ηBΔGij;（4）

ΔHij（k+1）=-ηBΔHij（k）。（5）

式中：G代表輸出層的權(quán)值和閾值，k表示第k步的學(xué)習(xí)過(guò)程，H代表隱含層的權(quán)值和閾值，學(xué)習(xí)速率為η，B表示網(wǎng)絡(luò)輸出值與目標(biāo)值z(mì)0的誤差，其誤差可以表示為B=∑m0，j=1（z0-yi）2/2。

在前向傳播過(guò)程中，訓(xùn)練信息從輸入層經(jīng)隱含層變換最終至輸出層，若輸出沒有達(dá)到期望目標(biāo)，則轉(zhuǎn)入誤差反向傳播，調(diào)整各神經(jīng)元間的聯(lián)系，循環(huán)重復(fù)上述過(guò)程直到誤差收斂。

1.2?ARIMA模型

自回歸積分滑動(dòng)平均模型（autoregressive integrated moving average model，簡(jiǎn)稱ARIMA），是對(duì)差分平穩(wěn)序列擬合的一種時(shí)間序列預(yù)測(cè)方法。ARIMA將時(shí)間序列表述為一個(gè)數(shù)學(xué)模型，該模型考察了時(shí)間序列的動(dòng)態(tài)和持續(xù)特征，揭示了時(shí)間序列內(nèi)在關(guān)系，適用于短期的時(shí)間序列預(yù)測(cè)[14]。

ARIMA（p，d，q）模型的一般形式：

ut=a+1ut-1+…+put-p+εt+θ1εt-1+…+θqεt-q。（6）

式中：p為自回歸模型的階數(shù)，d為差分次數(shù)，q為移動(dòng)平均階數(shù);ut為差分后的平穩(wěn)序列;a為常數(shù)，為自回歸模型系數(shù)，θ為移動(dòng)平均模型系數(shù)，εt為零均值白噪聲序列。

該方法建模步驟主要包括[15]：（1）時(shí)間序列平穩(wěn)性檢驗(yàn)。通過(guò)時(shí)序圖檢驗(yàn)時(shí)間序列的平穩(wěn)性，未通過(guò)平穩(wěn)性檢驗(yàn)即為非平穩(wěn)序列，采用差分運(yùn)算轉(zhuǎn)換為平穩(wěn)序列，并進(jìn)行平穩(wěn)性檢驗(yàn)。（2）模型初步識(shí)別。根據(jù)自相關(guān)系數(shù)圖和偏自相關(guān)系數(shù)圖，估計(jì)p、q的值，并利用AIC準(zhǔn)則確定模型的最佳階數(shù)。（3）參數(shù)估計(jì)和模型檢驗(yàn)。采用極大似然函數(shù)估計(jì)模型的參數(shù)，同時(shí)將估計(jì)的參數(shù)在模型中進(jìn)行檢驗(yàn)。（4）模型預(yù)測(cè)。根據(jù)已確定的最合適參數(shù)，進(jìn)行模型預(yù)測(cè)。

1.3?最小二乘支持向量機(jī)

20世紀(jì)90年代，俄羅斯數(shù)學(xué)家Vapnik提出了支持向量機(jī)（support vector machines，簡(jiǎn)稱SVM）的概念[16]。20世紀(jì)90年代末，Suykens等提出了最小二乘支持向量機(jī)（least squares support vector machine，簡(jiǎn)稱LSSVM）[17]，LSSVM是SVM的一種改進(jìn)算法，具有較強(qiáng)的泛化性，以結(jié)構(gòu)風(fēng)險(xiǎn)最小化為原則建立數(shù)據(jù)模型。該算法在很大程度上降低了樣本點(diǎn)在訓(xùn)練過(guò)程中的復(fù)雜度，優(yōu)于傳統(tǒng)支持向量機(jī)的運(yùn)算速度，是機(jī)器學(xué)習(xí)中應(yīng)用較為廣泛的一種建模方法。

LS-SVM模型將最小二乘線性系統(tǒng)作為損失函數(shù)，將凸二次規(guī)劃問(wèn)題轉(zhuǎn)化成線性方程組求解問(wèn)題，簡(jiǎn)化了計(jì)算的復(fù)雜性，其回歸過(guò)程可以表示為[18]

對(duì)給定的數(shù)據(jù)集S={（x1，y1），（x2，y2），…，（xN，yN），其線性回歸函數(shù)g（x）=wTφ（x）+k，其中xi表示第i個(gè)輸入向量，yi為與xi相對(duì)應(yīng)的輸出值，i=1，2，…，N，N為樣本容量，且w∈Rn，k∈R，φ（·）為解決非線性問(wèn)題的核函數(shù)。

回歸問(wèn)題對(duì)應(yīng)的優(yōu)化問(wèn)題目標(biāo)函數(shù)為公式（7），與之對(duì)應(yīng)的約束條件為公式（8）：

minw，k，eQ（w，k，e）=‖w‖2+γ∑Ni=1e2i2;（7）

yi=wTφ（xi）+k+ei。（8）

公式（7）中：w表示向量，k表示標(biāo)量，e表示誤差，γ表示正則化參數(shù)，ei是誤差變量。

引入一個(gè)拉格朗日乘子αi則LS-SVM優(yōu)化問(wèn)題對(duì)應(yīng)的拉格朗日函數(shù)為

L（w，k，e，α）=Q（w，k，e）-∑Ni=1αi[wTφ（xi）+k+ei-yi]。（9）

對(duì)公式（9）中每個(gè)參數(shù)求偏導(dǎo)并令其值為0，得到：

Lw=0w=∑Ni=1αiφ（xi）

Lk=0∑Ni=1αi=0

Lei=0αi=γei

Lαi=0wTφ（xi）+k+ei-yi=0。（10）

將公式（10）消去w和ei可以得到：

0ETED+I/γkα=0Y。（11）

式中：E=[1，1，…，1]，D為核函數(shù)矩陣，I為單位矩陣，α=[α1，α2，…，αN]，Y=[y1，y2，…，yN]。

令B=D+I/γ，可得到：

k=ETB-1yETB-1E

α=B-1（Y-kE）。（12）

因此，LS-SVM的預(yù)測(cè)函數(shù)為

y^（x）=∑Ni=1αiK（x，xi）+k。（13）

式中：K（x，xi）=exp-（x-xi）22σ2為徑向基核函數(shù)，其中σ為核參數(shù)。

LS-SVM模型常采用的核函數(shù)有高斯徑向基核函數(shù)（radial basis function，簡(jiǎn)稱RBF）、線性核函數(shù)和多項(xiàng)式核函數(shù)，本研究選用包含正則化參數(shù)γ、核參數(shù)σ的RBF核函數(shù)，這2個(gè)優(yōu)化參數(shù)對(duì)模型的泛化能力和預(yù)測(cè)精度影響極大[19]。

1.4?集成預(yù)測(cè)原理

集成預(yù)測(cè)模型可以克服單一模型的缺陷，優(yōu)化由于模型設(shè)定不合理而導(dǎo)致的預(yù)測(cè)偏誤，常用的有簡(jiǎn)單平均法和加權(quán)平均法[11]。筆者采用加權(quán)平均法對(duì)單一模型進(jìn)行集成，即根據(jù)3個(gè)模型訓(xùn)練數(shù)據(jù)的平均相對(duì)預(yù)測(cè)誤差確定各單一預(yù)測(cè)模型的權(quán)重，將預(yù)測(cè)值的加權(quán)結(jié)果作為最終的集成預(yù)測(cè)值。其單模型權(quán)重計(jì)算公式如下：

wj=S-1j/∑mj=1S-1j。（14）

式中：wj表示第j個(gè)模型的權(quán)重，Sj表示第j個(gè)模型的平均相對(duì)誤差的絕對(duì)值，m表示模型的數(shù)量。

2?結(jié)果與分析

2.1?數(shù)據(jù)描述

本研究數(shù)據(jù)來(lái)源于《中華人民共和國(guó)國(guó)家統(tǒng)計(jì)局》《陜西統(tǒng)計(jì)年鑒2018》，其詳細(xì)數(shù)據(jù)見表1（單位：萬(wàn)t），通過(guò)分析40年的數(shù)據(jù)，可看出陜西省蘋果產(chǎn)量1978—2017年整體保持增長(zhǎng)趨勢(shì)。40年產(chǎn)量數(shù)據(jù)的序列走勢(shì)見圖2。

2.2?1978—2017年蘋果產(chǎn)量數(shù)據(jù)的擬合與預(yù)測(cè)

2.2.1?基于BP神經(jīng)網(wǎng)絡(luò)模型的擬合與預(yù)測(cè)

采用BP神經(jīng)網(wǎng)絡(luò)對(duì)陜西省蘋果產(chǎn)量預(yù)測(cè)時(shí)，為提高預(yù)測(cè)結(jié)果的準(zhǔn)確性，將1978—2012年共35年蘋果產(chǎn)量數(shù)據(jù)作為輸入樣本和測(cè)試樣本數(shù)量?通過(guò)多次訓(xùn)練，最終將BP神經(jīng)網(wǎng)絡(luò)的隱含層設(shè)置為3層，輸入和輸出層各為1層。在運(yùn)算過(guò)程中對(duì)輸入層與隱含層和輸出層與隱含層之間的權(quán)值和閾值進(jìn)行多次修正，訓(xùn)練得到了預(yù)測(cè)精度高而且符合產(chǎn)量變化趨勢(shì)的網(wǎng)絡(luò)。

采用已訓(xùn)練好的網(wǎng)絡(luò)對(duì)1978—2013年蘋果產(chǎn)量進(jìn)行擬合（圖3-a），通過(guò)擬合效果可以看出訓(xùn)練的BP神經(jīng)網(wǎng)絡(luò)的輸出值和蘋果產(chǎn)量真實(shí)值基本吻合，可以用已訓(xùn)練好的網(wǎng)絡(luò)對(duì)2013—2017陜西蘋果產(chǎn)量進(jìn)行預(yù)測(cè)（圖3-b）。

2.2.2?基于ARIMA模型的擬合與預(yù)測(cè)

1978—2017年陜西省蘋果產(chǎn)量序列呈上升趨勢(shì)，可推斷該序列為非平穩(wěn)序列，并且根據(jù)該序列的自相關(guān)系數(shù)圖和偏自相關(guān)系數(shù)圖進(jìn)一步檢驗(yàn)了序列的平穩(wěn)性，因此判斷陜西省蘋果產(chǎn)量為非平穩(wěn)序列。對(duì)該序列進(jìn)行一階差分處理，一階差分處理后的自相關(guān)和偏自相關(guān)系數(shù)見圖4，可以看出一階差分后序列平穩(wěn)，且ADF檢驗(yàn)結(jié)果表明，一階差分后序列平穩(wěn)，因此可對(duì)陜西省蘋果產(chǎn)量序列構(gòu)建ARIMA（p，d，q）模型。

根據(jù)一階差分后的自相關(guān)系數(shù)和偏自相關(guān)系數(shù)的分布狀態(tài)以及AIC準(zhǔn)則對(duì)模型定階，通過(guò)反復(fù)檢驗(yàn)，最終將模型確定為ARIMA（1，1，0）。運(yùn)用MATLAB軟件對(duì)ARIMA（1，1，0）進(jìn)行建模預(yù)測(cè)，圖5-a為ARIMA模型對(duì)1978—2012蘋果產(chǎn)量擬合結(jié)果，圖5-b為該模型對(duì)2013—2017蘋果產(chǎn)量預(yù)測(cè)結(jié)果。

2.2.3?基于LS-SVM模型的擬合與預(yù)測(cè)

本研究利用最小二乘支持向量機(jī)對(duì)2013—2017年陜西省蘋果產(chǎn)量進(jìn)行預(yù)測(cè)，借助MATLAB軟件的LS-SVM工具箱進(jìn)行輔助建模（工具箱下載地址為：https：//www.esat.kuleuven.be/sista/lssvmlab/）。選用徑向基核函數(shù) RBF代替高維空間向量的內(nèi)積計(jì)算，回歸函數(shù)設(shè)為function estimation，用十倍交叉驗(yàn)證方法來(lái)不斷調(diào)整參數(shù)γ（正則參數(shù)，取決于對(duì)訓(xùn)練誤差最小化和估計(jì)函數(shù)平滑的權(quán)衡）和σ（內(nèi)核函數(shù)的參數(shù)），根據(jù)擬合值和預(yù)測(cè)效果，將參數(shù)γ確定為10，σ確定為5，使用函數(shù)simlssvm評(píng)估模型中新的點(diǎn)，最后得到基于蘋果產(chǎn)量數(shù)據(jù)的預(yù)測(cè)模型和預(yù)測(cè)值。圖6-a顯示了LS-SVM模型對(duì)1978—2012年陜西省蘋果產(chǎn)量的擬合效果，從圖中可以看出，LS-SVM對(duì)產(chǎn)量的擬合趨勢(shì)和真實(shí)值的變化趨勢(shì)基本吻合，即該模型可以用來(lái)預(yù)測(cè)蘋果產(chǎn)量，圖6-b為采用該模型預(yù)測(cè)的2013—2017年蘋果產(chǎn)量數(shù)值。

2.2.4?模型預(yù)測(cè)結(jié)果

BP神經(jīng)網(wǎng)絡(luò)、ARIMA、LS-SVM 模型對(duì)1978—2012年陜西省蘋果產(chǎn)量數(shù)據(jù)擬合值的平均相對(duì)誤差見表2，可以看出，BP神經(jīng)網(wǎng)絡(luò)對(duì)蘋果產(chǎn)量的擬合效果優(yōu)于ARIMA模型和LS-SVM模型。

3種模型對(duì)2013—2017年蘋果產(chǎn)量的預(yù)測(cè)值以及相對(duì)預(yù)測(cè)誤差見表3，可以看出，BP神經(jīng)網(wǎng)絡(luò)的相對(duì)誤差最大值為5.26，最小值為1.64，且通過(guò)計(jì)算該模型相對(duì)誤差的平均值為3.84;ARIMA模型的相對(duì)誤差最大值為6.25，最小值為1.31，其相對(duì)誤差的平均值為3.51;LS-SVM模型相對(duì)誤差最大值為9.05，最小值為1.73，相對(duì)誤差的平均值為5.02。

綜上所述，3種單一預(yù)測(cè)模型存在預(yù)測(cè)誤差波動(dòng)范圍比較大的問(wèn)題，無(wú)法更準(zhǔn)確地預(yù)測(cè)蘋果產(chǎn)量，因此將各模型訓(xùn)練數(shù)據(jù)擬合值的平均相對(duì)誤差值作為計(jì)算權(quán)重的基礎(chǔ)，根據(jù)權(quán)重計(jì)算公式，確定各單一模型的權(quán)重，再利用2013—2017年各模型預(yù)測(cè)的蘋果產(chǎn)量數(shù)值進(jìn)行加權(quán)?其加權(quán)結(jié)果作為最終的集成預(yù)測(cè)值。

2.3?集成預(yù)測(cè)結(jié)果

由于單一的BP神經(jīng)網(wǎng)絡(luò)、ARIMA和LS-SVM預(yù)測(cè)模型對(duì)陜西省蘋果產(chǎn)量2013—2017年的預(yù)測(cè)效果不佳，為提高預(yù)測(cè)精度，現(xiàn)將各單一預(yù)測(cè)模型的預(yù)測(cè)值進(jìn)行集成。首先根據(jù)模型訓(xùn)練數(shù)據(jù)擬合值的相對(duì)誤差計(jì)算出平均相對(duì)誤差;然后根據(jù)公式（14），計(jì)算出單一模型的權(quán)重（表4）;最后根據(jù)2013—2017年各模型的預(yù)測(cè)值以及相對(duì)應(yīng)的權(quán)重進(jìn)行加權(quán)，將所求得的結(jié)果作為最終的集成預(yù)測(cè)值。

利用表4預(yù)測(cè)模型的權(quán)重進(jìn)行加權(quán)計(jì)算，2013—2017年陜西省蘋果產(chǎn)量的集成預(yù)測(cè)值以及集成預(yù)測(cè)值相對(duì)誤差和平均相對(duì)誤差見表5。從表5可看出，集成預(yù)測(cè)的平均相對(duì)誤差為2.21%，小于各單一模型。

通過(guò)集成預(yù)測(cè)模型的預(yù)測(cè)值可以看出，集成后的效果更穩(wěn)定，預(yù)測(cè)誤差相比于單一各模型波動(dòng)較小，平均相對(duì)誤差在2.5%以內(nèi)，預(yù)測(cè)效果較好，該集成預(yù)測(cè)方法避免了單一預(yù)測(cè)模型預(yù)測(cè)的局限性，提高了預(yù)測(cè)精度，可用于今后蘋果產(chǎn)量的預(yù)測(cè)。

3?結(jié)論

本研究基于1978—2017年陜西省蘋果產(chǎn)量數(shù)據(jù)，分別利用BP神經(jīng)網(wǎng)絡(luò)、ARIMA、LS-SVM預(yù)測(cè)模型對(duì)蘋果產(chǎn)量進(jìn)行預(yù)測(cè)。結(jié)果表明，BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值的平均相對(duì)誤差較ARIMA和LS-SVM模型更小，但3種預(yù)測(cè)模型仍存在預(yù)測(cè)誤差波動(dòng)比較大的問(wèn)題，無(wú)法更準(zhǔn)確地預(yù)測(cè)蘋果產(chǎn)量。為規(guī)避單一預(yù)測(cè)模型的缺陷，使得預(yù)測(cè)結(jié)果更加準(zhǔn)確，本研究基于加權(quán)平均的思想，將3種模型的預(yù)測(cè)值依據(jù)權(quán)重的不同進(jìn)行集成，實(shí)證分析表明，集成預(yù)測(cè)的結(jié)果更準(zhǔn)確，預(yù)測(cè)性能更穩(wěn)健，相對(duì)誤差的波動(dòng)性更小，且平均相對(duì)誤差在2.5%以內(nèi)，可有效地實(shí)現(xiàn)蘋果產(chǎn)量的高精度預(yù)測(cè)。

參考文獻(xiàn)：

[1]李艷花，薛江偉. 基于SWOT分析的陜西蘋果國(guó)際競(jìng)爭(zhēng)力提升策略[J]. 干旱區(qū)資源與環(huán)境，2008，22（9）：30-36.

[2]王衛(wèi)東，李?凈，張福存，等. 基于BP神經(jīng)網(wǎng)絡(luò)的太陽(yáng)輻射預(yù)測(cè)——以蘭州市為例[J]. 干旱區(qū)資源與環(huán)境，2014，28（2）：185-189.

[3]吳玉霞，溫?欣. 基于ARIMA模型的短期股票價(jià)格預(yù)測(cè)[J]. 統(tǒng)計(jì)與決策，2016（23）：83-86.

[4]隗海林，包翠竹，李洪雪，等. 基于最小二乘支持向量機(jī)的怠速時(shí)間預(yù)測(cè)[J]. 吉林大學(xué)學(xué)報(bào)（工學(xué)版），2018，48（5）：1360-1365.

[5]翟?靜，曹?俊. 基于時(shí)間序列ARIMA與BP神經(jīng)網(wǎng)絡(luò)的組合預(yù)測(cè)模型[J]. 統(tǒng)計(jì)與決策，2016（4）：29-32.

[6]潘和平，張承釗. FEPA-金融時(shí)間序列自適應(yīng)組合預(yù)測(cè)模型[J]. 中國(guó)管理科學(xué)，2018，26（6）：26-38.

[7]龍?勇，蘇振宇，汪?於. 基于季節(jié)調(diào)整和BP神經(jīng)網(wǎng)絡(luò)的月度負(fù)荷預(yù)測(cè)[J]. 系統(tǒng)工程理論與實(shí)踐，2018，38（4）：1052-1060.

[8]張?侃，劉寶平，黃?棟. 基于EGA算法的小樣本非線性殘差灰色Verhulst計(jì)量組合預(yù)測(cè)模型[J]. 系統(tǒng)工程理論與實(shí)踐，2017，37（10）：2630-2639.

[9]金?淳，曹?迪，王?聰，等. 汽車零部件第三方物流倉(cāng)儲(chǔ)需求量集成預(yù)測(cè)模型[J]. 系統(tǒng)管理學(xué)報(bào)，2018，27（6）：1157-1165.

[10]蔡超敏，凌立文，牛?超，等. 國(guó)內(nèi)豬肉市場(chǎng)價(jià)格的EMD-SVM集成預(yù)測(cè)模型[J]. 中國(guó)管理科學(xué)，2016，24（增刊1）：845-851.

[11]李?蕾，李?超，丁雪辰，等. 我國(guó)執(zhí)業(yè)（助理）醫(yī)師需求集成預(yù)測(cè)——基于GM、ARIMA和VAR模型的實(shí)證研究[J]. 管理評(píng)論，2018，30（3）：171-178.

[12]張莉曼，張向先，李中梅，等. 基于BP神經(jīng)網(wǎng)絡(luò)的智庫(kù)微信公眾平臺(tái)信息傳播力評(píng)價(jià)研究[J]. 情報(bào)理論與實(shí)踐，2018，41（10）：93-99.

[13]周開利，康耀紅. 神經(jīng)網(wǎng)絡(luò)模型及其MATLAB仿真程序設(shè)計(jì)[M]. 北京：清華大學(xué)出版社，2005.

[14]潘?靜，張?穎，劉?璐. 基于ARIMA模型與GM（1，1）模型的居民消費(fèi)價(jià)格指數(shù)預(yù)測(cè)對(duì)比分析[J]. 統(tǒng)計(jì)與決策，2017（20）：110-112.

[15]陳?榮，梁昌勇，陸文星，等. 面向旅游突發(fā)事件的客流量混合預(yù)測(cè)方法研究[J]. 中國(guó)管理科學(xué)，2017，25（5）：167-174.

[16]Vapnik V N. The nature of statistical learning theory[M]. New York：Springer，2000：133-155.

[17]Suykens J，Vandewalle J. Least squares support vector machine classifiers[J]. Neural Processing Letters，1999，9（3）：293-300.

[18]張?慧，盧文冰，趙雄文，等. 基于最小二乘支持向量機(jī)和小波神經(jīng)網(wǎng)絡(luò)的電力線通信信道噪聲建模研究[J]. 電工技術(shù)學(xué)報(bào)，2018，33（16）：3879-3888.

[19]Suykens J K，Brabanter J D，Lukas L，et al. Weighted least squares support vector machines：robustness and sparse approximation[J]. Neurocomputing，2002，48（1/2/3/4）：85-105.