中國(guó)農(nóng)學(xué)通報(bào) 2018ChineseAgricultural Science Bulletin2018,34(16):115-125基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室溫濕度預(yù)測(cè)模型對(duì)比分析溫永菁1,2，李 春2，薛慶禹2，董朝陽2，黎貞發(fā)2，李秀芬1，陳思寧2（1沈陽農(nóng)業(yè)大學(xué)，沈陽110866；2天津市氣候中心，天津300074）摘 要：為構(gòu)建較準(zhǔn)確的日光溫室溫濕度預(yù)測(cè)模型，于20112014年冬季（1月、2月、12月）在天津市寶坻區(qū)開展溫室內(nèi)外環(huán)境監(jiān)測(cè)試驗(yàn)，并建立3種天氣類型（晴、多云、陰）下3個(gè)時(shí)段（08時(shí)、817時(shí)、1723時(shí)）逐步回歸與BP神經(jīng)網(wǎng)絡(luò)溫室內(nèi)溫濕度預(yù)測(cè)模型。結(jié)果表明；（1）溫室內(nèi)氣溫逐步回歸模型9種情況下模擬值與實(shí)際值的絕對(duì)誤差小于3的平均準(zhǔn)確率(3)為88%，平均均方根誤差(RMSE)為2；BP神經(jīng)網(wǎng)絡(luò)模型9種情況下模擬值與實(shí)際值的絕對(duì)誤差小于3的平均準(zhǔn)確率(3)為94%，平均均方根誤差(RMSE)為1.6。應(yīng)用BP神經(jīng)網(wǎng)絡(luò)建立的氣溫預(yù)測(cè)模型相對(duì)更為準(zhǔn)確穩(wěn)定；（2）相對(duì)濕度逐步回歸模型9種情況下模擬值與實(shí)際值的絕對(duì)誤差小于6%的平均準(zhǔn)確率(6%)為81%，平均均方根誤差(RMSE)為5.7%；BP神經(jīng)網(wǎng)絡(luò)模型9種情況下模擬值與實(shí)際值的絕對(duì)誤差小于6%的平均準(zhǔn)確率(6%)為80%，平均均方根誤差(RMSE)為6.7%。2類模型均不適宜預(yù)測(cè)817時(shí)日光溫室相對(duì)濕度，而1723時(shí)與08時(shí)應(yīng)用逐步回歸建立的濕度預(yù)測(cè)模型相對(duì)更準(zhǔn)確穩(wěn)定。關(guān)鍵詞：日光溫室；逐步回歸；BP神經(jīng)網(wǎng)絡(luò)；溫濕度模型中圖分類號(hào)：S625 文獻(xiàn)標(biāo)志碼：A 論文編號(hào)：casb17060118Temperature and Humidity Prediction Models in Solar Greenhouse:Comparative Analysis Based on Stepwise Regression and BP Neural NetworkWen Yongjing1,2, Li Chun2, Xue Qingyu2, Dong Chaoyang2, Li Zhenfa2, Li Xiufen1, Chen Sining2(1Shenyang Agricultural University, Shenyang 110866;2Tianjin Climate Center, Tianjin 300074)Abstract: To build a more accurate model for predicting temperature and humidity in solar greenhouse, amicroclimatic observing experiment was carried out at Baodi District of Tianjin from 2011 to 2014 in winter(January, February and December). The temperature and humidity prediction models in greenhouse by usingstepwise regression and BP neural network were established at 3 periods (0:00-8:00, 8:00-17:00, 17:00-23:00)of 3 kinds of weather types (sunny, cloudy, overcast). The results showed that: (1) the average accuracy rate ofthe absolute error of simulated and actual values less than 3 was 88%, and the root-mean-square error(REMS) was 2 under 9 conditions in greenhouse by using stepwise regression model of temperature; theaverage accuracy rate of the absolute error of simulated and actual values less than 3 was 94%, and the root-mean-square error (RMES) was 1.6 under 9 conditions in greenhouse by using BP neural network model oftemperature; the temperature prediction model established by BP neural network was more accurate and stable;(2) the average accuracy rate of the absolute error of simulated and actual values less than 6% was 81%, and基金項(xiàng)目：天津市農(nóng)業(yè)科技成果轉(zhuǎn)化與推廣項(xiàng)目“溫室小氣候資源高效利用及蔬菜茬口搭配技術(shù)集成與應(yīng)用”(201502150)；天津市科委青年基金“日光溫室黃瓜霜霉病監(jiān)測(cè)預(yù)警技術(shù)研究”(16JCQNJC14900)；天津市氣象局博士基金“日光溫室番茄低溫寡照影響評(píng)估與風(fēng)險(xiǎn)評(píng)價(jià)研究”(BSJJ201505)。第一作者簡(jiǎn)介：溫永菁，女，1993年出生，遼寧葫蘆島人，在讀碩士，主要從事農(nóng)業(yè)氣象研究。通信地址：300000天津市河西區(qū)氣象臺(tái)路100號(hào)天津市氣象局，E-mail：betterbaymax163.com。通訊作者：黎貞發(fā)，男，1964年出生，廣西橫縣人，正研級(jí)高級(jí)工程師，主要從事設(shè)施農(nóng)業(yè)氣象研究。通信地址：300000天津市河西區(qū)氣象臺(tái)路100號(hào)天津市氣象局，E-mail：lzfaaa126.com。收稿日期：2017-06-26，修回日期：2017-08-25。中國(guó)農(nóng)學(xué)通報(bào) http:/www.casb.org.cnthe root-mean-square error (REMS) was 5.7% under 9 conditions in greenhouse by using stepwise regressionmodel of relative humidity; and the average accuracy rate of the absolute error of simulated and actual valuesless than 6% was 80%, and the root-mean-square error (RMES) was 6.7% under 9 conditions in greenhouseby using BP neural network model of relative humidity. Both the 2 models are not suitable for predicting therelative humidity of solar greenhouse at 8:00- 17:00, while the humidity prediction model established bystepwise regression at 17:00-23:00 and 0:00-8:00 is more accurate and stable.Key words: solar greenhouse; stepwise regression; BP neural network; temperature and humidity model0 引言溫室小氣候是影響日光溫室作物生長(zhǎng)發(fā)育以及產(chǎn)量品質(zhì)的重要因素1，其中溫度條件是影響作物生命活動(dòng)的重要因素之一，溫室溫度適宜利于作物正常生長(zhǎng)發(fā)育，溫度過高過低都會(huì)對(duì)作物正常生長(zhǎng)發(fā)育產(chǎn)生影響，易導(dǎo)致減產(chǎn)；環(huán)境內(nèi)的水分情況對(duì)作物生長(zhǎng)發(fā)育影響亦較大，相對(duì)濕度是反映溫室內(nèi)水分情況的重要因子，溫室內(nèi)相對(duì)濕度的大小不僅直接影響作物蒸騰速率而且影響光合強(qiáng)度，溫室內(nèi)相對(duì)濕度過高易引發(fā)病蟲害的發(fā)生影響作物產(chǎn)量和品質(zhì)。研究探索準(zhǔn)確穩(wěn)定的日光溫室溫濕度預(yù)測(cè)模型對(duì)溫室作物的生產(chǎn)與發(fā)展至關(guān)重要。采用統(tǒng)計(jì)回歸法進(jìn)行線性建模一直被專家學(xué)者普遍使用，賀芳芳等2采用線性回歸分析法研究了上海地區(qū)荷蘭玻璃溫室溫濕度與外界氣象條件等影響因子的關(guān)系；李德等3采用逐步回歸法建立了秋季溫室最高與最低氣溫預(yù)報(bào)模型；張美玲等4采用逐步回歸法建立了日光溫室光合有效輻射預(yù)報(bào)模型；李寧等5采用主成分分析法建立了日光溫室內(nèi)低溫模型。近年來，BP神經(jīng)網(wǎng)絡(luò)在環(huán)境預(yù)測(cè)模擬方面應(yīng)用逐漸增多，初期Seginer等6利用神經(jīng)網(wǎng)絡(luò)建立了溫室氣候的模擬模型；Ferreira等7建立了溫室內(nèi)氣溫的神經(jīng)網(wǎng)絡(luò)模型；金志鳳等8利用BP神經(jīng)網(wǎng)絡(luò)法建立了楊梅大棚內(nèi)最高、最低氣溫預(yù)報(bào)模型；李倩等9利用BP神經(jīng)網(wǎng)絡(luò)法對(duì)南方塑料大棚不同季節(jié)溫濕度進(jìn)行了模擬；王春玲等10利用BP神經(jīng)網(wǎng)絡(luò)法建立不同天氣狀況下冬季日光溫室內(nèi)小氣候模型；劉淑梅等11利用BP神經(jīng)網(wǎng)絡(luò)法研究了日光溫室17天氣溫預(yù)報(bào)模型。多年來，建立日光溫室溫濕度預(yù)測(cè)模型普遍采用一種方法，且多為分季節(jié)建模。多數(shù)研究未能體現(xiàn)不同方法建模準(zhǔn)確性的偏差以及不同天氣狀況下不同時(shí)段日光溫室內(nèi)環(huán)境變化差異較大的情況。本研究選取統(tǒng)計(jì)回歸中經(jīng)典的逐步回歸法與BP神經(jīng)網(wǎng)絡(luò)為建模方法，根據(jù)20112013年冬季分天氣類型分時(shí)段數(shù)據(jù)建模，對(duì)比分析2種方法建立的日光溫室溫度濕度預(yù)測(cè)模型，以期探索建立更準(zhǔn)確的溫室溫濕度預(yù)測(cè)模型。1 資料與方法1.1 試驗(yàn)時(shí)間、地點(diǎn)試驗(yàn)于20112014年冬季溫室蓋膜期間在天津市寶坻區(qū)圣人莊村(3944N，11717E)進(jìn)行。1.2 試驗(yàn)材料1.2.1 日光溫室 天津地區(qū)典型的土墻式溫室，東西長(zhǎng)85 m，南北跨度9 m，脊高4.5 m，后墻高4.2 m，溫室內(nèi)種植番茄和黃瓜。1.2.2 資料來源 試驗(yàn)日光溫室內(nèi)的環(huán)境數(shù)據(jù)由天津市海洋氣象儀器廠生產(chǎn)的DZN1型小氣候儀觀測(cè)記錄，儀器安裝于溫室中部（溫室前沿及圍墻所構(gòu)成長(zhǎng)方形的中心點(diǎn)），采集間隔為10 min，觀測(cè)項(xiàng)目為溫室內(nèi)距地面1.5 m高度的空氣溫度、空氣相對(duì)濕度、總輻射，以及5、10、15、20、40 cm處土壤溫度。日光溫室外的氣象數(shù)據(jù)來自天津市寶坻區(qū)氣象局觀測(cè)站，要素包括空氣溫度、空氣濕度、平均風(fēng)速、不同深度土壤溫度。1.3 數(shù)據(jù)處理為保證數(shù)據(jù)的一致性與建模穩(wěn)定準(zhǔn)確性，試驗(yàn)選取20112013年冬季3個(gè)月數(shù)據(jù)為建模數(shù)據(jù)，2014年冬季3個(gè)月數(shù)據(jù)為檢驗(yàn)數(shù)據(jù)，并將觀測(cè)數(shù)據(jù)間隔處理為1 h，根據(jù)20112013年平均日變化情況（圖1）可知，08時(shí)白天溫度變化規(guī)律一致先升后降，817時(shí)傍晚溫度緩慢下降，08時(shí)夜晚溫度降為平穩(wěn)狀態(tài)，所以將全天分為08時(shí)、817時(shí)、1723時(shí)3個(gè)時(shí)05101520253001234567891011121314151617181920212223時(shí)刻氣溫/圖1 20112013年平均氣溫日變化情況116溫永菁等：基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室溫濕度預(yù)測(cè)模型對(duì)比分析段。根據(jù)氣象上的分類標(biāo)準(zhǔn)將日照百分率100%S60%、20%S60%、0S20%分為晴、多云、陰3種天氣類型12-14。1.4 研究方法1.4.1 逐步回歸 溫室外氣象因子直接或間接影響著溫室內(nèi)環(huán)境條件，因此建模之前篩選了25個(gè)溫室外氣象因子（室外1 h前氣溫、室外濕度、室外風(fēng)速、室外0 cm地溫、室外5 cm地溫、室外10 cm地溫、室外15 cm地溫、室外20 cm地溫、室外40 cm地溫、前1日室外最低氣溫、前1日室外最高氣溫、前1日室外平均氣溫、前1日室外相對(duì)濕度、前1日室外0 cm最高地溫、前1日室外0 cm最低地溫、前1日室外5 cm最高地溫、前1日室外5 cm最低地溫、前2日室外最低氣溫、前2日室外最高氣溫、前2日室外平均氣溫、前2日室外相對(duì)濕度、前2日室外0 cm最高地溫、前2日室外0 cm最低地溫、前2日室外5 cm最高地溫、前2日室外5 cm最低地溫）與溫室氣溫和溫室相對(duì)濕度進(jìn)行了相關(guān)性分析，保留了相關(guān)系數(shù)絕對(duì)值超過0.2的氣象因子進(jìn)行建模。表1中列舉了保留因子及其與模擬因子的相關(guān)系數(shù)。應(yīng)用經(jīng)典的統(tǒng)計(jì)回歸法逐步回歸建立日光溫室的溫濕度預(yù)測(cè)模型。采用Matlab 7.8軟件編程實(shí)現(xiàn)數(shù)據(jù)的相關(guān)性分析與逐步回歸。1.4.2 BP神經(jīng)網(wǎng)絡(luò) BP神經(jīng)網(wǎng)絡(luò)又稱誤差反向傳播(Back Propagation)神經(jīng)網(wǎng)絡(luò)，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一。其設(shè)計(jì)思想為信號(hào)正向傳播，誤差反向傳播。輸入信號(hào)正向傳播從輸入層通過隱藏層逐層處理，直至輸出誤差小于預(yù)先設(shè)定的目標(biāo)誤差或迭代次數(shù)超過最大訓(xùn)練次數(shù)停止，否則反向傳播15-18。設(shè)輸入層有n個(gè)神經(jīng)元，輸出層有m個(gè)神經(jīng)元，隱藏層有l(wèi)個(gè)神經(jīng)元；設(shè)i為輸入層神經(jīng)元的序號(hào)i=1,2,n；設(shè)j為隱藏層神經(jīng)元的序號(hào)j=1,2,l；設(shè)k為輸出層神經(jīng)元的序號(hào)k=1,2,m；設(shè)從輸入層神經(jīng)元i到隱藏層神經(jīng)元j的連接權(quán)值為wij，從隱藏層神經(jīng)元 j到輸出層神經(jīng)元k的連接權(quán)值為wjk；設(shè)第j個(gè)隱藏層神經(jīng)元的閾值為aj；設(shè)第k個(gè)輸出層神經(jīng)元閾值為bk；設(shè)第j個(gè)隱藏層神經(jīng)元的輸出為Hj；設(shè)輸入層第i個(gè)神經(jīng)元的值為xi，輸出層第k個(gè)神經(jīng)元的值為yk；設(shè)Q為調(diào)整激勵(lì)函數(shù)形勢(shì)的參數(shù)，則BP神經(jīng)網(wǎng)絡(luò)單胞元計(jì)算方法見公式(1)(8)。激勵(lì)函數(shù)：f(x)=11+e-xQ (1)計(jì)算第j個(gè)隱藏層神經(jīng)元輸出值：Hj=fi=1n( )wijxi-aj (2)計(jì)算第k個(gè)輸出層神經(jīng)元預(yù)測(cè)輸出值：Ok=j=1iHjwjk-bk (3)網(wǎng)絡(luò)預(yù)測(cè)誤差：ek=Yk-Qk (4)更新輸入層與隱藏層的連接權(quán)值：wij=wij+Hj( )1-Hjx( )i k=1m( )wjkek (5)更新隱藏層與輸出層的連接權(quán)值：wjk=wjk+Hjek (6)更新隱藏層神經(jīng)元閾值：aj=aj+Hj( )1-Hjk=1m( )wjkek (7)更新輸出層神經(jīng)元閾值：bk=bk+ek (8)反復(fù)迭代直至輸出誤差小于預(yù)先設(shè)定的目標(biāo)誤差或迭代次數(shù)超過最大訓(xùn)練次數(shù)為止19。本模型選定的相關(guān)參數(shù)值為：隱藏層神經(jīng)元數(shù)設(shè)為5，最大訓(xùn)練次數(shù)為40，初始學(xué)習(xí)速率為0.15，目標(biāo)誤差為0.0002。為提高訓(xùn)練效率和網(wǎng)絡(luò)泛化性能，采相對(duì)濕度室外1 h前氣溫OT室外濕度OH室外風(fēng)速OW室外0 cm地溫OGT0室外5 cm地溫OGT5前1日室外最低氣溫OTmin總輻射GR相關(guān)系數(shù)0.21-0.310.290.580.32-0.280.35氣溫室外1 h前氣溫OT室外濕度OH室外風(fēng)速OW室外0 cm地溫OGT0前1日室外相對(duì)濕度ORH1前2日室外相對(duì)濕度ORH2總輻射相關(guān)系數(shù)-0.360.46-0.41-0.670.210.24-0.37表1 日光溫室溫濕度預(yù)測(cè)模型入選因子相關(guān)系數(shù)注：OT為室外1 h前氣溫，OH為室外濕度，OW為室外風(fēng)速，OGT0為室外0 cm地溫，OGT5為室外5 cm地溫，OTmin為前1日室外最低氣溫，GR為總輻射，ORH1為前1日室外相對(duì)濕度，ORH2為前2日室外相對(duì)濕度。117中國(guó)農(nóng)學(xué)通報(bào) http:/www.casb.org.cn用歸一化方法將樣本數(shù)據(jù)進(jìn)行預(yù)處理。BP神經(jīng)網(wǎng)絡(luò)通過Matlab 7.8軟件編程實(shí)現(xiàn)。1.4.3 模型檢驗(yàn) 本試驗(yàn)采用模擬值與實(shí)際值的絕對(duì)誤差小于3的準(zhǔn)確率(3)和均方根誤差(RMSE)驗(yàn)證模型的精確度20，見公式(9)。RESE=i=1n( )SIMi-OBSi2n (9)式中：SIMi為模擬值，OBSi為實(shí)際值，n為樣本數(shù)。2 結(jié)果與分析由于結(jié)果中傍晚到夜間的預(yù)測(cè)效果相近，白天的預(yù)測(cè)效果差異較大，以下分為817時(shí)即白天，1723時(shí)和08時(shí)即傍晚到夜間2個(gè)部分進(jìn)行分析。2.1 基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室氣溫預(yù)測(cè)結(jié)果對(duì)比分析2.1.1 08時(shí)、1723時(shí)日光溫室氣溫預(yù)測(cè)模型 利用20112013年處理后數(shù)據(jù)應(yīng)用逐步回歸法建立不同天氣條件下傍晚到夜間日光溫室氣溫預(yù)測(cè)模型，見表2。利用2014年相應(yīng)時(shí)段數(shù)據(jù)進(jìn)行檢驗(yàn)，得到模擬值并進(jìn)行檢驗(yàn)分析。結(jié)果顯示，08時(shí)模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為57%、87.3%、98%，多云天為53%、82%、99%，晴天為52%、84%、95%，陰天為66%、96%、100%；模擬值與實(shí)際值的平均均方根誤差為1.3，多云天為1.4，晴天為1.5，陰天為1。1723時(shí)模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為44%、80.3%、95.3%，多云天為36%、70%、91%，晴天為45%、81%、95%，陰天為51%、90%、100%；模擬值與實(shí)際值的平均均方根誤差為1.5，多云天為1.8，晴天為1.5，陰天為1.2，見表3。結(jié)合預(yù)測(cè)值與實(shí)際值對(duì)比圖（見圖2）可以看出，3種天氣類型的預(yù)測(cè)結(jié)果基本一致，陰天的預(yù)測(cè)效果好于晴天和多云天，1723時(shí)、時(shí)段08時(shí)1723時(shí)天氣類型多云晴陰多云晴陰日光溫室氣溫預(yù)測(cè)逐步回歸模型y1=20.23+0.18x1-0.07x2-0.23x4+0.41x5y2=15.09+0.42x1-0.01x2-0.44x3-0.3x4+0.59x5-0.23x6+0.01x7y3=14.89+0.4x1-0.02x2-0.21x3-0.39x4+0.95x5-0.29x6y4=19.26+0.32x1-0.04x2-0.91x3-0.11x6y5=18.79+0.28x1-0.04x2-0.16x3+0.25x5-0.2x6y6=14.84+0.29x1-0.02x2-0.12x4+0.63x5-0.21x6不同方法逐步回歸BP神經(jīng)網(wǎng)絡(luò)時(shí)段08時(shí)1723時(shí)08時(shí)1723時(shí)天氣多云晴陰平均值多云晴陰平均值多云晴陰平均值多云晴陰平均值1準(zhǔn)確率/%535266573645514434495746.749453944.32準(zhǔn)確率/%82849687.370819080.388818785.383807077.73準(zhǔn)確率/%999510098919510095.398939896.394959093均方根誤差RMSE1.41.51.01.31.81.51.21.51.51.61.31.51.71.51.91.7樣本數(shù)96394189226.36925221017796394189226.369252210177注：x1為室外1 h前氣溫，x2為室外濕度，x3為室外風(fēng)速，x4為室外0 cm地溫，x5為室外5 cm地溫，x6為前1日室外最低氣溫，x7為總輻射。表2 不同天氣類型08時(shí)、1723時(shí)日光溫室氣溫預(yù)測(cè)逐步回歸模型表3 08時(shí)、1723時(shí)逐步回歸與BP神經(jīng)網(wǎng)絡(luò)日光溫室氣溫預(yù)測(cè)模型檢驗(yàn)結(jié)果118溫永菁等：基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室溫濕度預(yù)測(cè)模型對(duì)比分析08時(shí)為傍晚到夜間，日光溫室?guī)缀鯙槊荛]狀態(tài)，溫室氣溫變化規(guī)律基本穩(wěn)定，主要受溫室外氣象因素影響，其他因素干擾較少，因此光照影響較弱的陰天預(yù)測(cè)結(jié)果較其他天氣更好。BP神經(jīng)網(wǎng)絡(luò)建模數(shù)據(jù)和預(yù)測(cè)要素與逐步回歸保持一致，得到的模擬值進(jìn)行檢驗(yàn)分析。檢驗(yàn)結(jié)果為08時(shí)模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為46.7%、85.3%、96.3%，多云天為34%、88%、98%，晴天為49%、81%、93%，陰天為57%、87%、98%；模擬值與實(shí)際值的平均均方根誤差為1.5，多云天為1.5，晴天為1.6，陰天為1.3。1723時(shí)模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為44.3%、77.7%、93%，多云天為49%、83%、94%，晴天為45%、80%、95%，陰天為39%、70%、90%；模擬值與實(shí)際值的平均均方根誤差為1.7，多云天為1.7，晴天為1.5，陰天為1.9，見表3。結(jié)合圖2 不同天氣類型08時(shí)與1723時(shí)逐步回歸日光溫室氣溫預(yù)測(cè)模型模擬值與實(shí)際值對(duì)比(a)0null時(shí)多云1015202510152025模擬值/實(shí)際值/(b)0null時(shí)晴1015202510152025模擬值/實(shí)際值/(c)0null時(shí)陰1015202510152025模擬值/實(shí)際值/(d)17null3時(shí)多云1015202510152025模擬值/實(shí)際值/(e)17null3時(shí)晴1015202510152025模擬值/實(shí)際值/(f)17null3時(shí)陰1015202510152025模擬值/實(shí)際值/08時(shí)多云 08時(shí)晴08時(shí)陰 1723時(shí)多云1723時(shí)晴 1723時(shí)陰119中國(guó)農(nóng)學(xué)通報(bào) http:/www.casb.org.cn圖3可以看出，傍晚到夜間不同天氣類型的預(yù)測(cè)結(jié)果基本一致，陰天預(yù)測(cè)結(jié)果最好，多云和晴天較差，與逐步回歸預(yù)測(cè)模型進(jìn)行對(duì)比分析，BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型檢驗(yàn)結(jié)果略差于逐步回歸模型，1723時(shí)、08時(shí)為傍晚到夜間，日光溫室氣溫變化很穩(wěn)定，影響因素較少，溫室氣溫與溫室外氣象條件基本成線性關(guān)系，2種模型均可以實(shí)現(xiàn)日光溫室氣溫的預(yù)測(cè)。2.1.2 817時(shí)日光溫室氣溫預(yù)測(cè)模型 利用20112013年處理后數(shù)據(jù)應(yīng)用逐步回歸法建立不同天氣條件下白天日光溫室氣溫預(yù)測(cè)模型，見表4。利用2014年相應(yīng)時(shí)段數(shù)據(jù)進(jìn)行檢驗(yàn)，得到模擬值并進(jìn)行檢驗(yàn)分析。結(jié)果顯示，817時(shí)模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為27.7%、49.3%、68.7%，多云天為22%、45%、66%，晴天為18%、32%、50%，陰天為43%、71%、90%；模擬值與實(shí)際值的平均均方根誤差為3.3，多云天為圖3 不同天氣類型08時(shí)與1723時(shí)BP神經(jīng)網(wǎng)絡(luò)日光溫室氣溫預(yù)測(cè)模型模擬值與實(shí)際值對(duì)比(a)0null時(shí)多云1015202510152025模擬值/實(shí)際值/(b)0null時(shí)晴1015202510152025模擬值/實(shí)際值/(c)0null時(shí)陰1015202510152025模擬值/實(shí)際值/(d)17null3時(shí)多云1015202510152025模擬值/實(shí)際值/(e)17null3時(shí)晴1015202510152025模擬值/實(shí)際值/(f)17null3時(shí)陰1015202510152025模擬值/實(shí)際值/08時(shí)多云 08時(shí)晴08時(shí)陰1723時(shí)多云1723時(shí)晴1723時(shí)陰120溫永菁等：基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室溫濕度預(yù)測(cè)模型對(duì)比分析3.2，晴天為5，陰天為1.8，見表5。817時(shí)為白天，通風(fēng)、澆水、施肥、打叉等一系列農(nóng)事活動(dòng)都在該時(shí)段進(jìn)行，多方面因素影響溫室內(nèi)環(huán)境變化，單純的線性處理無法準(zhǔn)確地模擬溫室內(nèi)氣溫情況，晴天多云天預(yù)測(cè)結(jié)果較差，而陰天部分農(nóng)事活動(dòng)減少，光照影響較小，預(yù)測(cè)結(jié)果相對(duì)其他天氣較好。BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型817時(shí)檢驗(yàn)結(jié)果模擬值與實(shí)際值的絕對(duì)誤差小于1、2、3的準(zhǔn)確率(1、2、3)分別為46.7%、82.7%、94%，多云天為34%、75%、88%，晴天為53%、82%、94%，陰天為53%、91%、100%；模擬值與實(shí)際值的平均均方根誤差為1.7，多云天為2.3，晴天為1.6，陰天為1.2，見表5。白天3種天氣情況下的預(yù)測(cè)結(jié)果均較好，對(duì)比2個(gè)模型的檢驗(yàn)結(jié)果，BP神經(jīng)網(wǎng)絡(luò)遠(yuǎn)好于逐步回歸，白天溫室氣溫受多方面因素影響，變化規(guī)律不穩(wěn)定，與溫室外氣象條件成非線性關(guān)系，BP神經(jīng)網(wǎng)絡(luò)為非線性處理模型更適宜白天溫室內(nèi)氣溫預(yù)測(cè)。2.2 基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室相對(duì)濕度預(yù)測(cè)結(jié)果對(duì)比分析2.2.1 08時(shí)、1723時(shí)日光溫室相對(duì)濕度預(yù)測(cè)模型利用20112013年處理后數(shù)據(jù)應(yīng)用逐步回歸法建立不同天氣條件下傍晚到夜間日光溫室相對(duì)濕度預(yù)測(cè)模型，見表6。2014年數(shù)據(jù)作為檢驗(yàn)數(shù)據(jù)代入方程，得到模擬值并進(jìn)行檢驗(yàn)分析。BP神經(jīng)網(wǎng)絡(luò)也利用20112013年數(shù)據(jù)進(jìn)行訓(xùn)練建模，預(yù)測(cè)2014年溫室相對(duì)濕度并與實(shí)際值進(jìn)行檢驗(yàn)分析。08時(shí)、1723時(shí)逐步回歸模型晴天模擬值與實(shí)際值絕對(duì)誤差小于3的準(zhǔn)確率(6%)分別為91%和89%，均方根誤差為6.3和4.1；多云天模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)為93%和100%，均方根誤差為2.8和2.7；陰天模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)為97%和98%，均方根誤差為5和2.2。08時(shí)、1723時(shí)BP神經(jīng)網(wǎng)絡(luò)模型晴天模表4 不同天氣類型817時(shí)日光溫室氣溫預(yù)測(cè)逐步回歸模型注：x1為室外1 h前氣溫，x2為室外濕度，x3為室外風(fēng)速，x4為室外0 cm地溫，x5為室外5 cm地溫，x6為前1日室外最低氣溫，x7為總輻射。不同方法逐步回歸BP神經(jīng)網(wǎng)絡(luò)時(shí)段817時(shí)817時(shí)天氣多云晴陰平均值多云晴陰平均值1準(zhǔn)確率/%22184327.734535346.72準(zhǔn)確率/%45327149.375829182.73準(zhǔn)確率/%66509068.7889410094均方根誤差RMSE3.25.01.83.32.31.61.21.7樣本數(shù)8329914917783299149177時(shí)段08時(shí)1723時(shí)天氣類型多云晴陰多云晴陰日光溫室氣溫預(yù)測(cè)逐步回歸模型y1=98.46+0.02x5-0.58x7y2=99.48+0.15x1-0.08x7y3=97.7-0.2x1-0.39x3+0.12x4+0.03x5-0.11x7y4=96.26-0.42x1+0.13x4+0.03x5y5=94.97+0.16x1+0.07x2+0.08x5-0.11x6y6=92.26-0.44x1+0.05x2+0.03x5時(shí)段817時(shí)天氣類型多云晴陰日光溫室氣溫預(yù)測(cè)逐步回歸模型y7=15.06+0.17x1+0.39x4-0.19x6+0.02x7y8=19.2+0.22x1-0.04x2+0.23x4-0.22x6+0.02x7y9=13.64+0.22x1-0.01x2+0.2x3+0.26x4-0.18x6+0.02x7表5 817時(shí)逐步回歸與BP神經(jīng)網(wǎng)絡(luò)日光溫室氣溫預(yù)測(cè)模型檢驗(yàn)結(jié)果表6 不同天氣類型08時(shí)、1723時(shí)日光溫室相對(duì)濕度預(yù)測(cè)逐步回歸模型注：x1為室外1 h前氣溫，x2為室外濕度，x3為室外風(fēng)速，x4為室外0 cm地溫，x5為前1日室外相對(duì)濕度，x6為前2日室外相對(duì)濕度，x7為總輻射。121中國(guó)農(nóng)學(xué)通報(bào) http:/www.casb.org.cn擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)分別為92%和88%，均方根誤差為5.6和4.7；多云天模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)為99%和94%，均方根誤差為2.3和3.2；陰天模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)均為95%，均方根誤差為7.5和5.7，見表7。結(jié)合圖4和圖5可知，陰天模擬準(zhǔn)確率較高，擬合效果較好，好于其他2種天氣類型。2種模型不同精度下準(zhǔn)確率大小差異略小，但逐步回歸模型傍晚到夜間的均方根誤差幾乎都小于BP神經(jīng)網(wǎng)絡(luò)模型，模擬值與實(shí)際值的擬合效果也優(yōu)于BP神經(jīng)網(wǎng)絡(luò)。2.2.2 817時(shí)日光溫室相對(duì)濕度預(yù)測(cè)模型 表8列舉了白天20112013年處理數(shù)據(jù)建模的結(jié)果，采用以上方程預(yù)測(cè)2014年冬季1、2、12月同一天氣條件下相同時(shí)段溫室相對(duì)濕度值。利用相同時(shí)段相同室外氣象數(shù)據(jù)對(duì)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練建模，預(yù)測(cè)相同情況下的溫室相對(duì)濕度值。逐步回歸模型得到的模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)多云天、晴天、陰天分別為48%、圖4 不同天氣類型08時(shí)與1723時(shí)逐步回歸日光溫室相對(duì)濕度預(yù)測(cè)模型模擬值與實(shí)際值對(duì)比(a)0null時(shí)多云80901008090100模擬值/%實(shí)際值/%(b)0null時(shí)晴80901008090100模擬值/%實(shí)際值/%(c)0null時(shí)陰80901008090100模擬值/%實(shí)際值/%(d)17null3時(shí)多云80901008090100模擬值/%實(shí)際值/%(e)17null3時(shí)晴80901008090100模擬值/%實(shí)際值/%(f)17null3時(shí)陰80901008090100模擬值/%實(shí)際值/%08時(shí)多云 08時(shí)晴08時(shí)陰1723時(shí)多云1723時(shí)晴1723時(shí)陰122溫永菁等：基于逐步回歸與BP神經(jīng)網(wǎng)絡(luò)的日光溫室溫濕度預(yù)測(cè)模型對(duì)比分析34%、77%。模擬值與實(shí)際值的均方根誤差多云天、晴天、陰天分別為9.3、12.8、6.5。相應(yīng)的表9中BP神經(jīng)網(wǎng)絡(luò)模型模擬值與實(shí)際值絕對(duì)誤差小于6%的準(zhǔn)確率(6%)多云天、晴天、陰天分別為41%、43%、66%，均方根誤差多云天、晴天、陰天分別為10.6、12、8.6，見表9。2種模型的模擬準(zhǔn)確率均過低，均方根誤差也過大。說明817時(shí)日光溫室內(nèi)濕度變化規(guī)律易受澆水、施肥、通風(fēng)、打藥、開關(guān)棚門等一系列不確定因素影響，受室外環(huán)境條件影響較小，所建的2種模型均無法準(zhǔn)確預(yù)測(cè)溫室內(nèi)相對(duì)濕度。3 結(jié)論與討論（1）日光溫室氣溫受多種因素影響，要建立準(zhǔn)確且穩(wěn)定的預(yù)測(cè)模型有一定難度。本研究應(yīng)用逐步回歸和BP神經(jīng)網(wǎng)絡(luò)2種模型預(yù)測(cè)溫室內(nèi)氣溫，逐步回歸模型的白天晴好天氣預(yù)測(cè)準(zhǔn)確率較差，夜間陰雨天預(yù)測(cè)準(zhǔn)確率較高，預(yù)測(cè)模型穩(wěn)定性較差；BP神經(jīng)網(wǎng)絡(luò)模型具圖5 不同天氣類型08時(shí)與1723時(shí)BP神經(jīng)網(wǎng)絡(luò)日光溫室相對(duì)濕度預(yù)測(cè)模型模擬值與實(shí)際值對(duì)比(a)0null時(shí)多云80901008090100模擬值/%實(shí)際值/%(b)0null時(shí)晴80901008090100模擬值/%實(shí)際值/%(c)0null時(shí)陰80901008090100模擬值/%實(shí)際值/%(d)17null3時(shí)多云80901008090100模擬值/%實(shí)際值/%(e)17null3時(shí)晴80901008090100模擬值/%實(shí)際值/%(f)17null3時(shí)陰80901008090100模擬值/%實(shí)際值/%08時(shí)多云08時(shí)晴08時(shí)陰 1723時(shí)多云1723時(shí)晴1723時(shí)陰123中國(guó)農(nóng)學(xué)通報(bào) http:/www.casb.org.cn有處理非線性問題的能力，更適合溫室內(nèi)溫度預(yù)測(cè)建模，在不同天氣不同時(shí)段預(yù)測(cè)準(zhǔn)確率均較高且變化不大，其日光溫室氣溫預(yù)測(cè)模型準(zhǔn)確率較高穩(wěn)定性較好。（2）除去白天較晴好條件下，運(yùn)用BP神經(jīng)網(wǎng)絡(luò)與逐步回歸建立的日光溫室相對(duì)濕度預(yù)報(bào)模型模擬值與實(shí)際值絕對(duì)誤差小于等于6%的準(zhǔn)確率(6%)均達(dá)到94%，均方根誤差(RMSE)均在6%以內(nèi)。說明在傍晚到夜間陰天多云天，溫室內(nèi)濕度變化很規(guī)律，影響因子相對(duì)較少，運(yùn)用2種方法均可預(yù)測(cè)溫室內(nèi)相對(duì)濕度值。白天特別是晴好天氣，通風(fēng)時(shí)間較長(zhǎng)，農(nóng)事活動(dòng)較多，人為等其他因素干擾較大，因此外界氣象條件對(duì)溫室相對(duì)濕度模擬效果較差。（3）本研究選擇了冬季3個(gè)月溫室外數(shù)據(jù)進(jìn)行不同天氣類型不同時(shí)段的劃分，時(shí)段的劃分是主要考慮不同方法逐步回歸BP神經(jīng)網(wǎng)絡(luò)時(shí)段08時(shí)1723時(shí)08時(shí)1723時(shí)天氣多云晴陰平均值多云晴陰平均值多云晴陰平均值多云晴陰平均值2%準(zhǔn)確率/%72639275.735186840.380678075.746296747.34%準(zhǔn)確率/%9189969290559078.397869392877491846%準(zhǔn)確率/%93919793.7100899895.799929595.394889592.3均方根誤差RMSE2.86.35.04.72.74.12.232.35.67.55.13.24.75.74.5樣本數(shù)157753310406.7125577238313.3157753310406.7125577238313.3時(shí)段817時(shí)天氣類型多云晴陰日光溫室氣溫預(yù)測(cè)逐步回歸模型y7=91.58-x1-0.71x4-0.04x7y8=87.25-0.55x3-0.76x4-0.04x7y9=91.37+0.05x2-0.53x4+0.17x5-0.13x6-0.06x7不同方法逐步回歸BP神經(jīng)網(wǎng)絡(luò)時(shí)段817時(shí)817時(shí)天氣多云晴陰平均值多云晴陰平均值2%準(zhǔn)確率/%23114626.713153721.74%準(zhǔn)確率/%33226439.7323058406%準(zhǔn)確率/%4834775341436650均方根誤差RMSE9.312.86.59.510.612.08.610.4樣本數(shù)153710338400.3153710338400.3表7 08時(shí)、1723時(shí)逐步回歸日光溫室相對(duì)濕度預(yù)測(cè)模型檢驗(yàn)結(jié)果表8 不同天氣類型817時(shí)日光溫室相對(duì)濕度預(yù)測(cè)逐步回歸模型注：x1為室外1 h前氣溫，x2為室外濕度，x3為室外風(fēng)速，x4為室外0 cm地溫，x5為前一日室外相對(duì)濕度，x6為前2日室外相對(duì)濕度，x7為總輻射。表9 817時(shí)逐步回歸與BP神經(jīng)網(wǎng)絡(luò)日光溫室相對(duì)濕度預(yù)測(cè)模型檢驗(yàn)結(jié)果124