,,

(大連大學信息工程學院，遼寧 大連 116622)

0 引言

逆合成孔徑雷達(Inverse Synthetic Aperture Radar, ISAR)具有全天時遠距離目標成像的優(yōu)勢，已被廣泛應(yīng)用于軍事及民用等領(lǐng)域[1]。傳統(tǒng)的FFT成像方法，成像分辨率受發(fā)射信號帶寬和目標轉(zhuǎn)角的限制[2]，而發(fā)射大帶寬信號以及增加相干處理時間會產(chǎn)生龐大的數(shù)據(jù)，對雷達系統(tǒng)的A/D轉(zhuǎn)換速率、數(shù)據(jù)傳輸及存儲等有較高的要求[3]。因此，有必要研究一種以較少的觀測樣本實現(xiàn)高精度ISAR成像的算法。

近年來發(fā)展迅速的壓縮感知(Compressive Sensing, CS)理論[4]指出，當信號或信號在某個變換域是稀疏的，就可通過求解一個最優(yōu)化問題，以少量的信號測量值通過稀疏恢復算法高概率的恢復出原始信號。而對于ISAR成像而言，目標回波由目標上所有散射點回波疊加而成[5]，其中強散射點起主導作用，其數(shù)目遠小于回波采樣點數(shù)，滿足稀疏信號的特征[6]，因此可以將CS應(yīng)用于ISAR成像領(lǐng)域，緩解傳統(tǒng)成像方法帶來的超大數(shù)據(jù)采集、存儲和傳輸壓力。

將CS理論應(yīng)用于ISAR成像[7]的關(guān)鍵在于根據(jù)具體情況選擇或設(shè)計一種快速有效的重構(gòu)算法實現(xiàn)對原始信號的精確重構(gòu)。近年來眾多的綜述類文章都對稀疏重構(gòu)算法進行了歸納和分類[8]。凸優(yōu)化和貪婪追蹤算法是目前被廣泛使用的兩類算法。凸優(yōu)化算法重構(gòu)精度高，但計算比較復雜；相比之下，貪婪類算法在重構(gòu)精度和計算量之間取得了較好的平衡，迅速成為研究的熱點之一[9]。貪婪類算法的代表性算法為正交匹配追蹤算法(OMP)[10]。目前大多數(shù)的貪婪類算法都是對OMP算法加以改進，以提高算法的重構(gòu)精度與計算效率。例如StOMP算法是對原子識別步驟的改進[11]，而SP[12]和CoSaMp[13]等算法則是在原子識別步驟改進的基礎(chǔ)上，增加了冗余原子修剪的步驟，來遞歸的進行重構(gòu)支撐集選擇。對于雷達信號而言，由于不同目標散射點的散射強度不同，故信號幅度往往是隨機起伏的，對于這一類信號，若直接使用上述的OMP或SP等算法進行重構(gòu)，則會產(chǎn)生兩類重構(gòu)錯誤問題：一類是幅度相對小的信號分量不能被正確重構(gòu)；第二類是幅度相對較大的信號附近會產(chǎn)生“虛假”重構(gòu)信號[14]。產(chǎn)生這種現(xiàn)象的原因在于感知階段無法選出正確的原子，或是選出的原子不能很好地表達原信號，最終導致重構(gòu)失敗[15]?；诖?，本文提出基于改進的正則化子空間追蹤算法(MRSP)的ISAR成像方法。

1 壓縮感知ISAR成像模型

ISAR成像實質(zhì)上是用成像算法將回波中包含的散射點的位置和強度以圖像的方式表示。若雷達與目標間平動分量已經(jīng)被補償?shù)?，則ISAR目標成像等效為轉(zhuǎn)臺模型[16]。假設(shè)雷達發(fā)射線性調(diào)頻信號：

(1)

其中，tr是快時間，fc是載頻，γ是調(diào)頻率，Tp是脈沖寬度。距離脈沖壓縮后，忽略二次項與常數(shù)相位項，則目標上某散射點p(xp,yp)的回波信號表達式近似為：

(2)

