金家偉，阮懷林，孫 兵

（國防科技大學電子對抗學院，合肥 230037）

0 引言

物體或物體結(jié)構(gòu)的微動可能會對返回的信號產(chǎn)生頻率調(diào)制，從而產(chǎn)生多普勒頻移的邊帶，稱為微多普勒效應。彈頭目標為了保持飛行的穩(wěn)定性，不僅會繞自身對稱軸作自旋運動，還會受到?jīng)_擊力矩的作用，力矩消失后對稱軸將在平衡位置作圓錐運動，即進動。進動是目標保持穩(wěn)定的一種特殊微動，而誘餌一般不具備姿態(tài)控制裝置，因此，進動是中段彈頭特有的微動方式。通過提取具有進動錐體目標的幾何特征，可以區(qū)分它是彈頭還是誘餌。

至今為止，由于時頻分析可將一維信號映射到二維時頻平面，從而展現(xiàn)信號的時頻信息，因此時頻分析技術被廣泛應用于分析時變的微多普勒頻率特征，因為時頻分析可以通過將一維信號映射到二維時頻平面來同時產(chǎn)生信號的時頻信息。文獻［6］結(jié)合Gabor 時頻分布和變分模態(tài)分解，來估計目標的自旋頻率和錐旋頻率；文獻［7］利用短時分數(shù)階傅立葉變換分離肢體和軀干的微多普勒信號；文獻［8］采用基于加權迭代自適應的時頻分析方法結(jié)合逆Radon 變換，分離重構(gòu)不同散射點的微多普勒分量。然而，在實際應用中，大多數(shù)現(xiàn)有的時頻方法在提取微多普勒信息方面仍然存在一些不可避免的不足，包括邊界失真和干擾項。另外，當目標含有多個強散射中心時，微動目標返回的雷達信號是多分量的。因此，所需的時頻分析方法應該同時抑制交叉項和自動項。

近年來，一種新的分析非線性、非平穩(wěn)數(shù)據(jù)的方法——希爾伯特- 黃變換（hilbert-huang transform，HHT）越來越受到人們的重視。該方法的關鍵部分是經(jīng)驗模態(tài)分解（empirical mode decomposition，EMD），其主要思想是將需要進行處理的非線性非平穩(wěn)信號分解，得到一組頻率分量各有不同的本征模態(tài)函數(shù)（intrinsic mode functions，IMF）。通過希爾伯特變換（hilbert transform，HT），IMF 能夠產(chǎn)生隨時間變化的瞬時頻率，其最終表示形式是能量-頻率- 時間分布，稱為希爾伯特- 黃譜（hilberthuang spectrum，HHS）。希爾伯特-黃變換的優(yōu)點之一是希爾伯特- 黃譜不涉及頻率分辨率和時間分辨率的概念，而是瞬時頻率。因其良好的適應性，已被廣泛應用于微多普勒效應的研究，文獻［12］使用了希爾伯特-黃變換對振動目標進行雷達微多普勒特征分析；文獻［13］提出了一種基于經(jīng)驗模態(tài)分解的圓錐目標雷達微動特征提取方法。然而，在實際應用中，EMD 算法并不完美，同樣存在一些問題，如模態(tài)混疊、端點效應等。

針對基于希爾伯特- 黃變換的進動目標微多普勒特征提取問題，由于雷達信號不可避免地含有較高水平的隨機噪聲，其瞬時頻率和瞬時幅度非常分散，常規(guī)希爾伯特-黃變換難以將目標信號與噪聲信號分離。本文采用EMD 的衍生算法CEEMDAN先對雷達回波進行分解，然后再進行希爾伯特譜分析，以此提取進動目標的微多普勒周期。

1 目標模型

