利用自適應(yīng)子波變換提高對微弱運(yùn)動目標(biāo)的檢測性能

本文研究了長時間相參積累時微弱運(yùn)動目標(biāo)回波信號的特點，分析了常規(guī)檢測方法的局限性，針對多項式相位信號模型，結(jié)合復(fù)正交線性相位子波基函數(shù)的設(shè)計，提出了一種基于自適應(yīng)子波變換局域線性逼近的微弱運(yùn)動目標(biāo)檢測方法，該方法具有自適應(yīng)頻帶劃分，快速算法實現(xiàn)，計算量低，并對統(tǒng)計分布的雜散分量具有自動抑制的特點.理論分析和計算機(jī)仿真結(jié)果驗證了該方法的有效性.
　　關(guān)鍵詞:微弱信號檢測；相參積累；AWT檢測器

An Improved Coherent Integration Detector with Adaptive Wavelet Transform

WANG Jun,ZHANG Shou-hong,YANG Ke-hu
(Key Laboratory for Radar Signal Processing,Xidian University,Xi'an 710071,China)

　　Abstract:This paper proposes a novel long-term coherent integration method for detection of weak high order polynomial phase signal (PPS),such as weak echo signals from maneuvering targets in radar and sonar applications.The method exploits the time-frequency tiling property of adaptive wavelet transform (AWT),and is practically realzable.In this paper,the unitary and linear phase wavelet basis function together with one order local linear approaching procedure are presented to estimate the signal trace in time-frequency plane.Additionally,the Radon-Wigner transform (RWT) is exploited in subband outputs of AWT to achieve low computational complexity implementation.Analysis and simulation results confirm the efficacy of the method.
　　Key words:weak signal detection;long-term coherent integration;AWT detector

一、引　　言
　　在雷達(dá)、聲納等探測系統(tǒng)中對微弱運(yùn)動目標(biāo)的檢測一直是比較困難的，主要原因是這類目標(biāo)的回波強(qiáng)度小且多普勒頻率變化復(fù)雜.為了檢測強(qiáng)背景中的運(yùn)動目標(biāo)，除了常規(guī)雜波抑制、抗干擾和降低系統(tǒng)噪聲等措施外，一種比較有效的方法是利用相參積累技術(shù)來增強(qiáng)接收回波，即用時間換取能量.一般而言相參積累時間受目標(biāo)穿越波束和回波包絡(luò)移動等因素的影響，而積累的性能主要取決于目標(biāo)回波的相干性.通常，若目標(biāo)作機(jī)動飛行，則隨波束駐留時間的增加，則目標(biāo)回波相干性變差，此時傳統(tǒng)的譜分析方法是失效的.為此人們提出了基于線性調(diào)頻回波模型的各種檢測方法［1，2］，并建立了快速算法［3］.實際上，在相對較長的一段觀測時間內(nèi)由于目標(biāo)的復(fù)雜運(yùn)動而產(chǎn)生的回波信號已不能用線性調(diào)頻模型去近似.因此，有必要研究在這種情況下的信號檢測問題.本文基于推廣的多項式相位回波模型，提出了利用自適應(yīng)子波變換作線性逼近對目標(biāo)回波信號進(jìn)行長時間相參積累的檢測方法.由于回波信號的頻率變化呈現(xiàn)非線性特征，不同的時間上多普勒變化特征不同，利用自適應(yīng)子波劃分，可以將信號全局多普勒變化轉(zhuǎn)化為局域線性逼近，從而以區(qū)間時變的多普勒相參技術(shù)取代常規(guī)線性調(diào)頻模型下的RWT技術(shù).另外AWT方法存在快速變換，其分岔樹結(jié)構(gòu)易于實現(xiàn)，對回波的局域優(yōu)化分割還有助于包絡(luò)補(bǔ)償?shù)膶崿F(xiàn)，同時，多尺度多分辨性能使得AWT保持著對統(tǒng)計分布型雜散分量的自動抑制能力，結(jié)合子帶信號能量的RWT技術(shù)，可以在不增加太多運(yùn)算量的同時，得到對微弱運(yùn)動目標(biāo)良好的相參積累檢測性能.

二、長相參積累時目標(biāo)回波的相干性
　　為了檢測隱身等微弱運(yùn)動目標(biāo)信號，可以采用稀布陣綜合脈沖孔徑雷達(dá)的波束駐留工作方式或讓常規(guī)雷達(dá)工作在“燒穿”方式下，以增加回波脈沖數(shù).已知運(yùn)動目標(biāo)回波信號的多普勒頻率為

