基于DMFT的LFM信號參數估計

摘要：線(xiàn)性調頻信號是低截獲概率雷達常用的一種信號形式，如何在低信噪比情況下檢測線(xiàn)性調頻信號一直是人們
研究的焦點(diǎn)之一。在離散匹配傅里葉變換的基礎上對算法進(jìn)行改進(jìn)，并利用改進(jìn)后的算法分別對單分量和多分量線(xiàn)性調頻信號進(jìn)行仿真，仿真結果表明離散匹配傅里葉變換能夠在低信噪比情況下比較準確地估計出線(xiàn)性調頻信號的參數，不存在交叉項問(wèn)題。離散匹配傅里葉變換是一種針對線(xiàn)性調頻信號有效的參數估計方法。
關(guān)鍵詞：離散匹配傅里葉變換；線(xiàn)性調頻；參數估計；低信噪比

線(xiàn)性調頻(LFM)信號是低截獲概率雷達常用的一種信號形式，對LFM信號檢測和參數估計一直備受人們的關(guān)注。針對該信號的處理方法有短時(shí)Fourier變換、Wigner―Ville變換、分數階Fourier、Hough―Wigner等，都存在分辨率不夠高，交叉項嚴重或者運算量太大的問(wèn)題。而匹配傅里葉變換是一種線(xiàn)性變換，不存在多分量信號交叉項的影響，能在低信噪比條件下檢測信號，而且分辨率很高，是一種針對線(xiàn)性調頻信號有效地進(jìn)行參數估計的方法。

l LFM信號形式
LFM信號的復數形式表示為：

式中，A(t)為信號包絡(luò )函數，f0為中心頻率，k0=B／T為調頻斜率，B為調頻帶寬，T為信號持續時(shí)間。
對于實(shí)際需要處理的信號，都是經(jīng)過(guò)采樣的離散信號。LFM信號的離散形式為：

式中，Ts為采樣時(shí)間間隔，如果信號持續時(shí)間為T(mén)，那么采樣點(diǎn)數N=T／Ts。

2 DMFT基本原理
LFM信號t(t)的匹配傅里葉變換有如下兩種形式：

稱(chēng)式(3)和式(4)分別為二階匹配傅里葉變換和二步匹配傅里葉變換，對應其離散形式為：

由式(5)計算得到的譜圖可稱(chēng)為離散二階匹配傅里葉變換譜，其中k不為零，它表示了不同基條件下的匹配傅里葉變換；由式(6)計算得到的譜圖可稱(chēng)為離散二步匹配傅里葉變換譜，它表示在不同頻率補償條件下信號的匹配傅里葉變換。
無(wú)論對離散二階匹配傅里葉變換譜還是離散二步匹配傅里葉變換譜，在對應于信號(f0，k0)的位置上，信號能量會(huì )發(fā)生聚集，在譜上表現為尖峰。在匹配傅里葉變換譜分布圖上進(jìn)行二維搜索，尖峰的坐標(f0，k0)即為該LFM信號的線(xiàn)性頻率f0和線(xiàn)性調頻斜率k0。
由于離散二階匹配傅里葉變換和離散二步匹配傅里葉變換具有不同的分辨率，通過(guò)文獻表明二步匹配傅里葉變換總是有比二階匹配傅里葉變換更高的分辨率，因此下面的分析都采用離散二步匹配傅里葉變換進(jìn)行LFM信號的檢測和參數估計。

3 算法改進(jìn)
對離散匹配傅里葉變換的二維搜索求極大值可以在低信噪比條件下獲得較高精度的信號參數。但是當信號帶寬增加，采樣頻率提高時(shí)，采樣點(diǎn)數增加，運算量增大。下面從減少運算量的角度進(jìn)行算法改進(jìn)。對離散之后的信號進(jìn)行離散匹配傅里葉變換，借助傅里葉變換的快速算法思想，實(shí)現離散匹配傅里葉變換的快速算法。對于長(cháng)度為N的線(xiàn)性調頻信號序列x(n)，其N(xiāo)點(diǎn)離散匹配傅里葉變換定義如下：

其實(shí)質(zhì)是將一個(gè)輸入一維時(shí)間序列x(n)變換為關(guān)于線(xiàn)性頻率和調頻斜率的二維序列Xc(f，k)，其中f為線(xiàn)性調頻信號的初始頻率，k為調頻斜率。從式(7)可以看出，對于每一個(gè)固定的調頻斜率k來(lái)說(shuō)，{Xc(f，k)}0≤f，k≤N-1是信號x(n)WknN2的DFT；當調頻斜率k=0時(shí)式(7)就轉變?yōu)镈FT。對式(7)進(jìn)行改進(jìn)得：

{Xc(f，k)}0≤f，k≤N-1的計算可以通過(guò)x(n)WknN2的快速傅里葉變換得到。在式(7)中需要N3次復數運算，經(jīng)過(guò)式(8)變換，運算量減小為N2／2log2N，提高了運算速度。為了提高該算法的估計精度，還可以在搜索范圍內多次估計，分為粗估計和精估計。即首先在搜索范圍內選擇大步長(cháng)，估計出信號參數，然后再在估計值鄰近的區域內改變搜索步長(cháng)重新估計，從而達到需要的精度要求。

4 仿真實(shí)驗
4．1 單分量LFM信號仿真
先對單分量LFM信號s(t)進(jìn)行參數估計，s(t)=exp[j2π(f0t+1／2k0t2)]，經(jīng)過(guò)下變頻的信號線(xiàn)性頻率f0=200 MHz，信號時(shí)寬T=5μs，以帶寬200 MHz的信號進(jìn)行仿真，比較在不同信噪比條件下信號參數估計的結果，如圖1所示。