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

本文研究了長時間相參積累時微弱運動目標(biāo)回波信號的特點，分析了常規(guī)檢測方法的局限性，針對多項式相位信號模型，結(jié)合復(fù)正交線性相位子波基函數(shù)的設(shè)計，提出了一種基于自適應(yīng)子波變換局域線性逼近的微弱運動目標(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)中對微弱運動目標(biāo)的檢測一直是比較困難的，主要原因是這類目標(biāo)的回波強(qiáng)度小且多普勒頻率變化復(fù)雜.為了檢測強(qiáng)背景中的運動目標(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ù)雜運動而產(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ù)，可以在不增加太多運算量的同時，得到對微弱運動目標(biāo)良好的相參積累檢測性能.

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

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

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

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

式中θb為天線波束寬度，φ0為目標(biāo)勻速運動方向與波束軸線的夾角.式(2)僅僅假設(shè)目標(biāo)為勻速運動的情況，實際上，在數(shù)秒的相參時間內(nèi)，目標(biāo)的運動情況可能包括加減速及機(jī)動拐彎等復(fù)雜運動.若目標(biāo)作勻加速運動(機(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}的任