基于SARIMA-改進(jìn)RS-多步LSTM的農(nóng)機(jī)裝備需求組合預(yù)測(cè)模型

摘要：針對(duì)農(nóng)機(jī)裝備需求受實(shí)際農(nóng)業(yè)生產(chǎn)和其他多種因素影響，需求數(shù)據(jù)呈現(xiàn)周期性、非線性等特點(diǎn)，難以準(zhǔn)確預(yù)測(cè)農(nóng)機(jī)需求的問(wèn)題，提出一種集成SARIMA-改進(jìn)RS-多步LSTM的農(nóng)機(jī)裝備需求預(yù)測(cè)方法。通過(guò)確定參數(shù)組合，構(gòu)建季節(jié)性差分自回歸滑動(dòng)平均（SARIMA）模型。引入完全集合經(jīng)驗(yàn)?zāi)B(tài)分解（CEEMDAN）、改進(jìn)隨機(jī)搜索（RS）算法和多步長(zhǎng)短期記憶網(wǎng)絡(luò)（LSTM），構(gòu)建改進(jìn)RS-多步LSTM模型。將SARIMA模型和改進(jìn)RS-多步LSTM模型進(jìn)行最優(yōu)加權(quán)組合，得到組合預(yù)測(cè)模型。以某型號(hào)農(nóng)機(jī)裝備進(jìn)行實(shí)例驗(yàn)證，結(jié)果表明，所提方法能夠有效預(yù)測(cè)農(nóng)機(jī)裝備需求的時(shí)間序列，評(píng)價(jià)指標(biāo)MSE、MAE、R²分別為225.45、13.22和0.920 9。

關(guān)鍵詞：農(nóng)機(jī)裝備；需求預(yù)測(cè)；季節(jié)性；多步LSTM；模態(tài)分解

中圖分類(lèi)號(hào)：F323

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

文章編號(hào)：2095-5553（2024）12-0305-07收稿日期：2024年1月30日

修回日期：2024年4月15日

*基金項(xiàng)目：國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2020YFB1713500）

第一作者：呂鋒，男，1980年生，河南商丘人，博士，副教授；研究方向?yàn)楫a(chǎn)品服務(wù)系統(tǒng)設(shè)計(jì)。E-mail： lvfeng1980@haust.edu.cn

Combination prediction model of agricultural machinery equipment demand based on SARIMA-improved RS-multistep LSTM

Lü Feng¹， Wang Baosen¹， Chu Huili¹， Yang Cheng¹， Lü Ling²

（1. School of Mechatronics Engineering， Henan University of Science and Technology， Luoyang， 471003， China;2. Highway Business Development Centre， Liangyuan District， Shangqiu， 476000， China）

Abstract： In view of the fact that the demand for agricultural machinery equipment is affected by actual agricultural production and many other factors， and the demand data is cyclical and non-linear， making it difficult to accurately forecast the demand for agricultural machinery， a material demand forecasting method integrating SARIMA-improved RS-multistep LSTM was proposed. The seasonal differential autoregressive moving average （SARIMA） model was constructed by determining the parameter combination. The complete set Empirical Mode decomposition （CEEMDAN）， improved random search （RS） algorithm and multistep short term memory network （LSTM） were introduced to construct an improved RS-Multistep LSTM model. The optimal weighted combination of SARIMA model and improved RS multistep LSTM model was used to obtain a combined prediction model. Using a certain model of agricultural machinery equipment as an example for verification. The results have showed that the proposed method can effectively predict the time series of the demand for agricultural machinery equipment， the evaluation indicators MSE， MAE and R² are 225.45， 13.22 and 0.920 9 respectively.

Keywords： agricultural machinery equipment; demand prediction; seasonal; multistep LSTM; mode decomposition

0 引言

作為農(nóng)業(yè)生產(chǎn)的關(guān)鍵性投入要素，農(nóng)機(jī)裝備是現(xiàn)代農(nóng)業(yè)發(fā)展的重要支撐。隨著我國(guó)農(nóng)業(yè)發(fā)展方向和農(nóng)業(yè)生產(chǎn)方式轉(zhuǎn)型，農(nóng)機(jī)裝備需求也變得多樣化。為了及時(shí)滿(mǎn)足市場(chǎng)需求，農(nóng)機(jī)企業(yè)如何準(zhǔn)確預(yù)測(cè)農(nóng)機(jī)裝備需求變得至關(guān)重要，并將成為企業(yè)合理安排生產(chǎn)的依據(jù)。

