鄒建武，祝明波董 巍，李相平

(1.海軍航空工程學(xué)院，山東煙臺264001;2.海軍航空管制設(shè)備維修中心，北京100071)

1 引言

雷達(dá)方位超分辨就是在不改變雷達(dá)工作體制前提下，利用數(shù)字信號處理技術(shù)分辨同一雷達(dá)波束內(nèi)的幾個等距目標(biāo)，一直以來都是雷達(dá)探測領(lǐng)域研究的熱點(diǎn)。主要的方位超分辨方法有迭代反卷積法、RL算法、廣義逆濾波法和維納逆濾波算法等[1－9]。

Mark A.Richards等[1]提出了一種約束迭代反卷積法(CID)和一種快速迭代反卷積法 (FCID)，CID算法在特定情況下，方位分辨力改善了4倍，但算法收斂速度慢，計(jì)算量較大;FCID雖然提高了運(yùn)算速度，但是分辨性能下降。文獻(xiàn)[2－3]利用基于Bayesian準(zhǔn)則的RL算法實(shí)現(xiàn)了雷達(dá)方位超分辨，并對算法的超分辨性能進(jìn)行了深入研究，該算法在信噪比0 dB條件下，方位分辨力提高2倍。但是，上述兩種方法均是基于噪聲分布為泊松分布的情況下。丁義元等[4]提出利用廣義逆濾波方法提高實(shí)孔徑雷達(dá)角分辨力，在高信噪比的條件下，有效地改善了實(shí)孔徑雷達(dá)的角分辨力，但時域計(jì)算較為復(fù)雜。文獻(xiàn)[5－7]利用維納逆濾波算法對提高雷達(dá)方位角分辨力進(jìn)行了研究，在雷達(dá)回波信噪比大于30 dB的情況下，該方法可實(shí)現(xiàn)方位超分辨，證明了維納逆濾波算法實(shí)現(xiàn)方位超分辨是可行的，但不適用于低信噪比條件。

RL 算法由 Richardson[8]和 Lucy[9]先后獨(dú)立提出，主要用于射電天文學(xué)和醫(yī)學(xué)中圖像的復(fù)原，具有良好的頻譜外推能力，但其是基于噪聲分布為泊松分布的分辨算法。本文將RL算法應(yīng)用于噪聲服從高斯分布的雷達(dá)方位超分辨問題中，在噪聲存在的條件下，雖然RL算法也能分辨目標(biāo)，但噪聲使得分辨結(jié)果出現(xiàn)虛假目標(biāo)，且RL算法收斂速度慢。本文針對上述問題，將去噪和加速RL算法相結(jié)合，得到一種改進(jìn)RL算法，應(yīng)用于高斯噪聲條件下雷達(dá)超分辨問題中。

2 雷達(dá)方位回波信號統(tǒng)計(jì)觀測模型

為簡化分析，僅考慮方位角的變化。假設(shè)常規(guī)雷達(dá)在方位向上掃描，實(shí)孔徑雷達(dá)的方位回波信號表現(xiàn)為天線方向圖與目標(biāo)方位信息卷積，雷達(dá)目標(biāo)方位回波模型可以表示為

式中，x(θ)為目標(biāo)方位信息;h(θ)為雷達(dá)天線方向圖;n為加性噪聲，主要是接收機(jī)噪聲，本文主要考慮加性噪聲。將式(1)轉(zhuǎn)換為矩陣－向量表達(dá)式:

其中:

式(2)解釋為雷達(dá)方位回波信號統(tǒng)計(jì)觀測模型，當(dāng)n～pn(n)，假設(shè)y中每個元素的分布是相互獨(dú)立 ，則有

因此，雷達(dá)超分辨問題可以轉(zhuǎn)化為參數(shù)估計(jì)問題，常見的估計(jì)方法有極大似然估計(jì)，最大后驗(yàn)概率估計(jì)和矩陣估計(jì)方法。

