張海龍，陸廣華，于　慧，王東進(jìn)

(中國科學(xué)技術(shù)大學(xué) 中國科學(xué)院電磁空間信息重點實驗室,　合肥 230027)

?

基于壓縮感知的分布式MIMO成像雷達(dá)網(wǎng)格失配問題研究

張海龍，陸廣華，于慧，王東進(jìn)

(中國科學(xué)技術(shù)大學(xué) 中國科學(xué)院電磁空間信息重點實驗室,合肥 230027)

摘要：在基于壓縮感知的分布式多輸入多輸出(MIMO)成像雷達(dá)系統(tǒng)中，若成像目標(biāo)不在事先劃分的網(wǎng)格點上，即存在網(wǎng)格失配(Off-Grid)偏差時，成像結(jié)果將急劇惡化，文中針對該問題，分別在Off-Grid偏差一階和二階近似下提出了行之有效的解決方法。文中首先建立了系統(tǒng)存在Off-Grid偏差時的回波模型，在理想劃分網(wǎng)格點的基礎(chǔ)上，將Off-Grid偏差看作未知數(shù)，分別進(jìn)行一階和二階Taylor展開，建立擾動模型，分析Off-Grid偏差對成像結(jié)果的影響，最后給出了一種采用欠定系統(tǒng)局灶解法(FOCUSS)和約束總體最小二乘(CTLS)方法聯(lián)合估計目標(biāo)位置和散射系數(shù)的算法。與傳統(tǒng)的基于壓縮感知的成像方法相比，文中所提出的算法能夠很好地解決網(wǎng)格失配問題，提高目標(biāo)的反演精度。計算機(jī)仿真結(jié)果驗證了所提方法的有效性。

關(guān)鍵詞：分布式多輸入多輸出雷達(dá)；壓縮感知；網(wǎng)格失配；欠定系統(tǒng)局灶解法；約束總體最小二乘

0引言

多輸入多輸出(MIMO)雷達(dá)是近幾年提出的一種新體制雷達(dá)，正吸引著越來越多的學(xué)者的注意[1]。MIMO雷達(dá)采用多個發(fā)射天線和多個接收天線，并且不同于傳統(tǒng)的相控陣?yán)走_(dá)發(fā)射相干波形，MIMO雷達(dá)發(fā)射正交波形，接收端通過匹配處理來恢復(fù)各發(fā)射信號分量[2]。MIMO雷達(dá)有集中式和分布式兩大類：集中式MIMO雷達(dá)可以利用波形分集處理，提高雷達(dá)的抗截獲性能和雜波背景中探測低速、弱目標(biāo)的能力；分布式MIMO雷達(dá)可以通過大間距布陣形式獲得空間分

集增益，提高對起伏目標(biāo)的檢測性能[3]。在本文中，本文將重點研究分布式MIMO雷達(dá)。

當(dāng)發(fā)射機(jī)和接收機(jī)數(shù)目受限時，波數(shù)域的覆蓋范圍是不完整的，若仍采用基于匹配濾波的傳統(tǒng)成像方法將得不到好的成像效果[4]。在大多數(shù)的成像雷達(dá)應(yīng)用場景中，目標(biāo)散射點往往是稀疏分布的，也就是說，實際散射點的數(shù)目遠(yuǎn)小于潛在的散射點數(shù)目。目前很多學(xué)者致力于研究MIMO成像雷達(dá)的稀疏恢復(fù)技術(shù)[5-7]，比如正交匹配追蹤算法(OMP)和基追蹤算法(BP)。一般地，壓縮感知(CS)將稀疏恢復(fù)問題表示為y=Φx+e，其中x表示散射系數(shù)，y表示接收回波，Φ表示感知矩陣，e表示成像過程增加的噪聲。我們的

目的是得到x的最稀疏的解，即x中的非零元素個數(shù)最少的解。采用CS技術(shù)，利用目標(biāo)散射點的稀疏性，可以得到較好的成像效果[8]。

