劉 輝，田紅旗，李燕飛，張 雷

(1.中南大學(xué) 軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室，湖南 長沙 410075；2.德國Rostock大學(xué) 信息科學(xué)與電子工程學(xué)院，德國 羅斯托克 18109)

惡劣強(qiáng)風(fēng)是造成鐵路行車事故的主要自然災(zāi)害之一[1，2]。為應(yīng)對(duì)該災(zāi)害，國內(nèi)外學(xué)者在設(shè)計(jì)鐵路擋風(fēng)墻[3-5]、列車氣動(dòng)外形優(yōu)化[6-9]、強(qiáng)風(fēng)列車傾覆穩(wěn)定性建模[1，10-12]、氣動(dòng)數(shù)值計(jì)算[13-16]等方面開展了大量的研究工作。此外，近年來國際上開始致力于研建一種新型的鐵路強(qiáng)風(fēng)實(shí)時(shí)監(jiān)控與主動(dòng)預(yù)警系統(tǒng)。如德國國家鐵路公司的Nowcasting System大風(fēng)監(jiān)測預(yù)警系統(tǒng)[17]，東日本鐵路公司的Windas System大風(fēng)列車預(yù)警系統(tǒng)[18]以及我國自主開發(fā)的蘭新、青藏等強(qiáng)風(fēng)線路的大風(fēng)監(jiān)測與預(yù)警系統(tǒng)[19]。這類系統(tǒng)均根據(jù)強(qiáng)風(fēng)鐵路沿線測風(fēng)站的實(shí)測風(fēng)速融合列車和線路實(shí)現(xiàn)行車指揮。雖然這種基于實(shí)測風(fēng)速的模式目前也能用于指揮行車，但當(dāng)強(qiáng)風(fēng)消失后，鐵路部門由于無法預(yù)知該位置點(diǎn)的未來風(fēng)速變化趨勢，出于安全考慮，不得不對(duì)受控列車?yán)^續(xù)保持一段時(shí)間的“限速”或“停輪”狀態(tài)，降低了運(yùn)輸效率[17-19]。如果這類大風(fēng)監(jiān)測預(yù)警系統(tǒng)具備風(fēng)速預(yù)測功能，則可按如圖1所示的框架開展工作。即作為鐵路調(diào)度指揮系統(tǒng)TDCS的輔助系統(tǒng)，以列車行車計(jì)劃(包括車次、時(shí)間和到站)為主線，融合鐵路路況信息(如橋梁、隧道、大彎曲半徑等)和風(fēng)速傾覆穩(wěn)定性，計(jì)算獲得不同風(fēng)速等級(jí)的列車安全行車速度。同時(shí)，預(yù)警系統(tǒng)實(shí)時(shí)保持對(duì)沿線實(shí)測風(fēng)速的單步預(yù)測計(jì)算，當(dāng)發(fā)現(xiàn)某列車的前方將出現(xiàn)超等級(jí)的危險(xiǎn)大風(fēng)時(shí)，可以提前向列車發(fā)出“限速”或“停輪”命令。反之，系統(tǒng)通過預(yù)測發(fā)現(xiàn)某大風(fēng)將消失，則可以發(fā)出相應(yīng)的“恢復(fù)”命令。當(dāng)然，這涉及到系統(tǒng)集成、實(shí)時(shí)通信、數(shù)據(jù)融合與決策、行車優(yōu)化、風(fēng)速預(yù)測等眾多問題。本文重點(diǎn)研究鐵路風(fēng)速的高精度瞬時(shí)單步預(yù)測問題。

圖1 融合預(yù)測風(fēng)速的大風(fēng)監(jiān)測預(yù)警系統(tǒng)結(jié)構(gòu)圖