預(yù)測(cè)模型構(gòu)建是準(zhǔn)確預(yù)測(cè)農(nóng)機(jī)裝備需求的關(guān)鍵環(huán)節(jié)，常用的需求預(yù)測(cè)模型主要包括：回歸分析預(yù)測(cè)模型^［1］，建立回歸方程進(jìn)行預(yù)測(cè)；機(jī)器學(xué)習(xí)預(yù)測(cè)模型^{［2， 3］}，通過(guò)對(duì)歷史數(shù)據(jù)的學(xué)習(xí)，自動(dòng)發(fā)現(xiàn)數(shù)據(jù)中的規(guī)律和模式，做出預(yù)測(cè)；深度學(xué)習(xí)預(yù)測(cè)模型^［4］，用于處理復(fù)雜的時(shí)間序列數(shù)據(jù)，通過(guò)從原始數(shù)據(jù)中學(xué)習(xí)并提取有用的特征進(jìn)行預(yù)測(cè)。

眾多學(xué)者對(duì)此進(jìn)行大量的研究，取得了豐碩的研究成果。武慧榮等^［5］構(gòu)建GM（1，1）-MLP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行糧食產(chǎn)量預(yù)測(cè)。崔凱等^［6］提出了基于Group-BiLSTM-LightGBM的集成預(yù)測(cè)方法。王海程等^［7］研究構(gòu)建TCG-ODE模型對(duì)不同區(qū)域間的出租車(chē)需求量進(jìn)行精準(zhǔn)預(yù)測(cè)。黃國(guó)興等^［8］研究隨機(jī)森林回歸模型用于對(duì)艦船零部件備件的需求預(yù)測(cè)。付維方等^［9］提出自適應(yīng)預(yù)測(cè)框架，實(shí)現(xiàn)不同航材需求的多預(yù)測(cè)決策及同一航材多階段的動(dòng)態(tài)預(yù)測(cè)。張舜等^［10］提出PCA-SVR方法對(duì)農(nóng)機(jī)備件需求進(jìn)行預(yù)測(cè)。李瓊等^［11］考慮歷史需求數(shù)據(jù)在時(shí)間序列上的相鄰關(guān)聯(lián)性，研究GRU-BP組合神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)未來(lái)一段時(shí)間內(nèi)的產(chǎn)品需求量。Chen等^［12］針對(duì)通用飛機(jī)材料需求，建立PSO-BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

農(nóng)機(jī)裝備需求受多種因素的影響，需求數(shù)據(jù)呈現(xiàn)周期性、非線性等特點(diǎn)，增加了預(yù)測(cè)的難度。傳統(tǒng)的預(yù)測(cè)方法無(wú)法有效地捕捉和反饋外界影響因素，在處理具有非線性數(shù)據(jù)時(shí)效果有限。SARIMA模型^［13］具有統(tǒng)計(jì)性質(zhì)，有較好的穩(wěn)定性，能夠捕捉時(shí)間序列中的季節(jié)性和趨勢(shì)性規(guī)律。改進(jìn)RS-多步LSTM模型^{［14， 15］}具有動(dòng)態(tài)適應(yīng)性，學(xué)習(xí)能力較強(qiáng)，可以綜合考慮多個(gè)影響因素對(duì)農(nóng)機(jī)裝備需求的影響，在處理長(zhǎng)期依賴(lài)關(guān)系和非線性關(guān)系時(shí)表現(xiàn)出色?；诖?，通過(guò)構(gòu)建SARIMA-改進(jìn)RS-多步LSTM組合模型，利用SARIMA模型的穩(wěn)定性來(lái)約束改進(jìn)RS-多步LSTM模型的動(dòng)態(tài)行為，避免出現(xiàn)過(guò)度擬合或數(shù)據(jù)過(guò)適應(yīng)的問(wèn)題，以提高農(nóng)機(jī)裝備需求預(yù)測(cè)的準(zhǔn)確性。

1 農(nóng)機(jī)裝備需求影響因素辨識(shí)