其中，A表示回波信號幅度，Ta為觀測時間，λ為波長，c為光速。假設(shè)觀測時間內(nèi)目標轉(zhuǎn)角較小，則R(tm)≈R0+yp+xpωtm，R0為雷達與目標基準點的距離。根據(jù)式(2)所示脈沖壓縮后的信號，假設(shè)某距離單元包含K個較強的散射點，并考慮加性噪聲，則此單元的回波信號表示為：

(3)

其中，Bk和fk分別表示第k個散射點的散射系數(shù)和多普勒頻率。

ISAR數(shù)據(jù)中的強散射點數(shù)相對于全部數(shù)據(jù)來說具有稀疏特性，因此原始數(shù)據(jù)圖像可以通過CS技術(shù)進行感知重構(gòu)[17]。定義時間向量t=[1∶N]TΔt，N=Ta/Δt為脈沖數(shù)量，Δt=1/fr表示時間分辨率，fr為脈沖重復頻率。定義多普勒頻域的頻率分辨率為Δf，因此，離散的多普勒向量可以定義為fd=1：NTΔf。因此，構(gòu)建一個傅里葉基為：

此時，將式(3)某個距離單元的回波信號寫成矩陣形式為：

s=Ψα+n

(4)

其中，α為信號s在Ψ域的稀疏表示向量，其非零分量對應(yīng)于K個最大散射點的復幅度。

在基于壓縮感知的ISAR數(shù)據(jù)成像中，測量矩陣Φ∈RM×N(M

y=ΦΨα+n′=Θα+n′

(5)

考慮加性噪聲，ISAR信號重構(gòu)優(yōu)化問題可以寫為：

(6)

式中，ε為噪聲門限，Θ=φΨ為感知矩陣。如何設(shè)計有效的重構(gòu)算法重構(gòu)式(6)，避免重構(gòu)錯誤問題，提高重構(gòu)效率是本文的主要目的。

2 基于改進的正則化子空間追蹤算法的ISAR成像

2.1 改進的正則化子空間追蹤算法(MRSP)

OMP算法作為最早的貪婪算法之一，其對于之后的各種改良算法和新提出的算法都具有重要的指導意義。從現(xiàn)有的研究結(jié)果來看，對OMP計算量的改善通常以損失些許重構(gòu)性能為代價，而性能的改善則以較高的運算量或存儲量為成本[18]。對比之下，SP算法具有較好的計算速度和重構(gòu)概率，但是SP算法在每次回溯時，強信號分量總是得到保留，以至于在后續(xù)的回溯中，強信號分量遮擋弱信號分量，出現(xiàn)重構(gòu)錯誤現(xiàn)象[14]?？紤]實際應(yīng)用，本文提出MRSP算法，采用黃金比例分割法(0.618法)[19]確定正則化過程實現(xiàn)支撐集的二次篩選，以及回溯反饋的思路來更加準確地估計、擴充目標信號的真實支撐域，最終重構(gòu)出稀疏信號，使得壓縮感知重構(gòu)算法更具實用性。

MRSP算法在每次迭代過程中，首先計算殘差r，然后求殘差與感知矩陣Θ中各個原子(本文所提“原子”指組成感知矩陣Θ的各個列向量)之間的相關(guān)系數(shù)u，從u中挑選出絕對值最大的K(K為信號的稀疏度)個原子的索引值存入索引集J，然后采用正則化過程繼續(xù)篩選原子，根據(jù)式0.618|u(i)|≤|u(j)|，i，j∈J，將J中索引值所對應(yīng)的原子的相關(guān)系數(shù)分成若干組，選取能量最大的一組相關(guān)系數(shù)對應(yīng)的原子的索引值存入J0中，完成正則化過程；然后將該索引集J0并入當前支撐集索引Λ中，更新支撐集T并利用最小二乘法進行信號估計，再根據(jù)回溯思想，選取前m(m

值得注意的是，現(xiàn)有的重構(gòu)算法中涉及到正則化優(yōu)選原子的文獻都是根據(jù)式0.5|u(i)|≤|u(j)|,i,j∈J進行原子篩選，而本文所提算法采用黃金比例分割法確定正則化過程，根據(jù)式0.618|u(i)|≤|u(j)|,i,j∈J進行正則化識別候選集中能量最大的部分原子。黃金比例分割法是一種經(jīng)典的優(yōu)化算法，能夠以盡可能少的試驗次數(shù)盡快找到實驗中的最優(yōu)方案，故本文采用此種方法來確定正則化過程能更加精確快速的優(yōu)選出真實原子，提高支撐集的可靠性。此外，本文所提算法在回溯反饋時，只選取前m(m

輸入：測量向量y，感知矩陣Θ，稀疏度K；

1)初始化：迭代次數(shù)t=1；殘差ro=y；索引值集合J=?，Λ=?；支撐集T=?；

2)計算相關(guān)系數(shù)u=abs[ΘTrt-1](即計算〈rt-1，θj〉，1≤j≤N)，并將u中絕對值最大的K個值所對應(yīng)的索引值存入J；