fd=(2V/λ)cosφ(t)　(1)

其中λ為雷達(dá)波長，V為目標(biāo)的運(yùn)動速度，φ(t)為目標(biāo)運(yùn)動方向與雷達(dá)視線的夾角.常規(guī)波束掃描雷達(dá)在一次掃描中接收到的脈沖數(shù)由于受天線掃描的影響而較少，對這些回波，可以認(rèn)為其fd是不變的.但是在固定波束照射下，由于觀測時間的增加，即使對勻速直線運(yùn)動的目標(biāo)，由于φ(t)的變化而產(chǎn)生的多普勒頻率的變化已不能忽略，其大小為

Δfd=(4V/λ)sin(θb/2)sinφ0　(2)

式中θb為天線波束寬度，φ0為目標(biāo)勻速運(yùn)動方向與波束軸線的夾角.式(2)僅僅假設(shè)目標(biāo)為勻速運(yùn)動的情況，實際上，在數(shù)秒的相參時間內(nèi)，目標(biāo)的運(yùn)動情況可能包括加減速及機(jī)動拐彎等復(fù)雜運(yùn)動.若目標(biāo)作勻加速運(yùn)動(機(jī)動拐彎的情況類似)，即V(t)=V0+at,a為瞬時加速度，可以導(dǎo)出其fd的變化率為：

　(3)

其中R為目標(biāo)距離.由于在回波相參時間內(nèi)，φ(t)∈［0,θb］，且R很大，上式中忽略了(1/R)的高次項的影響.式(3)表明由于目標(biāo)相對雷達(dá)視角作勻加速機(jī)動而導(dǎo)致μf為高階時變的函數(shù).因此，為了分辨高階時變相位引入的非線性頻率調(diào)制，充分利用目標(biāo)相干性，則需要引入頻率局域分辨的思想，而利用傳統(tǒng)方法實現(xiàn)對非線性調(diào)頻信號和相參積累，實質(zhì)是一種全域分辨，得不到應(yīng)有的頻率分辨力.

三、復(fù)對稱正交子波基的設(shè)計
　　由于不同子波基函數(shù)的分辨性能不同，選擇有良好局域性能的子波基是重要的.為了有利于長時間雷達(dá)相參回波的積累，我們希望子波基函數(shù)具有線性相位、有限沖激響應(yīng)、正交分解的特點并存在快速變換實現(xiàn).
　　從Daubechies的FIR濾波器設(shè)計可知，其尺度函數(shù)滿足①緊支性，②正交性的約束關(guān)系.因此，只要能找到一個合適的尺度函數(shù)，在一定準(zhǔn)則的約束下，就能得到信號在多分辨空間的逼近，從而構(gòu)造出相應(yīng)的正交子波分解.Mallat給出了下面的關(guān)系，使得我們能夠利用濾波器組來產(chǎn)生相應(yīng)的尺度和子波函數(shù).
　　若設(shè)為由ck定義的離散濾波器的傳遞函數(shù)，則正交條件可表示為|H0(ω)|2+|H0(ω+π)|2=1,且|H0(0)|=1.若ck,k=1,…,N且|H0(ω)|≠0，則

　(4)

其中,H1(ω)=e-jω0(ω+π),“-”表示復(fù)共軛.并不是具有任意沖繳響應(yīng)的濾波器都可以用來產(chǎn)生子波，一般要求H0(ω)應(yīng)是一半帶低通濾波器，則H0(z)在z=-1處至少應(yīng)有一個零點存在.若H0(z)可以分解為下列形式：

　(5)

且滿足(1)F(0)=1,(2)B=sup|F(ejω)|,則Daubechies證明了若B＜2N-1，則由式(5)能夠生成緊支連續(xù)可導(dǎo)的尺度和子波函數(shù).
　　另外，為使子波函數(shù)具有較好的頻率分辨力，應(yīng)使其滿足一定的正則性要求，即

　(6)

Daubechies子波皆為滿足上述條件的實多項式的解，但不具備線性或廣義線性相位特性，這樣在處理中不可避免地帶來了新的相位分量，從而影響子帶信號的相參積累.
　　為了構(gòu)造線性相位或零相位復(fù)值正交子波變換，Lawton考慮了對稱性條件［4］，并由此得到了一定條件下(如N為偶數(shù))時的復(fù)值子波系數(shù).這里通過定義一種滿足式(5)，(6)以及Daubechies收斂條件的有效多項式，來獲得一般條件下復(fù)對稱正交子波函數(shù)的構(gòu)造.可以證明，常規(guī)的實對稱Haar子波，Daubiechies子波及Lawton復(fù)值子波函數(shù)，均是該多項式生成子波函數(shù)的特殊情況.定義：

