基于逐步回歸法和BP 神經(jīng)網(wǎng)絡(luò)模型的蘋果產(chǎn)量預(yù)測(cè)

田東霞 曹久才

（1 北京市昌平區(qū)氣象局，北京 102200；2 北京市門頭溝區(qū)氣象局，北京 102308）

果業(yè)生產(chǎn)是中國(guó)現(xiàn)代農(nóng)業(yè)發(fā)展的重要組成部分，氣候變化對(duì)果業(yè)發(fā)展有著重要影響。 目前，許多學(xué)者用不同方法開(kāi)展了產(chǎn)量預(yù)測(cè)工作， 并針對(duì)產(chǎn)量預(yù)測(cè)方法提出了建議：李蓬勃等[1]研究了BP 神經(jīng)網(wǎng)絡(luò)和多元線性回歸的產(chǎn)量預(yù)測(cè)模型， 并對(duì)比了不同統(tǒng)計(jì)方法下的模型預(yù)測(cè)誤差；莊 星等[2]運(yùn)用混合群智能算法優(yōu)化了神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型， 降低了產(chǎn)量預(yù)測(cè)誤差；李 環(huán)等[3]基于NARX 神經(jīng)網(wǎng)絡(luò)建立了糧食產(chǎn)量預(yù)測(cè)模型，取得了較好的效果；蘇 博等[4]運(yùn)用GM（1，N）灰色系統(tǒng)與BP 神經(jīng)網(wǎng)絡(luò)方法進(jìn)行了糧食產(chǎn)量預(yù)測(cè)比較研究，并分析了其預(yù)測(cè)能力；伍丹華等[5]以糧食種植面積、農(nóng)藥施用量、化肥施用量、農(nóng)業(yè)機(jī)械化總動(dòng)力、 農(nóng)用排灌動(dòng)力機(jī)械等5 個(gè)指標(biāo)為神經(jīng)網(wǎng)絡(luò)輸入值建立了預(yù)測(cè)模型，并進(jìn)行了產(chǎn)量預(yù)測(cè)。

本項(xiàng)目結(jié)合北京市昌平區(qū)近15 年的氣象數(shù)據(jù)對(duì)蘋果產(chǎn)量進(jìn)行預(yù)測(cè)， 旨在為蘋果防災(zāi)減災(zāi)和產(chǎn)業(yè)發(fā)展提供合理化建議，指導(dǎo)蘋果生產(chǎn)趨利避害，有效增加蘋果產(chǎn)量，提高蘋果產(chǎn)業(yè)氣象服務(wù)水平，為蘋果產(chǎn)業(yè)科學(xué)應(yīng)對(duì)氣候變化提供決策依據(jù)， 從而更好地為農(nóng)村穩(wěn)定、農(nóng)業(yè)增產(chǎn)、農(nóng)民增收提供針對(duì)性和專業(yè)性的氣象服務(wù)，全面助力鄉(xiāng)村振興。

1 資料與方法

1.1 資料來(lái)源

1.1.1產(chǎn)量資料。 昌平區(qū)蘋果產(chǎn)量數(shù)據(jù)資料來(lái)源于北京市昌平區(qū)統(tǒng)計(jì)局，包括蘋果種植面積和總產(chǎn)量。蘋果總產(chǎn)量受種植面積的影響， 因而采用單產(chǎn)數(shù)據(jù)（總產(chǎn)量/種植面積）分析產(chǎn)量與氣象因子的關(guān)系。

1.1.2氣象資料。 蘋果生育期分為萌芽至開(kāi)花期（3 月25 日至4 月30 日）、坐果期（5 月）、果實(shí)生長(zhǎng)前期（6 月）、果實(shí)生長(zhǎng)后期（7 月）、膨大期（8 月）、著色前期（9 月）、著色后期（10 月）。影響蘋果產(chǎn)量的因素有很多， 本文選取蘋果生長(zhǎng)過(guò)程中所需要的光、熱、水等氣象因素進(jìn)行分析。氣象資料來(lái)源于昌平國(guó)家地面氣象觀測(cè)站，包括2004—2018 年逐月平均溫度、最高氣溫、最低氣溫、氣溫日較差、極端最高氣溫、極端最低氣溫、日照時(shí)數(shù)、降水、濕度等數(shù)據(jù)。