實(shí)際上，鐵路風(fēng)速的高精度預(yù)測已成為國內(nèi)外爭先解決的核心技術(shù)之一。文獻(xiàn)[17]提出基于線性外推理論的鐵路風(fēng)速短期預(yù)測方法。文獻(xiàn)[18]提出基于卡爾曼濾波理論的鐵路風(fēng)速預(yù)測方法。文獻(xiàn)[20-22]開展了鐵路沿線風(fēng)速短期預(yù)測研究。在信號(hào)處理領(lǐng)域，小波分解理論被稱為“數(shù)學(xué)顯微鏡”，能夠?qū)⑷魏我环N非平穩(wěn)信號(hào)轉(zhuǎn)化為一組較為平穩(wěn)信號(hào)[23]。在人工智能領(lǐng)域，神經(jīng)網(wǎng)絡(luò)是信號(hào)非線性處理的優(yōu)秀代表[24]。因此，本文擬采用這兩種優(yōu)秀理論取長補(bǔ)短所形成的兩種不同混合算法(小波-神經(jīng)網(wǎng)絡(luò)法和小波型神經(jīng)網(wǎng)絡(luò)法)預(yù)測鐵路沿線的非平穩(wěn)風(fēng)速信號(hào)，并對(duì)其預(yù)測結(jié)果進(jìn)行對(duì)比分析。

1 鐵路原始風(fēng)速序列

本文運(yùn)用混合預(yù)測算法對(duì)我國青藏鐵路沿線3處監(jiān)控點(diǎn)的實(shí)測大風(fēng)風(fēng)速序列{X1(t)}、{X2(t)}和{X3(t)}進(jìn)行建模與預(yù)測。每處監(jiān)控點(diǎn)擁有400個(gè)采樣數(shù)據(jù)，所有采樣間隔均為3 min。取原始序列的前300個(gè)數(shù)據(jù)建立預(yù)測模型，后100個(gè)數(shù)據(jù)用于檢驗(yàn)?zāi)Ｐ汀?組原始風(fēng)速序列如圖2～圖4所示。

圖2 原始風(fēng)速序列{X1(t)}

圖3 原始風(fēng)速序列{X2(t)}

圖4 原始風(fēng)速序列{X3(t)}

2 小波-神經(jīng)網(wǎng)絡(luò)法預(yù)測鐵路風(fēng)速

2.1 建模原理

小波-神經(jīng)網(wǎng)絡(luò)法的建模流程如圖5所示，包括如下步驟：

(1)選擇小波法對(duì)擬預(yù)測鐵路風(fēng)速序列進(jìn)行多層分解與重構(gòu)計(jì)算，將原始非平穩(wěn)風(fēng)速序列轉(zhuǎn)化為多層較平穩(wěn)風(fēng)速序列，以減低后期分解層上建立高精度神經(jīng)網(wǎng)絡(luò)預(yù)測模型的難度。在本步驟中，為了確定最優(yōu)的小波分解參數(shù)，選擇對(duì)小波母函數(shù)和小波分解深度的不同組合進(jìn)行預(yù)測性能比較。

(2)利用神經(jīng)網(wǎng)絡(luò)法對(duì)各小波分解層分別建立合適的預(yù)測模型。并用所建的神經(jīng)網(wǎng)絡(luò)模型對(duì)各分解層風(fēng)速序列進(jìn)行超前單步預(yù)測計(jì)算。在本步驟中，為了選擇最優(yōu)的神經(jīng)網(wǎng)絡(luò)建模參數(shù)，利用時(shí)序ARIMA模型輔助確定神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)，同時(shí)對(duì)6種主流的神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法進(jìn)行性能篩選。主流學(xué)習(xí)算法包括帶動(dòng)量的梯度下降算法、自適應(yīng)LR動(dòng)量梯度下降算法、彈性梯度下降算法、Fletcher-Reeves共軛梯度算法、量化共軛梯度算法和擬牛頓算法。

(3)對(duì)獲得的各層風(fēng)速預(yù)測值進(jìn)行加權(quán)計(jì)算，輸出原始序列的預(yù)測值。

圖5 小波-神經(jīng)網(wǎng)絡(luò)法的建模和預(yù)測流程圖

2.2 計(jì)算步驟

2.2.1不同的小波分解參數(shù)對(duì)小波-神經(jīng)網(wǎng)絡(luò)法模型的預(yù)測性能影響

小波分解是將原始信號(hào)分解成一些呈現(xiàn)相對(duì)平穩(wěn)性的分解層數(shù)據(jù)的疊加[23]。從數(shù)學(xué)上解釋，小波分解是將基本小波函數(shù)φ(t)平移τ后，在不同尺度a下與待分析的原始信號(hào)x(t)做內(nèi)積，即

( 1 )

本文運(yùn)用小波法分解鐵路沿線非平穩(wěn)風(fēng)速序列的計(jì)算流程，如圖6所示。每進(jìn)行一次小波分解稱為一個(gè)分解深度。

