殷鵬遠(yuǎn)

(遼寧省錦州水文局，遼寧 錦州 121000)

灰色預(yù)測系統(tǒng)GM(1,1)模型及其Matlab實(shí)現(xiàn)

殷鵬遠(yuǎn)

(遼寧省錦州水文局，遼寧 錦州 121000)

灰色模型有嚴(yán)格的理論基礎(chǔ)，最大的優(yōu)點(diǎn)是實(shí)用，用灰色模型預(yù)測的結(jié)果比較穩(wěn)定，不僅適用于大數(shù)據(jù)量的預(yù)測。目前，灰色模型GM(1,1)已廣泛應(yīng)用于工程技術(shù)、社會、經(jīng)濟(jì)、農(nóng)業(yè)、生態(tài)、環(huán)境等各種系統(tǒng)的預(yù)測中。文章根據(jù)所建立的GM(1,1)模型以及模型分析，來預(yù)測未來安陽市旱災(zāi)發(fā)生年份，以期積極主動地采取措施進(jìn)行防旱抗旱工作提供科學(xué)依據(jù)。

灰色系統(tǒng)；GM(1,1)模型；Matlab；災(zāi)變預(yù)測

0 引 言

河南的水旱災(zāi)害經(jīng)過60a的大力治理，已經(jīng)初步形成了防御旱災(zāi)的工程體系，大面積的經(jīng)常性的旱災(zāi)已經(jīng)明顯減輕，對農(nóng)業(yè)的持續(xù)穩(wěn)定增產(chǎn)和城市、工業(yè)企業(yè)的安全起到了保障作用，社會效益、經(jīng)濟(jì)效益顯著。但是，在人口持續(xù)增長和經(jīng)濟(jì)迅速發(fā)展的情況下，對方于水旱災(zāi)害的要求將越來越高。展望2017年以后相當(dāng)長的時期內(nèi)，河南防御旱災(zāi)面臨的形式和問題將是嚴(yán)重的，需要繼續(xù)努力。

1 應(yīng)用 GM(1,1)模型對安陽旱災(zāi)發(fā)生年份的建模與預(yù)測

1.1 GM(1,1)模型的建立

GM(1,1)模型就是建立一個一階的微分方程模型，適合于一個變量的預(yù)測。建模過程中,首先將原始數(shù)據(jù)累加生成新的序列，使之有規(guī)律，成為一個遞增數(shù)列，然后利用生成后的數(shù)列進(jìn)行建模，預(yù)測時再還原成原數(shù)列，以恢復(fù)事物的原貌。

1)獲取原始數(shù)據(jù)矩陣:

X(0)=(Xi(0))N，n

(1)

2)對Xi(0)作1-AGO；

3)對X(0)作準(zhǔn)光滑性檢驗(yàn)，由

(2)

存在k0，使得ρ(k0)<0.5，則當(dāng)k>k0時準(zhǔn)光滑條件滿足；

4)檢驗(yàn)X(1)是否具有準(zhǔn)指數(shù)規(guī)律有

(3)

當(dāng)k>k0時，σ(1)(k)[1,1.5]，δ=0.5，準(zhǔn)指數(shù)規(guī)律滿足，故可對X(1)建立GM(1,N)模型。

5)對X(1)作緊鄰均值生成：

Z(1)=(z(1)(2),z(1)(3),…,z(1)(n))

(4)

(5)

7)確定模型：

(6)

及近似時間響應(yīng)式：

(7)

式中：x1(1)(0)取為x1(0)(1)。

8)累減還原式為：

(8)

9)檢驗(yàn)誤差：

殘差：ε(k)=x1(0)(k)-＾x0(k)

算出殘差平方和：s=εTε

1.2GM(1,1)模型的求解

根據(jù)上述建立的GM(1,1)模型，以豫西重旱區(qū)代表城市洛陽市為例，根據(jù)表1(洛陽市旱災(zāi)年份統(tǒng)計(jì)表)中的統(tǒng)計(jì)數(shù)據(jù)，通過對洛陽市旱災(zāi)年份進(jìn)行研究分析，經(jīng)Matlab軟件求解(具體程序見附錄：程序1)得到時間響應(yīng)序列函數(shù)：

(9)

做出殘差分析對照表(如下表4)，并求出殘差平方和s=215989.4

表4 殘差分析對照表

通過殘差檢驗(yàn)可以得到：

平均相對誤差為：

(10)

說明所建模型具有可靠性。