現(xiàn)有的大多數(shù)稀疏恢復(fù)算法都要求所有的散射點精確地位于事先劃分好的離散網(wǎng)格點上，但是當(dāng)散射點不在事先劃分的網(wǎng)格點上，即存在網(wǎng)格失配(Off-Grid)偏差時，稀疏恢復(fù)性能將會受到嚴(yán)重影響[9]。在MIMO雷達(dá)成像中，由于散射點分布在連續(xù)場景中，所以即使所劃分的離散網(wǎng)格點很密集，Off-Grid問題還是存在的。使用密集的網(wǎng)格點劃分可能會減輕失配的程度，但是對于Off-Grid目標(biāo)成像問題，這不是一個好的補救措施，因為密集的網(wǎng)格點劃分會顯著地增加感知矩陣的列相關(guān)性，違背可靠稀疏恢復(fù)所要遵循的限制等距條件(RIP)[10]?，F(xiàn)在已經(jīng)有一些學(xué)者致力于研究感知矩陣的失配問題，并將其研究成果應(yīng)用到波達(dá)角(DOA)的Off-Grid問題估算當(dāng)中。ZHU等[11]提出了一種稀疏總體最小二乘的方法(S-TLS)來緩和矩陣失配問題，但是該方法效率較低并且耗時較長。HAN等[12]提出了更快更穩(wěn)定的算法，總體最小欠定系統(tǒng)局灶解法(TLS-FOCUSS)和同步下降欠定系統(tǒng)局灶解法(SD-FOCUSS)，并將SD-FOCUSS算法應(yīng)用到DOA多次測量矢量(MMV)當(dāng)中。YANG等[13]建立了采用稀疏貝葉斯推理(SBI)對DOA的Off-Grid問題進(jìn)行估計的公式并通過期望最大化(EM)迭代得到源信號和矩陣偏差。

本文針對分布式MIMO雷達(dá)成像過程中所遇到的Off-Grid問題，提出了一種聯(lián)合估計目標(biāo)位置和散射系數(shù)的方法。

1網(wǎng)格失配下的成像方法

1.1網(wǎng)格失配一階近似下的成像方法

1.1.1網(wǎng)格失配一階近似下的成像模型

考慮分布式MIMO雷達(dá)，并且假設(shè)發(fā)射機(jī)、接收機(jī)和目標(biāo)位于同一個兩維平面上。系統(tǒng)共有M臺發(fā)射機(jī)和N臺接收機(jī)，第m臺發(fā)射機(jī)位于tm=[tmx,tmy]，第n臺接收機(jī)位于rn=[rnx,rny]。假設(shè)不同發(fā)射機(jī)所發(fā)射的信號是正交的且第m臺發(fā)射機(jī)所發(fā)射的基帶信號為sm(t)，則相應(yīng)的帶通信號可以表示為

sm=Re{sm(t)ej2πfct}

(1)

式中：Re{·}表示取實部；j表示虛數(shù)單位；fc表示載波頻率，并且在不同延時情況下不同發(fā)射信號之間的互相關(guān)為零。

假設(shè)成像目標(biāo)包含K個各向同性散射的散射點，σk和pk=[pkx,pky]分別表示第k個散射點的散射系數(shù)和位置。簡單起見，假設(shè)每個散射點的散射特性不隨觀察角度的改變而改變，并且系統(tǒng)發(fā)射窄帶信號，則第n臺接收機(jī)接收到的回波可以表示為

(2)

(3)

式中：c表示光速。

將所接收到的信號經(jīng)過一組匹配濾波器，由于發(fā)射信號是正交的，則在第q(q=1, 2, …,Q)個采樣時刻，第n臺接收機(jī)的第m個匹配濾波器的輸出可以表示為

(4)

式中：emn(q)為噪聲。

將成像區(qū)域按照最小分辨要求沿x方向劃分為U個離散單元、沿y方向劃分為V個離散單元，則在成像區(qū)域中共有R=UV個網(wǎng)格。假設(shè)所有的散射點都位于網(wǎng)格點上，則可以得到