P(z)=((1+z)/2)2N+2PN(z-1)　(7)

其中，且

　(8)

多項式PN(z)的2J個根給出了不同的解，由此可生成不同的子波基.通過對PN(z)作因式分解，可得：

H0(z)=((1+z)/2)1+J.q(z-1)　(9)

其中,xm和xm為PN(z)在單位圓內(nèi)的根，n為xn的復(fù)共軛，R為{1,2,…,N}的任一子集，U′和U的補(bǔ).利用式(7)～(9)生成的復(fù)值正交子波函數(shù)，具有最大平坦、有限沖激響應(yīng)、線性相位及快速分解特性.

四、利用AWT改善長時間相參積累性能
　　利用線性調(diào)頻模型的RWT變換實現(xiàn)積累，原理上看其時頻是平穩(wěn)相關(guān)，但由于目標(biāo)具有時變非平穩(wěn)特征，故需要進(jìn)行聯(lián)合時頻域分辨.利用上節(jié)的復(fù)正交線性相位子波函數(shù)，可以用AWT把對高階非線性多普勒變化用區(qū)間線性來近似.結(jié)合AWT的快速實現(xiàn)，這種逼近不需進(jìn)行多維搜索而能自適應(yīng)地分段匹配于信號.而且，不同子帶上頻率變化有相關(guān)性，無需作全程優(yōu)化，據(jù)此可進(jìn)一步降低算法的計算量，得到較好的檢測性能.
　　1.回波信號的自適應(yīng)子波變換
　　設(shè){x(n),n=1,2,…,N}是不同觀測期內(nèi)的回波信號，N為可能收到的總回波數(shù).對于自適應(yīng)子波的局域劃分，對每個分岔樹上的信號，可描述為：

　(10)

式中Wj,m(n)為第j個分解層上，m個子帶中第n個子波系數(shù)，Wj+1,2m(n)和Wj+1,2m+1(n)為Wj,m(n)的一個二叉樹，h(n)為酉尺度濾波器，g(n)為酉子波濾波器.快速離散自適應(yīng)子波分解可通過塔形算法實現(xiàn).復(fù)正交濾波器分解特性使得Wj,m(n)滿足能量守恒關(guān)系，即用任意的時頻瓦片劃分結(jié)構(gòu)覆蓋全頻率段，不同的瓦片間能量無冗余結(jié)構(gòu).這種分解有助于獲得優(yōu)化的局域匹配效果.
　　最佳子集的獲取與代價函數(shù)有關(guān)，通?？刹捎玫拇鷥r函數(shù)有多種，如門限法、最小lp范數(shù)法、最小方差和熵函數(shù)法等.其中最小熵意味著在一定的均方誤差條件下提取通帶信號能量的最集中表示.由于統(tǒng)計分布的雜散分量其熵函數(shù)較高，則最小熵準(zhǔn)則使得雜噪分量散布，利于子帶的局域相參積累和檢測處理.
　　2.AWT長時相干積累檢測方法
　　設(shè)回波信號模型由如下P階多項式相位信號表示，即x(t)=A.ej2πφ(t),其中,則檢測問題表示為：

Hm：r(t)=m.x(t)+n(t),m=0,1　(11)

這里r(t)為觀測信號，t∈［-Tc/2,Tc/2］.設(shè)n(t)為獨(dú)立隨機(jī)復(fù)高斯噪聲，Ω={P,ai,i=1,…,P}為信號參數(shù)空間.利用廣義最大似然估計可得到上述檢測問題的檢測統(tǒng)計量為

　(12)

式中和i是Ω的最大似然估計.式(12)實現(xiàn)的最優(yōu)解實際上是對回波非線性頻率變化的瞬時相位補(bǔ)償，因而能獲得檢測的性能上限.然而，最大似然估計是對Ω的(P+1)維聯(lián)合尋優(yōu)問題，具有很大運(yùn)算量，很難具體實現(xiàn)，而利用AWT多分辨，實質(zhì)上把最大似然估計對信號參數(shù)Ω的多維聯(lián)合搜索簡化為多個相關(guān)的一維搜索，同時保持檢測器性能基本不變，它是一種次最優(yōu)的實現(xiàn)方案.AWT長時間相干積累檢測器的檢測統(tǒng)計量為：

　(13)