崔姍姍 周建江 朱劼昊

(南京航空航天大學(xué)電子信息工程學(xué)院 南京 210016)

1 引言

高分辨距離像(High-Resolution Range Profile,HRRP)是寬帶雷達(dá)信號獲取的目標(biāo)散射點子回波在雷達(dá)視線方向上投影的矢量和，它能反映目標(biāo)的尺寸、強散射點位置、精細(xì)結(jié)構(gòu)和材料質(zhì)地等信息。HRRP相對于2維成像，不僅易于獲得與處理，同時還避免了復(fù)雜的運動補償運算，對目標(biāo)識別具有重要意義[1－5]。

HRRP具有目標(biāo)姿態(tài)敏感性，在一定的目標(biāo)劃分角域范圍內(nèi)，其回波幅值分布通常具有一定的統(tǒng)計特性[3－7]。HRRP幀內(nèi)回波幅值統(tǒng)計模型主要分為兩種，一種認(rèn)為HRRP各距離分辨單元回波幅值分布具有近似獨立性[8,9]，另一種則進(jìn)一步考慮了HRRP各距離分辨單元回波幅值的相關(guān)統(tǒng)計特性[10,11]。由于HRRP各距離分辨單元回波幅值特征之間具有弱相關(guān)性，本文主要考慮第1種統(tǒng)計模型，即假設(shè)HRRP各距離分辨單元回波幅值分布近似獨立。

文獻(xiàn)[9]基于第1種統(tǒng)計模型，提出一種非參數(shù)化方法-基于累計量的隨機學(xué)習(xí)算法(Stochastic Learning of the Cumulative,SLC)[12]來估計概率密度函數(shù)。在樣本量充足時，該方法不僅可以準(zhǔn)確地估計概率密度函數(shù)，而且不存在“窗寬”調(diào)節(jié)問題。然而，當(dāng)訓(xùn)練樣本量不足時，基于SLC的非參數(shù)化方法不能準(zhǔn)確地進(jìn)行概率密度估計，最終導(dǎo)致目標(biāo)平均識別率下降。

文獻(xiàn)[13]提出一種基于半?yún)?shù)化概率密度估計的雷達(dá)目標(biāo)識別方法，通過有效利用雷達(dá)目標(biāo)HRRP各距離分辨單元回波幅值分布的經(jīng)驗知識，并結(jié)合非參數(shù)化概率密度估計的優(yōu)點來估計雷達(dá)目標(biāo)HRRP各距離分辨單元回波幅值的概率密度，達(dá)到參數(shù)化方法和非參數(shù)化方法優(yōu)缺互補的目的，提高雷達(dá)目標(biāo)HRRP的識別率。

為了解決SLC非參數(shù)化概率密度估計方法在小樣本條件下估計精度不高的問題，本文結(jié)合半?yún)?shù)化概念，提出一種基于半?yún)?shù)化SLC的雷達(dá)目標(biāo)識別方法。該方法將通過SLC求解得到的非參數(shù)化概率密度估計表示為非參數(shù)化修正因子，再利用非參數(shù)化修正因子對Gamma分布假設(shè)條件下的參數(shù)化概率密度估計進(jìn)行修正，提高小樣本情況下SLC算法的雷達(dá)目標(biāo)HRRP識別率。下文主要安排如下：第2節(jié)和第3節(jié)分別介紹SLC概率密度估計方法以及半?yún)?shù)化SLC概率密度估計；第4節(jié)利用5種飛機目標(biāo)HRRP仿真數(shù)據(jù)對基于半?yún)?shù)化SLC的雷達(dá)目標(biāo)識別方法進(jìn)行分析和驗證，并與其它基于統(tǒng)計建模的方法進(jìn)行了比較；第5節(jié)給出結(jié)論。

2 SLC非參數(shù)化概率密度估計

已知1維隨機變量x的N個相互獨立的隨機采樣xn∈R,n=1,…,N,按升序排列，記為x1≤x2≤…≤xN。設(shè)x的分布函數(shù) u=F(x)，則xn對應(yīng)的分布函數(shù)值un=F(xn)(n=1,2,…,N)滿足u1≤u2≤…≤uN。當(dāng)N充分大時，可以證明 u n(n=1,2,…,N)應(yīng)服從[0,1]上的均勻分布[12]。