(5)

其中

(6)

y=Φσ+e

(7)

若存在Off-Grid偏差，即在第r個格點處存在x向偏差Δxr、y向偏差Δyr，則對于?r=1,2,…,R，均有|Δxr|≤0.5ρx，|Δγr|≤0.5ρy，其中，ρx和ρy分別為x和y方向上的網(wǎng)格大小。因此，對感知矩陣做一階Taylor近似，可得

(8)

式中：Φ=[Φx, Φy],Φ為考慮Off-Grid偏差以后計算得到的感知矩陣；Φ0=[Φx0，Φy0]，Φ0為不考慮Off-Grid偏差得到的感知矩陣；Φ′=[Φ′x, Φ′y]，Φ′為感知矩陣對Off-Grid偏差的一階微分；(g)x表示對x維的操作；(g)y表示對y維的操作；Λx1、Λy1分別為由Δxr、Δyr構(gòu)成的對角線矩陣。令Δx1=vec(Δxr)，Δy1=vec(Δyr)則Δ=[Δx1,Δy1]T為Off-Grid偏差。

經(jīng)過上述推導(dǎo)，可得考慮Off-Grid偏差以后將格點位置偏差進(jìn)行一階Taylor展開的成像模型為

y=(Φ+H1Λ1)σ+e

(9)

1.1.2網(wǎng)格失配一階近似下的求解方法

本節(jié)給出一種FOCUSS和CTLS相結(jié)合的CTLS-FOCUSS-1方法來求解式(9)所示問題。其中，F(xiàn)OCUSS用于稀疏恢復(fù)，CTLS用于求解目標(biāo)位置的格點偏差。

步驟1：稀疏恢復(fù)

在該步驟中，通過FOCUSS理論的迭代求解方法進(jìn)行稀疏恢復(fù)。設(shè)迭代標(biāo)號為l，在更新了校正后的目標(biāo)位置以及重新計算了感知矩陣以后，尋找最優(yōu)的估計值σ(l+1)以使得如下代價函數(shù)最小

(10)

令?F1/?σ(l+1)=0，可得

(11)

式(11)的左右兩邊都含有σ，無法直接獲得σ的顯示表達(dá)，因此這里采用迭代松弛方法求解σ

(12)

式中：s為內(nèi)循環(huán)次數(shù)；W(l,s)=W(σ(l))。

當(dāng)內(nèi)循環(huán)s收斂時，即可得σ(l+1)=σ(l+1,s)。

步驟2：格點校正

通過CTLS方法求解格點偏差，進(jìn)而更新格點位置。

定義擾動項μ1=[Bx1Δx1,By1Δy1，1/σee]T。則可以通過求解下式所示的優(yōu)化問題來求解格點偏差

(13)

(14)

(15)

(16)

式中：l為迭代次數(shù)。由原始劃分格點與Off-Grid偏差更新目標(biāo)格點位置

(17)

步驟3：參數(shù)更新

將步驟2中計算得到的目標(biāo)位置p(l+1)計算得到感知矩陣Φ(l+1)。

將迭代次數(shù)加1，重復(fù)步驟1到步驟3，直到σ變化很小或者達(dá)到預(yù)先設(shè)定的迭代次數(shù)。

1.2網(wǎng)格失配二階近似下的成像方法

1.2.1網(wǎng)格失配二階近似下的成像模型

前面所述方法都是基于1.1.1節(jié)中感知矩陣對Off-Grid偏差進(jìn)行泰勒展開并保留一階近似的模型上進(jìn)行的，接下來將推導(dǎo)Off-Grid偏差的二階近似模型，并給出相應(yīng)的CTLS-FOCUSS-2算法步驟。

當(dāng)存在Off-Grid偏差時，如式(18)所示，對感知矩陣做二階Taylor近似，可得

(18)