3 RL算法及其改進(jìn)

3.1 RL 算法

RL算法是一種具有頻譜外推能力的迭代算法，假設(shè)噪聲為泊松分布。采用最大似然法進(jìn)行估計(jì)，通過Picard迭代算法得到如下迭代公式:

3.2 改進(jìn) RL算法

RL算法是一種最大化可能性數(shù)據(jù)逼近算法，該算法在噪聲影響可忽略或較小時，隨著k不斷增大，會依概率收斂于x;但在噪聲不可忽略，尤其是在低信噪比條件下，經(jīng)過多次迭代，恢復(fù)的目標(biāo)中會出現(xiàn)虛假目標(biāo)，在應(yīng)用RL算法過程中，為消除噪聲的影響，在出現(xiàn)虛假目標(biāo)或者算法結(jié)果變差以前，終止迭代[2]，因而RL算法難以獲得較好的超分辨效果，且RL算法收斂速度較慢。基于上述問題，本文對RL算法進(jìn)行了改進(jìn)。

當(dāng)存在噪聲時，將式(2)代入式(8)得

從上式可知，若存在噪聲，則上述的收斂性難以保證，存在噪聲被放大的缺陷，消除噪聲對RL算法很重要。由于雷達(dá)目標(biāo)回波信息主要集中分布在低頻段，噪聲的信息主要集中在高頻段，因此我們將高頻段的噪聲去除，來降低噪聲對RL算法的影響。

綜上所述，本文中的改進(jìn)RL算法，首先利用低通濾波器[11]對雷達(dá)目標(biāo)信號進(jìn)行去噪處理，同時鑒于RL算法收斂速度慢的缺點(diǎn)，我們利用加速RL算法[12]來對濾波后的信號進(jìn)行超分辨運(yùn)算，改進(jìn)RL算法的迭代公式

式中，accRL(·)為加速RL算法，Y(ω)為y的傅里葉變換，G(ω)為低通濾波傳遞函數(shù)，F(xiàn)－1為傅里葉逆變換。加速RL算法是一種矢量外推加速方法，其數(shù)學(xué)表達(dá)如下:

3.3 迭代次數(shù)選擇

RL算法迭代次數(shù)直接影響方位超分辨結(jié)果和計(jì)算時間。在低信噪比條件下，用RL算法得到的結(jié)果會先接近真實(shí)目標(biāo)，再偏離真實(shí)目標(biāo)。本文選用基于均方誤差(MSE)準(zhǔn)則的迭代次數(shù)選擇方法:

在不存在噪聲時，隨著迭代的進(jìn)行，MSE的值會變小，最后趨于平穩(wěn)，考慮到超分辨效果和計(jì)算時間，此時選擇使MSE值趨向平穩(wěn)時的最小迭代次數(shù)即可。在存在噪聲的條件下，隨著迭代的進(jìn)行，MSE的值會先變小，后變大，此時可選擇MSE值最小時的迭代次數(shù)作為算法的總迭代次數(shù)。

4 RL算法與改進(jìn)RL算法對兩點(diǎn)目標(biāo)的分辨性能仿真

4.1 無噪聲條件下RL算法仿真

假設(shè)脈沖重復(fù)頻率為1 kHz，掃描速度為100°/s，掃描范圍為－15°～15°，天線方向圖采用辛克型，半功率波束寬度為2.7°。在不存在噪聲條件下，用RL算法對間隔1/2半功率波束寬度的兩個等幅點(diǎn)目標(biāo)(分別位于－0.675°和 0.675°)進(jìn)行了分辨，分辨結(jié)果與相應(yīng)的均方誤差如圖1和圖2所示。

圖1 無噪聲條件下的分辨結(jié)果Fig.1 Resolution result without noise

圖2 均方誤差Fig.2 MSE