為進(jìn)一步檢驗(yàn)?zāi)Ｐ偷恼_性，通過計(jì)算機(jī)仿真進(jìn)行深入研究，結(jié)果見圖1：

圖1 原始序列與模擬序列比照圖

通過上圖比照可以發(fā)現(xiàn)，實(shí)測值對模擬值具有較好的擬合，模擬效果良好。

綜上殘差分析以及仿真模擬結(jié)果，認(rèn)為所建立的GM(1,1)模型具有較高的精度。

(11)

(12)

(11)- (12)得：

(13)

(13)式可以用來預(yù)測安陽市旱災(zāi)年份。

1.3GM(1,1)模型的應(yīng)用

根據(jù)所建立的模型以及模型分析，來預(yù)測未來安陽市旱災(zāi)發(fā)生年份，以期積極主動地采取措施進(jìn)行防旱抗旱工作提供科學(xué)依據(jù)。用

(14)

進(jìn)行預(yù)測，可以得到：

2 結(jié) 語

1)在對干旱災(zāi)變問題研究中發(fā)現(xiàn)，此類問題的解法不止于文中提到的幾種，例如：馬爾科夫模型等。在實(shí)際的工作和學(xué)習(xí)中，我們可以根據(jù)實(shí)際情況，根據(jù)研究問題的性質(zhì)及要求，選用合理的模型分析。

2)在建立GM(1,1)模型分析問題的基礎(chǔ)上，而由模型對比部分分析得到結(jié)論：GM(1,1)在一定程度上具有局限性。而改進(jìn)的方法主要有：殘差GM(1 ,1) 模型、無偏灰色模型、參數(shù)優(yōu)化灰色模型、新陳代謝GM(1,1)模型等，這些模型都在不同場合下對GM(1,1)模型進(jìn)行了一定程度的改進(jìn)。

因此，可以進(jìn)一步地建立改進(jìn)的GM(1,1)模型，使得預(yù)測結(jié)果更加合理，這有待進(jìn)一步研究。

[1]河南省水利廳水旱災(zāi)害專著編輯委員會.河南水旱災(zāi)害[M].鄭州:黃河水利出版社,1999：19-21.

[2]河南省水文總站編.河南省歷代旱澇等水文氣候史料[M].鄭州:河南省水文總站，1982：27-33.

[3]河南省水文總站.河南省歷代大水、大旱年表[M].鄭州:河南省水文總站，1999：26-36.

[4]WashingtonRN.EstimatingandinterpretinghydrologicaldroughtindicesusingselectedcatchmentsinZimbabwe[J].PhysicsandChemistryoftheEarth, 2004, 29(02):1173-1180.

[5]張書濱.江西省農(nóng)業(yè)旱情預(yù)測模型的建立與應(yīng)用[J].南昌大學(xué)學(xué)報，2009(06):179-182.

[6]楊寶中.河南省干旱特點(diǎn)及水資源開發(fā)利用的研究[J].華北水利水電學(xué)院學(xué)報，2008，29(04)：1-4.

[7]房穩(wěn)靜.河南省冬小麥干旱災(zāi)害的空間分布研究[J].氣象與環(huán)境科學(xué)，2007(02):61-63.

[8]程炳巖．河南氣候概論[M]．北京:氣象出版社,1995：26-37.

[9]羅佑新.灰色系統(tǒng)理論及其在機(jī)械工程中的應(yīng)用[M].長沙:國防科技大學(xué)出版社, 2001：44-45.

[10]張雅君,劉全勝.城市需水量灰色預(yù)測的探討[J].中國給水排水,2002,18(03):28230.

GM(1,1) Model and Matlab Application of Gray Prediction System

YIN Peng-yuan

(Liaoning Provincial Jinzhou Hydrological Bureau, Jinzhou 121000, China)

Gray model owns strict theatrical basis, the biggest advantage is practical, the conclusions predicted by gray model is more stable, not only suitable for the current prediction of mass data quantities,but also applies widely in project technology, society, economy, agriculture, ecology and environment etc,.According to GM(1,1) model established and model analysis, this paper predicts the years of disasters will happen in Anyang City in the future, supply scientific accordance to adopt actively safe measures to control the draught disasters.

gray system; GM (1,1)model; Matlab; prediction of disasters

1007-7596(2017)07-0016-03

2017-06-20

殷鵬遠(yuǎn)(1990-)，男，河南許昌人，助理工程師，從事水情預(yù)報、水文勘測、水文資料整編工作。

S162.3

B