李會超 李 鴻 張 博
(長沙理工大學電氣與信息工程學院 湖南 長沙 410004)
在城市化高速推進的過程中,交通問題愈發成為制約城市發展的一個關鍵因素。作為一種能有效緩解城市交通問題的手段,智能交通系統ITS(Intelligent Transportation Systems)已經普遍應用于中大型城市交通系統中。短時交通流量的預測是交通控制與交通誘導等子系統的前提和關鍵,因此短時交通流量預測的準確性和實時性始終是學者研究ITS的重點[1]。
微觀層面上將間隔15 min以內的車流量預測稱之為短時交通流量預測。較早的預測模型如移動平均模型、歷史平均模型,其大量線性評估的理論基礎不需要建立精確的數學模型。但交通流量的時變性、非線性和高復雜性使得較早模型對交通流量的預測存在很大誤差。目前,廣泛應用的短時交通流量預測模型主要有以下三種:時間序列模型、非參數回歸模型和神經網絡模型[2]。其中神經網絡因其良好的非線性品質和自組織學習能力在預測領域得到廣泛應用,典型的代表就是小波神經網絡[3]。
在1992年,Zhang等[4]提出了小波神經網絡WNN(Wavelet Neural Network),隨后便成為了熱門的建模分析方法。小波神經網絡既繼承了BP神經網絡的網絡結構,又結合了小波分析理論,用小波基函數代替原有網絡中隱含層節點的Sigmoid函數進行函數映射,改變了神經網絡對非線性模型擬合的基函數[5]。小波神經網絡綜合了小波分析和神經網絡在處理非線性問題上的優勢,很好地解決了BP神經網絡結構設計的盲目性和易陷入局部最優的問題[6]。
小波神經網絡相對于其他網絡雖然具有更好的學習能力,但是小波參數初值設定直接影響網絡的逼近能力,如果網絡參數確定不好,將造成如下三種現象:1) 易陷入局部最優,誤差始終無法達到誤差限;2) 收斂速度慢,增加算法的搜索時間;3) 引起誤差曲線震蕩,不能保證誤差一直趨于最小。 因此,學者提出不少改進算法,文獻[7-8]采用粒子群算法優化WNN,但面臨粒子群算法易早熟和局部尋優能力差等缺點;文獻[9]采用蟻群算法優化WNN,但面臨蟻群算法全局搜索能力差,計算量大等缺點。鑒于遺傳算法有著較好的全局搜索能力,本文采用遺傳算法優化小波神經網絡的初始參數。
小波神經網絡分為松散型與融合型兩類,前者將小波分析作為神經網絡的前置手段,對網絡輸入初步優化,使得信息便于神經網絡的處理,但這類WNN并未改變隱含層的傳遞函數,依然存在神經網絡的固有缺陷;后者采用小波元替代神經元,用小波函數的伸縮、平移參數替代輸入層到隱含層的權值及隱含層閾值,是小波和神經網絡的直接融合。本文研究采用融合性WNN,其網絡拓撲結構如圖1所示。

圖1 融合性小波神經網絡拓撲結構圖
圖1中各個字母代表含義如下:
x1,x2,…,xM:輸入層數據;
wj,k:輸入層到隱含層的權值參數;
ψa,b(x):小波基函數,本文采用Morlet函數,其表達式如下:
(1)
a、b分別為小波基函數的伸縮參數和平移參數;
wi,j:隱含層到輸出層的權值參數;
σ(I-t):輸出層Sigmoid函數,其表達式如下:
(2)
I、t分別為輸出的輸入和閾值。
y1,y2,…,yN:輸出層數據。
為求得WNN的穩定結構,首先借助訓練集數據進行學習,其次借助測試集數據對網絡結構驗證。在學習過程的第一階段,將經隱含層處理的樣本數據作為輸出向前傳播,然后與期望輸出對比,若不匹配,則進入第二階段。此階段將輸出誤差以某種形式反向傳至各層神經元,從而獲取調整各個神經元權值的誤差信號。經過兩個階段重復迭代,修正WNN網絡的參數。在訓練中,以輸出誤差限定值或設定的迭代次數作為終止條件。整體思想如以下偽代碼所示:
While不滿足終止條件 {
遍歷訓練集中的每個樣本:
#1 將每個樣本輸入沿著網絡前向傳播計算
輸入層到隱含層:
隱含層輸入表達式:
(3)
隱含層輸出表達式:
(4)
隱含層到輸出層:
輸出層輸入表達式:
(5)
輸出層輸出表達式:
yn=sigmoid(In-tn)
(6)
#2 使用誤差項后項傳播
誤差函數:
(7)

輸出層權重更新表達式:
(8)

輸出層閾值調整表達式:
(9)

隱含層權值調整表達式:
(10)

伸縮參數調整表達式:
(11)

平移參數調整表達式:
(12)

式(8)-式(12)中η均為小波神經網絡的學習效率。
}
遺傳算法GA(Genetic Algorithm)是借鑒生物界優勝劣汰的進化規律演化而來的一種并行隨機搜索算法,其主要思想是把問題參數編碼為染色體,對染色體進行一系列運算。即由選擇算子獲得優良個體遺傳給下一代;利用交叉算子進行模式重組;利用變異算子進行模式突變。經過循環這些遺傳操作,種群向有利的方向進化,最終得到所求問題最優解。遺傳算法的隨機搜索能力,使其能夠避免陷入局部最優,且有著較好的收斂速度,可以在小波參數優化問題上充分發揮優勢。
遺傳算法優化的主要要素:
(1) 編碼設置 由于待優化參數較多,文中遺傳算法采用實數編碼,通過對種群中具有某種形式的個體進行遺傳操作,獲取群體中最優個體,然后通過解碼得到網絡最佳參數。
(2) 適應度函數 遺傳算法在迭代過程中將適應度函數作為個體優劣的評價依據,其評估結果可直接作為選擇算子的依據,不需要外部信息。本文中適應度函數確定為:
(13)
式中:E是小波神經網絡的誤差函數。
(3) 遺傳操作 遺傳操作是遺傳算法的核心,其目的是通過逐代優化獲得問題最優解,包括三個基本遺傳算子:選擇、交叉、變異。
選擇算子,在迭代過程中選擇優質的個體,淘汰劣質的個體,本文采用基于適應度比例選擇策略的輪盤賭法。
交叉算子,選擇操作后的新一代群體,按一定的概率在兩兩之間進行交叉操作。交叉概率PC的大小將影響算法的全局搜索能力,一般較大的交叉概率可以促進種群進化。
變異算子,隨機改變種群個體編碼串上的值,即某個個體的個別基因發生0-1反轉。變異操作按一定的概率Pm進行,其目的一方面是使算法具有局部搜索能力,另一方面是保持群體多樣性,防止早熟現象的發生。
利用遺傳算法對網絡權值優化的全局搜索能力與神經網絡良好的逼近能力,可以對復雜的短時交通流量建模,克服了傳統小波神經網絡訓練時易陷入局部極值的弊端。遺傳算法優化小波神經網絡的原理如圖2所示。