馮建英等^［16］通過(guò)研究發(fā)現(xiàn)影響農(nóng)戶(hù)購(gòu)買(mǎi)農(nóng)機(jī)的因素主要有受教育水平、家庭收入、政策因素等。陳旭等^［17］構(gòu)建向后逐步回歸模型，對(duì)農(nóng)機(jī)需求特征及影響因素進(jìn)行研究，發(fā)現(xiàn)農(nóng)戶(hù)人均純收入、農(nóng)業(yè)勞動(dòng)力等對(duì)農(nóng)機(jī)需求的影響統(tǒng)計(jì)上顯著。張標(biāo)等^［18］研究發(fā)現(xiàn)家庭年收入、農(nóng)業(yè)收入占比、示范戶(hù)等因素對(duì)農(nóng)戶(hù)購(gòu)買(mǎi)農(nóng)機(jī)具有顯著影響。參考現(xiàn)有研究成果，結(jié)合數(shù)據(jù)可得，提煉出主要農(nóng)作物價(jià)格、燃料價(jià)格、第一產(chǎn)業(yè)從業(yè)人員數(shù)量、農(nóng)村居民人均可支配收入4個(gè)影響因素預(yù)測(cè)農(nóng)機(jī)裝備需求。其中，主要農(nóng)作物價(jià)格和燃料價(jià)格取上月月初和月末的平均價(jià)格。第一產(chǎn)業(yè)從業(yè)人員數(shù)量為每年年末統(tǒng)計(jì)數(shù)據(jù)，農(nóng)村居民人均可支配收入為每年季度統(tǒng)計(jì)數(shù)據(jù)，均采用線性插值法充填為月度數(shù)據(jù)。數(shù)據(jù)來(lái)源為國(guó)家統(tǒng)計(jì)局、農(nóng)機(jī)360網(wǎng)、國(guó)家糧食和物資儲(chǔ)備局、上海石油天然氣交易中心。

2 基于SARIMA-改進(jìn)RS-多步LSTM的組合模型構(gòu)建

2.1 SARIMA模型構(gòu)建

SARIMA（p，d，q）（P，D，Q，s）模型適用于具有季節(jié)性時(shí)間序列的建模和預(yù)測(cè)。構(gòu)建SARIMA模型核心在于相關(guān)參數(shù)（p，d，q）（P，D，Q，s）的確定，這些參數(shù)決定了模型的預(yù)測(cè)性能。（1）對(duì)農(nóng)機(jī)裝備需求原始數(shù)據(jù)進(jìn)行時(shí)序分解，將數(shù)據(jù)分解為趨勢(shì)、季節(jié)性和殘差3個(gè)部分。（2）通過(guò)季節(jié)性曲線確定參數(shù)s。（3）通過(guò)趨勢(shì)曲線判斷原始數(shù)據(jù)是否需要差分，通過(guò)判斷差分次數(shù)確定參數(shù)d。（4）通過(guò)自相關(guān)圖（ACF）和偏自相關(guān)圖（PACF）確定參數(shù)p和參數(shù)q。（5）通過(guò)構(gòu)建嵌套循環(huán)算法，確定參數(shù)P、D、Q。遍歷p，d，q和P，D，Q，s中所有可能的參數(shù)組合，輸出不同參數(shù)組合下的AIC值（赤池信息量準(zhǔn)則）^［19］，比較AIC值，選擇最小AIC值對(duì)應(yīng)的組合作為模型的參數(shù)，構(gòu)建SARIMA模型。

2.2 改進(jìn)RS-多步LSTM模型構(gòu)建

2.2.1 CEEMDAN處理

農(nóng)機(jī)裝備需求數(shù)據(jù)包含周期性和非線性特性。CEEMDAN作為一種非線性信號(hào)處理方法，能夠有效地分解復(fù)雜信號(hào)，提取數(shù)據(jù)中的固有模式和特征^［20］。本文引入CEEMDAN對(duì)需求數(shù)據(jù)進(jìn)行處理。

CEEMDAN添加自適應(yīng)白噪聲序列

H_i（t）=H（t）+α×ω（t） （1）

式中：H（t）——原始時(shí)間序列；

H_i（t）——添加噪聲信號(hào)；

α——信噪比系數(shù)；

ω（t）——服從正態(tài)分布的白噪聲序列。

對(duì)H（t）進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解得到各個(gè)樣本的一階模態(tài)分量IMF₁（t）和余項(xiàng)r₁（t），對(duì)余項(xiàng)不斷分解得到各階的IMF（t）分量和余項(xiàng)r（t）。最終原始時(shí)間序列數(shù)據(jù)由一系列IMF（t）分量和余項(xiàng)r（t）構(gòu)成。

對(duì)農(nóng)機(jī)裝備需求、SARIMA殘差和影響因素等數(shù)據(jù)進(jìn)行CEEMDAN處理，每組數(shù)據(jù)經(jīng)過(guò)分解得到K_n個(gè)模態(tài)分量（IMF）和殘差（res）。其中，數(shù)據(jù)組數(shù)量為M，所有分解的IMF和res構(gòu)成矩陣X（t）。

提取具有周期性的IMF，與res合并為矩陣X′（t），作為預(yù)測(cè)模型的輸入端。

