光纖陀螺分形噪聲的預白化濾波方法

華春紅 任 章 張敏虎

(北京航空航天大學 自動化科學與電氣工程學院,北京 100191)

光纖陀螺分形噪聲的預白化濾波方法

華春紅 任 章 張敏虎

(北京航空航天大學 自動化科學與電氣工程學院,北京 100191)

光纖陀螺分形噪聲往往具有長相關特性,導致陀螺輸出信號慢漂移.為了濾除這類分形噪聲,提出了一種分數階預白化差分方法,通過該方法可以將分形噪聲轉化為高斯白噪聲;利用小波變換,在小波變換域中采用 Bayes軟閾值去噪方法去除該白化噪聲,從而形成了一種光纖陀螺分形噪聲的預白化濾波方法.對 VG949p型光纖陀螺實測數據處理結果表明,該方法能夠有效地去除其中的分形噪聲,抑制了陀螺低頻漂移,同時也濾掉了其高頻噪聲項.

光纖陀螺;隨機誤差;差分方程;小波分析

光纖陀螺隨機漂移噪聲主要包括白噪聲和1/fγ分形噪聲,其中分形噪聲往往具有長相關特性[1-2],采用自回歸-滑動平均模型描述不合適[3].文獻[4]對光纖陀螺分形噪聲采用多尺度白化濾波結構,其白化濾波參數是利用最大似然估計法估計的分形譜指數 γ.該估計算法復雜,且當存在其他噪聲干擾時,僅當譜指數 γ∈(0,1)時才是相合估計,否則以概率 1收斂到錯誤結果[5].文獻[5]提出的基于分形參數的小波軟閾值去噪方法難以很好地解決由長相關分形噪聲引起光纖陀螺慢漂移問題,因為該漂移處于低頻噪聲段,而小波閾值去噪方法只能濾掉高頻噪聲分量.因此如果能把長相關分形噪聲轉化為白噪聲,再采用小波閾值去噪可將其濾除.

針對以上問題,本文從信號長相關的角度來描述光纖陀螺隨機漂移噪聲,并提出了一種有效的預白化濾波方法.該方法僅通過分數階差分方法就能將光纖陀螺分形噪聲轉化為高斯白噪聲,有效地解決了其低頻漂移問題,為其低頻噪聲的濾波提供了一條新思路.

1 分形過程預白化差分方法

1.1 連續域分形過程

連續域內分形過程 x(t)描述為分數布朗運動 (FBM,Fractional Brownian Motion),其模型BH(t)為

將式(2)代入式(1),整理可得

由式(3)可以看出 BH(t)有高斯白噪聲的分數積分形式:

其中,q為分數因子.

將式(4)中 BH(t)表示的分形過程 x(t)改寫為分數微分形式:

由式(5)可以看出分形過程 x(t)的 q階微分結果為高斯白噪聲 ω(t).對式(5)進行傅里葉變換,得到功率譜密度為

1.2 離散域分形過程

用差分算子替代微分算子,式(5)的分數微分形式可以轉化為離散域的分數差分形式[7]:

式中,Δq為后向差分算子:

式中,L為滯后算子.

文獻[8-9]又將離散域分形序列{x(n)}描述為高斯白噪聲序列{ω(n)}的 -q階差分過程:

式(8)等價于下式:

由式(7)和式 (9)看出,1/fγ分形序列{x(n)}經過 q階差分后轉化為高斯白噪聲序列,即通過q階差分實現了 1/fγ分形序列的白化濾波.

1.3 分數因子 q的估計

由 1.1節可知,q與 1/fγ分形過程的長相關Hurst指數 H滿足下式:

通常估計 H指數的方法有:聚合時間序列方差法、周期圖法、重標度極差 R/S(the Rescaled adjusted range Statistic)分析法等[10],其中 R/S分析法是分析長相關指數的經典方法,算法如下.

假設給定的觀察序列為{xi,i=1,2,…,N},將其分為 k個長度為 n的區間,即 N=k×n.在可變長度為 n子區間內的累積極差 Rn和標準差Sn為