上述的仿真條件是針對小型實(shí)孔徑雷達(dá)平臺提出的，研究雷達(dá)對?；?qū)账阉髂繕?biāo)時的雷達(dá)目標(biāo)方位超分辨問題，此仿真不考慮真實(shí)場景的影響。圖1中的*表示點(diǎn)目標(biāo)的實(shí)際位置，由圖可看出，無噪聲時，利用RL算法可以分辨出間隔1/2半功率波束寬度的兩個點(diǎn)目標(biāo)，且隨著迭代次數(shù)的增加，MSE值趨于零。此時綜合考慮，我們選擇終止迭代次數(shù)為2 000。

4.2 噪聲存在條件下RL算法與改進(jìn)RL算法仿真

針對上述仿真參數(shù)，分別利用RL算法和改進(jìn)RL算法，在信噪比分別為30 dB和20 dB條件下對方位間隔分別為1/4和1/2半功率波束寬度的兩個等幅點(diǎn)目標(biāo)進(jìn)行了分辨，仿真結(jié)果如圖3～6所示。

圖3 30 dB條件下RL算法結(jié)果Fig.3 Results of RL algorithm with SNR=30 dB

圖4 30 dB條件下改進(jìn)RL算法結(jié)果Fig.4 Results of improved RL algorithm with SNR=30 dB

圖5 20 dB條件下RL算法結(jié)果Fig.5 Results of RL algorithm with SNR=20 dB

圖6 20 dB條件下改進(jìn)RL算法結(jié)果Fig.6 Results of improved RL algorithm with SNR=20 dB

由圖3和圖5可看出，在噪聲存在條件下，利用RL算法進(jìn)行分辨目標(biāo)，當(dāng)兩個目標(biāo)間隔為1/2半功率波束寬度時，噪聲對分辨結(jié)果沒造成太大影響，如圖3(b)和圖5(b)所示。但是，隨著兩目標(biāo)間隔的縮小，噪聲會使得分辨結(jié)果產(chǎn)生虛假目標(biāo)，如圖3(a)和圖5(a)所示。

由圖4(a)和圖6(a)可看出，在噪聲存在條件下，利用改進(jìn)RL算法能夠分辨兩個點(diǎn)目標(biāo)，而不產(chǎn)生虛假目標(biāo)。與圖3(a)和圖5(a)相比，改進(jìn)RL算法在一定程度上，可以降低噪聲對分辨結(jié)果的影響。由此可看出，改進(jìn)RL算法可用于高斯噪聲條件下的雷達(dá)方位超分辨問題的研究。

需要注意的是，由于RL算法受噪聲影響，產(chǎn)生虛假目標(biāo)，導(dǎo)致算法性能不穩(wěn)定，因此我們用維納逆濾波算法與改進(jìn)RL算法相比較，在不同信噪比條件下，兩種方法對兩個等幅點(diǎn)目標(biāo)所能分辨的最大分辨倍數(shù)(半功率波束寬度/改進(jìn)RL算法所能分辨的兩目標(biāo)間隔的最小角度)如表1所示。

表1 改進(jìn)RL和維納逆濾波算法的分辨性能Table1 Performance of improved RL algorithm and Wiener inverse filtering algorithm

由表1可知，隨著SNR值的降低，算法的分辨性能下降。與維納逆濾波相比，改進(jìn)RL算法可在低信噪比條件下分辨兩個點(diǎn)目標(biāo)，當(dāng)噪聲為0 dB時，雖然此時能分辨目標(biāo)，但是由于低頻段噪聲的影響，使得改進(jìn)RL算法的分辨精度變差，且有可能出現(xiàn)虛假目標(biāo)，如何降低低頻段噪聲的影響還有待于進(jìn)一步研究。

5 結(jié)束語