2.2.2 改進(jìn)RS-多步LSTM模型

構(gòu)建改進(jìn)RS-多步LSTM模型進(jìn)行未來(lái)多個(gè)時(shí)間步的預(yù)測(cè)。多步LSTM預(yù)測(cè)是將t時(shí)刻前d個(gè)歷史數(shù)據(jù)處理后的特征成分作為輸入，在一次預(yù)測(cè)中輸出L個(gè)農(nóng)機(jī)裝備需求預(yù)測(cè)值。輸入序列為［x_t-d+1，…，x_t-1，x_t］，輸出序列［Y_t+1，…，Y_t+L］為未來(lái)多個(gè)時(shí)間步的預(yù)測(cè)值。隨著預(yù)測(cè)窗口向右平移，每次預(yù)測(cè)輸出未來(lái)L個(gè)農(nóng)機(jī)裝備需求預(yù)測(cè)值，如圖1所示。多步LSTM模型的結(jié)構(gòu)包括輸入層、隱藏層和輸出層，如圖2所示。

輸入門(mén)

i_t=σ（W_ix_t+W_ih_t-1+b_i） （4）

遺忘門(mén)

f_t=σ（W_fx_t+W_fh_t-1+b_f） （5）

輸出門(mén)

o_t=σ（W_ox_t+W_oh_t-1+b_o） （6）

細(xì)胞狀態(tài)

c_t=f_t·c_t-1+i_t·tanh（W_cx_t+W_ch_t-1+b_c） （7）

輸出

h_t=o_t·tanh（c_t） （8）

式中：x_t——t時(shí)刻輸入向量；

h_t——t時(shí)刻輸出向量；

c_t——LSTM模型在t時(shí)刻的細(xì)胞狀態(tài)；

W_i、W_f、W_o、W_c——輸入門(mén)、遺忘門(mén)、輸出門(mén)、候選單元的權(quán)重矩陣；

b——偏置項(xiàng)。

超參數(shù)^［21］是用來(lái)控制模型學(xué)習(xí)過(guò)程的參數(shù)，對(duì)模型的性能有重要影響。為了確定多步LSTM模型最優(yōu)超參數(shù)組合，引入改進(jìn)RS進(jìn)行超參數(shù)優(yōu)化。

1） 設(shè)定一個(gè)超參數(shù)空間H，每個(gè)超參數(shù)h_i取值［a_i，b_i］，參數(shù)空間可以表示為式（9）。

H={h_i|a_i≤h_i≤b_i，i=1，2，…，n} （9）

2） 從超參數(shù)空間H中隨機(jī)抽樣，得到超參數(shù)組合

H_r=（h₁，h₂，…，h_n） （10）

使用H_r訓(xùn)練LSTM模型并評(píng)估性能P_r。

P_r=Evaluate（H_r） （11）

選擇性能最好的超參數(shù)為當(dāng)前最優(yōu)超參數(shù)

h_best=argmin_h∈HP_r（h） （12）

3） 進(jìn)行算法改進(jìn)，圍繞當(dāng)前最優(yōu)超參數(shù)，定義一個(gè)新的超參數(shù)空間H_new。從H_new中隨機(jī)抽樣得到新的超參數(shù)組合h_iter。

h_iter=（h₁′，h₂′，…，h_n′） （13）

h_i′～Uniform（a_i′，b_i′） （14）

使用h_iter訓(xùn)練LSTM模型并評(píng)估性能，得到新的性能指標(biāo)P_iter。

P_iter=Evaluate（h_iter） （15）

如果P_iter比之前的性能指標(biāo)更好，更新h_best。

h_best=argmin_h∈HP_iter（h_i′） （16）

4）

重復(fù)運(yùn)行上述步驟，直到達(dá)到設(shè)定的迭代次數(shù)。

2.3 基于最優(yōu)加權(quán)的組合預(yù)測(cè)模型構(gòu)建

在SARIMA模型和改進(jìn)RS-多步LSTM模型的基礎(chǔ)上，采用最優(yōu)加權(quán)組合算法構(gòu)建組合預(yù)測(cè)模型。

1） 構(gòu)建偏差矩陣。

式中：N——農(nóng)機(jī)裝備需求采樣數(shù)據(jù)總數(shù)；

e_1t、e_2t——

SARIMA模型和改進(jìn)RS-多步LSTM模型在t個(gè)采樣數(shù)據(jù)預(yù)測(cè)值與真實(shí)值的誤差。

2） 確定最優(yōu)權(quán)重。

使R=［1，1］^T，通過(guò)拉格朗日乘子法求解最優(yōu)權(quán)重

