李永梅 張立根 海云端 董越

摘 要：利用2000～2014年統(tǒng)計數(shù)據(jù)，分析了寧夏近15 a來的耕地動態(tài)變化特征，采用主成分分析法提取影響耕地變化的主要驅(qū)動因子，將驅(qū)動因子作為輸入數(shù)據(jù)，構(gòu)建了預(yù)測耕地資源變化趨勢的BP神經(jīng)網(wǎng)絡(luò)，并開展預(yù)測研究。結(jié)果表明：寧夏耕地資源變化經(jīng)歷了增加—迅速減少—波動增長—迅速增長4個階段，人均耕地經(jīng)歷了迅速減少—緩慢減少—迅速增加3個階段；其中，2002～2004年和2010～2012年的土地利用動態(tài)度K<0，2003年動態(tài)度最?。↘=-7.673%），2005～2010年和2012～2014年的土地利用動態(tài)度K>0，2014年動態(tài)度最大（K=11.424%）；寧夏耕地資源主要驅(qū)動因子概括為農(nóng)業(yè)發(fā)展、經(jīng)濟(jì)發(fā)展及退耕還林，經(jīng)濟(jì)發(fā)展對耕地產(chǎn)生巨大壓力，而農(nóng)業(yè)科技進(jìn)步某種程度減小了人口對耕地的壓力，退耕還林也是寧夏耕地減少的重要原因，三者相互作用共同對寧夏耕地產(chǎn)生影響；預(yù)測寧夏耕地資源的BP神經(jīng)網(wǎng)絡(luò)模型為3層（3×16×1），對2011～2014年的耕地資源預(yù)測取得較好效果，最小誤差僅為376 hm2。

關(guān)鍵詞：耕地資源；BP神經(jīng)網(wǎng)絡(luò)；主成分分析；寧夏

中圖分類號：F301.21 文獻(xiàn)標(biāo)識碼：A 文章編號：1006-060X（2017）01-0081-05

Abstract：Using the statistical data from 2000 to 2014， the paper analyzed the characteristics of the dynamic change of cultivated land in Ningxia in recent 15 years， extracted the main driving factors which affected cultivated land change using principal component analysis， built a BP neural network which is used to forecast cultivated land resources and carry out the prediction research. The results showed that： Cultivated land resource in Ningxia experienced four phases including increasing， rapidly reduced， fluctuations in growth， rapid growth， and per capita cultivated land experienced rapidly reduce - slowly - rapidly increasing stage. Dynamic analysis of land use， dynamic attitude was less than zero from 2002 to 2004 and from 2010 to 2012， minimum value was 7.673% in 2003， dynamic attitude was greater than zero from 2005 to 2010 and from 2012 to 2014， maximum value was 11.424% in 2014. Major driving factors of cultivated land resources in Ningxia are summarized as agricultural development， economic development and returning farmland to forest， economic development has a huge pressure on arable land， agricultural science and technology progress to some extent reduce the pressure， and returning farmland to forest is also the reason for decrease of cultivated land in Ningxia. BP neural network model which is used for prediction of cultivated land resource in Ningxia is three layers（3×16×1）， the cultivated land resource forecasting from 2011 to 2014 using the model achieve better effect， least error is only 376 hm2.

Key words： cultivated land resource； BP neural network； principal component analysis； Ningxia

土地利用/土地覆被（LUCC）研究是全球環(huán)境變化研究領(lǐng)域的核心內(nèi)容。耕地作為土地的精華部分，是人類寶貴的自然資源。保持一定數(shù)量的耕地資源是一個國家和地區(qū)土地安全乃至生態(tài)安全的核心。目前，耕地動態(tài)變化已成為LUCC重要的研究內(nèi)容之一[1]。我國是農(nóng)業(yè)大國，耕地資源直接影響著十幾億人口的糧食安全及農(nóng)業(yè)可持續(xù)發(fā)展[2]。隨著社會經(jīng)濟(jì)的發(fā)展，我國耕地資源面臨著巨大壓力，已引起政府部門的高度重視。國內(nèi)學(xué)者對不同空間尺度耕地資源的動態(tài)變化及驅(qū)動力進(jìn)行了相關(guān)研究[3]，旨在為合理利用及保護(hù)有限的耕地資源提供依據(jù)。

