王澤眾,劉鋒,黃宇,向崇文
離散周期Wigner-Hough變換及其檢測性能分析?
王澤眾1,劉鋒2,黃宇1,向崇文1
(1.海軍航空工程學院研究生管理大隊,山東煙臺264001;2.海軍航空工程學院電子信息工程系,山東煙臺264001)
分析了周期Wigner-Hough變換(PWHT)進行離散化計算時,對線性調(diào)頻連續(xù)波(LFMCW)信號的檢測性能。通過推導高斯白噪聲中LFMCW信號經(jīng)離散化PWHT后信噪比處理增益表達式,分析了離散PWHT的弱信號檢測性能。通過推導離散PWHT較離散Wigner-Hough變換(WHT)的檢測性能優(yōu)勢表達式,比較了離散PWHT與離散WHT算法的檢測性能。仿真實驗驗證了離散PWHT的信噪比處理增益具有準相干積累的效果,PWHT較WHT更優(yōu)的檢測性能。
時頻分析;周期Wigner-Hough變換;LFMCW;檢測性能分析
連續(xù)波(CW)雷達信號工作比為1,在相同的探測性能要求下,較脈沖雷達信號有更低的功率,雷達偵察系統(tǒng)對于連續(xù)波雷達信號的截獲比脈沖雷達信號更加困難[1]。單載頻連續(xù)波雷達信號可以實現(xiàn)速度測量,但是其距離分辨力較低,不適用于高性能的雷達系統(tǒng)。線性調(diào)頻連續(xù)波(LFMCW)信號既可以獲取遠程目標距離信息,同時又有較高的速度分辨率,又適合使用FFT處理器快速處理,成為高性能調(diào)頻連續(xù)波(FMCW)雷達[2-3]廣泛使用的雷達信號。并且,由于此類雷達信號輻射功率低,調(diào)制規(guī)律復雜,現(xiàn)有雷達偵察系統(tǒng)對此類信號的截獲性能亟待提高。對于LFMCW信號的有效截獲,是目前雷達偵察領(lǐng)域的研究熱點之一。
針對此問題,國內(nèi)外進行了廣泛研究。具有代表性的,文獻[4]采用鏡像濾波器組(QMFB)對LFMCW信號的檢測進行了研究,但是QMFB不具抗噪功能,需要首先對信號進行消噪處理,導致信號檢測與參數(shù)提取能力受消噪能力影響較大。文獻[5-8]在時頻平面上,通過Wigner-Hough變換(WHT),將LFMCW信號看成是多個LFM脈沖信號進行了信號檢測研究。WHT是LFM脈沖信號的廣義似然比檢測(GLRT)[9],即無任何先驗信息條件下的最優(yōu)檢測,但是在LFMCW信號持續(xù)時間內(nèi),其頻率調(diào)制的變化趨勢并不是持續(xù)的線性調(diào)頻(LFM),而是呈現(xiàn)出調(diào)頻率周期間斷的LFM特征,多個間斷LFM信號的WHT存在交叉項,將WHT應(yīng)用于LFMCW信號的檢測,并不是最優(yōu)檢測,并且其運算復雜度相對比較高。文獻[10-12]采用基于分數(shù)階Fourier變換(FRFT)的快速算法,在大大降低運算復雜度的情況下,取得了同WHT基本相似的LFMCW信號檢測能力,同樣無法克服多個間斷LFM信號的交叉項問題。
考慮偵察接收機的一次觀測時間,當觀測時間變長時,進入接收機的LFMCW信號有更多的調(diào)制周期,引起信號能量增大。理論上說,LFMCW信號的檢測能力應(yīng)隨進入的LFM調(diào)制周期數(shù)增加而相應(yīng)提高。但是,在LFMCW持續(xù)時間內(nèi),其頻率調(diào)制周期性變化并不是持續(xù)的LFM,基于QMFB、WHT或FrFT的LFMCW檢測的算法,并沒有隨著觀測時間的延長而產(chǎn)生輸出信噪比的改進。
周期Wigner-Hough變換(PWHT)是2010年由Francis G.Geroleo[13]提出的,將LFMCW信號的周期性LFM特征考慮在內(nèi),可以實現(xiàn)對LFMCW信號的廣義似然比檢測和極大似然估計。工程應(yīng)用和仿真分析中,往往通過計算機進行數(shù)值計算,需要研究PWHT的離散化算法。
對于離散化的PWHT,其性能直接影響工程應(yīng)用,本文以離散PWHT的LFMCW檢測性能為對象開展研究。首先,對LFMCW的離散PWHT進行了
研究,得出了其數(shù)值與采樣點數(shù)的關(guān)系。然后,對LFMCW經(jīng)PWHT處理后的輸入輸出信噪增益進行了推導,得出了關(guān)系表達式,分析了采樣頻率、采樣時間對輸入輸出信噪增益的影響。接下來,根據(jù)Wigner-Hough變換與周期Wigner-Hough變換的聯(lián)系,得到了周期Wigner-Hough變換較Wigner-Hough變換的檢測能力優(yōu)勢表示。最后,通過Matlab仿真驗證了離散化周期Wigner-Hough變換對LFMCW信號的檢測性能。
設(shè)采樣頻率為fs,采樣點數(shù)為N,LFMCW信號LFM周期為T,帶寬為BW,則調(diào)頻率Δ=BW/T,在采樣點數(shù)的觀測時間內(nèi),有M=N/(T·fs)個LFM調(diào)制周期。信號表示為
式中,x(n)表示LFMCW信號;n為采樣點,n=0,…,N-1。
將連續(xù)PWHT[13]的瞬時自相關(guān)函數(shù)和匹配函數(shù)進行離散化,得到:
式中,n和k表示時間采樣點。則任意離散信號y(n)的PWHT可以表示為
上式應(yīng)用于數(shù)值計算,需要確定n、k的求和范圍。由式(3)可知:
對此不等式方程組求解,可以得到:
當LFMCW雷達信號不含噪聲,則有:
LFMCW信號離散PWHT在參數(shù)集f?i=fi、?Δ= Δ、?τbias=τbias、?T=T處得到最大值,使LFMCW信號在坐標f?i=fi、?Δ=Δ、?τbias=τbias和?T=T處有最好的能量聚集性,并且峰值點幅度符合采樣點數(shù)N的二次方關(guān)系。隨著N的增大,峰值數(shù)值取得更大的增長。
由于隨機噪聲的存在,信號y(n)=x(n)+ ω(n)離散PWHT的峰值PWHT(fi,Δ,τbias,T)成為隨機變量,勢必影響對于LFMCW信號的檢測。存在噪聲的情況下,有必要分析離散PWHT算法對于LFMCW信號檢測性能。
假定噪聲為加性循環(huán)高斯白噪聲,且LFMCW信號與噪聲不相關(guān)。因為噪聲中LFMCW信號檢測是在離散PWHT處理后進行的,輸出信噪比SNRout與輸入信噪比SNRin有一定的關(guān)系,即處理增益。通過離散PWHT的處理增益對其檢測性能進行分析。
離散PWHT處理后噪聲污染的LFMCW信號輸出信噪比SNRout表示如下:
式中,分子表示LFMCW信號本身的離散PWHT,由式(9)可知其值為N2A2/22=N4A4/4;分母表示LFMCW信號加噪聲信號的離散PWHT的方差,其計算公式為
式中,E[·]表示數(shù)學期望。
加性循環(huán)高斯白噪聲,即噪聲的實部和虛部獨立同分布,均為高斯白噪聲。設(shè)ω(n)的實部和虛部均服從均值μ為0、方差為σ2ω的正態(tài)分布,則噪聲的離散自相關(guān)為
式中,δ(k)表示單位沖激函數(shù)。
由式(5)得:
將式(1)、(2)和(4)代入上式,得到:
其中:
綜合式(16)和(17),得到
則:
接下來求離散PWHT平方的期望,其表示式為
式中,x1=x( n+k),x2=x( n-k),x3=x( p+q),x4=x(p-q),ω1=ω(n+k),ω2=ω(n-k),ω3 =ω(p+q),ω4=ω(p-q);H(n,k,p,q)表示式(4)中離散PWHT匹配函數(shù)的平方,即:
由于LFMCW信號與噪聲不相關(guān),且零均值循環(huán)高斯白噪聲的3階矩為0,則式(21)表示如下:
將x1~x4、ω1~ω4和式(22)代入上式,得到
將式(20)和(24)代入式(11),可以得到
將式(9)和(25)代入式(10),可以得到輸出信噪比
由上式可以看出:當SNRin?1時,SNRout≈N/2 ·SNRin,采樣點數(shù)N成為SNRin的倍增器,說明在大信噪比條件下,離散PWHT可以大幅提高信噪比,有利于信號的檢測;但是,當0<SNRin?1時,SNRout?SNRin,離散PWHT處理后反而使SNRin變得更小。
為了描述SNRout與SNRin的關(guān)系,設(shè)計信噪比增益因子G,如下所示:
當SNRout大于SNRin時,G>0;SNRout小于SNRin時,G<0;SNRout等于SNRin時,G=0,即無增益。為了實現(xiàn)輸入信噪比SNRin的改善,需要采樣點數(shù)N的合理設(shè)置,使得到G>0。
圖1描述了SNRin和采樣點數(shù)N共同對G的影響,并在其中標出0 dB平面(即G=0),在0 dB平面以上表示SNRout大于SNRin,反之表示SNRout小于SNRin。由圖可以看出:關(guān)系曲面只要滿足一定條件,明顯穿越0 dB平面,其穿越臨界點由信噪比SNRin和采樣點數(shù)N共同決定。
取SNRin等于-20 dB、-15 dB、-10 dB、-5 dB、0 dB的切面,可以描述采樣點數(shù)N與SNRout的關(guān)系,如圖2所示。由圖可以看出:SNRout(dB形式)隨采樣點數(shù)N增大以指數(shù)級別提高,在-20 dB條件下,當N大于500時,SNRout大于SNRin。這說明:在低信噪比條件下,只要采樣點數(shù)N足夠大,可以滿足SNRout大于SNRin的要求,離散PWHT可以適用于低信噪比下的LFMCW信號檢測。
在觀測時間Tobs確定的情況下,采樣點數(shù)N主要由采樣頻率fs決定,有以下關(guān)系:
低信噪比下的LFMCW信號檢測,根本上是可以通過提高采樣頻率fs而不斷改善的。由于Tobs確定,N和fs為線性關(guān)系,采樣頻率fs與SNRout的關(guān)系與圖2相似。理論上說,LFMCW信號的低信噪比檢測可以通過提高fs而實現(xiàn)。
由于fs與SNRout的關(guān)系,當fs固定時即N固定時,SNRin與SNRout有確定的關(guān)系。設(shè)N等于8、16、32、64、128,同時加入SNR′out=SNRin(即輸出信噪比增益等于0 dB)作為對比,得到SNRin與G的關(guān)系,如圖3所示。
由圖可以看出:當N固定,SNRin小于某一特定值時,G<0 dB,有SNRout<SNR′out=SNRin;除此之外,SNRout均大于SNRin。此臨界點在N>16時,其值為SNRout均大于SNRin。而在實際接收機中,其采樣點數(shù)N一般遠大于此數(shù)值,低于此數(shù)值導致SNRout?SNRin的情況是不存在的。所以,離散WHT在實際信號環(huán)境中總能保證SNRout>SNRin,可以應(yīng)用于實際信號環(huán)境中(包括低信噪比環(huán)境)的LFMCW信號檢測。
綜上所述,可以得到如下結(jié)論:由于PWHT所處理的LFMCW信號其持續(xù)時間(或調(diào)頻帶寬)是有限的,在采樣頻率fs足夠高的條件下,采樣點數(shù)N也會足夠大以適應(yīng)低信噪比和參數(shù)估計精度要求,在運算量提高的代價下,無論在大信噪比還是低信噪比條件下,都可以實現(xiàn)LFMCW信號的檢測和參數(shù)估計。
設(shè)觀測時間Tobs內(nèi)一個LFMCW信號的LFM調(diào)制脈沖數(shù)為M,一個LFM調(diào)制周期的時間長度為T,采樣頻率為fs,則觀測時間內(nèi)的總采樣點數(shù)為
基于PWHT的處理增益為
由于WHT是基于一個LFM調(diào)制周期進行處理的,沒有PWHT的累積過程,其處理增益為一個LFM調(diào)制周期處理的增益,在此時間內(nèi)的采樣點數(shù)為Tfs,結(jié)合文獻[4]給出的WHT的處理增益,可以得到
由式(31)和式(32)可以得到PWHT較WHT的增益
取一個LFM周期內(nèi)的采樣點數(shù)Tfs分別為16、
32、128和256,觀測時間內(nèi)LFMCW信號包含LFM周期個數(shù)M=1,2,3,4,5,輸入信噪比SNRin取-50~10 dB,對PWHT較WHT的增益進行了Matlab仿真,結(jié)果如圖4所示。由圖4中任意一張圖可以看出,在一個LFM周期內(nèi)的采樣點數(shù)Tfs固定的情況下,觀測時間內(nèi)LFMCW信號包含LFM周期個數(shù)M越多,PWHT較WHT的增益越高。在M=5的條件下,PWHT較WHT的處理增益高達大約14 dB,充分證明了PWHT較WHT對LFMCW信號有準相干積累增益,PWHT較WHT處理更適合處理LFMCW信號。
值得注意的是,在M=1時,PWHT較WHT的處理增益為0 dB,說明PWHT處理單脈沖LFM信號時,其效果與WHT是相同的。
周期Wigner-Hough變換可以實現(xiàn)對LFMCW信號的廣義似然比檢測和極大似然估計。將其用于工程應(yīng)用和仿真時,往往進行離散化數(shù)值計算,離散化的PWHT其性能直接影響工程應(yīng)用,需要研究PWHT的離散化算法對該變換的影響。
本文以離散PWHT的LFMCW檢測性能為對象,研究了離散化周期Wigner-Hough變換的信噪比增益以及較WHT的檢測性能優(yōu)勢。得出結(jié)論為:離散PWHT的信噪比處理增益具有準相干積累的效果,PWHT較WHT更優(yōu)的檢測性能,更適合在噪聲環(huán)境中對LFMCW信號的檢測。
[1]Pace P E.Detecting and Classifying Low Probability of Intercept Radar[M].Norwood,MA:Horizon House Artech,2004.
[2]Levanon N,Mozeson E.Radar Signals[M].Hoboken,New Jersey:A John Wiley&Sons,Inc.,2004.
[3]楊建宇,凌太兵,賀峻.LFMCW雷達運動目標檢測與距離速度去耦合[J].電子與信息學報,2004,26(2):169-173.
YANG Jian-yu,LING Tai-bin,HE Jun.LFMCW radar moving target detection and range velocity decoupling[J].Journal of Electronics&Information Technology,2004,26(2):169-173.(in Chinese)
[4]戴幻堯,蔣鴻宇.基于濾波器組和高階累積量技術(shù)的LPI信號特征檢測的新方法[J].系統(tǒng)工程與電子技術(shù),2009,31(6):1336-1340.
DAI Huan-yao,JIANG Hong-yu.A new method of LPI signal feature detection based on filter bank and higher order cumulants[J].System Engineering and Electronics,2009,31(6):1336-1340.(in Chinese)
[5]Barbarossa S.Analysis of multicomponent LFM signals by a combined Wigner-Hough transform[J].IEEE Transactions on Signal Processing,1995,43(6):1511-1515.
[6]劉鋒,孫大鵬,黃宇,等.基于改進Wigner-Hough變換的多分量LFM信號特征提取[J].北京理工大學學報,2008,28(10):914-917.
LIU Feng,SUN Da-peng,HUANG Yu,et al.Feature extraction of multi-component LFM signals based on improved Wigner-Houghtransform[J].Journal of Beijing Institute of Technology,2008,28(10):914-917.(in Chinese)
[7]Liu Feng,Xu Huifa,Sun Dapeng,et al.Feature extraction of symmetrical triangular LFMCW signal using Wigner-Hough transform[J].Journal of Beijing Institute of Technology,2009,18(4):478-483.
[8]Milne P R,Pace P E.Wigner distribution detection and analysis of FMCW and P4 polyphase LPI waveforms[C]//Proceedings of 2002 IEEE International Conference on A-coustics,Speech and Signal Processing.Orlando,USA:IEEE,2002:3944-3947.
[9]劉建成,王雪松,劉忠,等.基于Wigner-Hough變換的LFM信號檢測性能分析[J].電子學報,2007,35(6):1212-1217.
LIU Jian-cheng,WANG Xue-song,LIU Zhong,et al.Analysis of LFM Signal Detection based on Wigner-Hough Transfrom[J].Acta Electronica Sinica,2007,35(6):1212-1217.(in Chinese)
[10]Qi Lin,Tao Ran,Zhou Siyong,et al.Detection and parameter estimation of multicomponent LFM signal based on the fractional Fourier transform[J].Science in China Ser. F,2004,47(2):184-198.
[11]陳凡,齊林,朱亞瓊.多分量微弱LFM信號的檢測[J].南京郵電大學學報(自然科學版),2007,27(5):29-33.
CHEN Fan,QI Lin,ZHU Ya-qiong.Detection of Weak Multi-Component LFM Signal Based on the Fractional Fourier Transform[J].Journal of Nanjing University of Posts and Telecommunications(Natural Science),2007,27(5):29 -33.(in Chinese)
[12]李家強,金榮洪,耿軍平,等.基于高斯短時分數(shù)階傅立葉變換的多分量LFM信號檢測與參數(shù)估計[J].電子與信息學報,2007,29(3):570-573.
LI Jia-qiang,JIN Rong-hong,GENG Jun-ping,et al.Detection and parameter estimation of multi-component LFM signal based on Gaussian STFT[J].Journal of Electronics&Information Technology,2007,29(3):570-573.(in Chinese)
[13]Geroleo F G,Pearce M B.Detection and Estimation of Multi-Pulse LFMCW Radar Signals[C]//Proceedings of 2010 IEEE International Radar Conference.Washington DC,USA:IEEE,2010:1009-1013.
WANG Ze-zhong was born in Jinan,Shandong Province,in 1983.He received the M.S.degree in 2008.He is currently working toward the Ph.D.degree.His research concerns interception and identification of complex modulation radar signal.
Emulate:godohmy2005@tom.com
劉鋒(1960—),男,陜西寶雞人,2010年獲博士學位,現(xiàn)為教授、博士生導師,主要研究方向為電子信息戰(zhàn)理論及應(yīng)用;
LIU Feng was born in Baoji,Shaanxi Province,in 1960.He received the Ph.D.degree in 2010.He is now a professor and also the Ph.D.supervisor.His research concerns the theory and application of electronic information war.
黃宇(1983—),男,湖南永州人,2009年獲碩士學位,現(xiàn)為博士研究生,主要研究方向為復雜調(diào)制雷達信號截獲與識別;
HUANG Yu was born in Yongzhou,Hunan Province,in 1983. He received the M.S.degree in 2009.He is currently working toward the Ph.D.degree.His research concerns interception and identification of complex modulation radar signal.
向崇文(1986—),男,湖南桑植人,2008年獲學士學位,現(xiàn)為碩士研究生,主要研究方向為雷達偵察信號處理。
XIANG Chong-wen was born in Sangzhi,Hunan Province,in 1986. He received the B.S.degree in 2008.He is now a graduate student. His research concerns radar reconnaissance signal processing.
Digitized Periodic Wigner-Hough Transform and its Performance Analysis
WANG Ze-zhong1,LIU Feng2,HUANG Yu1,XIANG Chong-wen1
(1.Graduate Students′Brigade,Naval Aeronautical and Astronautical Institute,Yantai 264001,China;2.Department of Electronic Engineering,Naval Aeronautical and Astronautical Institute,Yantai 264001,China)
Linear Frequency Modulation Continuous Wave(LFMCW)detection performance of the digitized Periodic Wigner-Hough Transform(PWHT)algorithm is studied.Firstly,the weak signal detection performance is analysed through the derivation of the gain of the input signal-to-noise ratio(SNR)of a noisy LFMCW signal to the output SNR after PWHT processing.Then,the PWHT and Wigner-Hough Transform(WHT)detection performances are compared by obtaining the promotion of the gain of the input SNR to the output SNR after PWHT processing and WHT processing.Finally,the quasi-coherent integration effect of the PWHT algorithm in LFMCW detection is demonstrated through simulations,and the fact that PWHT performs better in the detection of LFMCW signal than WHT is also verified.
time frequency analysis;digitized periodic Wigner-Hough transform;LFMCW;detection performance analysis
The National Natural Science Foundation of China(No.60902054);China Postdoctoral Science Foundation(No.20090460114)
TN971.1
A
10.3969/j.issn.1001-893x.2012.09.008
王澤眾(1983—),男,山東濟南人,2008年獲碩士學位,現(xiàn)為博士研究生,主要研究方向為復雜調(diào)制雷達信號截獲與識別;
1001-893X(2012)09-1452-07
2012-02-28;
2012-04-26
國家自然科學基金資助項目(60902054);中國博士后科學基金資助項目(20090460114)