［w₁ w₂］^T=E^-1R/R^TE^-1R （18）

s.t.{R^T［w₁ w₂］^T=1

w₁+w₂=1 （19）

式中：w₁、w₂——SARIMA模型和改進(jìn)RS-多步LSTM模型的權(quán)重系數(shù)。

3） 組合預(yù)測(cè)模型。模型組合后的農(nóng)機(jī)裝備需求計(jì)算公式為

Y=w₁y_1t+ w₂y_2t （20）

綜上，SARIMA-改進(jìn)RS-多步LSTM組合模型流程如圖3所示。

2.4 預(yù)測(cè)模型評(píng)估

為評(píng)價(jià)組合模型的性能，本文選取均方誤差MSE、平均絕對(duì)誤差MAE、決定系數(shù)R²作為評(píng)價(jià)指標(biāo)。

MSE=1/n∑n/i=1（y_i-?_i）² （21）

MAE=1/n∑n/i=1|y_i-?_i| （22）

R²=1-∑n/i=1（y_i-?_i）²/∑n/i=1（y_i-y_i-）² （23）

式中：n——觀測(cè)值的數(shù)量；

y_i——第i個(gè)觀測(cè)值的實(shí)際值；

?_i——第i個(gè)觀測(cè)值的預(yù)測(cè)值；

y_i-——觀測(cè)值的平均值。

3 實(shí)例驗(yàn)證

以LX1504型號(hào)農(nóng)機(jī)裝備60個(gè)月月度需求為時(shí)間序列樣本（數(shù)據(jù)來(lái)源于農(nóng)機(jī)360網(wǎng)）。以月度主要農(nóng)作物價(jià)格、燃料價(jià)格、第一產(chǎn)業(yè)從業(yè)人員數(shù)量、農(nóng)村居民人均可支配收入為影響因素時(shí)間序列樣本。

3.1 SARIMA參數(shù)選擇

繪制時(shí)間序列數(shù)據(jù)的歷史需求數(shù)據(jù)曲線，對(duì)時(shí)間序列進(jìn)行時(shí)序分解，如圖4所示，將原始時(shí)間序列分解為趨勢(shì)、季節(jié)性和殘差三個(gè)部分。圖4（b）曲線呈增長(zhǎng)趨勢(shì)，說(shuō)明原始時(shí)間序列不具有平穩(wěn)性，需要進(jìn)行差分處理。圖4（c）曲線呈現(xiàn)規(guī)則波動(dòng)，說(shuō)明原始時(shí)間序列季節(jié)性顯著，依據(jù)曲線波動(dòng)情況，確定參數(shù)s為12。對(duì)圖4（a）原始時(shí)間序列進(jìn)行一階差分得到平穩(wěn)序列，確定參數(shù)d為1。繪制序列自相關(guān)圖（ACF）和偏自相關(guān)圖（PACF）判斷參數(shù)階數(shù)，如圖5所示。

由圖5可知，ACF，PACF近似一階截尾，參數(shù)p和q分別可以取值0或1。使用兩個(gè)嵌套循環(huán)遍歷p，d，q和P，D，Q，s中的參數(shù)組合，得出不同參數(shù)組合下模型的AIC值，如表1所示。（1，1，1）（1，1，1，12）組合對(duì)應(yīng)的AIC值最小，選擇該組合作為SARIMA模型的參數(shù)組合，構(gòu)建模型進(jìn)行預(yù)測(cè)。

3.2 改進(jìn)RS-多步LSTM模型構(gòu)建

1） 數(shù)據(jù)CEEMDAN處理結(jié)果及分析。對(duì)農(nóng)機(jī)裝備需求、SARIMA殘差和主要農(nóng)作物價(jià)格、燃料價(jià)格、第一產(chǎn)業(yè)從業(yè)人員數(shù)量、農(nóng)村居民人均可支配收入4個(gè)影響因素?cái)?shù)據(jù)進(jìn)行CEEMDAN處理，結(jié)果如圖6、圖7所示。各組原始數(shù)據(jù)被分解成若干個(gè)不同頻率和波動(dòng)的IMF和一個(gè)res項(xiàng)，分解后的序列顯示各個(gè)序列的復(fù)雜性逐漸降低，表明算法成功地分離出了信號(hào)中不同尺度和頻率的成分。

圖6、圖7中，農(nóng)機(jī)裝備需求量、SARIMA殘差、主要農(nóng)作物價(jià)格、燃料價(jià)格、農(nóng)村居民人均可支配收入分解結(jié)果具有相似性，IMF₂、IMF₃、IMF₄表現(xiàn)出較規(guī)則的波動(dòng)，反映了原始數(shù)據(jù)的周期特性。res項(xiàng)反映了原始數(shù)據(jù)的長(zhǎng)期趨勢(shì)為上升趨勢(shì)。