1.2 研究方法

運(yùn)用逐步回歸法和BP 神經(jīng)網(wǎng)絡(luò)建立蘋果產(chǎn)量預(yù)測(cè)模型， 運(yùn)用2004—2018 年產(chǎn)量數(shù)據(jù)進(jìn)行運(yùn)算，運(yùn)用2019—2020 年產(chǎn)量數(shù)據(jù)進(jìn)行模型檢驗(yàn)和對(duì)比。

2 結(jié)果與分析

2.1 逐步回歸模型

多元逐步線性回歸是將全部自變量引入方程后，根據(jù)自變量對(duì)因變量影響程度的大小，利用相關(guān)顯著性檢驗(yàn)系數(shù)依次剔除無(wú)顯著意義的自變量，最終篩選出最具有統(tǒng)計(jì)學(xué)意義的自變量并建立回歸模型，以反映多個(gè)自變量與因變量之間的相關(guān)關(guān)系。

將可能影響昌平區(qū)蘋果產(chǎn)量的光、 熱等氣象因子作為自變量x、蘋果單產(chǎn)量為因變量y進(jìn)行逐步回歸分析。 通過(guò)逐步回歸法挑選出4 個(gè)影響蘋果產(chǎn)量的關(guān)鍵氣象因子，分別為4 月平均最低氣溫、8 月平均氣溫、10 月平均最低氣溫、10 月累計(jì)日照時(shí)數(shù)。根據(jù)2004—2018 年的關(guān)鍵氣候因子資料和產(chǎn)量數(shù)據(jù)建立蘋果產(chǎn)量預(yù)測(cè)模型， 得到蘋果產(chǎn)量的回歸模型及相關(guān)系數(shù)、F檢驗(yàn)結(jié)果。 方程通過(guò)了顯著性檢驗(yàn)（P<0.05），相關(guān)系數(shù)R為0.824，F(xiàn)為6.34，顯著性水平P為0.005 6。 結(jié)果表明，回歸模型對(duì)產(chǎn)量預(yù)測(cè)精度較高，預(yù)測(cè)模型的相關(guān)性和擬合效果好，顯著性水平值P＜0.05，故總體預(yù)測(cè)效果較好。四因子擬合方程為：

式中，y為蘋果單產(chǎn)（kg/hm2），x1、x2、x3、x4分別為8 月平均氣溫（℃）、10 月平均最低氣溫（℃）、10 月累計(jì)日照時(shí)數(shù)（h）、4 月平均最低氣溫（℃）。

2.2 BP 神經(jīng)網(wǎng)絡(luò)模型

BP 神經(jīng)網(wǎng)絡(luò)（back-propagation network）又稱反向傳播神經(jīng)網(wǎng)絡(luò)，即通過(guò)對(duì)樣本數(shù)據(jù)進(jìn)行訓(xùn)練，不斷修正網(wǎng)絡(luò)權(quán)值和閾值， 使誤差函數(shù)沿負(fù)梯度方向下降，逼近期望輸出，循環(huán)直到誤差達(dá)到允許的范圍之內(nèi)或訓(xùn)練次數(shù)達(dá)到預(yù)先設(shè)計(jì)的次數(shù)為止[6]。 BP 神經(jīng)網(wǎng)絡(luò)是一種多層的前饋神經(jīng)網(wǎng)絡(luò)，主要包含輸入層、隱含層和輸出層。每一層都由若干個(gè)節(jié)點(diǎn)組成，模型結(jié)構(gòu)如圖1 所示。

本文選用MATLAB 中的神經(jīng)網(wǎng)絡(luò)工具箱進(jìn)行網(wǎng)絡(luò)訓(xùn)練。 為了加快神經(jīng)網(wǎng)絡(luò)的收斂速度， 使用premnmx 函數(shù)對(duì)原始數(shù)據(jù)進(jìn)行歸一化處理。 該模型利用2004—2018 年數(shù)據(jù)進(jìn)行運(yùn)算，以逐步回歸法選出的影響蘋果產(chǎn)量的4 個(gè)氣象因子數(shù)據(jù)指標(biāo)作為輸入，以研究的預(yù)測(cè)產(chǎn)量指標(biāo)作為輸出，因而輸入層的節(jié)點(diǎn)數(shù)為4、輸出層的節(jié)點(diǎn)數(shù)為1，根據(jù)經(jīng)驗(yàn)公式確定隱含層的節(jié)點(diǎn)數(shù)為5， 網(wǎng)絡(luò)迭代次數(shù)為5 000 次，期望誤差為0.000 000 1。 經(jīng)驗(yàn)公式如下：式中，n1為隱含層節(jié)點(diǎn)數(shù)，n為輸入層節(jié)點(diǎn)數(shù)，m為輸出層節(jié)點(diǎn)數(shù)，k為1~10 之間的調(diào)節(jié)常數(shù)。