利用RL算法良好的頻譜外推能力對雷達(dá)方位超分辨問題進(jìn)行了研究，主要針對RL算法存在噪聲被放大的缺陷，提出一種改進(jìn)RL算法來降低噪聲的影響。初步仿真結(jié)果表明，改進(jìn)RL算法能降低噪聲的影響，可用于雷達(dá)方位超分辨，與維納逆濾波相比，噪聲適應(yīng)性更好，具有一定的實(shí)際應(yīng)用價值。但在低信噪比條件下，改進(jìn)RL算法的分辨性能下降，如何減小低信噪比條件下低頻段噪聲的影響尚待進(jìn)一步研究。

[1]Richards M A.Iterative Noncoherent Angular Superresolution[C]//Proceedings of 1988 IEEE National Radar Conference.Ann Arbor，MI:IEEE，1988:100－105.

[2]李惠.實(shí)孔徑雷達(dá)波束銳化算法研究[D].成都:電子科技大學(xué)，2010.LI Hui.Research on Beam Sharpening Algorithm of Real Aperture Radar[D].Chengdu:University of Electronic Science and Technology of China，2010.(in Chinese)

[3]Li D Y，Huang Y L，Yang J Y.Real Beam Radar Imaging Based on Adaptive Lucy－Richardson Algorithm[C]//Proceedings of 2011 IEEE Radar Conference.Chengdu:IEEE，2011:14－17.

[4]丁義元，楊建宇，張衛(wèi)華，等.改進(jìn)實(shí)孔徑雷達(dá)角分辨力的廣義逆濾波方法[J].電子學(xué)報(bào)，1993，21(9):15－19.DING Yi－yuan，YANG Jian－yu，ZHANG Wei－h(huán)ua，et al.Improvement of Angular Resolution of Real Aperture Radar Via Generalized Inverse Filtering[J].Acta Electronica Sinica，1993，21(9):15－19.(in Chinese)

[5]單榮光，李士國，朱力.去卷積實(shí)現(xiàn)雷達(dá)方位超分辨[J].現(xiàn)代雷達(dá)，1993(6):25－32.SHAN Rong－guang，LI Shi－guo，ZHU Li.Radar Azimuth Super－resolution Using Deconvolution[J].Modern Radar，1993(6):25－32.(in Chinese)

[6]Zhao K，Wang J G.Improved wiener filter superresolution algorithm for passive millimeter wave imaging[C]//Proceedings of 2011 IEEE Radar Conference.Chengdu:IEEE，2011:1768－1771.

[7]高明哲，祝明波.噪聲對維納濾波反卷積算法性能影響的分析[J].艦船電子工程，2012，32(12):35－37.GAO Ming－zhe，ZHU Ming－bo.Analysis on Noise's Impact on the Performance of Wiener Filter Deconvolution Algorithm[J].Ship Electronic Engineering，2012，32(12):35－37.(in Chinese)

[8]Richardson W H.Bayesian－based iterative method of image restoration[J].Journal of the Optical Society of A-merica，1972，62(1):55－59.

[9]Lucy L B.An iterative technique for the rectification of observed distributions[J].The Astronomical Journal，1974，79(6):745－754.

[10]周昌雄，蘇品剛，顏廷泰，等.應(yīng)用改進(jìn)自蛇模型和LR算法恢復(fù)毫米波圖像[J].激光與紅外，2012，42(4):468－472.ZHOU Chang－xiong，SU Pin－gang，YAN Ting－tai，et al.Millimeter wave image restoration based on the improved self－snake model and L－R algorithm[J].Laser＆ Infra-red，2012，42(4):468－472.(in Chinese)

[11]王志明，殷緒成.數(shù)字圖像處理與分析[M].北京:清華大學(xué)出版社，2012:72－73.WANG Zhi－ming，YIN Xu－cheng.Digital image processing and analysis[M].Beijing:Tsinghua University Press，2012:72－73.(in Chinese)

[12]David S，Biggs C，Andrews M.Acceleration of Iterative Image Restoration Algorithms [J].Applied Optics，1997，36(8):1766－1775.