將xn作為多層神經(jīng)網(wǎng)絡(luò)的輸入，隨機生成服從[0,1]上的均勻分布的N個點(n=1,2,…,N)，按升序排列，記為。比如，可取=n/ N(n=1,2,…,N)，將作為本次訓(xùn)練周期的期望輸出。期望輸出可以每次訓(xùn)練周期隨機生成，為了提高收斂速度，也可以每隔L個周期隨機生成一個期望輸出。

調(diào)整神經(jīng)網(wǎng)絡(luò)權(quán)值 w(t+1)=w(t)+η(t)? ε(w)/? w，增加循環(huán)次數(shù)，直至目標(biāo)函數(shù)ε(w)小于一個預(yù)先設(shè)定的閾值。其中，

網(wǎng)絡(luò)訓(xùn)練完畢，則輸入測試樣本x的網(wǎng)絡(luò)輸出H(w,x)即為其概率分布函數(shù)值。由于多層神經(jīng)網(wǎng)絡(luò)的激勵函數(shù)一般都是任意階可導(dǎo)函數(shù)，因此，概率密度函數(shù) f(x)=H'(w,x)。采用SLC估計概率密度，在樣本數(shù)充足的情況下，該算法不僅可以準(zhǔn)確地估計密度函數(shù)，而且不存在傳統(tǒng)的非參數(shù)化概率密度估計“窗寬”調(diào)節(jié)問題。但是隨著樣本數(shù)的減小，SLC算法表出概率密度函數(shù)的準(zhǔn)確度下降。對于HRRP，根據(jù)散射中心模型，單個距離單元回波幅值分布主要分為3種情況[8]：第1種情況是距離分辨單元內(nèi)沒有明顯的主散射中心，而只存在大量弱散射中心，其回波幅值服從瑞利分布；第2種情況是距離分辨單元內(nèi)同時存在一個主散射中心和大量弱散射中心，其回波幅值服從萊斯分布；第3種情況是距離分辨單元內(nèi)同時存在多個主散射中心(特別是包含2～3個主散射中心)

和大量弱散射中心，其回波幅值服從多峰分布，一般以雙峰分布為主。特別當(dāng)回波幅值呈多峰分布時，

SLC無法準(zhǔn)確描述出多峰情況。圖1給出了M2000

飛機模型在某方位角范圍內(nèi)兩種典型的HRRP距離單元幅值分布直方圖統(tǒng)計結(jié)果，實線表示SLC算法在樣本數(shù)充足情況下表出的概率密度估計值，虛線表示SLC算法在樣本數(shù)不足的情況下表出的概率密度估計值。由圖1可以看出，特別是對于多峰分布的情況，SLC算法在小樣本情況下表出的準(zhǔn)確度大大降低。

3 半?yún)?shù)化 SLC概率密度估計及 HRRP識別

常用概率密度估計方法大致可以分為兩類：一類是參數(shù)化概率密度估計，如高斯模型或 Gamma模型等，利用訓(xùn)練樣本對概率密度模型的參數(shù)進(jìn)行估計。參數(shù)化概率密度估計的優(yōu)點是在樣本數(shù)不足或者樣本特征維數(shù)較高時，也可以較為準(zhǔn)確地估計概率密度，其缺點是可能面臨“模型失配”問題，影響識別效果。另一類為非參數(shù)化概率密度估計。非參數(shù)化方法可以估計任意形式的概率分布，但是需要大量的訓(xùn)練要本，計算復(fù)雜度也較高。

半?yún)?shù)化概率密度估計[13－15]較好地結(jié)合了上述兩類方法的優(yōu)點，它在數(shù)據(jù)分布具有一定先驗知識的基礎(chǔ)上，通過非參數(shù)化方法對參數(shù)化概率密度估計進(jìn)行修正，達(dá)到提高概率密度估計精度的目的。

圖1 SLC算法概率密度估計結(jié)果