為研究不同小波母函數(shù)和小波分解深度對(duì)小波-神經(jīng)網(wǎng)絡(luò)模型預(yù)測性能的影響，選擇在小波法理論中最具代表性的三類母函數(shù)(Daubechies 函數(shù)、Haar函數(shù)和Symlets函數(shù))和兩種深度工況(3層和6層)的共計(jì)6種不同參數(shù)組合對(duì)圖2所示的鐵路原始風(fēng)速數(shù)據(jù)進(jìn)行小波分解與重構(gòu)計(jì)算。例如，Haar小波與3個(gè)分解深度的組合計(jì)算結(jié)果如圖7所示。對(duì)各個(gè)小波分解結(jié)果數(shù)據(jù)分別建立BP神經(jīng)網(wǎng)絡(luò)模型并進(jìn)行試預(yù)測計(jì)算，再利用式( 2 )～式( 4 )的精度評(píng)價(jià)指標(biāo)對(duì)試預(yù)測結(jié)果進(jìn)行對(duì)比分析，最終篩選出最優(yōu)的鐵路風(fēng)速小波分解參數(shù)組合。試預(yù)測的檢驗(yàn)樣本為原始風(fēng)速序列{X1(t)}的前300個(gè)數(shù)據(jù)。

圖6 三個(gè)分解深度的小波法計(jì)算流程

圖7 小波法分解結(jié)果示例

平均絕對(duì)誤差(MAE)

( 2 )

平均相對(duì)誤差(MPE)

( 3 )

均方根誤差(RMSE)

( 4 )

6種不同小波分解參數(shù)組合的預(yù)測評(píng)價(jià)結(jié)果見表1。為便于區(qū)分，將這6種不同組合分別命名為：H3(Haar小波+3層分解)、H6(Haar小波+6個(gè)分解深度)、D3(Daubechies小波+3個(gè)分解深度)、D6(Daubechies小波+6個(gè)分解深度)、S3(Symlets小波+3個(gè)分解深度)和S6(Symlets小波+6個(gè)分解深度)。

表1 不同小波分解參數(shù)的試預(yù)測結(jié)果分析

通過分析表1可知：

(1)對(duì)同一個(gè)小波母函數(shù)而言，增加分解深度并沒有提高后期神經(jīng)網(wǎng)絡(luò)模型的試預(yù)測精度。以Haar小波函數(shù)為例，其6層分解的平均相對(duì)誤差比3層分解的對(duì)應(yīng)指標(biāo)僅提高6.10%。但由于6層分解比3層分解多產(chǎn)生接近50%的分解層數(shù)據(jù)，因此6層分解在小波分解計(jì)算和神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練方面將多花費(fèi)接近50%的時(shí)間。

(2)對(duì)同一個(gè)分解深度而言，Symlets小波函數(shù)獲得最好的精度指標(biāo)，之后依次是Daubechies 小波函數(shù)和Haar小波函數(shù)。

(3)綜上，本文最終選擇“Symlets函數(shù)+3個(gè)分解深度”用于鐵路風(fēng)速預(yù)測。

2.2.2不同的神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法對(duì)小波-神經(jīng)網(wǎng)絡(luò)混合模型的預(yù)測性能影響

如圖5所示，小波-神經(jīng)網(wǎng)絡(luò)法除了需要確定合適的小波分解參數(shù)之外，還需要對(duì)每個(gè)分解層建立最優(yōu)的神經(jīng)網(wǎng)絡(luò)預(yù)測模型。在對(duì)分解層建立神經(jīng)網(wǎng)絡(luò)模型方面，參考文獻(xiàn)[2]提出了基于時(shí)間序列模型輔助確定神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)參數(shù)的做法，即先對(duì)各個(gè)分解層的數(shù)據(jù)建立時(shí)序ARIMA模型，然后將獲得的ARIMA差分方程的自變量認(rèn)定為神經(jīng)網(wǎng)絡(luò)的輸入節(jié)點(diǎn)，將其因變量認(rèn)定為神經(jīng)網(wǎng)絡(luò)的輸出節(jié)點(diǎn)，最后運(yùn)用嘗試法迭代篩選出最優(yōu)的神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點(diǎn)數(shù)。這種利用時(shí)間序列模型優(yōu)化神經(jīng)網(wǎng)絡(luò)模型的做法能夠有效避免網(wǎng)絡(luò)參數(shù)選定的隨機(jī)性和人為性，提高所建網(wǎng)絡(luò)的非線性擬合精度。