3)正則化：根據(jù)式0.618|u(i)|≤|u(j)|，i，j∈J，將J中索引值所對應(yīng)的原子的相關(guān)系數(shù)分成若干組，找出能量最大的子集J0；

4)更新Λt=Λt-1∪J0，Tt=Tt-1∪θj(j∈J0)；

6)更新殘差

7)t=t+1，如果t≤K則返回第(2)步繼續(xù)迭代；如果t>K或‖rt-rt-1‖2≤ε則停止迭代；

2.2 基于MRSP算法的CS-ISAR成像

基于本文所提MRSP算法的CS-ISAR成像算法步驟可以描述如下：

1)對雷達回波信號進行脈沖壓縮得到頻域上的稀疏信號；

2)設(shè)計測量矩陣Θ。為了能精確重建出原始信號，測量矩陣Θ與稀疏基Ψ需要滿足不相關(guān)的要求，故測量矩陣Θ采用隨機高斯矩陣，矩陣滿足均值為0，方差為1的正態(tài)分布。

3)稀疏采樣。測量矩陣在方位向總的N個脈沖中隨機抽取M個脈沖，即N維的回波信號經(jīng)測量矩陣Θ降維處理得到M維的觀測信號y，這個低維的觀測信號包含了重構(gòu)信號的絕大部分信息。

4)信號重構(gòu)。采用重構(gòu)算法MRSP對觀測信號y進行重構(gòu)。選定閾值ε，當殘差‖rt-rt-1‖2≤ε時停止迭代，否則繼續(xù)迭代過程。

5)雷達圖像后處理。將重構(gòu)信號進行適當組合即可得到雷達圖像。

3 仿真及實測數(shù)據(jù)實驗

3.1 MRSP算法重構(gòu)性能驗證

為驗證該算法的有效性與可行性，分別選取兩種信號對該算法進行仿真驗證。定義壓縮率γ=M/B，M為觀測數(shù)，N為原始信號長度。

1)隨機信號壓縮重構(gòu)

隨機產(chǎn)生一個長度N=256，稀疏度K=12的隨機信號。采用高斯隨機矩陣作為測量矩陣，利用MRSP算法對該隨機信號進行重構(gòu)。實驗中選擇壓縮率γ=0.25，即測量點數(shù)M=64。重構(gòu)信號與原信號的對比如圖1所示，重構(gòu)誤差如圖2所示。

2)線性調(diào)頻信號壓縮重構(gòu)

選擇線性調(diào)頻信號s(t)=exp(1.5×1012jπt2)，信號長度取N=1 024。采用離散傅里葉變換構(gòu)成的正交基作為稀疏基，采用高斯隨機矩陣作為測量矩陣，利用MRSP算法對調(diào)頻信號進行重構(gòu)。實驗中選擇壓縮率γ=0.5，即測量點數(shù)M=512，重構(gòu)信號與原信號的對比以及重構(gòu)誤差如圖3所示。

從以上仿真可以看出使用MRSP算法重建的信號能夠較好地反映原始信號，重構(gòu)誤差較小，故MRSP算法重構(gòu)性能較好，是一種有效可行的重構(gòu)算法。

3.2 重建性能比較

本部分利用仿真實驗驗證OMP、SP、MRSP三種算法的信號重建性能。仿真實驗選取一個長度為256的稀疏信號，稀疏度為 12，即包含 12個非零隨機數(shù)，其他元素全部為0。測量矩陣采用高斯隨機矩陣。仿真過程中，當殘差小于1×10-6時，表示信號重建成功，停止迭代。

