劉記紅， 徐少坤， 韓國(guó)強(qiáng)， 魏雁飛

(中國(guó)洛陽(yáng)電子裝備試驗(yàn)中心， 河南洛陽(yáng) 471003)

0 引言

彈道目標(biāo)識(shí)別是導(dǎo)彈攻防對(duì)抗勝負(fù)的關(guān)鍵環(huán)節(jié)。逆合成孔徑雷達(dá)(Inverse Synthetic Aperture Radar, ISAR)成像能夠反映目標(biāo)的尺寸、結(jié)構(gòu)、形狀等特征[1]，在雷達(dá)目標(biāo)識(shí)別領(lǐng)域發(fā)揮著極其重要的作用，已成為彈道目標(biāo)識(shí)別的主要途徑之一。對(duì)于彈道目標(biāo)而言，除高速軌道運(yùn)動(dòng)外，通常還伴隨有自旋、進(jìn)動(dòng)等復(fù)雜的運(yùn)動(dòng)形式，大大增加了ISAR成像的難度。進(jìn)動(dòng)是彈道中段彈頭目標(biāo)的重要物理特性之一，它使得目標(biāo)姿態(tài)角呈現(xiàn)近似正弦規(guī)律的非均勻變化，目標(biāo)上各散射點(diǎn)在距離上出現(xiàn)類(lèi)似正弦曲線的距離走動(dòng)，導(dǎo)致相應(yīng)的ISAR成像問(wèn)題較為復(fù)雜，常規(guī)成像算法難以獲得高質(zhì)量的ISAR圖像[2]。

彈道目標(biāo)成像問(wèn)題是ISAR成像領(lǐng)域的研究熱點(diǎn)，由于軍事敏感性，國(guó)外相關(guān)的研究成果很少見(jiàn)諸報(bào)道，但可從零星的報(bào)道中獲知某些導(dǎo)彈防御雷達(dá)具備了彈頭ISAR成像能力。國(guó)內(nèi)的研究起步較晚，直到近些年才有少量的文獻(xiàn)出現(xiàn)。目前，用于進(jìn)動(dòng)目標(biāo)的成像方法主要有基于時(shí)頻變換的方法[3-5]、基于復(fù)數(shù)逆投影的方法[6]、基于寬帶復(fù)數(shù)后向投影變換(Back Projection Transform, BPT)的方法[2]等。上述方法在一定條件下能夠獲得較好的成像結(jié)果，但對(duì)回波數(shù)據(jù)要求較高，成像性能敏感于目標(biāo)的進(jìn)動(dòng)參數(shù)，而在實(shí)際應(yīng)用過(guò)程中，獲取的數(shù)據(jù)段或目標(biāo)進(jìn)動(dòng)參數(shù)往往難以滿足要求，如彈道目標(biāo)的高速運(yùn)動(dòng)、反導(dǎo)系統(tǒng)雷達(dá)對(duì)多目標(biāo)的跟蹤需求、導(dǎo)彈突防中的電子干擾措施等因素都可能導(dǎo)致無(wú)法獲得充足的觀測(cè)數(shù)據(jù)，需要進(jìn)一步研究利用有限脈沖回波數(shù)據(jù)進(jìn)行彈道目標(biāo)ISAR成像的方法。

近年來(lái)，壓縮感知(Compressed Sensing, CS)理論在高分辨雷達(dá)成像領(lǐng)域受到了越來(lái)越多的關(guān)注[7-9]，其利用較少的測(cè)量數(shù)據(jù)即可重構(gòu)出雷達(dá)目標(biāo)圖像的特性，為進(jìn)一步完善傳統(tǒng)雷達(dá)成像技術(shù)提供了契機(jī)，也為彈道中段進(jìn)動(dòng)目標(biāo)成像提供了一種新思路。文獻(xiàn)[10]將CS理論引入到目標(biāo)微動(dòng)特性分析中，獲取了目標(biāo)回波的微多普勒時(shí)頻譜；文獻(xiàn)[11]研究了CS理論在含旋轉(zhuǎn)部件目標(biāo)成像中的應(yīng)用，獲取該類(lèi)目標(biāo)主體部分的ISAR圖像。本文針對(duì)彈道中段進(jìn)動(dòng)目標(biāo)的ISAR成像問(wèn)題，在分析旋轉(zhuǎn)對(duì)稱進(jìn)動(dòng)目標(biāo)回波模型的基礎(chǔ)上，充分利用雷達(dá)目標(biāo)散射率分布的稀疏性先驗(yàn)信息和目標(biāo)的進(jìn)動(dòng)信息，提出了基于CS的進(jìn)動(dòng)目標(biāo)成像算法，在降低脈沖重復(fù)頻率(Pulse Repetition Frequency, PRF)實(shí)際需求的同時(shí)，實(shí)現(xiàn)了進(jìn)動(dòng)目標(biāo)的高分辨成像。