除了神經(jīng)元數(shù)量，神經(jīng)網(wǎng)絡(luò)模型對(duì)學(xué)習(xí)算法也很敏感。但在文獻(xiàn)[2]中沒有這部分內(nèi)容。本文將補(bǔ)充開展這方面的工作。本文用圖2所示原始風(fēng)速序列{X1(t)}的前400個(gè)采樣點(diǎn)作為實(shí)驗(yàn)樣本，對(duì)6種主流的神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法進(jìn)行分類考核，其結(jié)果如圖8和表2所示。

圖8 不同神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法的擬合結(jié)果

學(xué)習(xí)算法指標(biāo)σ1/(m·s-1)σ2/%σ3/(m·s-1)帶動(dòng)量的梯度下降算法1.748813.932.2170自適應(yīng)LR動(dòng)量梯度下降算法1.06738.811.3416彈性梯度下降算法0.87448.281.1389Fletcher-Reeves共軛梯度算法0.91988.291.0481量化共軛梯度算法0.84096.471.0447擬牛頓算法0.34682.710.4525

分析圖8和表2可知：在相同的鐵路風(fēng)速學(xué)習(xí)樣本和條件下，擬牛頓算法獲得了最優(yōu)的學(xué)習(xí)效果，其擬合的平均相對(duì)誤差僅為2.71%，遠(yuǎn)高于其他的網(wǎng)絡(luò)學(xué)習(xí)算法。因此，本文最終選擇擬牛頓算法用于風(fēng)速預(yù)測計(jì)算。

3 小波型神經(jīng)網(wǎng)絡(luò)法預(yù)測鐵路風(fēng)速

小波型神經(jīng)網(wǎng)絡(luò)是以前饋型神經(jīng)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)為基礎(chǔ)，將小波基函數(shù)作為隱含層節(jié)點(diǎn)的傳遞函數(shù)，遵循誤差方向傳播的神經(jīng)網(wǎng)絡(luò)模型[24]。如圖9所示，小波基取代其他傳統(tǒng)函數(shù)(如Sigmoid函數(shù))作為隱含層節(jié)點(diǎn)的傳遞函數(shù)，圖中，Ni是輸入神經(jīng)元個(gè)數(shù)，Nj是隱含層小波基神經(jīng)元個(gè)數(shù)，Nk是輸出神經(jīng)元個(gè)數(shù)。

圖9 小波型神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

根據(jù)如圖9所示的小波型神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)，當(dāng)輸入風(fēng)速序列為{xi}時(shí)，網(wǎng)絡(luò)的隱含層輸出為

( 5 )

式中：h(j)為隱含層第j個(gè)節(jié)點(diǎn)輸出值；ωij為輸入層和隱含層的連接權(quán)值；hj為小波基母函數(shù)；τj為小波基函數(shù)hj的平移因子；aj為小波基函數(shù)hj的伸縮因子；Ni為輸入層節(jié)點(diǎn)數(shù)。

小波神經(jīng)網(wǎng)絡(luò)輸出層的計(jì)算公式為

( 6 )

式中：y(k)為輸出層的輸出值；ωjk為隱含層和輸出層的連接權(quán)值；h(j)為第j個(gè)隱含層節(jié)點(diǎn)的輸出；Nj為隱含層節(jié)點(diǎn)數(shù)；Nk為輸出層節(jié)點(diǎn)數(shù)。

本文采取梯度下降算法作為小波型神經(jīng)網(wǎng)絡(luò)模型的嵌入學(xué)習(xí)算法，具體計(jì)算公式如下。

(1)計(jì)算網(wǎng)絡(luò)預(yù)測誤差

( 7 )

(2)修正小波型神經(jīng)網(wǎng)絡(luò)權(quán)值和小波基函數(shù)系數(shù)

( 8 )

( 9 )

小波型神經(jīng)網(wǎng)絡(luò)模型也根據(jù)時(shí)序ARIMA模型確定網(wǎng)絡(luò)的輸出-輸出節(jié)點(diǎn)結(jié)構(gòu)。小波型神經(jīng)網(wǎng)絡(luò)模型將按照如圖10所示的計(jì)算框架完成鐵路風(fēng)速超前預(yù)測計(jì)算。