H即為所求的長相關 Hurst指數,當 0.5＜H＜1時,則該序列具有長相關性.但是當序列中同時存在長相關和短相關過程時,R/S分析法不能將二者區分開.為了解決這個問題,A.W.Lo于1991年提出了 MR/S(the Modified R/S)分析法.

MR/S分析法不再將時間序列分段,其累積極差 Rn和 R/S分析法的 Rn計算方法相同,只不過在標準差 Sn中增加了權重因子 λj(ρ):

1.4 q階預白化差分

采用 R/S分析法或 MR/S分析法估計得到 H后,根據式 (10)計算 q.假設差分后序列為{y(n)},則分數階差分過程如下:

對于任何 q＞-1的實數有

式中

當 q確定后,式(16)只是 k的函數[11].

因此 q階差分過程可以通過序列{x(n)}和序列{h(n)}的卷積實現.

2 小波軟閾值去噪

設光纖陀螺隨機噪聲采樣序列經過上述 q階預白化差分后,得到 f(n)可以表示為有用信號s(n)和差分后增強了的白噪聲 z(n)疊加:

式中,z(n)=σze(n)是預白化差分后方差為的廣義平穩白噪聲.

2.1 閾值估計

設預白化差分后的信號 f(n)小波系數服從零均值的廣義高斯分布 GGσs,β,如果求 Bayes風險最小的條件下理想閾值會使閾值去噪的誤差進一步減少,則 Bayes閾值估計值為[12]

σs為信號的標準差,對于式(18)的含噪信號模型有

2.2 軟閾值去噪算法

對含噪信號 f(n)進行軟閾值去噪,其過程如下:

1)將 f(n)進行多尺度小波分解,得到各尺度小波系數粗糙系數,j為分解尺度;

3 實驗結果與分析

光纖陀螺 VG949p靜止采樣 1h,采樣周期為1.0 s,取其中 2048 s數據進行分析,去除其均值并進行單位換算得到采樣序列標準差為 14.622,如圖 1所示.

分別采用 R/S和 MR/S分析方法計算采樣序列的 Hurst指數,結果為 0.8038和 0.7483,可見該序列具有長相關特性.取 MR/S分析法估計的 H指數,計算 q=H+0.5=1.2483,則預白化差分結果為,如圖 1所示,其標準差為 13.251.

圖1 采樣序列與預白化差分結果

將預白化差分結果 f(n)以 db5小波做 5尺度分解,采用小波域 Bayes軟閾值去噪方法去噪,見圖 2預白化濾波結果,其標準差為0.103.就同一組數據,采用文獻[5]基于分形參數的小波閾值去噪方法,簡稱為分形參數法,濾波結果見圖 2分形參數法結果,其標準差為11.844,顯然其低頻漂移仍然存在.

圖2 預白化濾波和分形參數法去噪結果

對本文預白化濾波結果和分形參數法去噪結果進行 Allan方差對比分析,如圖 3所示,從圖中可以看出,分形參數法能夠有效地濾除其中的高頻噪聲成分,但其低頻噪聲項沒有改變.而預白化濾波方法不但能很好地濾掉其高頻噪聲,而且也能有效地抑制光纖陀螺低頻噪聲項.

圖3 Allan方差綜合分析對比圖

4 結 論

分形噪聲是光纖陀螺主要噪聲源,其中長相關分形噪聲引起的慢漂移是光纖陀螺隨機漂移的主要成分,產生的誤差耦合到低頻慣性信號中,因此該類分形噪聲的濾波尤其重要.為了濾除這類噪聲,本文利用其長相關特征來分析其白化濾波的方法,并對陀螺輸出靜態隨機誤差進行處理,結果表明其各噪聲項都得到了有效抑制.

由于光纖陀螺隨機漂移具有時變性,在動態條件下采用本文方法進行在線濾波,需要就其全動態范圍進行隨機誤差測試,建立輸入角速率和長相關 Hurst指數的模型,并要估計最佳在線預白化差分濾波器的長度.

References)