1 彈道中段進(jìn)動(dòng)目標(biāo)回波模型

進(jìn)動(dòng)是一種自旋加錐旋的復(fù)合運(yùn)動(dòng)形式，彈道中段進(jìn)動(dòng)目標(biāo)的運(yùn)動(dòng)模型如圖1所示，其中3個(gè)坐標(biāo)系的建立方式同文獻(xiàn)[2]，O為目標(biāo)質(zhì)心，ωs,ωc分別為自旋和進(jìn)動(dòng)角頻率，目標(biāo)進(jìn)動(dòng)軸ON在YOZ平面內(nèi)，俯仰角為βN。假設(shè)初始時(shí)刻目標(biāo)自旋軸(z軸)也位于YOZ平面內(nèi)，Z軸與z軸之間的夾角為θ，則目標(biāo)的進(jìn)動(dòng)角為γ=π/2-θ-βN。

圖1 雷達(dá)與空間進(jìn)動(dòng)目標(biāo)幾何關(guān)系示意圖

彈道中段目標(biāo)的形體結(jié)構(gòu)比較簡(jiǎn)單，除少部分戰(zhàn)術(shù)彈道導(dǎo)彈外，彈頭及誘餌大都是旋轉(zhuǎn)對(duì)稱體。文獻(xiàn)[2]的分析表明，對(duì)于旋轉(zhuǎn)對(duì)稱體目標(biāo)而言，自旋不會(huì)帶來(lái)散射場(chǎng)的變化，目標(biāo)的三維姿態(tài)運(yùn)動(dòng)可以等效為在二維成像平面內(nèi)的旋轉(zhuǎn)，轉(zhuǎn)動(dòng)規(guī)律對(duì)應(yīng)目標(biāo)的姿態(tài)角變化規(guī)律。將目標(biāo)在三維空間的旋轉(zhuǎn)運(yùn)動(dòng)等效為在二維成像平面內(nèi)目標(biāo)對(duì)稱軸繞點(diǎn)O的轉(zhuǎn)動(dòng)，則等效旋轉(zhuǎn)角，即目標(biāo)的姿態(tài)角為

α(t)=arc cos(cosγsinβN+sinγcosβNcosωct)

(1)

由于進(jìn)動(dòng)角較小，通常為5°～15°，對(duì)應(yīng)α(t)的變化范圍較小，可用直線方程近似，即

α(t)≈arc cos(sinβN+γcosβNcosωct)≈

(2)

式中，θc=π/2-βN為雷達(dá)觀測(cè)目標(biāo)的平均視線角，它隨目標(biāo)的軌道運(yùn)動(dòng)而變化，相對(duì)于微動(dòng)帶來(lái)的姿態(tài)變化，θc的變化比較緩慢。

當(dāng)將目標(biāo)的進(jìn)動(dòng)等效為二維平面內(nèi)的轉(zhuǎn)動(dòng)時(shí)，目標(biāo)的等效散射中心模型也可等效到二維平面內(nèi)。設(shè)等效成像平面為Z′OY′，則在Z′OY′平面內(nèi)以α(t)為規(guī)律轉(zhuǎn)動(dòng)的目標(biāo)回波與進(jìn)動(dòng)狀態(tài)下的旋轉(zhuǎn)對(duì)稱目標(biāo)回波等效。假設(shè)雷達(dá)發(fā)射線性調(diào)頻信號(hào)，目標(biāo)共包含I個(gè)散射點(diǎn)，等效成像平面內(nèi)目標(biāo)上第i個(gè)散射點(diǎn)的坐標(biāo)為(zi,yi)，散射強(qiáng)度為δi，則消除了平動(dòng)影響的雷達(dá)目標(biāo)回波可表示為

