零相位數(shù)字濾波器在非平穩(wěn)信號(hào)處理中的應(yīng)用

數(shù)字濾波是數(shù)字?jǐn)?shù)據(jù)處理的一般方法。普通數(shù)字濾波在濾波過(guò)程中會(huì)發(fā)生一定的相位位移。為了解決這個(gè)問(wèn)題，引入了零相位數(shù)字濾波。兩個(gè)主要的零相位數(shù)字濾波算法（之前的序列、欺詐和結(jié)果序列）和rr（不安全、不安全、結(jié)果序列和結(jié)果序列）。首先，我們研究了采用合適的濾波初期模型、信號(hào)延伸度、frr和rr的方法，改善了零相位數(shù)字濾波的邊界問(wèn)題，并對(duì)穩(wěn)定信號(hào)濾波產(chǎn)生了良好的效果。在對(duì)非平穩(wěn)信號(hào)進(jìn)行濾波時(shí),與小波分解重構(gòu)和經(jīng)驗(yàn)?zāi)B(tài)分解等無(wú)相移的非平穩(wěn)信號(hào)處理方法相比,零相位數(shù)字濾波器具有截止頻率明確、計(jì)算量小的優(yōu)點(diǎn).尤其在計(jì)算能力較弱的嵌入式應(yīng)用系統(tǒng)中,零相位數(shù)字濾波器的這一優(yōu)點(diǎn)更加突出.調(diào)幅信號(hào)、調(diào)頻信號(hào)是兩類(lèi)典型的非平穩(wěn)信號(hào),很多工程實(shí)際信號(hào)可以借助這兩種信號(hào)進(jìn)行模擬.例如,電力系統(tǒng)的次同步振蕩、電壓暫態(tài)振蕩、開(kāi)關(guān)設(shè)備動(dòng)作時(shí)的振動(dòng)、發(fā)電機(jī)組軸系的扭振等就可用調(diào)幅信號(hào)模擬;旋轉(zhuǎn)機(jī)械調(diào)速過(guò)程中的轉(zhuǎn)速則多為調(diào)頻信號(hào).文獻(xiàn)利用線(xiàn)性調(diào)頻信號(hào)對(duì)零相位數(shù)字濾波器進(jìn)行了分析.因此,本文作者將主要針對(duì)調(diào)幅信號(hào)來(lái)討論零相位數(shù)字濾波器在非平穩(wěn)信號(hào)處理中的特點(diǎn).利用調(diào)幅信號(hào)對(duì)一種典型的帶通零相位數(shù)字濾波器進(jìn)行了仿真分析;對(duì)小波包分解重構(gòu)、經(jīng)驗(yàn)?zāi)B(tài)分解及零相位數(shù)字濾波器方法的濾波效果進(jìn)行了比較;探討了在調(diào)幅信號(hào)處理中存在的過(guò)渡過(guò)程,及其對(duì)零相位數(shù)字濾波的幅值特性以及頻率特性的影響;論述了影響濾波誤差的因素;提出了利用分段零相位濾波的方法改善零相位數(shù)字濾波器過(guò)渡過(guò)程的方法;用實(shí)測(cè)的次同步振蕩信號(hào)對(duì)兩種方法進(jìn)行了驗(yàn)證.1零相位濾波、ir數(shù)字濾波器對(duì)于一維時(shí)間序列,零相位數(shù)字濾波器的基本原理是,利用正向時(shí)間序列和翻轉(zhuǎn)時(shí)間序列通過(guò)濾波器時(shí)的相移相互抵消,從而實(shí)現(xiàn)濾波結(jié)果的零相移.利用傅里葉變換的時(shí)域位移性質(zhì)和時(shí)間翻轉(zhuǎn)性質(zhì),即可從理論上論證其正確性.以FRR方法為例,將長(zhǎng)度為N+1的單邊時(shí)間序列x[n]通過(guò)沖激響應(yīng)為h[n]的濾波器,得到一次濾波序列y1[n],再將一次濾波序列進(jìn)行時(shí)序反轉(zhuǎn)、平移得到數(shù)據(jù)序列y2[n],y2[n]通過(guò)濾波器h[n]得到二次濾波序列y3[n],最后將二次濾波序列進(jìn)行時(shí)序反轉(zhuǎn)、平移得到最終濾波輸出y[n].從而,有y1[n]=x[n]*h[n](1)y2[n]=y1[Ν-n](2)y3[n]=y2[n]*h[n](3)y[n]=y3[Ν-n](4)設(shè)x[n]?X(ω)、y1[n]?Y1(ω)、y2[n]?Y2(ω)、y3[n]?Y3(ω)、y[n]?Y(ω)、h[n]?H(ω)構(gòu)成傅里葉變換對(duì),且考慮到對(duì)實(shí)序列h[n]有H(-ω)=H*(ω).依式(1～4),FRR濾波輸出的頻域表達(dá)如下式推導(dǎo)Y(ω)=e-jωΝY3(-ω)=e-jωΝY2(-ω)Η(-ω)=e-jωΝejωΝY1(ω)Η(-ω)=Y1(ω)Η(-ω)=X(ω)Η(ω)Η(-ω)=X(ω)|Η(ω)|2(5)因而,FRR方法實(shí)現(xiàn)了零相移濾波,其傳遞函數(shù)頻域表示如下ΗFRR(ω)=Η(ω)Η(-ω)(6)同理可得,RRF方法實(shí)現(xiàn)零相移濾波時(shí),其傳遞函數(shù)頻域表示如下ΗRRF(ω)=Η(-ω)Η(ω)(7)選定沖激響應(yīng)序列為h[n]的普通數(shù)字濾波器,實(shí)際使用時(shí)可以通過(guò)兩個(gè)分立的濾波過(guò)程來(lái)實(shí)現(xiàn)零相位濾波,也可以直接構(gòu)造零相位數(shù)字濾波器,其沖激響應(yīng)序列分別為hFRR[n]和hRRF[n].hFRR[n]=h[n]*h[-n](8)hRRF[n]=h[-n]*h[n](9)其中h[-n]為h[n]的時(shí)間翻轉(zhuǎn)序列.h[n]的選擇可有FIR濾波器和IIR濾波器兩種類(lèi)型.在零相位濾波器的應(yīng)用中,IIR濾波器在相頻特性上的劣勢(shì)可以忽略,且在相同階數(shù)下IIR濾波器比FIR濾波器具有更好的截止特性,因而h[n]一般選用IIR數(shù)字濾波器.23khs.自適應(yīng)變化為比較小波包分解重構(gòu)、經(jīng)驗(yàn)?zāi)B(tài)分解和零相位數(shù)字濾波器3種方法的差異.選取含有3個(gè)衰減振蕩分量的輸入信號(hào)f(t),信號(hào)采樣率為1kHz.設(shè)定3個(gè)分量的振蕩頻率分別為10Hz、23Hz、39Hz,其計(jì)算公式如下f1(t)=e-3(t-1)u(t-1)sin(20π(t-1))(10)f2(t)=e-2(t-1)u(t-1)sin(46π(t-1))(11)f3(t)=e-(t-1)u(t-1)sin(78π(t-1))(12)輸入信號(hào)f(t)的計(jì)算公式如下f(t)=f1(t)+f2(t)+f3(t)(13)1小波包的濾波、并存在重要差異利用小波包可將信號(hào)分解到不同節(jié)點(diǎn),各節(jié)點(diǎn)對(duì)應(yīng)了不同頻帶的子波,對(duì)目標(biāo)節(jié)點(diǎn)系數(shù)進(jìn)行重構(gòu)即可得到所需的濾波輸出.考慮到重構(gòu)信號(hào)的不失真和光滑,采用了bior、sym和db小波基對(duì)多分量信號(hào)進(jìn)行處理,最終選定db8小波基對(duì)f(t)進(jìn)行5層分解重構(gòu),由于信號(hào)采樣率為1kHz,小波包5層分解的各子波頻帶寬度為15.625Hz,使10Hz、23Hz、39Hz3個(gè)頻率成分位于、、節(jié)點(diǎn)對(duì)應(yīng)的頻帶內(nèi),濾波結(jié)果如圖1所示(圖中虛線(xiàn)為濾波輸出,實(shí)線(xiàn)為輸入分量).小波包分解重構(gòu)的濾波輸出,突變點(diǎn)附近存在過(guò)渡過(guò)程,其余部分信號(hào)不光滑,在幅值和相位上都存在抖動(dòng),其中23Hz和39Hz分量抖動(dòng)比較嚴(yán)重.2mt經(jīng)驗(yàn)?zāi)B(tài)分解mt模式,m經(jīng)驗(yàn)?zāi)B(tài)分解方法是自適應(yīng)的將信號(hào)分解為不同時(shí)間尺度的本征模態(tài)函數(shù)(IMF)分量和趨勢(shì)項(xiàng),本征模態(tài)函數(shù)是一類(lèi)零均值的調(diào)幅調(diào)頻信號(hào).對(duì)含隨機(jī)寬帶噪聲的信號(hào),經(jīng)驗(yàn)?zāi)B(tài)分解是類(lèi)小波的二進(jìn)濾波器組,第一個(gè)本征模態(tài)函數(shù)的頻帶最高.由于m(t)由3個(gè)調(diào)幅正弦信號(hào)組成,經(jīng)驗(yàn)?zāi)B(tài)分解得到的前3個(gè)本征模態(tài)函數(shù)分量即對(duì)應(yīng)其濾波輸出.其中,第一個(gè)本征模態(tài)函數(shù)對(duì)應(yīng)39Hz分量、第二個(gè)本征模態(tài)函數(shù)對(duì)應(yīng)23Hz分量、第三個(gè)本征模態(tài)函數(shù)對(duì)應(yīng)10Hz分量,濾波結(jié)果如圖2所示(圖中虛線(xiàn)為濾波輸出,實(shí)線(xiàn)為輸入分量).經(jīng)驗(yàn)?zāi)B(tài)分解方法的濾波輸出,突變點(diǎn)附近存在過(guò)渡過(guò)程,但其余部分信號(hào)很光滑,幅值和相位的穩(wěn)定性強(qiáng)于小波包分解重構(gòu)方法,但是仍然存在抖動(dòng).我們又采用標(biāo)準(zhǔn)差為0.005、聚合次數(shù)為50次的聚合經(jīng)驗(yàn)?zāi)B(tài)分解對(duì)f(t)進(jìn)行濾波,發(fā)現(xiàn)能改善過(guò)渡過(guò)程.但是,同時(shí)造成了39Hz分量混疊到23Hz分量中,為減少混疊效應(yīng)需要增加聚合次數(shù),極大的增加了計(jì)算量.3過(guò)渡過(guò)程中的誤差檢驗(yàn)數(shù)字濾波器在對(duì)有限長(zhǎng)信號(hào)進(jìn)行濾波時(shí),相當(dāng)于給信號(hào)加了一個(gè)矩形窗,在濾波輸出信號(hào)上造成衰減振蕩.零相位數(shù)字濾波器進(jìn)行了兩次濾波,在輸出信號(hào)的首尾兩端形成失真.FRR方法第二次濾波的首端是輸入信號(hào)的尾端,濾波輸出的首端失真較小;RRF方法第二次濾波的首端是輸入信號(hào)的首端,濾波輸出的尾端失真較小.對(duì)于本文考慮的調(diào)幅正弦信號(hào),通過(guò)在輸入信號(hào)首尾兩端進(jìn)行信號(hào)延拓的方法,來(lái)消除端點(diǎn)效應(yīng).依據(jù)對(duì)端點(diǎn)效應(yīng)誤差的要求,延拓的數(shù)據(jù)長(zhǎng)度應(yīng)比擬于過(guò)渡過(guò)程時(shí)間,本文取1s長(zhǎng)度的數(shù)據(jù)點(diǎn).消除端點(diǎn)效應(yīng)的FRR方法和RRF,在使用中的差異可忽略不計(jì).本文的后續(xù)討論中,零相位數(shù)字濾波器均采用FRR方法實(shí)現(xiàn),且采用了首尾端的端點(diǎn)延拓,不再特別說(shuō)明.綜合考慮濾波器階躍響應(yīng)的過(guò)渡過(guò)程時(shí)間和振蕩次數(shù)等因素,選定用4階巴特沃斯帶通濾波器(直接型系數(shù)見(jiàn)表1)構(gòu)造零相位數(shù)字濾波器,中心頻率分別為10Hz、23Hz、39Hz,帶寬為3Hz,記為hbp10、hbp23、hbp39.濾波結(jié)果見(jiàn)圖3(圖中虛線(xiàn)為濾波輸出,實(shí)線(xiàn)為輸入分量).濾波輸出很光滑,除了突變點(diǎn)附近,幅值和相位與原始信號(hào)重合,效果好于小波包分解重構(gòu)和經(jīng)驗(yàn)?zāi)B(tài)分解方法.不考慮信號(hào)的突變,則零相位數(shù)字濾波器實(shí)現(xiàn)了較嚴(yán)格的零相移濾波.3種方法都存在過(guò)渡過(guò)程.以10Hz分量的誤差為例來(lái)考察過(guò)渡過(guò)程,將3種方法所得結(jié)果前3s的誤差示于圖4中.從上至下依次為小波包分解重構(gòu)、經(jīng)驗(yàn)?zāi)B(tài)分解、零相位濾波器誤差,實(shí)線(xiàn)為原始誤差,虛線(xiàn)為幅值放大100倍后的誤差.其中,經(jīng)驗(yàn)?zāi)B(tài)分解誤差極值最大,零相位濾波器與小波包分解重構(gòu)的誤差極值較小;零相位濾波與經(jīng)驗(yàn)?zāi)B(tài)分解誤差穩(wěn)定時(shí)間較短,小波包分解重構(gòu)的誤差穩(wěn)定時(shí)間較長(zhǎng);但零相位濾波器反時(shí)序的特性造成突變點(diǎn)前有信號(hào)的提前輸出.下面我們將對(duì)零相位濾波器的過(guò)渡過(guò)程詳細(xì)探討,并研究其幅值和頻率特性,從而對(duì)其過(guò)渡過(guò)程性能進(jìn)行改善.3過(guò)渡過(guò)程3.1同上轉(zhuǎn)臺(tái)濾波過(guò)濾算法零相位數(shù)字濾波器解決了濾波中存在的相移問(wèn)題,但是在對(duì)非平穩(wěn)信號(hào)進(jìn)行濾波時(shí),仍需注意濾波中可能存在的過(guò)渡過(guò)程問(wèn)題.尤其是在信號(hào)發(fā)生突變的位置影響較大,這是因?yàn)橥蛔凕c(diǎn)的存在類(lèi)似于將信號(hào)在該點(diǎn)截?cái)?引入了新的端點(diǎn).考察一個(gè)在第2秒時(shí)刻幅值由0躍變?yōu)?在第8秒再跳變?yōu)?的10Hz正弦信號(hào)x(t),其計(jì)算公式如下x(t)=sin(20π(t-2))[u(t-2)-u(t-8)](14)其中u(t)為單位階躍信號(hào).將x(t)通過(guò)hbp10進(jìn)行濾波,濾波結(jié)果如圖5所示,存在過(guò)渡過(guò)程.為了更清楚地表達(dá)過(guò)渡過(guò)程,將其幅值響應(yīng)曲線(xiàn)示于圖6.輸入信號(hào)的突變?cè)跒V波輸出的對(duì)應(yīng)位置上存在過(guò)渡過(guò)程,穩(wěn)定到±2%的時(shí)間約0.48s,最大超調(diào)量為3.4%;進(jìn)入平穩(wěn)階段后,濾波輸出與輸入信號(hào)能很好的吻合,誤差均值在10-5量級(jí),而誤差的方差在10-7量級(jí),平穩(wěn)階段時(shí)間越長(zhǎng),誤差越小.對(duì)這類(lèi)幅值發(fā)生階躍的正弦信號(hào),可以采用濾波輸出信號(hào)最大幅值的中值時(shí)刻來(lái)近似定位突變點(diǎn),從而去除突變點(diǎn)前的虛假信號(hào).采用該方法估計(jì)x(t)的突變點(diǎn)的估計(jì)誤差只有0.005s,實(shí)際應(yīng)用中受波形形態(tài)及噪聲的影響,誤差變大.事實(shí)上,不同特點(diǎn)的信號(hào)通過(guò)同一零相位數(shù)字濾波器,濾波信號(hào)幅值和頻率表現(xiàn)出的特性也會(huì)不同.為進(jìn)行具體的仿真分析,選取衰減振蕩、發(fā)散振蕩、發(fā)散振蕩經(jīng)采取措施后又衰減振蕩這3種情況進(jìn)行研究,利用幅值衰減、幅值發(fā)散、幅值先發(fā)散后衰減這3種不同調(diào)幅正弦信號(hào)進(jìn)行模擬.3.23濾波器hbp10濾波特性先構(gòu)造調(diào)幅正弦信號(hào)的幅值包絡(luò)a1(t)、a2(t)、a3(t),其中a1(t)為衰減幅值、a2(t)為發(fā)散幅值、a3(t)為先發(fā)散后衰減幅值,其計(jì)算公式如下a1(t)=e-(t-1)u(t-1)(15)a2(t)=e(t-9)[1-u(t-9)](16)a3(t)=e(t-5)[1-u(t-5)]+e-(t-5)u(t-5)(17)對(duì)應(yīng)的調(diào)幅正弦信號(hào)分別為m1(t)、m2(t)、m3(t),其計(jì)算公式如下m1(t)=a1(t)sin(20πt)(18)m2(t)=a2(t)sin(20πt)(19)m3(t)=a3(t)sin(20πt)(20)采用一個(gè)高斯白噪聲信號(hào)wn(t)作為噪聲干擾信號(hào),wn(t)的統(tǒng)計(jì)分布為N(0,0.01).對(duì)應(yīng)的含噪聲的輸入信號(hào)分別為s1(t)、s2(t)、s3(t),其表達(dá)式分別為下列公式所示s1(t)=m1(t)+wn(t)(21)s2(t)=m2(t)+wn(t)(22)s3(t)=m3(t)+wn(t)(23)s1(t)、s2(t)、s3(t)的波形如圖7所示,經(jīng)過(guò)零相位帶通數(shù)字濾波器hbp10濾波后,濾波輸出信號(hào)如圖8所示.可明顯看出,這3種信號(hào)的濾波輸出的最大幅值小于的x(t)的濾波輸出的最大幅值;過(guò)渡過(guò)程時(shí)間也小于x(t)的過(guò)渡過(guò)程時(shí)間;其中s3(t)的濾波結(jié)果最好,誤差均值和方差都在10-5量級(jí).可見(jiàn)輸入信號(hào)的形態(tài)特點(diǎn)對(duì)濾波結(jié)果有影響.4振幅值和頻率特性4.1濾波結(jié)果的頻率方差濾波輸出的幅值和頻率可利用希爾伯特變換求得.其基本原理如下.對(duì)信號(hào)ν(t),其希爾伯特變換由下式定義ν?(t)=1π∫-∞+∞ν(τ)t-τdτ(24)對(duì)應(yīng)的幅值包絡(luò)a(t)可按下式計(jì)算a(t)=ν2(t)+ν?2(t)(25)對(duì)應(yīng)的相位θ(t)可按下式計(jì)算θ(t)=arctanν?(t)ν(t)(26)對(duì)應(yīng)的瞬時(shí)頻率定義為f(t)=12πdθ(t)d(t)(27)下面利用希爾伯特變換獲取濾波結(jié)果的幅值包絡(luò)和頻率,對(duì)零相位數(shù)字濾波器的幅值響應(yīng)和頻率響應(yīng)特點(diǎn)進(jìn)行詳細(xì)分析.4.2放大電路一般特性s1(t)、s2(t)、s3(t)經(jīng)過(guò)零相位帶通數(shù)字濾波器hbp10濾波后,理想的幅值應(yīng)為a1(t)、a2(t)、a3(t).實(shí)際按式(25)計(jì)算出的濾波輸出信號(hào)幅值分別如圖9～圖11所示.可明顯看出,s1(t)、s2(t)、s3(t)的濾波輸出的最大幅值分別只有0.81、0.80和0.93,明顯小于的x(t)的濾波輸出的最大幅值1.034;s3(t)的最大值沒(méi)有明顯延遲,s1(t)、s2(t)過(guò)渡過(guò)程時(shí)間為0.22s、0.22s,也明顯小于對(duì)x(t)濾波時(shí)的過(guò)渡過(guò)程時(shí)間0.48s;信號(hào)幅值較小的部分噪聲影響比較大.由于不像x(t)有穩(wěn)定的信號(hào)階段可以使濾波輸出完成過(guò)渡達(dá)到穩(wěn)定,上述3類(lèi)信號(hào)衰減系數(shù)越大幅值損失越大.將不同衰減系數(shù)的10Hz分量濾波輸出的幅值示于圖12中,其結(jié)果印證了該結(jié)論.對(duì)于零相位帶通數(shù)字濾波器而言,與普通帶通數(shù)字濾波器相似,其幅值響應(yīng)具有如下特點(diǎn):1)濾波器類(lèi)型一定時(shí),較寬的頻帶可以獲得較大的幅值和較小的時(shí)延;提高濾波器階數(shù)可加速阻頻帶衰減,減少臨近頻帶的干擾,但也會(huì)帶來(lái)通帶幅值的振蕩.2)輸入信號(hào)幅值包絡(luò)變化時(shí),濾波輸出的時(shí)延從幅值包絡(luò)變化率大的部分偏向變化率小的部分;幅值包絡(luò)變化越平緩,濾波后的時(shí)延也越小;濾波器階數(shù)大于等于4時(shí),信號(hào)幅值包絡(luò)越接近階躍,超調(diào)量越大.3)因?yàn)榘自肼暉o(wú)法完全濾除,幅值較小的部分噪聲影響比較明顯.4.3濾波、濾波效率及信號(hào)背景分析按式(27)式定義對(duì)m1(t)、m2(t)、m3(t)的濾波輸出信號(hào)進(jìn)行計(jì)算,得到的瞬時(shí)頻率分別示于圖13～圖15中.其中,首尾兩端的振蕩是由于希爾伯特變換的邊界效應(yīng)造成的.另外需要注意的地方是,在信號(hào)突變點(diǎn)附近(約為一個(gè)過(guò)渡過(guò)程時(shí)間),濾波輸出存在較大的頻率突變,說(shuō)明濾波輸出信號(hào)在該位置存在畸變,在圖13～圖14中分別位于0.755s、9.246s處.取中間頻率穩(wěn)定的3～8s時(shí)間段進(jìn)行分析,可得頻率均值為10.000Hz,方差小于10-5量級(jí).其中,m3(t)在中間突變點(diǎn)處的頻率處于9.99Hz與10.01Hz之間,明顯好于m1(t)、m2(t)在突變點(diǎn)處的濾波結(jié)果.為考察噪聲對(duì)頻率響應(yīng)的影響,取s1(t)信號(hào)為例.圖16給出了s1(t)經(jīng)濾波后的瞬時(shí)頻率曲線(xiàn),可見(jiàn)濾波結(jié)果變差,信號(hào)幅值較小的部位上噪聲造成相位的波動(dòng),使波形相對(duì)發(fā)生了較大的畸變.仍取3～8s時(shí)間段進(jìn)行分析,可得頻率均值為9.962Hz,方差在10-1量級(jí).零相位帶通數(shù)字濾波器的頻率響應(yīng)具有如下特點(diǎn):1)中心頻率和頻帶寬度變化對(duì)穩(wěn)定的頻率段影響很小,主要對(duì)信號(hào)突變點(diǎn)附近的頻率有影響.2)當(dāng)輸入信號(hào)的幅值包絡(luò)變化時(shí),濾波輸出的頻率畸變點(diǎn)位置與幅值時(shí)延的方向相反,位于突變點(diǎn)偏幅值包絡(luò)變化率較大的一側(cè).3)噪聲的存在,造成信號(hào)的畸變,濾波輸出頻率存在波動(dòng),幅值較小的部分噪聲影響比較明顯.5分段零相位濾波由上述仿真結(jié)果可知,m3(t)信號(hào)的濾波結(jié)果明顯好于m1(t)與m2(t).m3(t)與m1(t)、m2(t)的不同之處在于,其突變點(diǎn)兩邊的幅值是對(duì)稱(chēng)的.基于此,作者提出分段零相位濾波方法,用以改善零相位濾波器在處理存在突變的調(diào)幅信號(hào)時(shí)的過(guò)渡過(guò)程性能,其具體步驟如下:1)對(duì)信號(hào)進(jìn)行零相位濾波.2)在濾波幅值包絡(luò)的中值點(diǎn)附近,搜索濾波輸出中絕對(duì)值最小的點(diǎn),進(jìn)行突變點(diǎn)的近似定位.3)將輸入信號(hào)按定位結(jié)果分段,對(duì)分段后的信號(hào)進(jìn)行端點(diǎn)延拓,再分別進(jìn)行零相位濾波.4)將分段濾波結(jié)果組合,形成濾波輸出.零相位濾波器采用hbp10,利用零相位濾波和分段零相位濾波方法分別對(duì)s1(t)進(jìn)行濾波,濾波結(jié)果示于圖17中.從圖中可以明顯觀察到,分段零相位濾波的結(jié)果不但比原始的零相位濾波結(jié)果過(guò)渡過(guò)程縮短,而且也解決了反時(shí)序提前輸出的問(wèn)題.可見(jiàn),分段零相位濾波能縮短過(guò)渡過(guò)程時(shí)間,消除反時(shí)序提前輸出,極大改善零相位濾波器的濾波效果.其缺點(diǎn)是,增加了濾波次數(shù)和端點(diǎn)延拓次數(shù),增加了計(jì)算量.6零相位數(shù)字濾波器和分段零相位數(shù)字濾波器為從次同步振蕩中分離出主要的振蕩分量,將該方法應(yīng)用到次同步振蕩信號(hào)的濾波中.在某次試驗(yàn)中測(cè)得振蕩信號(hào),采樣率為1kHz,數(shù)據(jù)長(zhǎng)度為95s,如圖18所示.通過(guò)時(shí)頻分析知該振蕩信號(hào)中主要包含有13Hz、50Hz、25Hz等幾個(gè)較大的頻率成分,其中在次同步頻域內(nèi)13Hz振幅最大.以13Hz分量為例,仍采用4階巴特沃斯帶通濾波器(直接型的中心頻率為13Hz,帶寬3Hz),構(gòu)造零相位