圖2 遺傳算法優化小波神經網絡流程圖
算法的基本步驟如下:
Step1設種群規模為M。隨機生成初始種群N=(N1,N2,…,NM),采用實數編碼對個體Ni編碼。
Step2用Step1中的種群N進行訓練,此時的WNN參數由初始化獲得。
Step3計算種群N中每個個體的適應度值。若滿足終止條件則跳至Step6,不滿足則執行Step4。
Step4對于適應度大的個體,選擇其中一部分直接繼承給下一代,同時將遺傳操作后產生的最差個體用父代適應度最大的個體替代。
Step5對當前一代群體進行交叉和變異等遺傳操作,產生下一代群體,轉Step3。
Step6循環結束,將最后一代群體中最優個體的解碼還原值作為WNN的初始參數。
Step7采用Step6獲取的最優參數建立短時交通流量預測模型。
采集長沙市某路口連續5天的地磁車流量數據作為實驗樣本,每15 min統計一次交通流量數據,共得到480個實驗數據,并取其中384個數據作為訓練集,剩余96個數據作為測試集。
為評價GA-WNN算法的性能,仿真實驗以未經優化的WNN算法作為對比模型,以E、MAE、RMSE、MAPE作為評價指標,計算公式如下:
(14)
(15)
(16)
(17)

遺傳算法迭代次數為120次,交叉概率Pc=0.9,變異概率Pm=0.05;GA-WNN及WNN網絡模型參數設置:學習率為0.08,訓練次數為100次,尋優次數為20次。圖3-圖6為本實驗在MATLAB環境下的仿真結果,表1為兩種預測模型各項評價指標。

圖3 遺傳算法迭代誤差進化曲線圖

圖4 預測模型迭代訓練誤差曲線圖

圖5 WNN與GA-WNN預測誤差圖

圖6 WNN與GA-WNN預測結果

評價指標算法模型MAERMSEMAPEWNN5.458 36.265 08.66%GA-WNN4.791 75.588 37.79%
如圖3所示,隨著遺傳迭代的進行,誤差在逐漸縮小,當進行至72次時,網絡趨于穩定。迭代完成后,保存網絡,此時獲得小波神經網絡初始參數值。如圖4所示,WNN網絡在迭代60次時趨于穩定,而GA-WNN網絡在迭代進行至15次時已趨于穩定,表明GA-WNN預測模型收斂速度明顯優于WNN。由圖5可知,優化后的算法可將預測誤差控制在一個相對較小的范圍內波動,并且WNN模型Emax=15,而GA-WNN模型Emax=12。由圖6可知,WNN與GA-WNN雖整體上都能反映短時交通流量變化趨勢,但后者明顯具有更好的擬合效果。如表1所示,后者各項評價指標均小于前者,結合圖6,表明GA-WNN模型對短時交通流量具有更高的預測精度。
綜合分析可知,本文提出的基于遺傳算法優化小波神經網絡的預測模型對于復雜的短時交通流量具有良好的預測精度,能較好地反映交通流量的走勢,有著良好的實用價值。
本文針對短時交通流量的特性,建立了小波神經網絡預測模型,并就小波神經網絡收斂速度慢,易陷入局部極值的缺點,提出了用遺傳算法優化小波神經網絡參數。仿真結果表明,該方法具有學習速率快、逼近能力強和全局尋優等優點,對短時交通流量的預測具有更高的精度。