為了提高在小樣本情況下的平均識別率，結(jié)合半?yún)?shù)化概率密度估計概念，本文提出一種半?yún)?shù)化SLC概率密度估計方法，用于HRRP識別。該方法將通過 SLC求解得到的非參數(shù)化概率密度估計表示為非參數(shù)化修正因子，再利用非參數(shù)化修正因子對Gamma分布假設(shè)條件下的參數(shù)化概率密度估計進(jìn)行修正，達(dá)到參數(shù)化方法和非參數(shù)化方法優(yōu)缺互補的目的。最終提高在小樣本情況下的目標(biāo)識別率。

本文采用的非參數(shù)化修正因子的形式[14]為

上式中，P(x,r)表示樣本x處于以自身為球心，半徑為r的球體B中的概率，即

其中，P1為SLC算法對隨機變量x的概率估計值，P0為參數(shù)化方法對隨機變量x的概率估計值。因此，與式(1)對應(yīng)的半?yún)?shù)化SLC概率密度估計函數(shù)可以表示為

在非參數(shù)化修正因子式(1)中，球半徑r實際上充當(dāng)了光滑因子的作用。當(dāng)r增加時，非參數(shù)化修正因子曲線變得更加光滑。當(dāng)球半徑r趨近于0時，式(1)趨近于非參數(shù)化 SLC概率密度估計值和參數(shù)化概率密度估計值之比，此時，式(3)退化為非參數(shù)化SLC概率密度估計。而當(dāng)球半徑r增加到包含整個樣本空間時，P0,P1以及式(1)剛好等于1，此時，式(3)又退化為參數(shù)化概率密度估計。因此，在實際使用半?yún)?shù)化SLC概率密度估計時，需要對球半徑r進(jìn)行設(shè)定，以滿足最佳概率密度估計需求。

圖2給出了M2000飛機模型在以0°俯仰角、60°方位角為中心的姿態(tài)下，15°角域劃分內(nèi)兩種典型的HRRP距離分辨單元回波幅值分布的直方圖、Gamma概率密度模型和半?yún)?shù)化SLC概率密度模型統(tǒng)計結(jié)果。由圖 2可見，在回波幅值呈單峰分布的情況下，參數(shù)化方法和半?yún)?shù)化 SLC方法都能比較好地估計其概率密度；但在回波幅值呈雙峰分布的情況下，參數(shù)化方法出現(xiàn)了“模型失配”問題，半?yún)?shù)化 SLC方法仍能比較準(zhǔn)確地估計其分布情況。

圖2 M2000飛機模型HRRP參數(shù)化與半?yún)?shù)化SLC概率密度估計

本節(jié)最后對基于半?yún)?shù)化 SLC的雷達(dá)目標(biāo)HRRP識別過程進(jìn)行小結(jié)，該方法主要包括訓(xùn)練和測試兩個階段。

訓(xùn)練階段：

(1) 將待識別目標(biāo)的訓(xùn)練樣本按角域分幀，各幀內(nèi)的HRRP對齊并幅度歸一化；

(2) 按式(3)統(tǒng)計每幀 HRRP 各距離分辨單元回波幅值的半?yún)?shù)化概率密度估計。由于雷達(dá)目標(biāo)HRRP各距離分辨單元回波幅值分布具有近似獨立性，因此，可以認(rèn)為每幀HRRP的概率密度就是該幀HRRP各距離分辨單元回波幅值概率密度的乘積。

識別階段：

(1) 假設(shè)各種飛機目標(biāo)先驗概率均相等，計算測試樣本在每幀HRRP半?yún)?shù)化概率密度估計模型下的概率密度值；

(2) 根據(jù)貝葉斯決策規(guī)則把測試樣本判別為最大概率密度值對應(yīng)的目標(biāo)。

4 實驗與分析

4.1 實驗數(shù)據(jù)及其預(yù)處理

實驗選用南京航空航天大學(xué)目標(biāo)特性研究中心提供的Su27,F16,M2000,J8II和J6等5種飛機全方位角轉(zhuǎn)臺仿真數(shù)據(jù)。雷達(dá)發(fā)射帶寬500 MHz，并且在計算HRRP之前，分別對所有雷達(dá)I/Q兩路通道數(shù)據(jù)添加高斯白噪聲，主要分析信噪比為 20 dB條件下的HRRP數(shù)據(jù)。由于實驗采用轉(zhuǎn)臺數(shù)據(jù)，因此不存在HRRP時移敏感性問題，針對幅度敏感性，所有HRRP均已實現(xiàn)2-范數(shù)歸一化處理。

