徐健,邵山
(西安工程大學(xué) 電信學(xué)院,陜西 西安 710048)
小波變換在電壓閃變信號(hào)中的檢測(cè)分析
徐健,邵山
(西安工程大學(xué) 電信學(xué)院,陜西 西安 710048)
為了實(shí)現(xiàn)對(duì)電力系統(tǒng)中電壓閃變信號(hào)更為精確地檢測(cè),提出了小波變換的多分辨率分析的原理。由于小波變換具有較好的時(shí)頻局部化特性,進(jìn)而提取出電壓閃變信號(hào)中的頻率和幅值等參數(shù)。最后通過MATLAB仿真證明小波變換是一種更為有效的電壓閃變分析方法。
電力系統(tǒng);電壓閃變;小波變換;多分辨率分析;MATLAB
電壓閃變作為電能質(zhì)量的重要指標(biāo)之一,通常是由大功率沖擊性負(fù)載的投入所引起的。電壓閃變不僅會(huì)引起人的視覺疲勞也會(huì)影響到電子儀器、自動(dòng)控制設(shè)備等正常工作[1]。如果能對(duì)電力系統(tǒng)中的電壓閃變的信號(hào)進(jìn)行有效的檢測(cè),將會(huì)給電壓閃變的治理帶來很大的幫助,對(duì)整個(gè)電力系統(tǒng)的安全穩(wěn)定的運(yùn)行也有十分重要的意義。
隨著電力系統(tǒng)的不斷發(fā)展,也逐漸形成了多種對(duì)電壓閃變的檢測(cè)方法。如傅里葉變換以及 Kalam濾波、最小誤差二乘法、小波變換、S變換等。文中主要介紹前兩種方法及其應(yīng)用,最后通過對(duì)比,小波變換具有良好的時(shí)頻局部化特性[2],可以更加方便運(yùn)用到電壓閃變信號(hào)檢測(cè)中去。
公式(1)中:V(t)是電網(wǎng)電壓的額定值;A(t)為包絡(luò)信號(hào);ω代表信號(hào)的工頻角頻率;M是調(diào)幅系數(shù),M通常為 0.01,最大可取到 0.1;a(t)為調(diào)制信號(hào)。調(diào)制信號(hào)表示為
式(2)中,為調(diào)制信號(hào)的角頻率,代表諧波的次數(shù)。
2.1 傅里葉變換
傅里葉變換是信號(hào)分析的最重要的方法,他架起了時(shí)域和頻域之間的橋梁。如果周期信號(hào)f(t)滿足狄里赫利條件,則的傅里葉級(jí)數(shù)可表示為:
任何周期函數(shù)都可以表示為正弦函和余弦函數(shù)的疊加,并且這些正余弦函都是一組正交的函數(shù)基[3]。
式(3)中a0為常數(shù)分量,an、bn代表頻率分量系數(shù)。其中:
傅里葉變換可以將函數(shù)f(t)分解為不同頻率的正余弦函數(shù)的組合,也說明了傅里葉變換能夠準(zhǔn)確地提取出信號(hào)的幅頻特性。
2.2 小波變換及其多分辨特性
小波變換的含義是把基本小波的函數(shù)ψ(t)做位移τ后,在不同的尺度a下和待分析信號(hào)x(t)做內(nèi)積[4]:
等效的頻域表示為
式(7)中,X(ω),ψ(t)分別是x(t),ψ(t)的傅里葉變換。
將任意L2(R)空間中的函數(shù)f(t)在小波基下展開,稱這種展開為函數(shù)f(t)的連續(xù)小波變換(Continue Wavelet Transform),其表達(dá)式為:
小波變換也是一種積分變換,Wf(a,b)是小波變換系數(shù),函數(shù) ψ(t)為基小波或母小波[5],a為尺度參數(shù),b為平移參數(shù),基小波波形表現(xiàn)為正負(fù)交替兩端振蕩衰減,均值為零?;〔?ψ(t)的傅里葉變換 ψ(w)滿足容許性條件[6]:
小波變換的多分辨率分析將信號(hào)在高頻和低頻頻帶上細(xì)化,得到越來越精細(xì)的時(shí)頻描述,可以精確地再現(xiàn)信號(hào)在不同層次上的信息,多分辨率分析的方法就是將信號(hào)的信息在不同頻帶上體現(xiàn)出來的過程,多分辨率子空間V0可以用有限子空間來逼近[7],即:
Vj是尺度為j的尺度空間,也是反映空間Vj-1內(nèi)信號(hào)近似部分的低頻子空間;Wj是小波空間,也是反映空間Vj-1內(nèi)信號(hào)細(xì)節(jié)部分的高頻子空間[8]。Vj和Wj相互正交,從公式(7)又可以看出Wj是Vj在空間中Vj-1的補(bǔ)子空間。
對(duì)任意信號(hào)f(t)∈L2(R),將其進(jìn)行多分辨率分解為:
式(12)中:
φj,k(t),ψj,k(t)分別代表尺度函數(shù)φ(t)和小波函數(shù)ψ(t)經(jīng)過尺度和位移變化后的函數(shù)集合。cj(k)為尺度系數(shù),dj(k)為小波系數(shù)。
含有突變信號(hào)的閃變模型為:
V=1V(標(biāo)幺值),Ω=14 πrad/s,M=0.04,ω=100 πrad/s采樣頻率為800 Hz,采樣點(diǎn)數(shù)為2 000,采樣時(shí)間為2.5 s,電壓閃變波形如圖1所示
圖1 含突變信號(hào)的閃變波形Fig.1 Flicker mutation signal waveform
經(jīng)過快速傅里葉(FFT)變換得到的圖形如圖2所示。
圖2 含突變信號(hào)的FFT圖形Fig.2 FFT transform graphic of mutation signals
由圖2可以看出FFT只能讀出閃變信號(hào)v(t)的幅頻特征,并沒有對(duì)閃變信號(hào)發(fā)生突變的時(shí)刻給以精確地區(qū)分,因而存在了某種程度的局限性。
根據(jù)多分辨率概念,采樣頻率為800 Hz,選用db20小波函數(shù)對(duì)閃變信號(hào)經(jīng)過第一次分解后得到a1(0~200 Hz)和d1 (200~400 Hz)兩個(gè)頻帶。再將a1分解得到a2(0~100 Hz)和d2(100~200 Hz)兩個(gè)頻帶。以此類推經(jīng)過六次分解后得到了a6(0~6.25 Hz)和d6(6.25~12.5 Hz)兩個(gè)頻帶,此時(shí)7 Hz的閃變頻率就落到了的頻帶上,低頻的頻帶中已經(jīng)沒有了閃變信號(hào)分量,因此多分辨率分解到此為止。a1到a6為信號(hào)分解后的各層低頻信號(hào)(逼近信號(hào))如圖3,d1到a6為信號(hào)分解后的各層高頻信號(hào)(細(xì)節(jié)信號(hào))如圖4。另外在圖4的d1和d2中可以觀察到閃變信號(hào)波形有兩個(gè)明顯的突變點(diǎn),分別對(duì)應(yīng)著突變的開始時(shí)刻和結(jié)束時(shí)刻。這就是小波變換優(yōu)越于傅里葉變換的突出所在。小波變換多分辨率分析能夠?qū)?fù)雜的電壓閃變信號(hào)頻率劃分到各子頻帶中,特別適合于調(diào)幅波中混合多種頻率成分的情況。這是小波的多分辨率特性在電壓閃變信號(hào)分析中的價(jià)值的最要的體現(xiàn)。
通過對(duì)閃變信號(hào)的仿真將傅里葉變換和小波變換的不同特點(diǎn)體現(xiàn)出來,仿真分析的結(jié)果可以看出小波變換在電壓閃變信號(hào)的時(shí)頻分析方面比傅里葉變換具有更強(qiáng)的分析能力。而快速傅里葉變換分析信號(hào)時(shí)只能反應(yīng)頻域上的信息,對(duì)時(shí)域不能做出任何的識(shí)別分析能力。小波的多分辨特性能將電壓閃變不同頻率的信號(hào)分離,從而能觀察到信號(hào)的細(xì)微變化,也可以觀察出突變的開始和結(jié)束時(shí)刻,因此小波變換換更適用于對(duì)電壓閃變信號(hào)的分析。小波變換為我們的電壓閃變檢測(cè)工作提供很有價(jià)值的信息,對(duì)后續(xù)閃變問題的治理也具有十分重要的實(shí)際意義。
圖3 信號(hào)v(t)小波分解后各層逼近信號(hào)Fig.3 Approximate components on each floor after the wavelet transform
圖4 信號(hào)v(t)小波分解后各層細(xì)節(jié)信號(hào)Fig.4 Detailed components on each floor after the wavelet transform
[1]劉會(huì)金,吳玉蓉,熊玲玲.基于Daubechies小波的多分辨分解在電壓閃變信號(hào)分析中的應(yīng)用[J].繼電器,2004,32(13): 6-10.
[2]張勇,管聲啟.基于小波提升分解的帶鋼表面缺陷檢測(cè)[J].西安工程大學(xué)學(xué)報(bào),2013(4):1-2.
[3]謝小磊.電壓波動(dòng)與閃變的小波分析[D].南京:河海大學(xué),2006.
[4]郭曉麗,朱志松,許小梅.短時(shí)間閃變嚴(yán)重度Pst的計(jì)算方法的研究[J].南通工學(xué)院學(xué)報(bào):自然科學(xué)版,2004,3(3): 66-69.
[5]郝曉平,黨幼云.基于FPGA的信號(hào)奇異性檢測(cè)[J].西安工程大學(xué)學(xué)報(bào),2011(6):1-2.
[6]Dwight F M,Kraig J.OLEJNIZAK.小波基礎(chǔ)及應(yīng)用教程[M].楊志華,楊力華譯.北京:機(jī)械工業(yè)出版社,2006.
[7]張志剛.電能質(zhì)量講座 第八講電壓波動(dòng)與閃變[J].低壓電器,2007(16):56-60.
[8]桑松,柴玉華,孫影.基于小波變換和快速傅里葉變換的諧波檢測(cè)[J].電測(cè)與儀表,2012,49(559):29-32.
[9]任震,黃群古,黃雯瑩.多層多分辨分析算法及在電力系統(tǒng)中的應(yīng)用[J].電網(wǎng)技術(shù),2012(6):5-7.
Detection and analysis of wavelet transform in voltage flicker signal
XU Jian,SHAO Shan
(College of Electronic Information of Xi’an Polytechnic University,Xi’an 710048,China)
In order to more accurately detect the power system voltage flicker signal,,this paper proposes the principles of multi-resolution wavelet analysis.Wavelet transform has good time-frequency localization properties and then extract the voltage flicker signal parameters such as frequency and amplitude.Finally through MATLAB simulation,it is proved that the wavelet transform is more effective in analysis of voltage flicker.
power system;Voltage Flicker;Wavelet transform;Multi-resolution analysis;MATLAB
TP
A
1674-6236(2016)04-0155-03
2015-04-14 稿件編號(hào):201504139
陜西省教育廳專項(xiàng)科研計(jì)劃(2013JK1015)
徐 ?。?963—),男,陜西西安人,碩士,副教授。研究方向:電能質(zhì)量檢測(cè)、信號(hào)與系統(tǒng)等。