基于小波分析AR(P)-SVR組合模型在大壩變形預(yù)測中的應(yīng)用

沈哲輝，黃騰，邱偉，鄭浩

(河海大學(xué)地球科學(xué)與工程學(xué)院，江蘇南京210098)

摘要：監(jiān)測序列經(jīng)小波分解后，得到低頻分量和高頻分量。對低頻分量采用自回歸AR(P)模型預(yù)測，對高頻分量采用支持向量回歸機(jī)SVR模型預(yù)測，最后將各分量進(jìn)行小波重構(gòu)，得到監(jiān)測序列的預(yù)測值。結(jié)果表明，此種預(yù)測方法比直接使用SVR模型或經(jīng)小波分解后再采用SVR模型預(yù)測精度高。

關(guān)鍵詞：小波分解；AR(P)模型；SVR模型；小波重構(gòu)；預(yù)測

中圖分類號：TU196

收稿日期：2014-04-30;修回日期：2014-11-23

作者簡介：沈哲輝(1991-)，男，碩士研究生.

ApplicationofAR(P)-SVRcombinationmodelbasedonwaveletanalysisindamdeformationpredictionSHENZhe-hui,HUANGTeng,QIUWei,ZHENGHao

(SchoolofEarthSciencesandEngineering,HehaiUniversity,Nanjing210098,China)

Abstract：Low frequency and high frequency components are obtained through wavelet decomposition.The low frequency components are adopted in the AR(P) model to make predictions, while the high frequency components make predictions with SVR model. Then the predicted data after reconstructing them are obtained.Results show this model has higher prediction accuracy than SVR model without wavelet decomposition and the model that only uses SVR to predict each component after wavelet decomposition.

Keywords：waveletdecomposition;AR(P)model;SVRmodel;waveletreconstruction;prediction

大壩變形預(yù)測的方法有灰色模型法、自回歸模型、神經(jīng)網(wǎng)絡(luò)等。但由于大壩變形受影響因素多，組合模型的應(yīng)用使大壩變形預(yù)測精度有了提高[1-4]。

小波變換是一個時間和頻率的局域變換，能有效地提取信號，通過伸縮和平移對函數(shù)或信號進(jìn)行多尺度分析[5]，并且小波分解后的平穩(wěn)性比原信號好的多。支持向量機(jī)是一種新興機(jī)器學(xué)習(xí)方法，克服了神經(jīng)網(wǎng)絡(luò)的過擬合、收斂速度慢、容易陷入局部極值等缺點，被公認(rèn)為是人工神經(jīng)網(wǎng)絡(luò)的替代方法[6]。小波分解后的低頻分量，用線性自回歸模型預(yù)測，而高頻分量仍用SVR預(yù)測模型，重構(gòu)后相比低、高頻分量使用SVR預(yù)測精度頗為理想。

1建立組合模型

1.1小波分解

利用Mallat算法，首先將離散的大壩監(jiān)測數(shù)據(jù)分成高頻部分d1和低頻部分c1，然后再把低頻部分c1進(jìn)一步分解，重復(fù)以上過程可以得到任意分辨率下的高頻和低頻部分。具體算法[7]：

(1)

(2)

式中:j=1,2,…，J-1為分解尺度;k，n為平移系數(shù)；Cj+1,k，Dj+1,k分別為在2j分辨率下的低頻近似分量系數(shù)和高頻細(xì)節(jié)分量系數(shù)；gn-2k，hn-2k分別為高通濾波系數(shù)和低通濾波系數(shù)。

式中:ψ(x)為小波函數(shù)，φ(x)為尺度函數(shù)。

利用分解后的小波系數(shù)可重構(gòu)原來的序列小波重構(gòu)的Mallat算法為

(3)

監(jiān)測序列經(jīng)小波分解后，其低頻分量用自回歸AR(P)模型進(jìn)行見面預(yù)測，高頻分量用支持向量回歸機(jī)SVR模型進(jìn)行建模預(yù)測。

1.2自回歸AR(P)模型

AR(P)模型稱為自回歸模型，由于該模型在建模時只需解線性方程組，而又不涉及白噪聲序列值，計算簡便，所以在時序分析中，AR(P)是采用最多的一種模型[8]。AR(P)模型表示為

(4)

本文采用FPE定階準(zhǔn)則，即最小最終預(yù)報誤差準(zhǔn)則對p定階，回歸系數(shù)采用最小二乘法求得。

1.3支持向量機(jī)模型