為驗證測量數(shù)M變化對三種算法重建性能的影響，設(shè)定量測數(shù)從24以等間隔5增加到 256，每個測量數(shù)下進行500次蒙特卡洛實驗，三種算法重建成功概率結(jié)果如圖4所示。

可見，當測量數(shù)較少時，三種算法的重建概率都較低；隨著測量數(shù)的增加，三種算法重建概率都逐漸增大。對比來講，隨著測量數(shù)的增多，MRSP算法重建成功概率要高于OMP、SP算法。

接下來驗證信號稀疏度K變化對三種算法重建性能的影響。設(shè)定稀疏度從1以等間隔5增加到70，每個稀疏度下進行500次蒙特卡洛實驗，重建成功概率結(jié)果如圖5所示。

由圖5可得，當稀疏度K>16時，OMP算法重建概率開始下降；當K>36時，SP、MRSP算法的重建概率開始下降，但SP算法重建概率下降速度較快，而MRSP算法下降速度較平緩。由此可知，隨著信號稀疏度的增加，MRSP算法性能更加穩(wěn)定、準確重建效果較OMP算法、SP算法更好。

3.3 實測數(shù)據(jù)成像

本節(jié)選取YaK-42飛機的一組實測數(shù)據(jù)驗證本文所提算法成像性能。其相關(guān)的雷達參數(shù)為：載頻5.52 GHz，脈寬25.6 μs，脈沖重復頻率400 Hz。實驗中選取256個脈沖序列作為全孔徑數(shù)據(jù)量。

為驗證本文所提算法的成像性能，對回波信號分別進行傳統(tǒng)R-D成像、OMP算法、SP算法以及本文所提MRSP算法成像。仿真時，傳統(tǒng)R-D算法利用256個全回波數(shù)據(jù)進行成像；基于壓縮感知的OMP算法、SP算法以及本文所提算法成像時，采用離散傅里葉變換構(gòu)成的正交基作為稀疏基，采用高斯隨機矩陣作為測量矩陣，選擇壓縮率γ=2/3，即測量回波脈沖個數(shù)為M=170。仿真結(jié)果如圖6所示。

圖6(a)是利用256個回波的R-D成像結(jié)果，圖6(b)—(d)是基于CS算法的成像結(jié)果。比較可見，利用CS成像算法，僅需170個回波脈沖即可得到相應(yīng)的成像效果，甚至更少的回波都能得到有效成像，說明壓縮感知方法是有效可行的，能夠利用較少的回波數(shù)據(jù)進行有效成像，在數(shù)據(jù)完整性、傳輸以及存儲空間等方面具有很高的優(yōu)越性。

比較圖6(b)—(d)成像結(jié)果，相同運行環(huán)境，相同回波脈沖個數(shù)條件下，MRSP算法較OMP算法、SP算法得到的成像結(jié)果聚焦性更好、更加清晰。對比這三種成像算法的峰值信噪比PSNR和運行時間：OMP算法的PSNR為5.474 3 dB，運行時間為12.246 1 s；SP算法的PSNR為4.407 8 dB，運行時間為10.420 9 s；MRSP算法的PSNR為5.662 9 dB，運行時間為7.238 4 s；與SP算法相比，MRSP算法成像結(jié)果峰值信噪比PSNR提高了約1.2 dB，運行時間提高約30%；與OMP算法相比，MRSP算法成像結(jié)果峰值信噪比PSNR提高了約0.2 dB，運行時間提高約40%。由此可以得出，本文提出的MRSP算法與OMP、SP算法相比，成像效果更好、運行時間更短，具有更好的優(yōu)越性。

4 結(jié)論