式中：Φ″為感知矩陣對Off-Grid偏差的二階微分；Λx2、Λy2分別為由(Δxr)2、(Δyr)2構(gòu)成的對角線矩陣，并令Δx2=(vec((Δxr)2)，Δy2=vec((Δyr)2)。

則經(jīng)過上述推導(dǎo)，可得考慮Off-Grid偏差以后，將格點位置偏差進(jìn)行二階Taylor展開的成像模型為

y=(Φ+H1Λ1+H2Λ2)σ+e

(19)

式中：H2=[Φ″x， Φ″y]；Λ2=[Λx2,Λy2]T。

1.2.2網(wǎng)格失配二階近似下的求解方法

對于該問題的求解，與Off-Grid偏差一階近似的過程相似，僅需在1.1節(jié)的求解過程中，替換第二步即可，下面僅給出在Off-Grid偏差二階近似的情況下采用CTLS方法求解Off-Grid偏差的過程。

定義擾動項μ2=[Bx1Δx1, By1Δy1,Bx2Δx2,Bx2Δy2,1/σee]T，則可以通過求解下式的優(yōu)化問題來求解格點偏差

(20)

(21)

(22)

其中，式(21)仍可采用牛頓法進(jìn)行求解，設(shè)置的初值為該問題的最小二乘解。由μ2進(jìn)而可以反解出Off-Grid偏差Δx1和Δy1

(23)

按照式(17)由原始劃分格點與Off-Grid偏差得到更新后的目標(biāo)格點位置p(l+1)。

2仿真實驗及分析

在本節(jié)中，我們通過給出一些仿真結(jié)果來證明所提算法CTLS-FOCUSS-1和CTLS-FOCUSS-2的有效性，所做的仿真是在Matlab中完成的，并且CTLS-FOCUSS-1和CTLS-FOCUSS-2的迭代終止條件均設(shè)置為 或者達(dá)到最大迭代次數(shù)50。

仿真參數(shù)設(shè)置如表1所示。

表1　仿真參數(shù)

目標(biāo)的散射點分布情況如圖1所示。在成像場景中，沿x方向和y方向均劃分為40個離散單元，并且相鄰單元之間的間隔為1 m。從圖1可以看出，這里共有10個散射點，并且散射點均不在事先離散劃分好的網(wǎng)格上。

圖1　目標(biāo)散射點分布情況

在下面的仿真當(dāng)中，除了我們所提算法CTLS-FOCUSS-1和CTLS-FOCUSS-2，還包括其他的算法：MF，OMP，F(xiàn)OCUSS，和TLS-FOCUSS。TLS-FOCUSS算法是在對Off-Grid偏差進(jìn)行一階Taylor展開的基礎(chǔ)上，假設(shè)所加噪聲為高斯白噪聲并且Off-Grid偏差服從高斯分布的情形下，通過對式(24)的求解得到的最大后驗解(MAP)。

(24)

采用MF，OMP，F(xiàn)OCUSS，TLS-FOCUSS和我們所提算法(CTLS-FOCUSS-1，CTLS-FOCUSS-2)的恢復(fù)結(jié)果如圖2a)~f)所示，其中，圓圈表示的是真實散射點的位置。如前所述，由于波數(shù)域覆蓋不是完整的，所以采用匹配濾波的方法得不到好的恢復(fù)結(jié)果。由于Off-Grid散射點的存在，OMP和FOCUSS算法也是失效的。與OMP和FOCUSS算法相比，TLS-FOCUSS算法能夠得到改善的恢復(fù)結(jié)果，但是恢復(fù)結(jié)果中存在虛假散射點。CTLS-FOCUSS-1是在Off-Grid偏差一階近似的情況下對網(wǎng)格失配偏差做校正，所以其恢復(fù)精度有限。相比之下，CTLS-FOCUSS-2在當(dāng)前仿真條件下能

圖2　Off-Grid目標(biāo)的恢復(fù)結(jié)果

夠?qū)ff-Grid散射點成像并得到好的恢復(fù)結(jié)果。