ΔRi(tm)=zicosα(tm)+yisinα(tm)

(4)

根據(jù)前述分析，目標(biāo)的姿態(tài)角主要取決于其進(jìn)動(dòng)參數(shù)，目前關(guān)于進(jìn)動(dòng)參數(shù)的估計(jì)算法有很多[12]，若已經(jīng)得到目標(biāo)的平均視線角和進(jìn)動(dòng)參數(shù)，則可獲得目標(biāo)的姿態(tài)角變化規(guī)律，進(jìn)而基于式(5)重構(gòu)目標(biāo)的二維圖像。由式(2)可知，進(jìn)動(dòng)使得中段目標(biāo)的成像面臨如下問(wèn)題： 1) 進(jìn)動(dòng)使得方位向的采樣間隔呈類(lèi)正弦規(guī)律變化，期間目標(biāo)姿態(tài)角在兩倍的進(jìn)動(dòng)角范圍內(nèi)非均勻變化，且不能用勻加速轉(zhuǎn)動(dòng)模型進(jìn)行等效； 2) 與高速自旋目標(biāo)不同，由于進(jìn)動(dòng)角一般較小，進(jìn)動(dòng)過(guò)程中目標(biāo)的姿態(tài)角變化范圍不大，基于窄帶信息的自旋目標(biāo)成像方法[13]不適用于進(jìn)動(dòng)目標(biāo)。

2 基于壓縮感知的進(jìn)動(dòng)目標(biāo)成像

鑒于CS利用少量測(cè)量數(shù)據(jù)即可實(shí)現(xiàn)高分辨雷達(dá)成像的特點(diǎn)，這里將CS引入到旋轉(zhuǎn)對(duì)稱進(jìn)動(dòng)目標(biāo)的ISAR成像中，其原理是基于目標(biāo)回波信號(hào)的稀疏性，利用少量非相干測(cè)量通過(guò)非線性優(yōu)化重構(gòu)雷達(dá)目標(biāo)圖像。

2.1 回波模型的線性化表征

雷達(dá)回波的稀疏性表征和模型的線性化處理是CS成像的先決條件。將成像區(qū)間離散化為P×Q的二維網(wǎng)格，每個(gè)網(wǎng)格點(diǎn)代表可能存在散射點(diǎn)的位置坐標(biāo)(zp,yq)，相應(yīng)位置處的散射強(qiáng)度表示為δpq(p=0,1,…,P-1,q=0,1,…,Q-1)。離散間隔大小決定了重構(gòu)圖像的分辨率，當(dāng)某一個(gè)網(wǎng)格交點(diǎn)的坐標(biāo)(zp,yq)上存在等效散射中心時(shí)，此網(wǎng)格點(diǎn)上的散射強(qiáng)度δpq≠0；反之，δpq=0。由于目標(biāo)所涵蓋的僅是一部分網(wǎng)格交點(diǎn)，且目標(biāo)僅包含有限個(gè)等效散射中心，PQ個(gè)網(wǎng)格交點(diǎn)中絕大多數(shù)位置上都不存在散射中心，因此目標(biāo)散射率分布具有很強(qiáng)的稀疏性。

假設(shè)距離向的采樣點(diǎn)數(shù)為N(N≤P)，對(duì)應(yīng)的基頻采樣點(diǎn)為fn(n=0,1,…,N-1)，方位向用于成像的脈沖數(shù)為M(M≤Q)，則根據(jù)式(5)可構(gòu)造進(jìn)動(dòng)目標(biāo)的觀測(cè)矩陣如下：

Φ=[Φ(0),Φ(1),…,Φ(M-1)]

(6)

設(shè)Y為測(cè)得的N×M維的回波數(shù)據(jù)矩陣，δ=[δpq]P×Q為目標(biāo)的二維散射率分布矩陣，按列堆疊后的矢量化表示形式分別為y=vec(Y)和σ=vec(δ)，則式(5)可表示為如下矩陣形式：

y=Φσ

(7)

2.2 雷達(dá)圖像形成