空間進動錐體目標模型如圖1 所示，坐標系為雷達坐標系（U，V，W），雷達靜止于原點Q。O 為目標質(zhì)心，以O 為原點、目標對稱軸Oz 為z 軸建立目標本體坐標系O-xyz。以O 為原點建立參考坐標系O-XYZ，以初始時刻與目標對稱軸Oz、進動軸OZ共面且垂直于OZ 的方向為Y 軸，X 軸根據(jù)右手準則確定。目標在平動的同時，以角速度ω繞對稱軸z軸作自旋運動，同時以角速度ω繞OZ 軸錐旋運動（ω和ω均采用參考坐標系中的表達式），自旋軸和錐旋軸之間的夾角為進動角θ。

圖1 目標模型

在光學區(qū)，雷達目標的整體散射特性通?？梢缘刃槿舾蓚€散射中心的疊加。不失一般性，假設目標的微動是周期性進動，目標等效為K 個散射中心，各散射中心各向同性，雷達發(fā)射的電磁波為連續(xù)單頻波，載頻為f。雷達向由i 個散射點組成的目標發(fā)射電磁波，不考慮目標平動帶來的影響，返回的基帶信號可表示為

2 改進的希爾伯特-黃變換

2.1 希爾伯特-黃變換

2.2 完備總體經(jīng)驗模態(tài)分解

Torres 等提出了一種對噪聲有自適應分解特性的完備總體經(jīng)驗模態(tài)分解（CEEMDAN），在每一階段的分解中都加入一個特定的白噪聲，再通過計算殘差得到該階段的模態(tài)。CEEMDAN 是一種噪聲輔助數(shù)據(jù)分析方法，所加入的白噪聲在整個時頻空間分布相對均勻，以提供時頻空間中濾波的參考尺度，而不同尺度的信號自動投影到由白噪聲建立的合適的參考尺度上。雖然每次試驗產(chǎn)生的結(jié)果都會包含噪聲，但是當試驗的數(shù)量足夠時，其結(jié)果的總體平均當中的噪聲將消失，即總體平均當中所添加的白噪聲往往會相互抵消，從而只包含所分解得到的信號。

3 仿真與分析

為驗證希爾伯特- 黃變換在進動目標微多普勒特征提取方面的效果，本文分別以EMD 和CEEMDAN 來進行信號分解，然后進行希爾伯特譜分析，并對比兩種方法的提取效果。

仿真條件：假設雷達工作在10 GHz，且雷達坐標系中本體坐標系原點O 的坐標為（400，500，100）km，本地坐標系和參考坐標系之間的初始歐拉角（x-y-z 序列）為（30°，60°，45°）。假設目標繞z 軸旋轉(zhuǎn)，目標上有兩個散射點：第1 散射點A 位于錐頂，在本體坐標系中的坐標為（0，0，1）m；第2 個散射點B 位于錐底的尾翼，在本體坐標系中的坐標為（0.5，0，-0.5）m。雷達照射時間為8 s。旋轉(zhuǎn)頻率為f=1 Hz，圓錐運動頻率為f=0.5 Hz。假設環(huán)境噪聲為高斯白噪聲。

當不存在噪聲時，基于EMD 和CEEMDAN 的進動周期提取結(jié)果分別如圖2 和圖3 所示，圖2（a）和圖3（a）分別是兩種算法獲取的回波信號希爾伯特譜，圖2（b）和圖3（b）分別是兩種算法獲取的IMF前3 個分量的希爾伯特譜。

圖2 基于EMD 的估計結(jié)果（無噪聲）

希爾伯特譜是時間、頻率和能量的三維譜，當沒有噪聲時，圖2（a）和圖3（a）都展現(xiàn)了回波信號時頻分布的細節(jié)，其中的頻率都是取模后的頻率，也就是正頻率。微多普勒信號經(jīng)過EMD 后都被分解在第1 個IMF 分量中，如圖2（b）所示；而微多普勒信號經(jīng)過CEEMDAN 后被分解在各個IMF 分量中，如圖3（b）所示。由圖2 和圖3 可看出，兩種算法在不存在噪聲時，都能夠準確估計出進動目標的微多普勒周期為2 s。

圖3 基于CEEMDAN 的估計結(jié)果（無噪聲）

