關(guān)鍵詞:大氣環(huán)境質(zhì)量評價;軟測量;神經(jīng)網(wǎng)絡(luò);倍斜率聚類法;等效數(shù)值法 余丙榮(1971-)
男,安徽宿松人,工學(xué)學(xué)士、安徽大學(xué)電子與通信工程碩士生,研究方向?yàn)殡娏﹄娮印⒖刂评碚摗?BR>
1 引言
目前,環(huán)境污染已經(jīng)成為了一個嚴(yán)峻的問題,特別是城市、工業(yè)區(qū),環(huán)境污染甚是嚴(yán)重。大氣污染也是環(huán)境污染問題的重要部分,現(xiàn)在幾乎每個城市每天都要向環(huán)境部門報告其空氣污染情況和空氣污染指數(shù),以確定該城市的工業(yè)生產(chǎn)和環(huán)境治理力度是否合格。所以對空氣環(huán)境的質(zhì)量評價即空氣污染指數(shù)的計算是不可缺少的部分,不僅要求準(zhǔn)確,而且還要快速。近年來,在環(huán)境質(zhì)量綜合評價方面已提出了很多評價模型,各有各的優(yōu)缺點(diǎn)。本文將用神經(jīng)網(wǎng)絡(luò)的方法來對空氣環(huán)境質(zhì)量作評價,通過應(yīng)用于大氣環(huán)境質(zhì)量評價的具體實(shí)例來與兩個經(jīng)典評價模型倍斜率聚類法、等效數(shù)值法作比較并分析其評價的結(jié)果。
2 軟測量方法概述
軟測量技術(shù)是通過選擇一些容易測量且有關(guān)的變量(輔助變量),構(gòu)造一個與輸出建立聯(lián)系的數(shù)學(xué)模型,且用計算機(jī)軟件實(shí)現(xiàn)估計出輸出變量(主導(dǎo)變量)的方法。建立軟測量的方法有很多,在經(jīng)驗(yàn)建模法中,介紹一下人工智能的方法——人工神經(jīng)網(wǎng)絡(luò)。
圖1 徑向基神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)
污染物 Ⅰ Ⅱ Ⅲ Ⅳ SO2 0-0.05 0.05-0.15 0.15-0.25 0.25-∞ NOx 0-0.10 0 -0.10 0.10-0.15 0.15-∞ TSP 0-0.12 0.12 -0.30 0.30-0.50 0.50- ∞ 分級 SO2 NOx TSP Ⅱ 0.15 0.10 0.12 Ⅱ 0.15 0.10 0.30 Ⅳ 0.25 0.15 0.50
2.1 RBF神經(jīng)網(wǎng)絡(luò)模型
RBF神經(jīng)網(wǎng)絡(luò)也叫徑向基神經(jīng)網(wǎng)絡(luò),它是一種局部逼近的神經(jīng)網(wǎng)絡(luò)。其結(jié)構(gòu)由三層組成:輸入層、隱含層和輸出層,如圖1所示,隱含層第i個節(jié)點(diǎn)的輸出表示為:
(1)
式中,ui是第個i節(jié)點(diǎn)的輸出,是第i個節(jié)點(diǎn)的標(biāo)準(zhǔn)化常數(shù),q是隱含層的節(jié)點(diǎn)數(shù),x=(x1,x2,...,xM)T是樣本輸入;Ci是第i個節(jié)點(diǎn)高斯函數(shù)的中心向量,此向量是一個與輸入樣本x的維數(shù)相同的列向量,即Ci=(Ci1,Ci2,...,CiM)T 由式(1)可知,節(jié)點(diǎn)的輸出范圍在0和1之間,且輸入樣本越靠近節(jié)點(diǎn)中心,輸出值愈大。當(dāng)x=Ci時, Ui=1。
從隱含層到輸出層即Ui(x)→yk為線性映射,表示式為:
(2)
式中,ui是隱含層第i個節(jié)點(diǎn)的輸出;yk是輸出層第k個節(jié)點(diǎn)的輸出;wki是隱含層到輸出層的加權(quán)系數(shù);θk是輸出層的閾值;q是隱含層節(jié)點(diǎn)數(shù)。
2.2 RBF神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)過程
設(shè)訓(xùn)練樣本有N個,則系統(tǒng)的總誤差表示式為:
(3)
式中,N為模式樣本對數(shù);L為網(wǎng)絡(luò)輸出節(jié)點(diǎn)數(shù);表示在樣本p的作用下的第k個神經(jīng)元的期望輸出;
表示在樣本p的作用下的第k個神經(jīng)元的實(shí)際輸出。
其學(xué)習(xí)過程分為兩個階段:無教師學(xué)習(xí)階段和有教師學(xué)習(xí)階段。
2.2.1 無教師學(xué)習(xí)階段
可以通過k-均值聚類算法來對所有樣本輸入進(jìn)行聚類,求得各隱含層節(jié)點(diǎn)的RBF的中心向量,其算法如下:
(1)給定各隱節(jié)點(diǎn)初始中心向量Ci(0)和判定停止計算的;
(2)計算歐式距離,并求最小距離節(jié)點(diǎn);
(4)
式中,k為樣本序號,r為中心向量Ci(k-1)與輸入樣本x(k)距離最近的隱節(jié)點(diǎn)序號;
(3)調(diào)整中心向量
(5)
式中,是學(xué)習(xí)速率。
為取整運(yùn)算。每經(jīng)過q個樣本之后,調(diào)小一次學(xué)習(xí)速率,逐漸減至零。
(4)判斷聚類質(zhì)量
(6)
對于全部樣本反復(fù)進(jìn)行以上運(yùn)算,直到以上條件或聚類結(jié)束。
2.2.2 有教師學(xué)習(xí)階段
當(dāng)Ci確定以后,求隱含層至輸出層之間的權(quán)值Wki(k=1,2,...,L;i=1,2,...,q)就成了線性優(yōu)化問題,其權(quán)值的學(xué)習(xí)算法為:
(7)
式中,,ui(x)為高斯函數(shù);
為學(xué)習(xí)速率,且0<
<1;tk和yk分別表示第k個輸出分量的期望值和實(shí)際值。由于向量u只要少量幾個元素為1,故需要調(diào)整量少,學(xué)習(xí)速度比較快。
3 軟測量技術(shù)在大氣環(huán)境質(zhì)量評價上的應(yīng)用
表1和表2為國家大氣環(huán)境質(zhì)量分級標(biāo)準(zhǔn)(GB3095-1996)以及本文要用到的分級標(biāo)準(zhǔn)濃度限值。
表1 大氣環(huán)境質(zhì)量分級標(biāo)準(zhǔn) (mg/m3)
表2 分級標(biāo)準(zhǔn)濃度限值 (mg/m3)
這里假設(shè)環(huán)境質(zhì)量標(biāo)準(zhǔn)為連續(xù)的數(shù),表2表示當(dāng)污染物開始達(dá)到限度值時,分級標(biāo)準(zhǔn)開始達(dá)到標(biāo)準(zhǔn)值。分級指標(biāo)=(int)環(huán)境質(zhì)量標(biāo)準(zhǔn),即取整。下面以各污染物濃度作為該軟測量技術(shù)方法的輔助變量,以分級標(biāo)準(zhǔn)作為主導(dǎo)變量通過MATLAB軟件來建立評價模型,各設(shè)計參數(shù)取默認(rèn)值,其結(jié)果如圖2所示。
圖2 環(huán)境質(zhì)量模型的擬合評價效果
其輸出層網(wǎng)絡(luò)訓(xùn)練的誤差曲面如圖3所示。
圖3 訓(xùn)練誤差曲面圖
表3為某市大氣污染指標(biāo)實(shí)測值。
表3 某市大氣污染指標(biāo)實(shí)測值
監(jiān)測點(diǎn) |
A1 |
A2 |
A3 |
A4 |
A5 |
A6 |
A 7 |
A8 |
A9 |
A10 |
SO2 |
0.125 |
0.022 |
0.131 |
0.097 |
0.067 |
0.014 |
0.176 |
0.029 |
0.1 |
0.038 |
NOx |
0.086 |
0.013 |
0.016 |
0.027 |
0.035 |
0.018 |
0.123 |
0.026 |
0.075 |
0.038 |
TSP |
0.239 |
0.188 |
0.101 |
0.417 |
0.657 |
0.409 |
0.415 |
0.274 |
0.225 |
0.462 |
根據(jù)上面已建立好的評價模型,對此市大氣污染指標(biāo)進(jìn)行評價,并與倍斜率聚類法、等效數(shù)值法列表作比較,其結(jié)果如表4所示。
表4 各軟測量方法的結(jié)果比較
監(jiān)測點(diǎn) |
A1 |
A2 |
A3 |
A4 |
A5 |
A6 |
A7 |
A8 |
A9 |
A10 |
評價數(shù)值結(jié)果 |
2.6702 |
2.1333 |
1.9955 |
3.2580 |
3.9199 |
3.0081 |
3.5391 |
2.5266 |
2.5428 |
3.2877 |
評價數(shù)值指標(biāo) |
Ⅱ |
Ⅱ |
Ⅰ |
Ⅲ |
Ⅲ |
Ⅲ |
Ⅲ |
Ⅱ |
Ⅱ |
Ⅲ |
倍余率聚類法 |
Ⅱ |
Ⅱ |
Ⅱ |
Ⅱ |
Ⅲ |
Ⅱ |
Ⅲ |
Ⅱ |
Ⅱ |
Ⅱ |
等效數(shù)值法 |
Ⅱ |
Ⅱ |
Ⅱ |
Ⅲ |
Ⅲ |
Ⅲ |
Ⅲ |
Ⅱ |
Ⅱ |
Ⅲ |
從表4可以看出,用軟測量技術(shù)的神經(jīng)網(wǎng)絡(luò)技術(shù)作評價,其結(jié)果與倍斜率聚類法、等效數(shù)值法的結(jié)果基本上一致,特別是與等效數(shù)值法,基本一樣。從計算過程來看,軟測量技術(shù)計算非常簡捷,只需要輸入數(shù)據(jù),通過計算機(jī)軟件,立即得到結(jié)果,并從數(shù)據(jù)中可以看出其偏向前一級或后一級的趨向。而倍斜率聚類法等中隸屬函數(shù)的建立以及計算相對復(fù)雜。
根據(jù)神經(jīng)網(wǎng)絡(luò)這一軟測量方法所建立的評價模型,可以確定其輔助變量與主導(dǎo)變量的關(guān)系,即y=f(x), x=[x1,x2,...,xp]T 為輔助變量,y=[y1,y2,...,yt]T為主導(dǎo)變量。本例中它們之間的函數(shù)關(guān)系為:tx=simurb(px,w1,b1,w2,b2),px為輔助變量,tx為主導(dǎo)變量。
4 結(jié)束語
用神經(jīng)網(wǎng)絡(luò)的軟測量方法對大氣環(huán)境質(zhì)量評價,因?yàn)橐?jīng)過多次訓(xùn)練,其評價誤差可達(dá)到很小,本例為10-16以下,而且通過神經(jīng)網(wǎng)絡(luò)的方法,可以很方便地得到輸入與輸出的函數(shù)關(guān)系。使用方便,在環(huán)境評價中的應(yīng)用前景廣闊。
參考文獻(xiàn):
[1] 李澤應(yīng),曹菁菁等. RBF 神經(jīng)網(wǎng)絡(luò)及其在水質(zhì)評價中的應(yīng)用.兵工自動化,2006.7.
[2] 阮沈勇,王永利. MATLAB程序設(shè)計.北京:電子工業(yè)出版社,2004.1.
[3] 李國勇.智能控制及其MATLAB實(shí)現(xiàn).北京:電子工業(yè)出版社,2005.5.
[4] 喬敬萍.幾種環(huán)境質(zhì)量評價方法的應(yīng)用.山西大學(xué)學(xué)報(自然科學(xué)版),2004.