[1]Skaloud J,Bruton A M,Schwarz K P.Detection and filtering of short-term(1/fγ)noise in inertial sensors[J].Journal of the Institute of Navigation,1999,46(2):97-107

[2]Han Junliang,Ge Shengmin,Shen Yi,et al.Modeling and simulation of digital closed-loop fiber optic gyroscope[C]//Proceedings of the 6th World Congress on Intelligent Control and Automation.Dalian:IEEE,2006:1659-1663

[3]Torre K,Delignières D,Lemoine L.Detection of long-range dependence and estimation of fractal exponents through ARFIMA modeling[J].British Journal of Mathematical and Statistical Psychology,2007,60:85-106

[4]Gao Yukai,Deng Zhenglong.Extraction of FOG signal from fractal noise[C]//Proceedings of the 2003 International Conference on Neural Networks and Signal Processing.Nanjing:IEEE,2003:768-771

[5]武華伍,任章.光纖陀螺中分形噪聲的參數估計和去除[J].北京航空航天大學學報,2008,34(6):630-633 Wu Huawu,Ren Zhang.Parameter estimation and elimination to fractal noise of fiber optic gyro[J].Journal of Beijing University of Aeronautics and Astronautics,2008,34(6):630-633(in Chinese)

[6]楊福生.小波變換的工程分析與應用[M].北京:科學出版社,2003:217-241 Yang Fusheng.Engineering analysis and application of wavelet transform[M].Beijing:Science Industry Press,2003:217-241(in Chinese)

[7]陳彥輝,謝維信.隨機分形信號參數的分數差分估計[J].電子與信息學報,2001,23(1):9-15.Chen Yanhui,Xie Weixin.Fractional difference estimation for the parameters of random fractal signal[J].Journal of Electronics and Information Technology,2001,23(1):9-15(in Chinese)

[8]Deriche M,Tewfik A H.Maximum likelihood estimation of the parameters of discrete fractionally differenced Gaussian noise process[J].IEEE Transactions on Signal Processing,1993,41(10):2977-2989

[9]Rao R,Lee S,Bayraktar E,et al.Signal processing models for discrete-time self-similar and multifractal processes[C]//The thirty-seventh Asilomar Conference on Signals Systems and Computers.California:IEEE,2003,1:8-12

[10]Xie Chi,Yang Ni.Analysis and empirical research on the long memory property of exchange rate time series[C]//International conference on wireless communications,networking and mobile computing.Shanghai:IEEE,2007:5911-5914

[11]Jin Xiu,Yao Jin.Empirical study of ARFIMA model based on fractional differencing[J].Physica A,2007,377:138-154

[12]成禮智,王紅霞,羅永.小波的理論與應用[M].北京:科學出版社,2004:269-286 Cheng Lizhi,Wang Hongxia,Luo Yong.The theory and application of wavelet[M].Beijing:Science Press,2004:269-286(in Chinese)

(編 輯:趙海容)

Pre-whitening filter of fractalnoise in fiber op tic gyroscope

Hua Chunhong Ren Zhang Zhang Minhu

(School of Automation Science and Electrical Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)

Fractal noise in fiber optic gyroscope(FOG)usually is long-term dependent,which causes FOG a rather slow drift.In order to filter thisnoise,a pre-whitening filter based on fractional difference(FD)was proposed.By this filter,fractal noise could be transformed into Gaussian white noise(GWN).And through wavelet analysis,the pre-whitening noise could be eliminated by Bayes soft-threshold denoising in waveletdomain.On the two aforementioned steps,it came into a pre-whitening filter method of fractal noise.After the test data processing of FOG-VG949p,the results show that its fractal noise was removed efficiently,so the slow drift is restrained.The high frequency noise terms of FOG also greatly decrease.

fiber optic gyroscope(FOG);random errors;difference equations;wavelet analysis

V 241.5

A

1001-5965(2011)02-0216-04

2009-12-24

航空科學基金資助項目(20070851011)

華春紅(1976-),女,江西崇仁人,博士生,huachh@asee.buaa.edu.cn.