下面我們將研究成像誤差與信噪比(SNR)之間的關(guān)系。SNR在-5 dB~30 dB變化，變化間隔為5 dB，目標(biāo)恢復(fù)誤差和Off-Grid偏差均是通過對超過30次的仿真試驗取平均得到的結(jié)果。在每次仿真試驗中，Off-Grid偏差在沿x方向或者y方向的一個離散單元內(nèi)均勻分布。

采用OMP，F(xiàn)OCUSS，TLS-FOCUSS，CTLS-FOCUSS-1和CTLS-FOCUSS-2算法的目標(biāo)恢復(fù)結(jié)果的歸一化均方誤差(NMSE)隨SNR的變化曲線如圖3a)所示。從圖3a)中可以看出，隨著信噪比的增大，采用OMP和FOCUSS算法的目標(biāo)恢復(fù)誤差仍然維持在一個較高水平。采用TLS-FOCUSS和CTLS-FOCUSS-1算法將得到改善的恢復(fù)結(jié)果，并且后者的目標(biāo)恢復(fù)效果和Off-Grid恢復(fù)誤差均優(yōu)于前者，這是因為在仿真過程中，Off-Grid偏差是均勻分布的，TLS-FOCUSS算法中有關(guān)Off-Grid偏差Λ的假設(shè)與之不符。從圖3a)和圖3b)可以看出，采用CTLS-FOCUSS-2算法可以得到好的目標(biāo)恢復(fù)結(jié)果，并且Off-Grid恢復(fù)誤差得到明顯提高，可見在對Off-Grid偏差進(jìn)行二階Taylor展開的情況下采用該方法求解目標(biāo)位置和散射系數(shù)的有效性。

圖3　不同信噪比下幾種算法恢復(fù)結(jié)果對比

3結(jié)束語

針對在分布式MIMO雷達(dá)中使用壓縮感知成像所面臨的Off-Grid問題，分別在網(wǎng)格偏差一階和二階近似下提出了行之有效的CTLS-FOCUSS-1和CTLS-FOCUSS-2算法。該算法通過循環(huán)迭代的過程，將非凸優(yōu)化問題轉(zhuǎn)化為三個主要步驟：稀疏恢復(fù)、格點校正和參數(shù)更新，實現(xiàn)Off-Grid偏差的自適應(yīng)調(diào)整，具有很好的抗噪性，相比傳統(tǒng)的基于壓縮感知的成像方法能夠得到更高的恢復(fù)精度。本文最后通過數(shù)字仿真驗證了所提算法的有效性，顯示了其在實際系統(tǒng)中的應(yīng)用前景。

參 考 文 獻(xiàn)

[1]陳浩文，黎湘，莊釗文. 一種新興的雷達(dá)體制——MIMO雷達(dá)[J]. 電子學(xué)報，2012, 40(6): 1190-1198.

CHEN Haowen, LI Xiang, ZHUANG Zhaowen. A rising radar system——MIMO radar[J]. Actc Electronica Sinica, 2012, 40(6): 1190-1198.

[2]FISHLER E, HAIMOVICH A, BLUM R，et al. MIMO radar: an idea whose time has come[C]// IEEE Radar Conference. Boston, MA: IEEE Press, USA, 2007: 71-78.

[3]裴英，肖文書. 相干MIMO雷達(dá)關(guān)鍵技術(shù)研究[J]. 現(xiàn)代雷達(dá)，2013, 35(3): 9-12.

PEI Ying, XIAO Wenshu. A study on key technology of the coherent MIMO radar[J]. Modern Radar, 2013, 35(3): 9-12.

[4]LIU C, CHEN W. Sparse frequency diverse mimo radar imaging[C]//2012 Conference on Signals, Systems and Computers. [S.l.]: IEEE Press, 2012: 853-857.

[5]CHEN C Y, VAIDYANATHAN P P. Compressive sensing in MIMO radar[C]// In Proceedings of 42nd Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA: IEEE Press, 2008: 41-44.