第一產(chǎn)業(yè)從業(yè)人員數(shù)量原始曲線存在微小波動(dòng)。分解后，IMF₂、IMF₄呈現(xiàn)出較規(guī)則的波動(dòng)，反映了原始數(shù)據(jù)的周期特性。res項(xiàng)和原始曲線高度相同，反映了原始數(shù)據(jù)的趨勢(shì)。

將以上具有周期性的IMF及res集成為一個(gè)輸入矩陣，作為預(yù)測(cè)模型的輸入端。

2）改進(jìn)RS超參數(shù)優(yōu)化。設(shè)置多步LSTM的隱藏單元數(shù)量、dropout比例、學(xué)習(xí)率、LSTM層數(shù)、是否返回每個(gè)時(shí)間步的輸出以及是否保存每個(gè)時(shí)間步的狀態(tài)等超參數(shù)組合為優(yōu)化對(duì)象。對(duì)優(yōu)化前后的模型損失率進(jìn)行比較分析，如圖8所示。由圖8可知，改進(jìn)RS-多步LSTM模型相比未進(jìn)行優(yōu)化的平均多步LSTM模型，模型的損失率更低，表示模型預(yù)測(cè)性能更好。

3.3 預(yù)測(cè)結(jié)果分析及評(píng)價(jià)

基于SARIMA-改進(jìn)RS-多步LSTM組合模型進(jìn)行預(yù)測(cè)，驗(yàn)證模型的有效性。圖9為各預(yù)測(cè)模型預(yù)測(cè)結(jié)果，表2為預(yù)測(cè)模型效果評(píng)價(jià)。

由表2可知，SARIMA-改進(jìn)RS-多步LSTM組合模型的MSE、MAE均小于其他4種模型，決定系數(shù)R²大于其他4種模型，達(dá)0.920 9。表明對(duì)于受季節(jié)性、多因素影響的農(nóng)機(jī)裝備需求預(yù)測(cè)問(wèn)題，本文構(gòu)建的組合模型能夠提供更加可靠、準(zhǔn)確的預(yù)測(cè)結(jié)果，證明了SARIMA-改進(jìn)RS-多步LSTM組合模型的可靠性和有效性。

4 結(jié)論

構(gòu)建SARIMA-改進(jìn)RS-多步LSTM組合模型，對(duì)周期性、非線性的某農(nóng)機(jī)裝備制造商LX1504型號(hào)農(nóng)機(jī)裝備需求時(shí)間序列進(jìn)行預(yù)測(cè)。

1） 引入CEEMDAN對(duì)數(shù)據(jù)特征信息進(jìn)行分解和提取，經(jīng)CEEMDAN處理后多步LSTM模型的MSE、MAE分別降低15.35%和4.29%，R²由0.828 4增加到0.854 7，預(yù)測(cè)精度提高3.17%，降低數(shù)據(jù)非平穩(wěn)性對(duì)預(yù)測(cè)結(jié)果的影響。

2） 設(shè)計(jì)改進(jìn)RS算法對(duì)多步LSTM模型超參數(shù)進(jìn)行優(yōu)化，降低模型的損失率。改進(jìn)RS算法對(duì)多步LSTM超參數(shù)優(yōu)化后，MSE、MAE分別降低21.00%和10.03%，R²由0.854 7增加到0.885 2，預(yù)測(cè)精度提高3.57%。

3） 基于最優(yōu)加權(quán)組合，構(gòu)建組合預(yù)測(cè)模型。相比單一SARIMA模型和改進(jìn)RS-多步LSTM模型，組合模型的MSE、MAE分別降低59.37%、31.06%和32.14%、12.28，R²分別由0.805 3和0.885 2增加到0.920 9，預(yù)測(cè)精度分別提高14.35%和4.03%。

4） 綜上，SARIMA-改進(jìn)RS-多步LSTM組合模型相比其他模型MSE、MAE指標(biāo)均小于其他4種模型，在預(yù)測(cè)準(zhǔn)確性和誤差控制方面表現(xiàn)更好，決定系數(shù)R²大于其他4種模型，達(dá)0.920 9，有更高的預(yù)測(cè)精度，能夠更好地?cái)M合數(shù)據(jù)并反映原始數(shù)據(jù)的特征，驗(yàn)證本文所提方法的有效性，能為農(nóng)機(jī)企業(yè)的農(nóng)機(jī)裝備需求預(yù)測(cè)提供參考。