當雷達回波信號中混雜著噪聲時，同樣的仿真條件下，假設信噪比SNR=10 dB，圖4（a）為基于EMD 的回波信號希爾伯特譜，圖4（b）是基于EMD得到的IMF 前3 個分量的希爾伯特譜。

圖4 基于EMD 的估計結(jié)果（SNR=10 dB）

從圖4 可以看出，隨著信噪比的降低，信號的希爾伯特譜中的微多普勒頻率逐漸淹沒在噪聲中，無法用以提取進動目標微多普勒周期，如圖4（a）所示。由于EMD 未能將微多普勒信號和噪聲分離開，而是都混在第1 個IMF 分量中，如圖4（b）所示，隨著信噪比SNR 的降低，將同樣無法通過第1 個IMF分量的希爾伯特譜來提取目標信息。但在其他條件不變的情況下，改用CEEMDAN 對雷達回波進行處理，然后進行希爾伯特譜分析。圖5（a）為回波信號希爾伯特譜，圖5（b）為采用CEEMDAN 所獲取的IMF 前3 個分量的希爾伯特譜。

圖5 基于CEEMDAN 的估計結(jié)果（SNR=10 dB）

從圖5（a）中可以看出，當信噪比持續(xù)降低時，與圖4（a）相似，其中的微多普勒信號伴隨著大量的噪聲，難以直接通過回波信號地希爾伯特譜來提取微多普勒信息。但通過CEEMAN 分解后所得到的一組IMF，噪聲頻率主要集中于第1 個IMF 分量中，雖然第1 個IMF 分量中也包含了一部分的微多普勒頻率信號，但更多的微多普勒信號被分解到了第2 個和第3 個分量中，從而達到將噪聲和目標信號分離的目的，如圖5（b）所示，則從該圖中的第2 個和第3 個IMF 分量的希爾伯特譜中，可以準確提取出目標的微多普勒周期。

進一步降低信噪比進行對比分析，圖6（a）和圖6（b）分別給出了信噪比分別為5 dB 和1 dB 時，基于EMD 得到的IMF 前3 個分量的希爾伯特譜。

圖6 EMD 獲取的IMF 前3 個分量的希爾伯特譜

由于EMD 算法難以從噪聲分離出微多普勒信號，在其IMF 分量的希爾伯特譜中，即使信噪比為5 dB，仍然無法估計出進動周期，如圖6（a）所示。但采用CEEMAN 時，圖7（a）和圖7（b）分別給出了信噪比分別為5 dB 和1 dB 時，基于CEEMDAN 得到的IMF 前3 個分量的希爾伯特譜。

圖7 CEEMDAN 獲取的IMF 前3 個分量的希爾伯特譜

將圖7 與圖6 對比可知，隨著信噪比的進一步降低，即使信噪比低至1 dB，仍能通過其IMF 分量的希爾伯特譜提取出進動周期，如圖7（b）所示。這也進一步驗證了在低信噪比環(huán)境中，以CEEMDAN改進的希爾伯特-黃變換具備更加良好的性能。仿真中假設環(huán)境噪聲為高斯白噪聲，旨在驗證CEEMDAN 具備比EMD 更好的抗噪性能，但在實際工程中的噪聲不一定是高斯白噪聲，算法具體的抗噪性能需要以實際雷達回波數(shù)據(jù)進行驗證，這是本文算法存在的局限性。

4 結(jié)論

本文將CEEMDAN 算法替換EMD 算法，以改進希爾伯特-黃變換，用以提取進動目標的微多普勒周期。仿真實驗表明，改進的希爾伯特-黃變換能夠獲取進動目標微多普勒信號的希爾伯特譜，從而提取出目標的進動周期；相對于EMD 算法，CEEMDAN 算法能夠?qū)⑽⒍嗥绽招盘柡驮肼曅盘栠M行分離，從而有效克服了常規(guī)希爾伯特-黃變換在噪聲環(huán)境中性能較差的不足。