□文/劉尚俊 劉君
(塔里木大學(xué)經(jīng)濟(jì)與管理學(xué)院新疆·阿拉爾)
我國是紡織服裝業(yè)大國,棉花的穩(wěn)定供給對(duì)我國棉紡業(yè)意義重大。2011年新疆棉花種植面積2,393.9萬畝,產(chǎn)量289.8萬噸,連續(xù)19年保持面積、單產(chǎn)、總產(chǎn)、調(diào)出量全國第一。棉花產(chǎn)業(yè)發(fā)展的穩(wěn)定與否,不僅關(guān)系到國家棉花安全和棉紡工業(yè)穩(wěn)定發(fā)展,也關(guān)系到新疆農(nóng)村經(jīng)濟(jì)發(fā)展和社會(huì)的穩(wěn)定。而棉花產(chǎn)業(yè)的穩(wěn)定發(fā)展,離不開棉花產(chǎn)量的預(yù)測(cè)。分析與預(yù)測(cè)新疆棉花產(chǎn)量,不僅可以了解新疆棉花產(chǎn)量的走勢(shì),而且有助于統(tǒng)籌安排新疆棉花的種植、消費(fèi)、出口等相關(guān)事宜,從而穩(wěn)定棉花市場(chǎng)的供求,同時(shí)兼顧棉農(nóng)的利益。
棉花屬于純經(jīng)濟(jì)作物,受市場(chǎng)價(jià)格變化影響很大,因此棉花種植面積具有很大波動(dòng)性,同時(shí)氣候變化對(duì)于棉花產(chǎn)量的影響至關(guān)重要。因此,相比糧食作物來說,棉花產(chǎn)量的預(yù)測(cè)具有較大難度。
時(shí)間序列預(yù)測(cè)和灰色系統(tǒng)GM(1,1)等模型均是假設(shè)所有的影響因素都蘊(yùn)含在單一歷史序列中,主要依靠總產(chǎn)量數(shù)據(jù)建立預(yù)測(cè)模型,需要的數(shù)據(jù)較少,比較容易操作;適合于具有長(zhǎng)期趨勢(shì)的序列,對(duì)于波動(dòng)比較劇烈的序列預(yù)測(cè)效果較差。目前,對(duì)新疆棉花產(chǎn)量預(yù)測(cè)的研究大部分都是此類。
回歸分析預(yù)測(cè)和神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)屬于因果關(guān)系預(yù)測(cè),假定一個(gè)因素的變動(dòng)是由另一個(gè)或幾個(gè)變量引起的,通過掌握自變量的變動(dòng)可以知道因變量的變動(dòng)趨勢(shì)。要求占有盡可能多的資料,而對(duì)數(shù)列的波動(dòng)趨勢(shì)沒有特別要求?;貧w分析需要假設(shè)關(guān)系的數(shù)量模型形式,然后用最小二乘法擬合,而神經(jīng)網(wǎng)絡(luò)則不需要假設(shè)數(shù)量關(guān)系的形式,通過反復(fù)多次的學(xué)習(xí)和訓(xùn)練達(dá)到模擬變量關(guān)系的目的。本文擬選用BP神經(jīng)網(wǎng)絡(luò)對(duì)棉花總產(chǎn)量進(jìn)行預(yù)測(cè)。
(一)BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。BP網(wǎng)絡(luò)是一種按誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多層前饋網(wǎng)絡(luò),是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一。BP網(wǎng)絡(luò)能學(xué)習(xí)和存貯大量的輸入-輸出模式映射關(guān)系,而無需事前揭示描述這種映射關(guān)系的數(shù)學(xué)方程。BP神經(jīng)網(wǎng)絡(luò)模型拓?fù)浣Y(jié)構(gòu)包括輸入層、隱層和輸出層。當(dāng)有信息向網(wǎng)絡(luò)輸入時(shí),信息先由輸入層傳至隱層節(jié)點(diǎn),再傳至輸出層。
(二)BP神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法。以只包含一層隱層節(jié)點(diǎn)神經(jīng)網(wǎng)絡(luò)為例,其算法的步驟如下:
(1)設(shè)置變量和參量:
為輸入向量,M代表輸入向量維數(shù),N為訓(xùn)練樣本個(gè)數(shù)。
WMI(n)=(wij)M×I為第n次迭代時(shí)輸入層與隱層I間的權(quán)值向量;WIP(n)=(wij)I×P為第n次迭代時(shí)隱層I與輸出層之間的權(quán)值向量。
(k=1,2,…,N)為第n次迭代時(shí)網(wǎng)絡(luò)的實(shí)際輸出,P為輸出向量的維數(shù)。dk(n)=[dk1(n),dk2(n),…,dkP(n)],(k=1,2,…,N)為期望輸出。
(2)初始化。賦較小的隨機(jī)非零值于WMI(0),WIP(0)。
(3)輸入樣本Xk,令n=0。
(4)對(duì)輸入樣本Xk,前向計(jì)算BP網(wǎng)絡(luò)每層神經(jīng)元的輸入信號(hào)u和輸出信號(hào)v,將得出的實(shí)際輸出Yk(n)與期望輸出dk(n)相比較,判斷是否滿足要求,若滿足轉(zhuǎn)至(7),不滿足轉(zhuǎn)至(5)。
(5)判斷n是否達(dá)到最大迭代次數(shù),若大于轉(zhuǎn)至(7),否則對(duì)輸入樣本Xk,反向計(jì)算每層神經(jīng)元局部梯度δ。其中:δPp=
(6)按下式計(jì)算權(quán)值修正△W,并修正權(quán)值,η為學(xué)習(xí)速率。n=n+1,轉(zhuǎn)至(4)。
(7)判斷是否學(xué)完所有的訓(xùn)練樣本,是則結(jié)束,否則返回(3)。
(一)模型設(shè)定。以新疆棉花種植面積X1,有效灌溉面積X2,化肥使用量為輸入神經(jīng)元X3,以新疆棉花總產(chǎn)量y為輸出神經(jīng)元,利用以上算法建立BP神經(jīng)網(wǎng)絡(luò)模型,根據(jù)1995年到2010年的歷史數(shù)據(jù)進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練和預(yù)測(cè)模擬實(shí)驗(yàn)。
模型輸入變量:棉花種植面積(萬畝)X1;有效灌溉面積(萬畝)X2;機(jī)械擁有量(萬千瓦)X3。期望輸出:棉花總產(chǎn)量(萬噸)Y。
由于每個(gè)指標(biāo)的單位和變動(dòng)范圍不同,需要對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化(標(biāo)志值落在全區(qū)間[–1,1])。標(biāo)準(zhǔn)化方法:X'=
表1 神經(jīng)元網(wǎng)絡(luò)訓(xùn)練結(jié)果與回歸分析結(jié)果比較表
模型設(shè)定一個(gè)隱層。隱層節(jié)點(diǎn)數(shù)量根據(jù)經(jīng)驗(yàn)設(shè)定為5,學(xué)習(xí)速率設(shè)定0.01,目標(biāo)精度0.0001,最大訓(xùn)練次數(shù)50000。
(二)模型運(yùn)行與分析。選用多元線性回歸模型與神經(jīng)網(wǎng)絡(luò)模型進(jìn)行比較,運(yùn)行結(jié)果見表1。(表1)其中,1995~2009年數(shù)據(jù)為樣本擬合值,2010年數(shù)據(jù)為樣本預(yù)測(cè)值。將神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)值與多元線性回歸預(yù)測(cè)值進(jìn)行比較,其中:相對(duì)誤差=(預(yù)測(cè)值–實(shí)際值Y)/實(shí)際值Y。
人工神經(jīng)網(wǎng)絡(luò)BP模型能自主地從訓(xùn)練實(shí)例中學(xué)習(xí),建立自變量與因變量之間的非線性映射關(guān)系。與線性回歸預(yù)測(cè)一樣,指標(biāo)的選擇對(duì)神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)和預(yù)測(cè)有較大的影響,應(yīng)結(jié)合實(shí)際問題對(duì)預(yù)測(cè)指標(biāo)進(jìn)行適當(dāng)?shù)娜∩?,力求所建立的神?jīng)網(wǎng)絡(luò)預(yù)測(cè)模型有較高的預(yù)測(cè)精度。鑒于棉花種植成本連年提高,加之新疆尤其是南疆水資源匱乏,新疆維吾爾自治區(qū)政府2008年提出減棉增果號(hào)召,因此2009年和2010年新疆棉花種植面積和產(chǎn)量明顯萎縮。所以,進(jìn)行新疆棉花產(chǎn)量預(yù)測(cè)時(shí)應(yīng)綜合考慮宏觀政策、棉花價(jià)格等諸因素,從而進(jìn)一步提高預(yù)測(cè)精度。
[1]2011年我區(qū)農(nóng)業(yè)農(nóng)村工作大事盤點(diǎn):新疆農(nóng)業(yè)信息網(wǎng).http://www.xj-agri.gov.cn/.
[2]肖雙喜,劉小和.棉花種植面積影響因素分析——基于新疆、河南、江蘇和山東四省的調(diào)查[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2008.4.
[3]楊偉靜.基于BP神經(jīng)網(wǎng)絡(luò)的鐵路客流量預(yù)測(cè)研究[J].合作經(jīng)濟(jì)與科技,2010.13.
[4]朱輝.GM(1,1)模型在新疆生產(chǎn)建設(shè)兵團(tuán)棉花總產(chǎn)量預(yù)測(cè)中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué),2008.20.