多元統(tǒng)計分析[4-7]是目前耕地資源變化分析的常用方法。人工神經(jīng)網(wǎng)絡(luò)（Artifical Neutral Network，簡稱 ANN） 方法由于具有很強(qiáng)的聯(lián)想記憶及泛化推理能力，在預(yù)報非線性演變特征問題方面具有更好的能力[8]。筆者以寧夏為研究區(qū)域，將傳統(tǒng)的多元統(tǒng)計分析法和BP神經(jīng)網(wǎng)絡(luò)方法相結(jié)合，分析了近15 a來寧夏自治區(qū)的耕地動態(tài)變化特征，提取影響耕地變化的主要驅(qū)動因子，構(gòu)建了BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型，并對寧夏耕地資源的未來變化進(jìn)行了預(yù)測，以期為寧夏的社會經(jīng)濟(jì)與耕地資源協(xié)調(diào)發(fā)展提供參考，為農(nóng)業(yè)規(guī)劃及土地管理提供依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

寧夏位于東經(jīng)104°17′～107°39′，北緯35°14′～

39°23′之間，地處西北地區(qū)東北內(nèi)陸，北倚賀蘭山、南靠六盤山，與甘肅、內(nèi)蒙古和陜西省毗鄰。寧夏的氣候、地形地貌、土地資源等自然條件呈明顯帶狀分布，自北向南分為北部引黃灌區(qū)、中部干旱帶和南部山區(qū)。北部引黃灌區(qū)素有“塞上江南”之美譽(yù)，以灌溉農(nóng)業(yè)為主；中部干旱帶晝夜溫差大、光照充足、工業(yè)污染少，適宜發(fā)展有機(jī)特色農(nóng)業(yè)；南部黃土丘陵區(qū)溝壑縱橫， 氣候溫涼爽，適宜發(fā)展生態(tài)農(nóng)業(yè)[9]。寧夏的農(nóng)業(yè)資源獨(dú)具特色，發(fā)展?jié)摿薮?，全區(qū)有土地面積5.19×106 hm2，占全國土地面積的0.54%；2014年耕地面積為1.289×106 hm2，其中旱地約占總面積的64%，水澆地只占36%[10]。寧夏耕地總體質(zhì)量較差，限制因素較多，北部受土壤鹽漬化制約，中部受土壤貧瘠、風(fēng)大沙多制約，南部受土壤侵蝕制約[9]，因此科學(xué)的土地管理及種植業(yè)規(guī)劃顯得尤為重要。

1.2 數(shù)據(jù)來源

耕地資源短期內(nèi)主要受人類社會經(jīng)濟(jì)活動的影響。為了全面科學(xué)地對寧夏的耕地資源進(jìn)行分析預(yù)測，需要選用多個社會經(jīng)濟(jì)因子構(gòu)建一個指標(biāo)體系。只有指標(biāo)選取合理全面，分析結(jié)果才會準(zhǔn)確有效。根據(jù)研究區(qū)實(shí)際情況，選取了2000～2014年的社會經(jīng)濟(jì)指標(biāo)來分析寧夏耕地資源動態(tài)變化（表1）。社會經(jīng)濟(jì)數(shù)據(jù)和耕地數(shù)據(jù)均來源于《寧夏統(tǒng)計年鑒》（2001～2015年）。

1.3 研究方法

1.3.1 土地利用動態(tài)度 土地利用動態(tài)度可以定量地描述一定時間段內(nèi)，某種土地利用類型數(shù)量相對于前一時期土地利用數(shù)量的變化速度。模型為[11]：

式中，K為某時間段內(nèi)某種土地利用類型的變化率；Uib和Uia分別表示研究末期和研究初期時某種土地利用類型的數(shù)量；△T為研究時段長，一般以年為單位。

1.3.2 主成分分析法 主成分分析法就是在盡可能不損失信息或者少損失信息的情況下，將多個變量減少為少數(shù)幾個公共因子，這幾個因子可以揭示大部分?jǐn)?shù)據(jù)中的信息，這樣既減少了變量個數(shù)，避免了信息重疊，又能再現(xiàn)變量間的內(nèi)在聯(lián)系[11]。研究采用該方法從影響耕地資源的眾多社會經(jīng)濟(jì)因子中提取主要驅(qū)動因子。