圖10 小波型神經(jīng)網(wǎng)絡(luò)模型的預(yù)測流程

4 預(yù)測結(jié)果分析與比較

為充分比較小波型神經(jīng)網(wǎng)絡(luò)法和小波-神經(jīng)網(wǎng)絡(luò)法兩種不同算法的預(yù)測性能，本文分別用小波型神經(jīng)網(wǎng)絡(luò)模型、小波-神經(jīng)網(wǎng)絡(luò)模型和單種BP神經(jīng)網(wǎng)絡(luò)模型對(duì)如圖2～圖4所示的3組鐵路原始風(fēng)速數(shù)據(jù)進(jìn)行超前多步預(yù)測計(jì)算。在預(yù)測過程中，將建模數(shù)據(jù)和預(yù)測數(shù)據(jù)分開，即用風(fēng)速序列的前300個(gè)數(shù)據(jù)建立預(yù)測模型，后100個(gè)數(shù)據(jù)檢驗(yàn)?zāi)Ｐ汀?/p>

(1)小波型神經(jīng)網(wǎng)絡(luò)模型對(duì)不同測風(fēng)站風(fēng)速的預(yù)測結(jié)果

小波型神經(jīng)網(wǎng)絡(luò)混合模型對(duì)3組鐵路原始風(fēng)速序列的預(yù)測結(jié)果如圖11～圖13所示。利用式( 2 )～式( 4 )預(yù)測精度評(píng)價(jià)指標(biāo)評(píng)估所獲得的預(yù)測值結(jié)果見表3。

圖11 {X1(t)}序列的小波型神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖12 {X2(t)}序列的小波型神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖13 {X3(t)}序列的小波型神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

算例預(yù)測精度指標(biāo)σ1/(m·s-1)σ2/%σ3/(m·s-1)圖11算例1.703710.812.0733圖12算例1.48128.001.7849圖13算例1.65849.441.9256

(2)小波-神經(jīng)網(wǎng)絡(luò)模型對(duì)不同測風(fēng)站風(fēng)速的預(yù)測結(jié)果

小波-神經(jīng)網(wǎng)絡(luò)混合模型對(duì)3組鐵路原始風(fēng)速序列的預(yù)測結(jié)果如圖14～圖16所示。其對(duì)應(yīng)的預(yù)測評(píng)價(jià)指標(biāo)見表4。

圖14 {X1(t)}序列的小波-神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖15 {X2(t)}序列的小波-神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖16 {X3(t)}序列的小波-神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

算例預(yù)測精度指標(biāo)σ1/(m·s-1)σ2/%σ3/(m·s-1)圖14算例0.45973.310.6726圖15算例0.22201.560.2749圖16算例0.41272.450.5649

(3)單種神經(jīng)網(wǎng)絡(luò)模型對(duì)不同測風(fēng)站風(fēng)速的預(yù)測結(jié)果

單種神經(jīng)網(wǎng)絡(luò)模型對(duì)3組鐵路原始風(fēng)速序列的預(yù)測結(jié)果如圖17～圖19所示。其對(duì)應(yīng)的預(yù)測評(píng)價(jià)指標(biāo)見表5。

圖17 {X1(t)}序列的神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖18 {X2(t)}序列的神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

圖19 {X3(t)}序列的神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果

算例預(yù)測精度指標(biāo)σ1/(m·s-1)σ2/%σ3/(m·s-1)圖17算例1.778512.532.1000圖18算例1.54349.951.8363圖19算例1.669011.131.9668

(4)三種不同算法的性能對(duì)比

為了揭示小波分解在小波型神經(jīng)網(wǎng)絡(luò)法和小波-神經(jīng)網(wǎng)絡(luò)法兩者不同混合算法中的貢獻(xiàn)度，將兩種混合算法對(duì)單種神經(jīng)網(wǎng)絡(luò)模型的性能優(yōu)化百分比進(jìn)行統(tǒng)計(jì)，見表6、表7。