在20世紀(jì)90年代，Vapnik提出支持向量機(jī)是數(shù)據(jù)挖掘中的一項新技術(shù)，是借助于最優(yōu)化的方法解決機(jī)器學(xué)習(xí)問題的新工具，它是一種以結(jié)構(gòu)風(fēng)險最小化原則為基礎(chǔ)的機(jī)器學(xué)習(xí)算法，通過尋找一個線性約束的二次規(guī)劃問題得到全局最優(yōu)解，因而不存在局部極小值的問題，快速算法保證了收斂速度[9]。

支持向量回歸機(jī)(SupportVectorRegression，SVR)是支持向量在函數(shù)回歸領(lǐng)域的應(yīng)用。支持向量回歸機(jī)首先考慮用線性回歸函數(shù)

(5)

擬合(xi,yi),i=1,2,…,n；xi∈Rn為輸入量，yi∈R為輸出量，即需要確定ω和b。在非線性的情況下，通過一個非線性映射將數(shù)據(jù)映射到高維特征空間中，并在這個空間進(jìn)行線性回歸。需要引入核函數(shù)K(x,xi)代替內(nèi)積運算。常用的核函數(shù)為

1)多項式核：K(x,xi)=(x,xi+1)d，d為階數(shù)；

3)感知器核：K(x,xi)=tanh(βxi+b)。

本文采用的核函數(shù)為徑向基RBF核函數(shù)(高斯核函數(shù))。RBF核函數(shù)是一個普通的核函數(shù)，通過參數(shù)的選擇，適用于任意分布的樣本，因此RBF核

函數(shù)是目前支持向量機(jī)中被應(yīng)用最廣泛的一種核函數(shù)[10]。選擇合適的損失參數(shù)ε和懲罰參數(shù)C，構(gòu)造高維特征空間中求回歸問題：

(6)

(7)

得到非線性回歸函數(shù)為

(8)

本文采用網(wǎng)格法尋求最優(yōu)的損失參數(shù)ε和懲罰參數(shù)C。

1.4小波分解下的AR(P)-SVR組合模型

本文構(gòu)建小波分解的AR(P)-SVR組合模型。

1)將監(jiān)測序列進(jìn)行小波分解，小波函數(shù)選用Daubechies系列小波函數(shù)db3。將監(jiān)測序列分解為2層，并進(jìn)行單支重構(gòu)，得到低頻分量a2和高頻分量d1，d2。

2)用自回歸AR(P)模型對低頻分量進(jìn)行預(yù)測，得到擬合值和預(yù)測值。

3)用支持向量回歸機(jī)SVR對各層高頻分量進(jìn)行建模并預(yù)測。對每個高頻分量，分別構(gòu)建一個多輸入、單輸出的支持向量回歸機(jī)預(yù)測模型[11]。

4)小波重構(gòu)，將用AR(P)模型預(yù)測后的低頻分量和用SVR建模預(yù)測后的高頻分量進(jìn)行小波重構(gòu)，得到監(jiān)測序列的擬合值和預(yù)測值，建模流程如圖1所示。

圖1　建模流程

2實例分析

本文采用福建省某大壩的50期順河向水平位移觀測數(shù)據(jù)，利用前45期數(shù)據(jù)建模，并預(yù)測后5期，用后5期數(shù)據(jù)檢驗預(yù)測的準(zhǔn)確性。

本文采用緊支撐正交小波db3小波進(jìn)兩層小波分解。在對低頻分量a2進(jìn)行AR(P)建模時，采用FPE定階準(zhǔn)則，最后確定階數(shù)為4階。對高頻分量d1和d2進(jìn)行SVR建模時，根據(jù)經(jīng)驗，取輸入量維數(shù)為5，采用網(wǎng)格法尋找最優(yōu)參數(shù)。圖2為a2的自回歸模型擬合及預(yù)測曲線，圖3為d2支持向量回歸機(jī)擬合及預(yù)測曲線，圖4為d1支持向量回歸機(jī)擬合及預(yù)測曲線。

圖2　a 2的自回歸模型擬合及預(yù)測曲線

圖3　d 2的支持向量回歸機(jī)模型擬合及預(yù)測曲線

圖4　d 1的支持向量回歸機(jī)模型擬合預(yù)測曲線

得到各個分量的擬合值和預(yù)測值后，再進(jìn)行小波重構(gòu)，根據(jù)S=a2+d2+d1進(jìn)行重構(gòu)，得到最終的擬合值及預(yù)測值，并將5期的預(yù)測值與后5期的實測值進(jìn)行對比，用相對誤差衡量預(yù)測精度。為驗證本文模型(模型1)的預(yù)測效果，對比使用SVR模型預(yù)測(模型2)的精度和經(jīng)小波分解的各分量都采用SVR模型預(yù)測(模型3)的精度。各模型預(yù)測結(jié)果比較見表1。