1.3.3 BP神經(jīng)網(wǎng)絡(luò) BP神經(jīng)網(wǎng)絡(luò)采用誤差反向傳播（Back-Propagation）算法對網(wǎng)絡(luò)進(jìn)行訓(xùn)練，是目前應(yīng)用較廣較成功的人工神經(jīng)網(wǎng)絡(luò)，其中80%～90%的模型是采用BP網(wǎng)絡(luò)或它的變化形式[12]。其基本原理是，通過“信息順傳播”和“誤差逆?zhèn)鞑ァ眱蓚€過程交替反復(fù)進(jìn)行，不斷調(diào)整網(wǎng)絡(luò)的權(quán)值和閾值，使得誤差信號最小，最終使網(wǎng)絡(luò)的實(shí)際輸出逐漸向各自所對應(yīng)的期望輸出逼近，直到網(wǎng)絡(luò)輸出的誤差達(dá)到可接受的程度[12]。網(wǎng)絡(luò)中神經(jīng)元按輸入層、隱層及輸出層分層排列，輸入層和輸出層均為一層，隱層可以是一層或多層。層內(nèi)各神經(jīng)元不連接，層間各神經(jīng)元完全連接，其拓?fù)浣Y(jié)構(gòu)如圖1所示。

圖1中，Aj是輸入樣本，j=1，2，…m，m為輸入樣本個數(shù)；Qi是隱層的輸出向量，i==1，2，…s，s為隱層神經(jīng)元個數(shù)；Cp是輸出層，p=1，2，…n，n是輸出層神經(jīng)元個數(shù)。Qi、Cp分別以式和表示，其中，Wij為輸入層與隱層的連接權(quán)，Vpi為隱層與輸出層的連接權(quán)[13-14]。

2 結(jié)果與分析

2.1 耕地資源動態(tài)變化分析

由圖2可見，21世紀(jì)以來寧夏耕地經(jīng)歷了先增加—迅速減少—波動增長—迅速增長4個變化階段，人均耕地變化經(jīng)歷了迅速減少—緩慢減少—迅速增加3個階段。耕地面積最大值出現(xiàn)在2001年，為1.299×106 hm2；最小值出現(xiàn)在2004年，為1.105 ×106 hm2。人均耕地最大值為2000年的0.233 26×106 hm2，最小值為2012年的0.174 91×106 hm2。

2000～2001年，耕地面積由1.293×106 hm2增加至 1.299×106 hm2 ，凈增加了6 000 hm2；該時期盡管耕地在增加，但由于人口的增加，人均耕地呈現(xiàn)迅速減少的趨勢，該趨勢一直持續(xù)到2004年。2001～2004年，該時期耕地面積迅速減少，僅3 a時間耕地面積減少了1.94×105 hm2，年均減少量為4.85×104 hm2；2004～2014年，耕地面積由1.105×106 hm2增加至1.289×106 hm2，10 a時間增加耕地為1.84×105 hm2，年均增加量為1.84×104 hm2。2004～2013年間，由于耕地面積的小幅度增長，而且人口增長減慢，所以該時期人均耕地面積進(jìn)入緩慢減少階段；并且2013年之后，由于耕地面積的大幅度增長，使得人均耕地也出現(xiàn)迅速增加的趨勢。

由圖3可知，2002～2004年和2010～2012年耕地動態(tài)度為負(fù)數(shù)，說明以前一年的耕地面積為基準(zhǔn)，該時間段耕地面積呈現(xiàn)減少趨勢，而且2002～2004年間耕地面積減少速度更快，其中2003年動態(tài)度最?。↘= -7.673%），表明該年耕地面積減少最快。2005～2010年和2012～2014年的耕地動態(tài)度為正數(shù)，即與上一年耕地面積相比，其面積是增長的；其中，2014年的增長速度最快（K= 11.424 %）。

2.2 耕地資源驅(qū)動因子分析

以2000～2014年序列資料為基礎(chǔ)，以選取的13個社會經(jīng)濟(jì)因子作為樣本，對其進(jìn)行因子分析前的擬合檢驗(yàn)，結(jié)果顯示KMO=0.615，Bartlett=392.917，P<0.000 1，變量之間具有較強(qiáng)相關(guān)性，且相關(guān)矩陣不是單位矩陣，各指標(biāo)并非獨(dú)立，適合做主成分分析。