通過分析表3～表5可知：小波-神經(jīng)網(wǎng)絡(luò)法的預(yù)測性能在三種預(yù)測算法中最優(yōu)，然后依次是小波型神經(jīng)網(wǎng)絡(luò)法和單種神經(jīng)網(wǎng)絡(luò)法。以評(píng)價(jià)指標(biāo)體系中最為重要的平均相對(duì)誤差為例，在本文的3組預(yù)測實(shí)例中，小波-神經(jīng)網(wǎng)絡(luò)法的誤差分別為3.31%、1.56%和2.45%，屬于高精度5%的預(yù)測范疇。而小波型神經(jīng)網(wǎng)絡(luò)法的同類指標(biāo)分別為10.81%、8.00%和9.44%；單種神經(jīng)網(wǎng)絡(luò)法的同類指標(biāo)分別為12.53%、9.95%和11.13%。特別需要指出的是小波-神經(jīng)網(wǎng)絡(luò)法實(shí)現(xiàn)了對(duì)鐵路瞬間突變風(fēng)速數(shù)據(jù)的有效追蹤。

表6 小波型神經(jīng)網(wǎng)絡(luò)法提高神經(jīng)網(wǎng)絡(luò)預(yù)測精度 %

表7 小波-神經(jīng)網(wǎng)絡(luò)法提高神經(jīng)網(wǎng)絡(luò)預(yù)測精度 %

通過分析表6和表7可知，在3組不同的預(yù)測實(shí)例中，小波型神經(jīng)網(wǎng)絡(luò)混合模型對(duì)單種神經(jīng)網(wǎng)絡(luò)模型的平均相對(duì)誤差優(yōu)化幅度分別為13.73%、19.60%和15.18%，而小波-神經(jīng)網(wǎng)絡(luò)混合模型對(duì)單種神經(jīng)網(wǎng)絡(luò)的平均相對(duì)誤差優(yōu)化幅度分別為73.58%、84.32%和77.99%。后者的優(yōu)化幅度遠(yuǎn)優(yōu)于前者。這說明小波法用于前端分解非平穩(wěn)風(fēng)速所獲得的貢獻(xiàn)度高于小波法用于后期優(yōu)化神經(jīng)網(wǎng)絡(luò)傳遞函數(shù)所獲得的貢獻(xiàn)度。究其原因，小波-神經(jīng)網(wǎng)絡(luò)法由于采用了小波法分解原始風(fēng)速序列，不僅直接降低了后期神經(jīng)網(wǎng)絡(luò)模型的高精度擬合難度，而且提供了用于建立更多神經(jīng)網(wǎng)絡(luò)所需要的分解層數(shù)據(jù)。這使得小波-神經(jīng)網(wǎng)絡(luò)混合模型實(shí)現(xiàn)了從單個(gè)神經(jīng)網(wǎng)絡(luò)到多個(gè)神經(jīng)網(wǎng)絡(luò)群的進(jìn)化，因此可以獲得較為理想的預(yù)測結(jié)果。而小波型神經(jīng)網(wǎng)絡(luò)法雖然利用小波函數(shù)作為網(wǎng)絡(luò)隱含層傳遞函數(shù)可以優(yōu)化傳統(tǒng)神經(jīng)網(wǎng)絡(luò)的非線性搜索性能，但究其根本它始終還是只有一個(gè)神經(jīng)網(wǎng)絡(luò)模型參與實(shí)際預(yù)測計(jì)算，因此小波傳遞函數(shù)所帶來的優(yōu)化幅度無法與小波前端分解所帶來的效益相提并論。

5 結(jié)論

(1)文章的三組預(yù)測實(shí)例表明，小波-神經(jīng)網(wǎng)絡(luò)法和小波型神經(jīng)網(wǎng)絡(luò)法這兩種混合算法在不同程度上都吸收了小波法的信號(hào)細(xì)分功能，具有對(duì)鐵路風(fēng)速跳躍點(diǎn)出色的辨識(shí)和分解能力。同時(shí)，兩者也繼承了神經(jīng)網(wǎng)絡(luò)的非線性隱射能力，能夠應(yīng)對(duì)鐵路風(fēng)速非線性變化趨勢的預(yù)測要求。由于兩者融合了小波分解的算法特性，因此他們都獲得了比單種神經(jīng)網(wǎng)絡(luò)更加優(yōu)秀的預(yù)測算法性能。但兩者相比較，前者通過運(yùn)用小波法分解原始風(fēng)速數(shù)據(jù)從而獲得了對(duì)多分解層建立多神經(jīng)網(wǎng)絡(luò)預(yù)測“群”的機(jī)會(huì)，因此達(dá)到5%的高精度預(yù)測級(jí)別。后者通過將小波函數(shù)作為網(wǎng)絡(luò)傳遞函數(shù)的手段從而提高了傳統(tǒng)單個(gè)神經(jīng)網(wǎng)絡(luò)的全局搜索性能，因此達(dá)到10%的中等精度預(yù)測級(jí)別。前者比后者獲得了更好的優(yōu)化性能。