根據(jù)CS理論，當(dāng)式(7)中的Φ滿足一定條件[8]時(shí)，通過(guò)非線性優(yōu)化即可重構(gòu)σ，進(jìn)而通過(guò)重排得到目標(biāo)的二維圖像。該方法的優(yōu)勢(shì)在于，能夠克服方位向非均勻采樣的影響，保證角度變化范圍較小情況下的成像分辨率，且當(dāng)PRF不滿足Nyquist采樣準(zhǔn)則而使回波產(chǎn)生欠采樣時(shí)依然有效。值得注意的是，此時(shí)測(cè)量矩陣Φ的維數(shù)較高(MN×PQ)，將在計(jì)算中占據(jù)大量的存儲(chǔ)空間，而且大矩陣的乘法和求逆運(yùn)算計(jì)算量較大，使得算法整體復(fù)雜度較高。事實(shí)上， CS重構(gòu)過(guò)程中無(wú)需采用所有的回波數(shù)據(jù)，可對(duì)雷達(dá)回波進(jìn)行降采樣處理。假設(shè)距離向隨機(jī)采樣K點(diǎn)，方位向獲得的脈沖回波數(shù)為L(zhǎng)(K

上述方法通過(guò)CS算法實(shí)現(xiàn)了進(jìn)動(dòng)目標(biāo)雷達(dá)圖像的二維聯(lián)合重構(gòu)，所需的存儲(chǔ)量和計(jì)算量較大。根據(jù)觀測(cè)矩陣元素的構(gòu)造過(guò)程可知，目標(biāo)散射率在方位向和距離向的分布均與目標(biāo)姿態(tài)角(慢時(shí)間)耦合在一起，無(wú)法得到觀測(cè)模型的解耦形式，不能通過(guò)距離向和方位向分別重構(gòu)的方式降低算法的復(fù)雜度。為進(jìn)一步降低所提成像算法的存儲(chǔ)需求，可利用每個(gè)測(cè)量頻點(diǎn)的信息分別進(jìn)行重構(gòu)，然后對(duì)各個(gè)頻點(diǎn)的重構(gòu)結(jié)果進(jìn)行相干疊加，進(jìn)而通過(guò)重排得到目標(biāo)的高分辨ISAR圖像，這里稱之為“分頻處理方法”。

yn=Φnσ

(8)

通過(guò)CS算法求解式(8)可得到第n個(gè)測(cè)量頻點(diǎn)處的重構(gòu)結(jié)果，對(duì)各個(gè)頻點(diǎn)數(shù)據(jù)執(zhí)行相同的操作，然后疊加合成重構(gòu)結(jié)果，即可得到清晰的目標(biāo)ISAR圖像。

3 算法性能分析

從理論上講，所提基于CS的進(jìn)動(dòng)目標(biāo)成像算法充分利用了雷達(dá)目標(biāo)散射率分布的稀疏性先驗(yàn)信息和目標(biāo)的進(jìn)動(dòng)信息，通過(guò)少量回波數(shù)據(jù)即可實(shí)現(xiàn)進(jìn)動(dòng)目標(biāo)的高分辨成像，有效降低了進(jìn)動(dòng)目標(biāo)成像所需的PRF，提高了數(shù)據(jù)利用率，付出的代價(jià)是重構(gòu)過(guò)程復(fù)雜度的增加。從成像質(zhì)量上來(lái)看，二維數(shù)據(jù)聯(lián)合處理的重構(gòu)模型與雷達(dá)觀測(cè)模型較匹配，性能相對(duì)穩(wěn)定，能獲得較理想的成像結(jié)果；分頻處理方法將各個(gè)頻點(diǎn)測(cè)量數(shù)據(jù)分別處理，忽略了各頻點(diǎn)回波數(shù)據(jù)間的耦合效應(yīng)，在相同測(cè)量條件下圖像質(zhì)量將有所損失。

從資源需求和運(yùn)算復(fù)雜度上看，以運(yùn)算效率較高的SL0算法[14]為例，二維聯(lián)合成像算法需要的總存儲(chǔ)量為O(KLPQ)，每次迭代運(yùn)行時(shí)間為O(KLPQ)，總的計(jì)算復(fù)雜度為O(JL0KLPQ)，其中J=O(1),L0=O(1)分別表示外層和內(nèi)層循環(huán)次數(shù)；分頻處理方法需要的存儲(chǔ)量為O(LPQ)，每個(gè)頻點(diǎn)每次迭代重構(gòu)的時(shí)間復(fù)雜度為O(LPQ)，共計(jì)重復(fù)K次，故總的運(yùn)行時(shí)間為O(JL0KLPQ)，與二維聯(lián)合成像算法在一個(gè)量級(jí)上。可見(jiàn)，在同等測(cè)量條件下，二者的運(yùn)算開(kāi)銷(xiāo)大致相同，但二維聯(lián)合處理方式需要更多的存儲(chǔ)空間。實(shí)際中，彈道中段進(jìn)動(dòng)目標(biāo)的形體較小且結(jié)構(gòu)簡(jiǎn)單，上述兩種成像處理方法均是可行的。

4 仿真實(shí)驗(yàn)及分析

實(shí)驗(yàn)中利用仿真數(shù)據(jù)和暗室合成數(shù)據(jù)驗(yàn)證了所提成像算法的有效性，由于篇幅限制，這里僅給出暗室數(shù)據(jù)結(jié)果。暗室測(cè)量目標(biāo)的外形尺寸如圖2所示，測(cè)量掃頻范圍為8～12GHz，掃頻間隔20MHz，目標(biāo)姿態(tài)為俯仰角0°，方位角0°～90°變化(鼻翼方向?yàn)?°)，變化間隔為0.2°。選用8～12GHz頻段和相應(yīng)姿態(tài)下的測(cè)量數(shù)據(jù)合成進(jìn)動(dòng)目標(biāo)的回波，從而驗(yàn)證所提CS成像算法的性能，CS方法重構(gòu)過(guò)程中采用SL0算法。合成參數(shù)設(shè)置如下：目標(biāo)起始姿態(tài)角為0°，進(jìn)動(dòng)周期為1s，雷達(dá)信號(hào)帶寬B=4 GHz，測(cè)量頻點(diǎn)數(shù)N=201，PRF為500 Hz，先后合成進(jìn)動(dòng)角為5°和10°時(shí)的目標(biāo)回波，成像過(guò)程中采用一個(gè)進(jìn)動(dòng)周期內(nèi)的回波數(shù)據(jù)。