根據(jù)主成分分析法步驟，利用SPSS17.0數(shù)據(jù)分析軟件對樣本數(shù)據(jù)進(jìn)行主成分分析，得出標(biāo)準(zhǔn)化矩陣、相關(guān)系數(shù)矩陣、特征值及貢獻(xiàn)率（表2）。結(jié)果表明，前2個主成分累積貢獻(xiàn)率達(dá)90.309%，符合主成分分析中最小Em>85%的選取標(biāo)準(zhǔn)，前2個成分基本代表了所有驅(qū)動因子對耕地的影響，損失的信息僅占9.691%，這樣既剔除了冗余因子，使數(shù)據(jù)得到壓縮，同時又保留了最大的信息量。因此，提取了2個主成分，且進(jìn)一步計算出主成分載荷，如表2所示。

由表2可知，X2、X12與第一主成分之間具有很強(qiáng)的相關(guān)性，其載荷分別為0.993和0.991；X9與第二主成分相關(guān)性較強(qiáng)，載荷為0.805，明顯高于其他因子的貢獻(xiàn)。因此，寧夏耕地資源主要驅(qū)動因子可以概括為農(nóng)業(yè)發(fā)展因素、經(jīng)濟(jì)發(fā)展及退耕還林3類。

農(nóng)業(yè)發(fā)展、經(jīng)濟(jì)發(fā)展及退耕還林三者相互作用，共同影響著寧夏耕地資源的變化趨勢。經(jīng)濟(jì)發(fā)展促進(jìn)了工業(yè)化進(jìn)程和城市化水平，使寧夏的基本建設(shè)項(xiàng)目持續(xù)增加，進(jìn)而土地資源的需求量隨之增加，耕地資源面臨著巨大的壓力。同時，經(jīng)濟(jì)發(fā)展帶動農(nóng)業(yè)科技和耕作技術(shù)的進(jìn)步，減輕了人口對耕地資源的壓力。寧夏的土地墾殖率很高，坡地被不同程度地開墾，退耕還林面廣，也是寧夏耕地資源減少的原因之一。

2.3 BP神經(jīng)網(wǎng)絡(luò)模型設(shè)計與訓(xùn)練

在設(shè)計BP神經(jīng)網(wǎng)絡(luò)模型時，主要考慮網(wǎng)絡(luò)層數(shù)和各層神經(jīng)元個數(shù)[8]。由于至今未有指導(dǎo)神經(jīng)網(wǎng)絡(luò)構(gòu)建的完善理論，所以該研究在借鑒前人經(jīng)驗(yàn)的基礎(chǔ)上，采用試探法，通過不斷的試驗(yàn)和改進(jìn)，最終確定BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。

Kolmogorov定理中任一連續(xù)函數(shù)f：[0，1]I→RJ，f可以精確地用一個三層的BP神經(jīng)網(wǎng)絡(luò)來實(shí)現(xiàn)（其中， 為輸入層的神經(jīng)元個數(shù)， 為輸出層的神經(jīng)元個數(shù)）。該定理表明：單隱層的BP神經(jīng)網(wǎng)絡(luò)能映射所有的連續(xù)函數(shù)，只有當(dāng)隱層的神經(jīng)元個數(shù)已經(jīng)很多，網(wǎng)絡(luò)性能仍不能得到提高時，才會考慮增加隱層數(shù)量。同時，目前還沒有一個很好的解析式能用來計算隱層神經(jīng)元的個數(shù)，它的確定一般是根據(jù)經(jīng)驗(yàn)先設(shè)置較小數(shù)目對網(wǎng)絡(luò)進(jìn)行訓(xùn)練，之后逐漸增加隱節(jié)點(diǎn)數(shù)，通過反復(fù)的試驗(yàn)，最后選擇網(wǎng)絡(luò)誤差最小時對應(yīng)的隱層節(jié)點(diǎn)數(shù)。