4.2 雷達(dá)目標(biāo)HRRP識別結(jié)果及分析

實驗主要比較分析了半?yún)?shù)化SLC算法應(yīng)用于所選HRRP實驗數(shù)據(jù)的識別性能，尤其在小樣本情況下。在仿真實驗中，SLC所對應(yīng)的神經(jīng)網(wǎng)絡(luò)只有一個隱層，包含3個神經(jīng)元。其中輸入層到隱層的激勵函數(shù)為tanh(x)，隱層到輸出層的激勵函數(shù)為logsig(x),η=0.0001,λ=1,Δ=0.001。

表1 是樣本數(shù)充足時3種方法對5種飛機目標(biāo)的詳細(xì)識別結(jié)果。其中，參數(shù)化概率密度估計選用Gamma模型。由表1可見，當(dāng)訓(xùn)練樣本數(shù)充足時，非參數(shù)化SLC方法直接從訓(xùn)練數(shù)據(jù)本身出發(fā)，估計其分布情況，具有良好的靈活性，識別率優(yōu)于參數(shù)化方法。半?yún)?shù)化SLC方法的平均識別效果介于參數(shù)化與非參數(shù)化SLC方法之間，而且在Su27和J6飛機目標(biāo)的識別率上高于非參數(shù)化方法。

表1 信噪比為20 dB條件下訓(xùn)練樣本數(shù)充足時3種方法的目標(biāo)識別率(%)

但是當(dāng)訓(xùn)練樣本數(shù)減小，非參數(shù)化SLC方法的優(yōu)勢便逐漸下降，實驗中分別選取訓(xùn)練樣本數(shù)為60、55、50、45、40、35、30、25、20、15，圖 3 是 3種方法的平均識別率隨著樣本量減少時的變化。由圖3可以看到，在小樣本情況下，非參數(shù)化SLC方法的平均識別率下降迅速，低于參數(shù)化方法。半?yún)?shù)化SLC方法有效利用了 HRRP 各距離單元幅值分布的經(jīng)驗知識，結(jié)合了參數(shù)化方法的優(yōu)點，因而顯示出良好的識別性能，在3種方法中平均識別率處于最高。

表2為訓(xùn)練樣本數(shù)為25時，3種方法的詳細(xì)識別效果。雖然從5種飛機目標(biāo)個體的識別結(jié)果來看，半?yún)?shù)化SLC方法未必是最優(yōu)的，但是就平均識別率來講，半?yún)?shù)化SLC方法高于非參數(shù)化SLC方法。通過結(jié)合半?yún)?shù)化概念，提高了SLC算法在小樣本雷達(dá)目標(biāo)識別中的識別效果。在利用半?yún)?shù)化SLC概率密度估計時，面臨最優(yōu)光滑因子r的選擇。圖4為訓(xùn)練樣本數(shù)為25時，半?yún)?shù)化SLC方法在不同的光滑因子作用下所取得的不同的平均識別率。由圖4可以很清楚地看到，當(dāng)光滑因子取值為1.2時，識別效果最好。當(dāng)光滑因子大于1.2時，識別率迅速下降。

圖3 小樣本情況下3種方法隨樣本數(shù)變化的識別效果

圖4 半?yún)?shù)化SLC方法識別率隨光滑因子的變化

表2 信噪比為20 dB條件下訓(xùn)練樣本數(shù)不足時3種方法的目標(biāo)識別率(%)

5 結(jié)束語