圖2 暗室測(cè)量目標(biāo)外形尺寸圖

圖3給出了進(jìn)動(dòng)角為10°時(shí)合成進(jìn)動(dòng)目標(biāo)回波的距離壓縮結(jié)果和常規(guī)距離-多普勒(Range-Doppler, RD)算法成像結(jié)果?？梢?jiàn)，進(jìn)動(dòng)目標(biāo)的距離-慢時(shí)間域回波具有一定的周期性，合成數(shù)據(jù)反映了進(jìn)動(dòng)目標(biāo)的姿態(tài)變化規(guī)律，常規(guī)的RD算法無(wú)法獲得清晰的ISAR圖像。

(b)RD成像結(jié)果圖3 暗室合成進(jìn)動(dòng)目標(biāo)距離-慢時(shí)間域回波及RD成像結(jié)果

為驗(yàn)證所提CS成像方法的有效性，并與現(xiàn)有的BPT方法進(jìn)行比較，對(duì)目標(biāo)進(jìn)動(dòng)角γ=5°，10°的情況，分別先后采用BPT方法、二維數(shù)據(jù)聯(lián)合處理的CS方法和分頻處理的CS方法進(jìn)行成像，結(jié)果如圖4、圖5所示。其中，圖示結(jié)果的動(dòng)態(tài)范圍取為30 dB，為便于存儲(chǔ)和計(jì)算， CS方法重構(gòu)過(guò)程中對(duì)回波數(shù)據(jù)進(jìn)行了降采樣，距離向隨機(jī)采樣31個(gè)頻點(diǎn)，方位向進(jìn)行19倍的均勻降采樣，分頻處理方式未對(duì)距離向進(jìn)行降采樣。

(a)γ=5°

(b)γ=10°圖4基于BPT的進(jìn)動(dòng)目標(biāo)成像結(jié)果