本文提出了基于改進的正則化子空間追蹤算法的ISAR成像方法。該方法采用壓縮感知理論，將SP算法與正則化思想結(jié)合，采用回溯反饋和貪婪精選的方法進行支撐集選擇，并利用黃金比例分割法確定正則化過程實現(xiàn)支撐集的二次篩選，提高了對稀疏信號支撐集的估計精度；減少了傳統(tǒng)ISAR成像中對大數(shù)據(jù)的要求，降低了雷達成像系統(tǒng)的復雜性，以較少的觀測信息獲得較高質(zhì)量的目標圖像，在一定程度上解決了大量數(shù)據(jù)傳輸及存儲的難題。仿真實驗結(jié)果表明該算法對較少的回波數(shù)據(jù)進行重構(gòu)，既能獲得更加清晰的成像結(jié)果，又能進一步提高運算速度，是一種有效可行的重構(gòu)算法。

參考文獻:

[1]江東, 童寧寧, 馮為可, 等. 基于欠定系統(tǒng)局灶解法的逆合成孔徑雷達成像[J]. 探測與控制學報, 2016, 38(1)：99-103.

[2]吳敏, 邢孟道, 張磊. 基于壓縮感知的二維聯(lián)合超分辨ISAR成像算法[J]. 電子與信息學報, 2014, 36(1)：187-193.

[3]李文靜, 陳紅衛(wèi). 一種基于壓縮感知的ISAR成像方法[J].計算機仿真, 2015, 32(8)：10-13.

[4]Donoho D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4)：1289-1306.

[5]黃培康, 殷紅成, 許小劍. 雷達目標特性[M]. 北京： 電子工業(yè)出版社, 2005.

[6]李民. 基于壓縮感知的雷達高分辨成像技術(shù)研究[D]. 哈爾濱：哈爾濱工業(yè)大學, 2014.

[7]HuChenxi, Liu Yimin, Li Gang, et al. Improved FOCUSS method for reconstruction of cluster structured sparse signals in radar imaging[J]. Science China Information Sciences, 2012, 55(8)：1776-1788.

[8]Qaisar S, Bilal R M, Iqbal W, et al. Compressive sensing： From theory to applications, a survey[J]. Journal of Communications & Networks, 2013,15：443-456.

[9]Chandran K S, Mishra A, Shirhatti V, et al. Comparison of matching pursuit algorithm with other signal processing techniques for computation of the time-frequency power spectrum of brain signals[J]. Journal of Neuroscience the Official Journal of the Society for Neuroscience, 2016, 36(12)：3399.

[10]Tropp J A, Gilbert A C. Signal recovery from random measurements via orthogonal matching pursuit[J]. IEEE Transactions on Information Theory, 2007, 53(12)：4655-4666.

[11]Donoho D L, Tsaig Y, Drori I, et al. Sparse solution of underdetermined systems of linear equations by stagewise orthogonal matching pursuit[J]. IEEE Transactions on Information Theory, 2012, 58(2)：1094-1121.

[12]Dai W,Milenkovic O. Subspace pursuit for compressive sensing signal reconstruction[J]. IEEE Transactions on Information Theory, 2009, 55(5)：2230-2249.

[13]Needell D,Tropp J A. CoSaMP： Iterative signal recovery from incomplete and inaccurate samples[J]. Applied & Computational Harmonic Analysis, 2008, 26(3)：301-321.

[14]丁函,王毅,袁磊,等.一種新的貪婪回溯子空間追蹤算法研究[J/OL]. [2016-10-12]. http：// www. arocmag.com/article/02-2017-10-036.html.

[15]梁丹亞, 李宏偉, 李琦,等. 一種改進的壓縮感知重構(gòu)算法[J]. 計算機仿真, 2014, 31(9)：310-313.

[16]李少東,楊軍,馬曉巖.基于壓縮感知的ISAR高分辨成像算法[J].通信學報, 2013(9)：150-157.

[17]Zhang Lei, Xing Mengdao, Qiu Chengwei, et al. Achieving higher resolution ISAR imaging with limited pulses via compressed sampling[J]. IEEE Geoscience & Remote Sensing Letters, 2009, 6(3)：567-571.

[18]Majumdar A, Krishnan N, Pillai S R B, et al. Extensions to orthogonal matching pursuit for compressed sensing[C]//2011 National Conference on Communications (NCC). Bangalore: India Institute of Science,2011： 1-5.

[19]華羅庚.優(yōu)選法平話及其補充[M].北京：國防工業(yè)出版社, 1972.