參 考 文 獻(xiàn)

［1］ 何登旭， 李艷芳， 劉向虎， 等. 基于泛函網(wǎng)絡(luò)的非線性回歸預(yù)測(cè)模型及學(xué)習(xí)算法［J］. 計(jì)算機(jī)工程與應(yīng)用， 2008， 44（24）： 74-77.

He Dengxu， Li Yanfang， Liu Xianghu， et al. Non-linear regression forecast model based on functional networks and learning algorithm ［J］. Computer Engineering and Applications， 2008， 44（24）： 74-77.

［2］ 李之紅， 申天宇， 文琰杰， 等. 基于混合機(jī)器學(xué)習(xí)框架的網(wǎng)約車(chē)訂單需求預(yù)測(cè)與異常點(diǎn)識(shí)別［J］. 交通信息與安全， 2023， 41（3）： 157-165， 174.

Li Zhihong， Shen Tianyu， Wen Yanjie， et al. Order demand prediction and anomaly-point identification for online car-hailing orders based on hybrid machine learning framework ［J］. Journal of Transport Information and Safety， 2023， 41（3）： 157-165， 174.

［3］ Maulud D， Abdulazeez A M. A review on linear regression comprehensive in machine learning ［J］. Journal of Applied Science and Technology Trends， 2020， 1（4）： 140-147.

［4］ 吳庚奇， 牛東曉， 耿世平， 等. 多價(jià)值鏈視角下基于深度學(xué)習(xí)算法的制造企業(yè)產(chǎn)品需求預(yù)測(cè)［J］. 科學(xué)技術(shù)與工程， 2021， 21（31）： 13413-13420.

Wu Gengqi， Niu Dongxiao， Geng Shiping， et al. Product demand forecasting of manufacturing enterprises based on deep learning algorithm from the perspective of multi-value chains ［J］. Science Technology and Engineering， 2021， 21（31）： 13413-13420.

［5］ 武慧榮， 陳少陽(yáng)， 崔淑華. 基于GM（1，1）-MLP神經(jīng)網(wǎng)絡(luò)模型的大宗貨物運(yùn)輸需求預(yù)測(cè)［J］. 公路交通科技， 2023， 40（10）： 233-240.

Wu Huirong， Chen Shaoyang， Cui Shuhua. Transportation demand forecast of bulk cargo based on GM（1，1）-MLP neural network model ［J］. Journal of Highway and Transportation Research and Development， 2023， 40（10）： 233-240.

［6］ 崔凱， 郭宇， 錢(qián)偉偉， 等. 基于Group-BiLSTM-LightGBM的物料配送需求量預(yù)測(cè)［J］. 組合機(jī)床與自動(dòng)化加工術(shù)， 2022（8）： 164-168.

Cui Kai， Guo Yu， Qian Weiwei， et al. Material distribution demand forecasting method based on Group-BiLSTM-LightGBM ［J］. Modular Machine Tool amp; Automatic Manufacturing Technique， 2022（8）： 164-168.

［7］ 王海程， 馬紀(jì)穎， 張?jiān)锋拢?等. 融合時(shí)序關(guān)聯(lián)動(dòng)態(tài)圖與常微分方程的區(qū)域間出租車(chē)需求預(yù)測(cè)［J］. 計(jì)算機(jī)應(yīng)用研究， 2024， 41（3）： 794-798， 860.

Wang Haicheng， Ma Jiying， Zhang Yuanyuan， et al. Inter-regional taxi demand forecasting based on time series correlation dynamic graph and ordinary differential equation ［J］. Application Research of Computers， 2024， 41（3）： 794-798， 860.

［8］ 黃國(guó)興， 曹先懷， 錢(qián)曉飛. 一種基于隨機(jī)森林的備件預(yù)測(cè)模型研究［J］. 運(yùn)籌與管理， 2021， 30（10）： 165-168.

Huang Guoxing， Cao Xianhuai， Qian Xiaofei. Research on spare parts prediction model based on random forest ［J］. Operations Research and Management Science， 2021， 30（10）： 165-168.

［9］ 付維方， 穆彩虹， 劉英杰. 基于機(jī)器學(xué)習(xí)方法的航空消耗件需求自適應(yīng)預(yù)測(cè)［J］. 科學(xué)技術(shù)與工程， 2022， 22（11）： 4609-4617.

Fu Weifang， Mu Caihong， Liu Yingjie. Adaptive demand prediction of aviation consumable spare parts based on machine learning method ［J］. Science Technology and Engineering， 2022， 22（11）： 4609-4617.

