基于自適應權重GPSR算法的ISAR成像

史潤佳，黃一飛，蔣忠進

（東南大學毫米波國家重點實驗室，江蘇南京 210096）

0 引 言

逆合成孔徑雷達（ISAR）成像在軍事和民用領域得到大量研究和應用，對ISAR 成像的精度與速度的要求也越來越高。壓縮感知理論自提出以來，在可壓縮信號處理方面體現出突出性能，得到了大量的關注和研究。在ISAR 成像中，雷達散射中心理論符合壓縮感知對于未知參數稀疏分布的要求，使得壓縮感知理論被成功應用于ISAR 成像，同時降低了對數據采樣率的要求［1-3］。

壓縮感知理論主要包括3個部分：信號的稀疏表示、測量矩陣以及重構算法。重構算法作為模型參數估計的重要步驟，一直是壓縮感知理論的研究重點。根據重構模型直接尋求l0范數最小化是一個NP 難問題，因此將l0范數的最優化問題轉化為l1范數正則化最小二乘問題［4］。為解決此問題，已經被提出的算法有基追蹤（BP）算法［5］、內點法［6］、稀疏重構（SpaRSA）算法［7］、閾值迭代（ISA）算法［8］、梯度投影稀疏重構（GPSR,Gradient Projection Sparse Reconstruction）算法［9］等。

為改進算法的性能，文獻［10］對SpaRSA 算法進行了改進，使用一種非單調策略來確定迭代步長以加快收斂速度。文獻［11］提出了迭代近似梯度投影（IAGP）算法，用近似的梯度來代替GPSR算法中的梯度模型，降低了計算量。文獻［12］在常規GPSR 算法的基礎上，通過Hessian 矩陣的近似得到迭代的二次近似模型，并結合延遲策略得出一種新的迭代步長，提高了算法的重構效率。

基于正則化參數在凸優化稀疏重構算法中的重要性［13-16］，本文提出一種自適應權重GPSR（AW GPSR）算法，并將其用于ISAR成像，此處的權重即為正則化參數。該快速算法中，在常規GPSR 算法的基礎上，給ISAR 圖像中每個散射點配以不同的正則化參數，并在迭代過程中自適應調整此正則化參數，以提升ISAR成像的迭代效率。

1 ISAR成像模型

在很短的相干累積時間內，ISAR 成像基本原理可等效為圖1所示的二維轉臺模型。由于目標尺寸遠遠小于目標到雷達的距離，所以在對回波數據完成運動補償及相位校正以后，目標上任一散射點p(a,b)到雷達的距離可近似為

圖1 ISAR成像二維轉臺模型

式中，a和b為散射點的坐標，R0為雷達到目標中心的距離，θ為目標轉動的角度。

根據理想散射中心模型，目標的譜域后向散射場可表示為

式中，f為入射頻率，xd為第d個散射點的散射幅度，d=1，2，…，D。其中exp(-j4πfR0/c)項不影響各個散射點之間的相位差，故在后續推導中將其忽略。

將雷達向目標發射信號的方位角離散為M個角點，第m個角點為θm，m=1,…,M。每個角點發射N個頻率步進的脈沖，第n個頻點為fn，n=1,…,N。由此，在第m個方位角的第n個頻點的散射場經推導可得

將二維成像場景離散采樣P行Q列，在壓縮感知里，P通常取M的倍數，Q通常取N的倍數。為使散射中心不發生越距離單元徙動，雷達的相對轉動角度一般很小，因此可近似認為sinθm≈θm，cosθm≈1。由于在ISAR 成像中，雷達脈沖帶寬相對于中心頻率非常小，則可以認為在頻率變化時，波長近似不變。基于以上條件，通過推導可得散射場表達式：

將上式改寫為矩陣形式可得

式中Y∈CM×N為散射回波信號的觀測矩陣，X∈CP×Q為散射幅度矩陣，也是需要估計的未知參數。Φr∈CM×P與Φa∈CN×Q分別為距離向和方位向的信息矩陣，其具體表達式為

由于強散射點數量遠小于散射幅度矩陣X的尺寸，所以矩陣X具有稀疏性，可以利用壓縮感知理論的重構算法，基于觀測信號Y予以重構。