(2)相對(duì)于其他算法，如時(shí)間序列分析法、單種自適應(yīng)神經(jīng)模糊推理法、卡爾曼濾波法等而言，雖然小波-神經(jīng)網(wǎng)絡(luò)法和小波型神經(jīng)網(wǎng)絡(luò)法這兩種混合算法所包含的建模步驟更多，但通過借助MATLAB神經(jīng)網(wǎng)絡(luò)和小波分析工具箱可以方便地實(shí)現(xiàn)其計(jì)算步驟的編程和集成。因此他們可以在相關(guān)鐵路風(fēng)速預(yù)警系統(tǒng)中推廣使用。

參考文獻(xiàn)：

[1]高廣軍,田紅旗,姚松,等.蘭新線強(qiáng)橫風(fēng)對(duì)車輛傾覆穩(wěn)定性的影響[J].鐵道學(xué)報(bào),2004,26(4):36-40.

GAO Guangjun,TIAN Hongqi,YAO Song,et al.Effect of Strong Cross-wind on the Stability of Trains Running on the Lanzhou-Xinjiang Railway Line[J].Journal of the China Railway Society,2004,26(4):36-40.

[2]劉輝,田紅旗,李燕飛.基于小波分析法與神經(jīng)網(wǎng)絡(luò)法的非平穩(wěn)風(fēng)速信號(hào)短期預(yù)測優(yōu)化算法[J].中南大學(xué)學(xué)報(bào):自然科學(xué)版,2011,42(9):2 704-2 711.

LIU Hui,TIAN Hongqi,LI Yanfei.Short-term Forecasting Optimization Algorithm for Wind Speed from Wind Farms Based on Wavelet Analysis Method and Rolling Time Series Method[J].Journal of Central South University:Science and Technology,2011,42(9)：2 704-2 711.

[3]李田,張繼業(yè),張衛(wèi)華.橫風(fēng)下高速列車通過擋風(fēng)墻動(dòng)力學(xué)性能[J].鐵道學(xué)報(bào),2012,34(7):30-35.

LI Tian,ZHANG Jiye,ZHANG Weihua.Dynamic Performance of High-speed Train Passing Windbreak in Crosswind[J].Journal of the China Railway Society,2012,34(7):30-35.