2.3 模型檢驗(yàn)

利用2019—2020 年數(shù)據(jù)進(jìn)行檢驗(yàn)：逐步回歸預(yù)測(cè)模型計(jì)算得出的2019 年預(yù)測(cè)產(chǎn)量平均誤差為21%，2020 年預(yù)測(cè)產(chǎn)量平均誤差為15.5%；BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型計(jì)算得出的2019 年預(yù)測(cè)產(chǎn)量平均誤差為3.8%，2020 年預(yù)測(cè)產(chǎn)量平均誤差為12.4%（表1）。 通過(guò)2 年數(shù)據(jù)可以看出，利用BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型反演得到的蘋果產(chǎn)量與實(shí)際產(chǎn)量吻合度更高，較逐步回歸法具有更高的預(yù)測(cè)精度。

表1 逐步回歸法和BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型檢驗(yàn)結(jié)果

利用2004—2020 年數(shù)據(jù)分析對(duì)比2 種預(yù)測(cè)模型（圖2），可以看出，2 種預(yù)測(cè)模型擬合效果均較好，均能較好地預(yù)測(cè)今后蘋果產(chǎn)量趨勢(shì)，且BP 神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)更接近實(shí)際產(chǎn)量。 通過(guò)實(shí)際產(chǎn)量與預(yù)測(cè)產(chǎn)量的相關(guān)性分析可看出： 逐步回歸模型的預(yù)測(cè)產(chǎn)量與實(shí)際產(chǎn)量相關(guān)性分析方程斜率為0.715 9，擬合度（R2）為0.657 2（圖3），預(yù)測(cè)結(jié)果的平均相對(duì)誤差為11.6%；神經(jīng)網(wǎng)絡(luò)模型具有高度非線性和較強(qiáng)的泛化能力，BP 神經(jīng)網(wǎng)絡(luò)的模型預(yù)測(cè)結(jié)果優(yōu)于逐步回歸法模型，神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型相關(guān)性分析方程斜率較大（為1.016 9）、截距最小，擬合度（R2）為0.971 6（圖4），預(yù)測(cè)產(chǎn)量平均相對(duì)誤差較?。?.2%）。

但是，BP 神經(jīng)網(wǎng)絡(luò)模型存在一定的局限性，其缺少對(duì)主導(dǎo)預(yù)測(cè)因子進(jìn)行篩選的功能[4]，只能依照逐步回歸法得出的結(jié)果， 將挑選出的影響蘋果產(chǎn)量的關(guān)鍵氣象因子（4 月平均最低氣溫、8 月平均氣溫、10 月平均最低氣溫、10 月累計(jì)日照時(shí)數(shù)）作為神經(jīng)網(wǎng)絡(luò)輸入層參數(shù)，也存在收斂速度慢、迭代步數(shù)多、易于陷入局部極小等缺點(diǎn)。

3 結(jié)論

本文運(yùn)用逐步回歸法和BP 神經(jīng)網(wǎng)絡(luò)建立蘋果產(chǎn)量預(yù)測(cè)模型。通過(guò)檢驗(yàn)發(fā)現(xiàn)，2 種模型擬合效果均較好，均能較好地預(yù)測(cè)今后蘋果產(chǎn)量趨勢(shì)。 BP 神經(jīng)網(wǎng)絡(luò)模型有較高的預(yù)測(cè)精度， 但存在一定的局限性。 產(chǎn)量預(yù)測(cè)能有效指導(dǎo)蘋果生長(zhǎng)發(fā)育過(guò)程中針對(duì)氣候進(jìn)行趨利避害，為果農(nóng)進(jìn)行蘋果生產(chǎn)管理提供參考。