(a)二維聯(lián)合處理(γ=5°)

(b)二維聯(lián)合處理(γ=10°)

(c)分頻處理(γ=5°)

(d)分頻處理(γ=10°)圖5基于CS的進(jìn)動(dòng)目標(biāo)成像結(jié)果

從重構(gòu)結(jié)果可以看出， BPT方法和CS方法均能夠有效地對(duì)進(jìn)動(dòng)目標(biāo)成像。相比之下，BPT方法的成像結(jié)果具有較高的旁瓣和較多的背景雜波，且對(duì)進(jìn)動(dòng)角大小較為敏感；CS方法特別是二維聯(lián)合處理的成像結(jié)果清晰地反映了目標(biāo)的散射點(diǎn)分布，具有較好的旁瓣抑制能力，成像性能受進(jìn)動(dòng)角大小的影響較小，且適用于獲得數(shù)據(jù)十分有限的情況；分頻處理的CS成像質(zhì)量略低于二維聯(lián)合處理方法，這是因?yàn)樗盍蚜烁鞑蓸宇l點(diǎn)數(shù)據(jù)間的相關(guān)性，與理論分析一致。

為分析帶寬對(duì)成像性能的影響，圖6給出了合成回波數(shù)據(jù)帶寬為2 GHz(9～11 GHz)、目標(biāo)進(jìn)動(dòng)角為10°情況下的成像結(jié)果?？梢钥闯觯?BPT方法的成像性能受帶寬影響較大，這是由于它在重構(gòu)過(guò)程中利用的是目標(biāo)的一維距離像信息，受sinc函數(shù)主瓣寬度的影響，帶寬越小影響越大。相比之下， CS方法利用的是目標(biāo)的回波數(shù)據(jù)，盡管成像質(zhì)量也隨著帶寬的降低有所下降，但仍較為準(zhǔn)確清晰地反映了目標(biāo)的二維散射率分布。

(a)BPT方法

(b)CS方法(二維聯(lián)合處理)圖6帶寬2 GHz條件下的進(jìn)動(dòng)目標(biāo)成像結(jié)果

此外，實(shí)驗(yàn)過(guò)程中發(fā)現(xiàn)，一味地增加CS重構(gòu)所用的樣本數(shù)并不一定能改善成像質(zhì)量，相反可能導(dǎo)致成像結(jié)果的嚴(yán)重失真。分析原因可知，與高速自旋目標(biāo)不同，由于進(jìn)動(dòng)角較小，進(jìn)動(dòng)目標(biāo)的姿態(tài)角變化范圍相對(duì)較小，當(dāng)采樣較密時(shí)，稀疏字典的原子間可能具有較高的相似性，使得重構(gòu)性能對(duì)噪聲非常敏感，只有在高信噪比下才能獲得較好的重構(gòu)結(jié)果；當(dāng)采樣較疏時(shí)，稀疏字典原子間的相似性較低，重構(gòu)性能對(duì)噪聲更為穩(wěn)健。

5 結(jié)束語(yǔ)

以彈道中段彈頭目標(biāo)為代表的進(jìn)動(dòng)目標(biāo)成像是ISAR成像研究的一個(gè)重要方面。本文在詳細(xì)分析彈道中段進(jìn)動(dòng)目標(biāo)電磁散射特性的基礎(chǔ)上，研究了基于壓縮感知的進(jìn)動(dòng)目標(biāo)成像算法。首先提出了基于二維數(shù)據(jù)聯(lián)合處理的CS成像方法，然后針對(duì)該方法重構(gòu)過(guò)程中存在的存儲(chǔ)空間需求較大的問(wèn)題，提出了基于分頻處理的進(jìn)動(dòng)目標(biāo)成像方法，對(duì)兩種算法的性能進(jìn)行了分析比較，通過(guò)暗室合成數(shù)據(jù)驗(yàn)證了所提方法的有效性。與現(xiàn)有方法相比，所提方法不僅可以在少量數(shù)據(jù)條件下實(shí)現(xiàn)進(jìn)動(dòng)目標(biāo)的高分辨成像，充分改善了成像質(zhì)量，而且提高了對(duì)目標(biāo)進(jìn)動(dòng)參數(shù)的穩(wěn)定性。