為適應梯度投影稀疏重構（GPSR）算法的需要，將式（5）作矢量化處理，改為一維表達式：

式中y∈CMN×1、x∈CPQ×1分別為觀測矩陣Y、散射幅度矩陣X矢量化后的一維數組，Θ=Θa?ΘTr∈CMN×PQ為矢量化后的感知矩陣，符號?為Kronecker積。

2 GPSR參數重構算法

2.1 常規GPSR算法

在壓縮感知參數重構中，可基于l1范數最小化原則求解：

針對上式，可建立參數重構的目標函數為

τ為正則化參數。為了將其轉化成一個二次問題，將x分為正負兩部分，引入置換式：

其中(x)+=max{0,x}，ui=(xi)+，vi=(-xi)+，i=1,2,3,…,I，I=PQ是向量x的長度。然后x的l1范數可表示為

1I=[1,1,…,1]T為長度為I的單位列向量。

然后可將式（10）改寫為下邊界受限的二次方程：

將上式再次改寫為

式中，

常規GPSR 算法利用梯度投影算法來求解式（14）［17］。在第k次迭代中，將參數估計值zk沿著負梯度-?F(z(k))方向進行修正，修正公式如下：

式中，λ(k)∈[0,1]，α(k)需要確定。

然后定義函數g(jk)為

式中，j=1,2,…,J，J=2PQ也是矢量z的長度，且有g(k)={g(1k),g(2k),…,g(Jk)}。

確定第k次迭代中的初始迭代步長因子α(0k)：

常規GPSR算法步驟如下：

Step0：給定初始化的z(0)，選擇常數β∈(0,1),μ∈(0,1/2)，令k=0。

Step1：通過式（18）計算初始步長因子α(0k)。

Step2：通過回溯線搜索，選擇α(k)為序列α(0k),βα(0k),β2α(0k)…中第一個滿足下式的數字：

Step3：更新z值：z(k+1)=(z(k)-α(k)?F(z(k)))+。

Step4：測試z(k+1)是否滿足終止條件。如滿足條件則終止迭代；否則令k=k+1，并轉回Step1。

2.2 自適應權重GPSR算法

自適應權重GPSR（AW GPSR）算法中，對于x中不同幅值的元素賦予不同的正則化參數，以加快參數重構的收斂速度。其參數優化的目標函數為

式中，wI∈RI為權重向量，表示為

式中，wi是對應于參數x中的第i個元素的正則化參數。

得到參數優化的目標函數為

式中，wu為對應于u的權重矢量，wv為對應于v的權重矢量。將目標函數進一步寫為

式中，

令w=表示對應于z的權重矢量。

在AW GPSR算法里，也是利用梯度投影算法來求解未知參數z，其原理和常規GPSR算法相同。只是在迭代過程中，還存在一個加權系數w的更新問題。

為參數z加權以提高迭代效率的思想是，在迭代過程中，對于一個值很小的zi，可以為其添加一個較大的正則化參數wi，從而使它在后面的迭代中快速地收斂至0；對于一個值較大的zi，可以為其添加一個較小的正則化參數wi，以使其在下次的迭代中保持不變。

在該算法中，通過前一次迭代的z(k)和w(k)對正則化參數w進行更新，得到w(k+1)。設z(k)的支撐集為Γ(k)，也就是z(k)中非零元素的下標集合，則正則化參數w更新為

式中，η=，L=2PQ是待重構參數z的長度。

AW GPSR算法步驟如下：

Step0：給定初始化的z(0)，w(0)，選擇常數β∈(0,1),μ∈(0,1/2)，令k=0。

Step2：通過式（18）計算初始迭代步長因子α(0k)。

Step3：通過回溯線搜索，選擇α(k)為序列α(0k),βα(0k),β2α(0k)…中第一個滿足下式的數字：

Step4：更新z(k+1)=(z(k)-α(k)?F(z(k)))+；根據式（24），更新wi(k+1)，得w(k+1)。

Step5：測試z(k+1)是否滿足終止條件。如滿足條件則終止迭代；否則令k=k+1并轉回Step1。

3 實驗結果及分析

本文對AW GPSR 算法進行了Matlab 編程和ISAR 成像測試，使用仿真數據和實測數據進行ISAR 成像實驗，并與常規GPSR 算法、距離-多普勒（R-D）算法、ESPRIT 算法在計算時間與成像效果方面進行了比較。