[4]李鯤.大風(fēng)區(qū)高速鐵路新型防風(fēng)設(shè)施研究[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版,2012,43(2):756-762.

LI Kun.Research on New Anti-wind Facility of High-speed Train in Strong Wind Area[J].Journal of Central South University:Science and Technology,2012,43(2):756-762.

[5]CHU C R,CHANG C Y,HUANG C J,et al.Windbreak Protection for Road Vehicles Against Crosswind[J].Journal of Wind Engineering and Industrial Aerodynamics,2013,116(5):61-69.

[6]CHELI F,GIAPPINO S,ROSA L,et al.Experimental Study on the Aerodynamic Forces on Railway Vehicles in Presence of Turbulence[J].Journal of Wind Engineering and Industrial Aerodynamics,2013,123(6):311-316.

[7]SIMONOVIC A,SVORCAN J,STUPAR S.Aerodynamic Characteristics of High Speed Train under Turbulent Cross Winds:A Numerical Investigation Using Unsteady-RANS Method[J].FME Transactions,2014,42(1):10-18.

[8]張雷,楊明智,張輝, 等.高速鐵路隧道洞門對(duì)隧道空氣動(dòng)力效應(yīng)的影響[J].鐵道學(xué)報(bào),2013,35(11):76-79.

ZHANG Lei,YANG Mingzhi,ZHANG Hui,et al.Influence of Tunnel Portals on Tunnel Aerodynamic Effect in Operation of High-speed Railways[J].Journal of the China Railway Society,2013,35(11):76-79.

[9]WAGNER S,EHRENFRIED K,DILLMANN A.Numerical Simulation of Train-Tunnel Entry Using a BEM in Time Domain[C]//New Results in Numerical and Experimental Fluid Mechanics VIII.Berlin:Springer Berlin Heidelberg,2013:739-746.

[10]周丹,田紅旗,楊明智,等.強(qiáng)側(cè)風(fēng)下客車在不同路況運(yùn)行的氣動(dòng)性能比較[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版,2008,39(3):554-559.

ZHOU Dan,TIAN Hongqi,YANG Mingzhi,et al.Comparison of Aerodynamic Performance of Passenger Train Traveling on Different Railway Conditions up Strong Cross-wind[J].Journal of Central South University:Science and Technology,2008,39(3):554-559.

[11]熊小慧,梁習(xí)鋒.橫風(fēng)作用下貨車篷布結(jié)構(gòu)強(qiáng)度計(jì)算[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版,2012,43(8):3 280-3 286.

XIONG Xiaohui,LIANG Xifeng.Numerical Simulation of Tarpaulin Structural Strength of Railway Vehicle under Cross Wind Condition[J].Journal of Central South University:Science and Technology,2012,43(8):3 280-3 286.

[12]ALLAIN E,PARADOT N.Aerodynamics in Train Cross Wind Studies[J].International Journal of Aerodynamics,2014,4(1):10-23.

[13]LI T,ZHANG J Y,ZHANG W H.Co-simulation of High-speed Train Fluid-structure Interaction Dynamics in Crosswinds[J].Journal of Vibration Engineering,2012,25(2):138-145.

[14]GILBERT T,BAKER C J,QUINN A.Gusts Caused by High-speed Trains in Confined Spaces and Tunnels[J].Journal of Wind Engineering and Industrial Aerodynamics,2013,121(5):39-48.

[15]BAKER C.A Framework for the Consideration of the Effects of Crosswinds on Trains[J].Journal of Wind Engineering and Industrial Aerodynamics,2013,123:130-142.

[16]BAKER C,CHELI F,ORELLANO A,et al.Cross-wind Effects on Road and Rail Vehicles[J].Vehicle System Dynamics,2009,47(8):983-1 022.

[17]HOPPMANN U,KOENIG S,TIELKES T,et al.A Short-term Strong Wind Prediction Model for Railway Application:Design and Verification[J].Journal of Wind Engineering and Industrial Aerodynamics,2002,90(10):1 127-1 134.

[18]KOBAYASHI N,SHIMAMURA M.Study of a Strong Wind Warning System[J].JR East Technical Review,2003(2)：61-65.

[19]賈國裕.蘭新鐵路大風(fēng)災(zāi)害及其對(duì)策[J].路基工程,2008(2):195-197.

JIA Guoyu.Wind Hazards and Countermeasures in Lanzhou-Xinjiang Railway[J].Subgrade Engineering,2008(2):195-197.

[20]潘迪夫,劉輝,李燕飛,等.青藏鐵路格拉段沿線風(fēng)速短時(shí)預(yù)測方法[J].中國鐵道科學(xué),2009,30(5):129-133.

PAN Difu,LIU Hui,LI Yanfei,et al.A Short-term Forecast Method for Wind Speed Along Golmud-Lhasa Section of Qinghai-Tibet Railway[J].China Railway Science,2009,30(5):129-133.

[21]王瑞,史天運(yùn),王彤.基于RBF神經(jīng)網(wǎng)絡(luò)的鐵路沿線短時(shí)風(fēng)速預(yù)測方法[J].中國鐵道科學(xué),2011,32(5)：132-134.

WANG Rui,SHI Tianyun,WANG Tong.Prediction Method for Short-time Wind Speed Along Railway Based on RBF Neural Network[J].China Railway Science,2011,32(5)：132-134.

[22]苗秀娟.瞬態(tài)風(fēng)荷載下的列車運(yùn)行安全性研究[D].長沙:中南大學(xué),2011:20-35.

[23]楊淑平,易國棟,袁修貴,等.一種基于分塊小波的人臉識(shí)別算法[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版,2013,44(5)：1 902-1 909.

YANG Shuping,YI Guodong,YUAN Xiugui,et al.A Face Recognition Algorithm Based on Blocking Wavelet Transforms[J].Journal of Central South University:Science and Technology,2013,44(5)：1 902-1 909.

[24]史峰,王小川,郁磊,等.MATLAB神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析[M].北京：北京航空航天大學(xué)出版社,2010：208-210.