[1]金勝,朱天林. ISAR高分辨率成像方法綜述[J]. 雷達(dá)科學(xué)與技術(shù), 2016, 14(3):251-260. JIN Sheng, ZHU Tianlin. A Review of High-Resolution Inverse Synthetic Aperture Radar Imaging Methods[J]. Radar Science and Technology, 2016, 14(3):251-260.(in Chinese)

[2]胡杰民. 復(fù)雜運(yùn)動(dòng)目標(biāo)高分辨雷達(dá)成像技術(shù)研究[D]. 長(zhǎng)沙：國(guó)防科學(xué)技術(shù)大學(xué), 2011.

[3]DING Xiaofeng, FAN Meimei, WEI Xizhang, et al. Narrowband Imaging Method for Spatial Precession Cone-Shaped Targets[J]. Science China Technological Sciences, 2010, 53(4):942-949.

[4]WANG Tao, WANG Xuesong, CHANG Yuliang, et al. Estimation of Precession Parameters and Generation of ISAR Images of Ballistic Missile Targets[J]. IEEE Trans on Aerospace and Electronic Systems, 2010, 46(4):1983-1995.

[5]PAN Xiaoyi, WANG Wei, LIU Jin, et al. Modulation Effect and Inverse Synthetic Aperture Radar Imaging of Rotationally Symmetric Ballistic Targets with Precession[J]. IET Radar, Sonar & Navigation, 2013, 7(9):950-958.

[6]胡曉偉,童寧寧,王建業(yè),等. 基于組網(wǎng)雷達(dá)的空間旋轉(zhuǎn)對(duì)稱進(jìn)動(dòng)目標(biāo)三維重構(gòu)[J]. 系統(tǒng)工程與電子技術(shù), 2016, 38(10):2281-2286. HU Xiaowei, TONG Ningning, WANG Jianye, et al. Three-Dimensional Reconstruction for Spatial Precession Targets with Rotational Symmetry in Radar Networks[J]. Systems Engineering and Electronics, 2016, 38(10):2281-2286.(in Chinese)

[7]ZHANG Shunsheng, Zhang Wei, ZONG Zhulin, et al. High-Resolution Bistatic ISAR Imaging Based on Two-Dimensional Compressed Sensing[J]. IEEE Trans on Antennas and Propagation, 2015, 63(5):2098-2111.

[8]李少東,楊軍,陳文峰,等. 基于壓縮感知理論的雷達(dá)成像技術(shù)與應(yīng)用研究進(jìn)展[J]. 電子與信息學(xué)報(bào), 2016, 38(2):495-508.

[9]SUN Chao, WANG Baoping, FANG Yang, et al. High-Resolution ISAR Imaging of Maneuvering Targets Based on Sparse Reconstruction[J]. Signal Processing, 2015, 108(1):535-548.

[10]葉淋美. 基于壓縮感知的雷達(dá)信號(hào)處理應(yīng)用研究[D]. 廈門(mén)：廈門(mén)大學(xué), 2014.

[11]LIU Hongchao, JIU Bo, LIU Hongwei, et al. A Novel ISAR Imaging Algorithm for Micromotion Targets Based on Multiple Sparse Bayesian Learning[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(10):1772-1776.

[12]洪靈,戴奉周,劉宏偉. 基于三維重構(gòu)的空間目標(biāo)進(jìn)動(dòng)參數(shù)估計(jì)方法[J]. 電波科學(xué)學(xué)報(bào), 2015, 30(2):237-243.

[13]王保平,方陽(yáng),孫超,等. 壓縮感知窄帶自旋目標(biāo)雷達(dá)成像[J]. 遙感學(xué)報(bào), 2015, 19(2):254-262. WANG Baoping, FANG Yang, SUN Chao, et al. Narrowband Radar Imaging for Spinning Targets Based on Compressed Sensing[J]. Journal of Remote Sensing, 2015, 19(2):254-262.(in Chinese)

[14]MOHIMANI H, BABAIE-ZADEH M, JUTTEN C. A Fast Approach for Overcomplete Sparse Decomposition Based on Smoothed0Norm[J]. IEEE Trans on Signal Processing, 2009, 57(1):289-301.