仿真實驗使用的目標模型為Su27 戰機，其CAD 模型如圖2所示，電磁仿真的基本參數為：發射信號中心頻率為10 GHz，帶寬為0.5 GHz，步進頻采樣，采樣頻點為64；俯仰角固定為85°；方位角范圍為［-1.432 4°，1.432 4°］，采樣64 個角點；所以仿真得到的回波二維譜矩陣的尺寸為64×64。對應的徑向和橫向分辨率均為0.3 m。實驗環境為CPU 主頻3.20 Hz、64 位Windows10 操作系統、Matlab R2018a。

圖2 Su27飛機CAD模型

為證明AW GPSR 算法在目標函數收斂速度方面的效果，本實驗記錄并給出了AW GPSR 算法與常規GPSR 算法的目標函數迭代曲線對比圖，如圖3所示，其橫坐標為迭代次數，縱坐標代表目標函數值，從圖中可以看出，AW GPSR 算法的目標函數收斂速度優于常規的GPSR算法。

圖3 目標函數收斂趨勢比較

在對于仿真數據的成像效果方面，實驗比較了4種算法在信噪比分別為0，5，10 dB條件下進行ISAR成像的結果，如圖4所示。

圖4 基于仿真數據的ISAR成像結果

從圖中可以看出，R-D 算法和ESPRIT 算法成像的尺寸與回波二維譜矩陣的尺寸相同，都為64×64。而常規GPSR 算法與AW GPSR 算法屬于壓縮感知類算法，散射幅度矩陣的尺寸可以大于回波二維譜矩陣，本實驗設定的散射幅度矩陣X的尺寸為128×128。從ISAR 成像效果來看，AW GPSR算法要好于R-D算法和ESPRIT算法，與常規GPSR算法相當。

為更準確地比較AW GPSR 算法與常規GPSR算法的ISAR成像效果，本文引入了ISAR圖像的重構均方誤差指標，即以常規GPSR 算法在無噪聲條件下所成ISAR圖像為參考圖像，計算每幅ISAR圖像相對于參考圖像的均方誤差。表1給出了兩種算法在不同信噪比條件下的重構均方誤差，從表中可以看出，兩種算法在不同信噪比下的重構均方誤差都十分接近，證明了AW GPSR 算法與常規GPSR成像效果相當。

表1 兩種算法誤差對比

表2為4 種算法對于Su27 戰機仿真數據的ISAR 成像時間對比，其中R-D 算法成像時間最短，ESPRIT 算法成像時間最長；AW GPSR 算法成像時間要短于常規的GPSR 算法，在本次實驗中，AW GPSR 算法的成像時間相較于常規GPSR 算法大概減少了23%。

表2 4種算法成像時間對比

本實驗采用實測數據進行了ISAR 成像，實測數據來自于Yak-42 飛機。發射信號的中心頻率為5.52 GHz，帶寬為0.4 GHz，脈沖重復頻率為400 Hz，徑向分辨率為0.375 m。上述4 種算法對于實測數據的ISAR成像結果如圖5所示。

圖5 ISAR成像實測數據

由圖5可知，AW GPSR 算法的ISAR 成像效果與常規GPSR 算法相當，但明顯好于傳統的R-D 算法和ESPRIT算法。關于上述4種算法采用Yak-42飛機實測數據進行ISAR 成像的計算時間對比，其情況與仿真數據類似，此處不再列出。

4 結束語

本文在GPSR 算法的基礎上，研究了正則化參數對算法迭代效率的影響，并提出一種AW GPSR算法，該算法針對ISAR 圖像中的不同散射點使用不同正則化參數，以提高ISAR 成像迭代效率。本文采用仿真數據和實測數據，對比了AW GPSR 算法、常規GPSR 算法、傳統的R-D 算法以及ESPRIT算法的計算時間和ISAR 成像效果。相較于常規GPSR算法，AW GPSR算法的參數重構迭代收斂更快，因此具有更短的ISAR 成像計算時間，其計算時間縮短了約23%。在ISAR 成像效果方面，AW GPSR 算法和常規GPSR 算法相當，但明顯優于傳統的R-D算法和ESPRIT算法。