[6]YU Y, PETROPULU A P, POOR H V. CSSF MIMO radar: Low-complexity compressive sensing based MIMO radar that uses step frequency[J]. IEEE Transactions on Aerospace Electronic System, 2012(48): 1490-1504.

[7]TAN X, ROBERTS W, LI J, et al. Sparse learning via iterative minimization with application to MIMO radar imaging[J]. IEEE Transactions on Signal Processing, 2011, 59(3): 1088-1101.

[8]YU Y, PETROPULU A P, POOR H V. MIMO radar using compressive sampling[J]. IEEE Journal of Selected Topics in Signal Processing, 2010, 4(1): 146-163.

[9]CHI Y, PEZESHKI A, SCHARF L, et al. Sensitivity to basis mismatch in compressed sensing[J]. IEEE Transactions on Signal Processing, 2011, 59(3): 2182-2195.

[10]HE X, LIU C, LIU B, et al. Sparse frequency diverse MIMO radar imaging for off-grid target based on adaptive iterative MAP[J]. Remote Sensing, 2013, 5(2): 631-647.

[11]ZHU H, LEUS G, GIANNAKIS G B. Sparsity-cognizant total least-squares for perturbed compressive sampling[J]. IEEE Transactions on Signal Processing, 2011, 59(5): 2002-2016.

[12]HAN X, ZHANG H, LI G. Fast algorithms for sparse recovery with perturbed dictionary[J]. arXiv preprint arXiv:1111.6237, 2011.

[13]YANG Z, XIE L, ZHANG C. Off-grid direction of arrival estimation using sparse Bayesian inference[J]. IEEE Transactions on Signal Processing, 2013, 61(1): 38-43.

張海龍男，1989年生，在讀碩士。研究方向為微波成像、陣列信號處理和壓縮感知。

陸廣華男，1970年生，副教授。研究方向為微波與毫米波技術(shù)及其在精確制導(dǎo)、通訊系統(tǒng)、雷達(dá)系統(tǒng)等領(lǐng)域中的應(yīng)用。

于慧女，1988年生，碩士。研究方向為微波成像、陣列信號處理和壓縮感知。

王東進(jìn)男，1955年生，教授。研究方向為微波與毫米波技術(shù)及其在精確制導(dǎo)、通信系統(tǒng)、雷達(dá)系統(tǒng)等領(lǐng)域中的應(yīng)用。

A Study on Off-Grid Problem for Distributed MIMO

Imaging Radar Based on Compressive Sensing

ZHANG Hailong，LU Guanghua，YU Hui，WANG Dongjin

(Key Laboratory of Electromagnetic Space Insformation of CAS, USTC,Hefei 230027, China)

Abstract：In distributed MIMO radar system based on compressive sensing, if the imaging target is not positioned at the pre-discretized grid locations, that is to say, there exists Off-Grid deviation in the system, the imaging results will be deteriorated sharply. In order to solve this problem, an effective solution under one order and two order approximation of Off-Grid deviation respectively is proposed in this paper. The echo model when there exists Off-Grid deviation in the system is established first and then, based on the ideal discretized grid, the perturbation model is established by regarding the Off-Grid deviation as unknown and carrying out one order and two order Taylor expansion of it. After that, the influence of Off-Grid deviation on imaging result is analysed. Finally, the algorithm based on FOCUSS and CTLS is proposed to jointly estimate the target location and its scattering coefficient. Compared with the traditional imaging method based on compressive sensing, the algorithm proposed in this paper can solve the Off-Grid problem effectively and improve the inversion accuracy of target recovery. The effectiveness of the proposed algorithm is confirmed by several experimental results.

Key words：distributed MIMO radar; compressive sensing; Off-Grid; FOCUSS; CTLS

收稿日期：2015-07-26

修訂日期：2015-09-28

通信作者：張海龍Email：hailong@mail.ustc.edu.cn

中圖分類號：TN957.52

文獻(xiàn)標(biāo)志碼：A

文章編號：1004-7859(2015)12-0049-05