表1　預(yù)測結(jié)果比較

由表1計算可得，模型2(SVR)的預(yù)測值相對誤差的平均值為1.45%，模型3(W-SVR)的預(yù)測值相對誤差的平均值為1.44%，可見模型3預(yù)測精度比模型2的預(yù)測精度有提高，但不明顯。本文建立的組合模型，即模型1的預(yù)測值相對誤差的平均值為1.24%，可見本文的組合模型預(yù)測效果明顯比其他兩種模型的預(yù)測效果好。

3結(jié)束語

本文將小波分析、自回歸AR(P)模型和支持向量回歸機(jī)SVR結(jié)合在一起，構(gòu)建用于預(yù)測大壩變形的組合模型。通過小波分解，將低頻分量和高頻分量分別用自回歸模型和支持向量回歸機(jī)模型建模預(yù)測、重構(gòu)。經(jīng)實驗發(fā)現(xiàn)，監(jiān)測序列經(jīng)小波分解后，相比SVR,其低頻分量更適合用自回歸模型建模預(yù)測。小波重構(gòu)后，本文構(gòu)建的組合模型的預(yù)測精度相對于經(jīng)小波分解的各分量都采用SVR模型預(yù)測的精度高。

參考文獻(xiàn)：

[1]王新洲,范千,許承權(quán),等.基于小波變換和支持向量機(jī)的大壩變形預(yù)測[J].武漢大學(xué)學(xué)報：信息科學(xué)版,2008,33(5)：469-471.

[2]郭鍇，陳偉清，文鴻雁，等.非等時距灰色線性組合模型在變形監(jiān)測中的應(yīng)用[J].測繪科學(xué)，2014,39(7)：131-135.

[3]簡伶俐，羅志清，趙福洪.基于二次平滑和指數(shù)曲線的串聯(lián)式組合模型在沉降預(yù)測中的應(yīng)用[J].測繪工程，2014,23(9)：45-47.新洲,范千,許承權(quán),等.基于小波變換和支持向量機(jī)的大壩變形預(yù)測[J].武漢大學(xué)學(xué)報：信息科學(xué)版,2008,33(5)：469-471.

[2]郭鍇，陳偉清，文鴻雁，等.非等時距灰色線性組合模型在變形監(jiān)測中的應(yīng)用[J].測繪科學(xué)，2014,39(7)：131-135.

[3]簡伶俐，羅志清，趙福洪.基于二次平滑和指數(shù)曲線的串聯(lián)式組合模型在沉降預(yù)測中的應(yīng)用[J].測繪工程，2014,23(9)：45-47.

[4]唐佑輝，黃騰，趙兵帥.GOM和時間序列組合模型在地表沉降預(yù)測中的應(yīng)用[J].測繪工程，2014,23(2)：70-72.

[5]胡曉虎.基于小波-BP神經(jīng)網(wǎng)絡(luò)的風(fēng)電場短期風(fēng)速預(yù)測[J].銅陵學(xué)院學(xué)報,2012,11(4)：107-109.

[6]崔艷,程躍華.小波支持向量機(jī)在交通流量預(yù)測中的應(yīng)用[J].計算機(jī)仿真,2011,28(7)：353-356.

[7]劉新東,陳煥遠(yuǎn),佘彩綺.基于小波閥值降噪和BP神經(jīng)網(wǎng)絡(luò)的超短期風(fēng)電功率預(yù)測[J].世界科技研究與發(fā)展,2011,33(6)：1006r1010.

[8]郝少峰，方源敏，楊建文，等.基于熵權(quán)法的組合模型在滑坡變形預(yù)測中的應(yīng)用[J].測繪工程，2014,23(7)：62-64.

[9]唐亮貴,程代杰.基于小波的支持向量機(jī)預(yù)測模型及應(yīng)用[J].計算機(jī)科學(xué),2006,33(3)：202-204.

[10]唐爭氣，譚志強(qiáng).灰色時序組合模型在基坑監(jiān)測中的運用[J].測繪工程，2014,23(2)：49-53.

[11]黃惠峰，張獻(xiàn)州.高速鐵路沉降變形的組合預(yù)測方法應(yīng)用[J].測繪工程，2014,23(9)：48-51.

[責(zé)任編輯：李銘娜]