［10］ 張舜， 楊欣， 魏津瑜. 基于預(yù)防性維修的農(nóng)機(jī)備件需求預(yù)測(cè)問(wèn)題研究［J］. 中國(guó)農(nóng)機(jī)化學(xué)報(bào)， 2020， 41（6）： 225-230.

Zhang Shun， Yang Xin， Wei Jinyu. Research on forecasting of agricultural machinery spare parts demand based on preventive maintenance ［J］. Journal of Chinese Agricultural Mechanization， 2020， 41（6）： 225-230.

［11］ 李瓊， 陳佳， 徐斌. 多產(chǎn)品相互制約條件下的制造業(yè)產(chǎn)品需求預(yù)測(cè)研究［J］. 計(jì)算機(jī)應(yīng)用與軟件， 2021， 38（3）： 59-69.

Li Qiong， Chen Jia， Xu Bin. Demand forecasting of manufacturing products under the condition of multi product mutual restriction ［J］. Computer Applications and Software， 2021， 38（3）： 59-69.

［12］ Chen X， Wang T， Liang W. General aircraft material demand forecast based on PSO-BP neural network ［J］. International Journal of Control and Automation， 2016， 9（5）： 407-418.

［13］ 李潔， 彭其淵， 楊宇翔. 基于SARIMA模型的廣珠城際鐵路客流量預(yù)測(cè)［J］. 西南交通大學(xué)學(xué)報(bào)， 2020， 55（1）： 41-51.

LiJie， Peng Qiyuan， Yang Yuxiang. Passenger flow prediction for Guangzhou-Zhuhai intercity railway based on SARIMA model ［J］. Journal of Southwest Jiaotong University， 2020， 55（1）： 41-51.

［14］ Hochreiter S， Schmidhuber J. Long short-term memory ［J］.Neural Computation， 1997， 9（8）： 1735-1780.

［15］ Yu Y， Si X， Hu C， et al. A review of recurrent neural networks： LSTM cells and network architectures ［J］. Neural Computation， 2019， 31（7）： 1235-1270.

［16］ 馮建英， 穆維松， 張領(lǐng)先， 等. 基于消費(fèi)者購(gòu)買(mǎi)意愿的農(nóng)機(jī)市場(chǎng)需求分析［J］. 商業(yè)研究， 2008（2）： 191-194.

Feng Jianying， Mu Weisong， Zhang Lingxian， et al. An empirical study on consumers’ purchase intention in agricultural machinery Market ［J］. Commercial Research， 2008（2）： 191-194.

［17］ 陳旭， 楊印生， 魏思琳. 基于向后逐步回歸模型的我國(guó)農(nóng)機(jī)需求特征及影響因素研究［J］. 數(shù)理統(tǒng)計(jì)與管理， 2017， 36（5）： 774-782.

Chen Xu， Yang Yinsheng， Wei Silin. A study on characteristics and influencing factors of agricultural machinery demand in China based on backward stepwise regression model ［J］. Journal of Applied Statistics and Management， 2017， 36（5）： 774-782.

［18］ 張標(biāo)， 張領(lǐng)先， 傅澤田， 等. 農(nóng)戶(hù)農(nóng)機(jī)需求及購(gòu)買(mǎi)行為分析——基于18省的微觀調(diào)查數(shù)據(jù)實(shí)證［J］. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2017， 22（11）： 208-223.

Zhang Biao， Zhang Lingxian， Fu Zetian， et al. Analysis on farmers’ requirements for agricultural machinery and purchasing behavior： Based on the survey data of 18 provinces ［J］. Journal of China Agricultural University， 2017， 22（11）： 208-223.

［19］ Kuiper R. AIC-type theory-based model selection for structural equation models ［J］. Structural Equation Modeling， 2022， 29（1）： 151-158.

［20］ Torres M E， Colominas M A， Schlotthauer G， et al. A complete ensemble empirical mode decomposition with adaptive noise ［C］. Proceedings of 2011 IEEE International Conference on Acoustics， Speech and Signal Processing， IEEE， 2011： 4144-4147.

［21］ 韓惠， 于巧， 祝義. 超參數(shù)優(yōu)化對(duì)跨版本缺陷預(yù)測(cè)影響的實(shí)證研究［J］. 計(jì)算機(jī)科學(xué)與探索， 2023， 17（12）：3052-3064.

Han Hui， Yu Qiao， Zhu Yi. Impact of hyperparameter optimization on cross-version defect prediction： An empirical study ［J］. Journal of Frontiers of Computer Science and Technology， 2023， 17（12）： 3052-3064.