研究中構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)實(shí)際上是通過輸入數(shù)據(jù)對輸出進(jìn)行預(yù)測，是一種具有函數(shù)逼近功能的網(wǎng)絡(luò)，因此使用單個隱層即可。隱層節(jié)點(diǎn)數(shù)首先根據(jù)Kolmogorov定理設(shè)置為2I+1，采用2000～2010年的農(nóng)業(yè)生產(chǎn)總值、國民生產(chǎn)總值和造林面積作為網(wǎng)絡(luò)輸入，耕地面積作為輸出，利用net=newff（minmax（p），[s1，s2]，{‘tansig，‘logsig}）構(gòu)建BP神經(jīng)網(wǎng)絡(luò)，其他都采用了默認(rèn)函數(shù)；用函數(shù)init（）初始化權(quán)值和閾值；網(wǎng)絡(luò)訓(xùn)練調(diào)用函數(shù)net=train（net，p，t）；性能函數(shù)為均方誤差“mse”；最大訓(xùn)練次數(shù)設(shè)置為

8 000次，學(xué)習(xí)速率設(shè)置為0.01和目標(biāo)誤差設(shè)置為0.01。采用同一樣本數(shù)據(jù)反復(fù)試驗(yàn)，最終確定用于預(yù)測寧夏耕地資源的BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為3層，共由20神經(jīng)元組成（3×16×1），隱層的神經(jīng)元個數(shù)為16。

2.4 BP神經(jīng)網(wǎng)絡(luò)的耕地預(yù)測

訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)不僅“儲存記憶”了網(wǎng)絡(luò)輸入輸出間的“函數(shù)”關(guān)系，而且具備了很強(qiáng)的“泛化”能力。因此，當(dāng)向網(wǎng)絡(luò)輸入訓(xùn)練時未見過的非樣本數(shù)據(jù)時，它也能輸出合適的映射，用于對寧夏耕地面積的預(yù)報。將寧夏2010～2014年共計5 a的農(nóng)業(yè)生產(chǎn)總值、國民生產(chǎn)總值和造林面積作為輸入數(shù)據(jù)，對耕地面積的變化進(jìn)行預(yù)測，結(jié)果見表3。從表3中可以看出，預(yù)測結(jié)果較為理性，最小誤差為376 hm2；同時，發(fā)現(xiàn)隨著預(yù)測時間的推移，預(yù)測誤差隨之增加。

3 結(jié)論與討論

3.1 結(jié) 論

利用2000～2014年統(tǒng)計數(shù)據(jù)，從耕地面積、人均耕地及動態(tài)度對寧夏近15 a來耕地面積動態(tài)變化特征進(jìn)行了分析，表明近15 a來寧夏耕地經(jīng)歷了先增加—迅速減少—波動增長—迅速增長4個變化階段，人均耕地變化經(jīng)歷了迅速減少—緩慢減少—迅速增加3個階段。 2002～2004年和2010～2012年耕地動態(tài)度K<0，2005～2010年和2012～2014年的K>0。

采用主成分分析法，從許多社會經(jīng)濟(jì)指標(biāo)中提取了2個主成分，其累計貢獻(xiàn)率為90.309%，結(jié)合主成分載荷得分確定農(nóng)業(yè)生產(chǎn)總值、國民生產(chǎn)總值和造林面積為主要影響因子，其載荷得分分別為0.993、0.991和0.805。因此，影響寧夏耕地資源變化的驅(qū)動因子概括為農(nóng)業(yè)發(fā)展因素、經(jīng)濟(jì)發(fā)展及退耕還林，三者相互作用共同影響著寧夏耕地資源的變化趨勢。

根據(jù)人工神經(jīng)網(wǎng)絡(luò)原理，采用BP算法，構(gòu)建的BP神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)為3層（3×16×1），共有20個神經(jīng)元，其中隱層神經(jīng)元為16個。采用2011～2014年統(tǒng)計數(shù)據(jù)對耕地進(jìn)行預(yù)測，取得較好的效果，其最小誤差僅為376 hm2。

3.2 討 論

耕地資源演變是多種驅(qū)動因素綜合作用的結(jié)果，由于其具有極強(qiáng)的經(jīng)濟(jì)屬性，該研究在構(gòu)建驅(qū)動因子指標(biāo)時只考慮了社會經(jīng)濟(jì)因素。在選取社會經(jīng)濟(jì)因子時雖然遵循了完備性和非重疊性原則，但因子間仍存在某種程度的信息重疊。采用主成分分析法，精選出有代表性的主要驅(qū)動因子，這樣抓住主要矛盾，保證分析的科學(xué)準(zhǔn)確，同時降低了訓(xùn)練樣本的維數(shù)，提高了網(wǎng)絡(luò)收斂速度。