非參數(shù)化SLC方法直接從訓(xùn)練數(shù)據(jù)本身出發(fā)估計其分布情況，不存在“窗寬”調(diào)節(jié)問題，具有良好的靈活性。但是在訓(xùn)練樣本不充足的情況下，其概率密度估計的準(zhǔn)確度大大降低，導(dǎo)致將其應(yīng)用于雷達(dá)目標(biāo)識別時識別率下降。為此，本文提出將SLC算法結(jié)合半?yún)?shù)化概念，通過SLC求解得到的非參數(shù)化修正因子對Gamma分布假設(shè)條件下的參數(shù)化概率密度估計進(jìn)行修正，從而提高小樣本情況下雷達(dá)目標(biāo)的識別率?；?5 種飛機目標(biāo) HRRP 的仿真實驗表明，在小樣本情況下，半?yún)?shù)化SLC概率密度估計下的雷達(dá)目標(biāo)識別方法相比非參數(shù)化SLC方法具有更好的識別效果。

[1]王鵬輝,劉宏偉,杜蘭,等.基于線性動態(tài)模型的雷達(dá)高分辨距離像小樣本目標(biāo)識別方法[J].電子與信息學(xué)報,2012,34(2):305-311.Wang Peng-hui,Liu Hong-wei,Du Lan,et al..Linear dynamic model based radar HRRP target recognition under small training set conditions[J].Journal of Electronics &Information Technology,2012,34(2): 305-311.

[2]Shi Lei,Wang Peng-hui,Liu Hong-wei,et al..Radar HRRP statistical recognition with local factor analysis by automatic Bayesian Ying-Yang harmony learning[J].IEEE Transactions on Signal Processing,2011,59(2): 610-617.

[3]Du L,Wang P,Liu H,et al..Bayesian spatiotemporal multitask learning for radar HRRP target recognition[J].IEEE Transactions on Signal Processing,2011,59(7):3182-3196.

[4]Zheng Tang-hong,Wang Yang,and Li Shi-guo.An HRRP preprocessing method and its application in radar targetrecognition[C].IEEE CIE International Conference on Radar,Chengdu,China,Oct.24-27,2011,1: 646-649.

[5]Wang P,Dai F,Pan M,et al..Radar HRRP target recognition in frequency domain based on autoregressive model[C].IEEE Radar Conference,Kansas City,USA,2011:714-717.

[6]Webb A R.Gamma mixture models for target recognition[J].Pattern Recognition,2000,33(12): 2045-2054.

[7]Copsey K and Webb A.Bayesian gamma mixture model approach to radar target recognition[J].IEEE Transactions on Aerospace and Electronic Systems,2003,39(4):1201-1217.

[8]Du Lan,Liu Hong-wei,Bao Zheng,et al..A two-distribution compounded statistical model for radar HRRP target recognition[J].IEEE Transactions on Signal Processing,2006,54(6): 2226-2238.

[9]趙峰,張軍英,劉敬,等.基于非參數(shù)化概率密度估計的雷達(dá)目標(biāo)識別[J].電子與信息學(xué)報,2008,30(7): 1740-1743.Zhao Feng,Zhang Jun-ying,Liu Jing,et al..Radar target recognition based on nonparmetric density estimation[J].Journal of Electronics & Information Technology,2008,30(7):1740-1743.

[10]Du Lan,Liu Hong-wei,and Bao Zheng.Radar HRRP statistical recognition based on hypersphere model[J].Signal Processing,2008,88(5): 1176-1190.

[11]Du Lan,Liu Hong-wei,and Bao Zheng.Radar HRRP statistical recognition: parametric model and model selection[J].IEEE Transactions on Signal Processing,2008,56(5): 1931-1944.

[12]Malik Magdon-Ismail and Amir Atiya.Density estimation and random variate generation using multilayer networks[J].IEEE Transactions on Neural Networks,2002,13(3):497-520.

[13]朱劼昊,周建江,吳杰,等.基于半?yún)?shù)化概率密度估計的雷達(dá)目標(biāo)識別[J].電子與信息學(xué)報,2010,32(9): 2161-2166.Zhu Jie-hao,Zhou Jian-jiang,Wu Jie,et al..Radar target recognition based on semi-parametric density-estimation[J].Journal of Electronics & Information Technology,2010,32(9):2161-2166.

[14]Chaudhuri P,Ghosh A K,and Oja H.Classification based on hybridization of parametric and nonparametric classifiers[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2009,31(7): 1153-1164.

[15]Hjort N L and Glad I K.Nonparametric density estimation with a parametric start[J].Annals of Statistics,1995,23(3):882-904.