耕地資源預(yù)測是典型的非線性問題，預(yù)測因子和預(yù)報量之間沒有確定的物理關(guān)系。神經(jīng)網(wǎng)絡(luò)模型無需建立數(shù)學(xué)模型，沒有把輸入因子和預(yù)測量間的非線性關(guān)系限制在多項(xiàng)式等函數(shù)上，進(jìn)而增強(qiáng)了預(yù)測的可靠性。同時，由于該方法無需人為確定權(quán)重，于是提高了預(yù)測結(jié)果的客觀性。

研究發(fā)現(xiàn)，隨著樣本數(shù)量的增加，網(wǎng)絡(luò)神經(jīng)學(xué)習(xí)能力和預(yù)報精度也隨之提高。今后在構(gòu)建和訓(xùn)練網(wǎng)絡(luò)時應(yīng)盡可能選取較多的樣本數(shù)據(jù)，同時采用一些智能算法對BP神經(jīng)網(wǎng)絡(luò)進(jìn)行改進(jìn)，加強(qiáng)預(yù)測水平的研究，使其在地區(qū)耕地資源預(yù)測方面得到更廣泛實(shí)際的應(yīng)用。

參考文獻(xiàn)：

[1]劉 夢. 寶雞市耕地面積變化及其驅(qū)動力分析[J]. 安慶師范學(xué)院學(xué)報（自然科學(xué)版），2015，21（4）：91-95.

[2]白小梅，李壁成，白均堂，等. 陜西省旬邑縣耕地面積變化及驅(qū)動力分析[J]. 水土保持研究，2010，17（5）：59-62，67.

[3]吳 業(yè)，楊桂山，萬榮榮. 耕地變化與社會經(jīng)濟(jì)發(fā)展關(guān)系研究進(jìn)展[J]. 地理科學(xué)進(jìn)展，2008，27（1）：90-98.

[4]李 博，何騰兵，陶福然. 貴州省耕地面積的時空變化及驅(qū)動力分析[J]. 貴州農(nóng)業(yè)科學(xué)，2013，41（3）：140-145.

[5]張 梅，吳 迪. 河南省耕地面積動態(tài)變化趨勢及驅(qū)動力分析[J]. 中州學(xué)刊，2012， 192（6）：51-54.

[6]郭雅雯，趙敏娟. 基于STIRPAT模型的陜西省耕地面積變化社會經(jīng)濟(jì)驅(qū)動力分析[J]. 水土保持研究，2012，19（6）：198-202.

[7]吳美瓊，陳秀貴. 基于主成分分析法的欽州市耕地面積變化及其驅(qū)動力分析[J]. 地理科學(xué)，2014，34（1）：54-59.

[8]金 龍，吳建生，林開平，等. 基于遺傳算法的神經(jīng)網(wǎng)絡(luò)短期氣候預(yù)測模型[J]. 高原氣候，2005，24（6）：980-986.

[9]趙亞峰，賈科利，文 琦. 寧夏耕地壓力動態(tài)變化規(guī)律分析[J]. 干旱區(qū)資源與環(huán)境，2014，28（1）：62-65.

[10]郁光磊，和 瑋，璩向?qū)? 寧夏耕地面積變化及驅(qū)動力分析[J]. 農(nóng)業(yè)科學(xué)研究，2007，28（4）：1-3.

[11]李永梅. 基于RS的杭州市城區(qū)擴(kuò)展及驅(qū)動力研究[D]. 金華：浙江師范大學(xué)，2011.

[12]張圣楠，郭文義，肖力墉. 基于MATLAB的BP神經(jīng)網(wǎng)絡(luò)的設(shè)計與訓(xùn)練[J]. 內(nèi)蒙古科技與經(jīng)濟(jì)，2005，（17）：95-98.

[13]韓力群. 人工神經(jīng)網(wǎng)絡(luò)理論、設(shè)計及應(yīng)用（第2版）[M]. 北京：化學(xué)工業(yè)出版社，2007.47-62.

[14]李永梅，張學(xué)儉，張立根. 基于BP神經(jīng)網(wǎng)絡(luò)的寧夏水資源需求量預(yù)測[J]. 水資源與水工程學(xué)報，2014，25（6）：98-101.

（責(zé)任編輯：成 平）