本發(fā)明涉及一種部分均勻海雜波背景下的平滑廣義似然比檢測(cè)方法，屬于雷達(dá)目標(biāo)檢測(cè)技術(shù)領(lǐng)域。

背景技術(shù)：
在海面目標(biāo)檢測(cè)中，采用匹配于海雜波統(tǒng)計(jì)及相關(guān)特性的自適應(yīng)目標(biāo)檢測(cè)算法是一種普遍采用的技術(shù)手段。因而，待檢測(cè)單元雜波的特性與自適應(yīng)檢測(cè)器的設(shè)計(jì)、檢測(cè)性能密切相關(guān)。隨著雷達(dá)距離分辨率的進(jìn)一步提高，雷達(dá)接收回波由先前的均勻雜波變?yōu)榉蔷鶆螂s波環(huán)境。部分均勻的海雜波具有瞬時(shí)功率波動(dòng)較大的特點(diǎn)，這種較大的頻率波動(dòng)會(huì)直接影響著檢測(cè)器的性能。目前，海雜波背景下的自適應(yīng)檢測(cè)器設(shè)計(jì)往往為了簡(jiǎn)化計(jì)算，減少檢測(cè)器的復(fù)雜度，而假設(shè)所處理的雷達(dá)接收回波為均勻雜波特性，又或是為了應(yīng)對(duì)部分均勻海雜波，而引入了較為復(fù)雜的信號(hào)處理算法。例如，北京環(huán)境特性研究所申請(qǐng)的發(fā)明專利：海雜波抑制以及海雜波背景中目標(biāo)檢測(cè)的方法和系統(tǒng)(專利申請(qǐng)?zhí)枺篊N201310556638.9，公開號(hào)：CN103645467A)。該專利申請(qǐng)從實(shí)測(cè)時(shí)空色散關(guān)系中提取出實(shí)測(cè)速度項(xiàng)參數(shù)，確定出基于所述本征速度項(xiàng)參數(shù)的海雜波的本征時(shí)空色散關(guān)系，進(jìn)而重構(gòu)得到估計(jì)的海雜波一維距離像的歷程圖。最終通過(guò)實(shí)測(cè)所得歷程圖與估計(jì)的歷程圖相減來(lái)獲得抑制海雜波的一維距離像的歷程圖的圖像數(shù)據(jù)，從而達(dá)到消除多普勒頻移的目的。該專利申請(qǐng)的不足之處在于：在考慮消除雷達(dá)移動(dòng)所導(dǎo)致的多普勒頻移這一事實(shí)時(shí)，忽略了部分均勻海雜波這另一客觀事實(shí)。所以，其最終得到抑制海雜波的圖像數(shù)據(jù)引入了非均勻海雜波的干擾特性。又例如，西安電子科技大學(xué)申請(qǐng)的發(fā)明專利：海雜波背景下子帶自適應(yīng)GLRT-LTD檢測(cè)方法(專利申請(qǐng)?zhí)枺篊N201510030360.0，公開號(hào)：CN104569948A)，該專利申請(qǐng)通過(guò)構(gòu)造子帶濾波器的方式來(lái)應(yīng)對(duì)部分均勻海雜波，從而實(shí)現(xiàn)準(zhǔn)確的檢測(cè)判決，提高了檢測(cè)性能。但該專利的主要不足是：引入的子代濾波器過(guò)于復(fù)雜，計(jì)算量較大。這勢(shì)必會(huì)影響檢測(cè)器的信號(hào)處理速度。而本發(fā)明能夠很好地解決上面的問(wèn)題。

技術(shù)實(shí)現(xiàn)要素：
本發(fā)明目的在于解決了上述現(xiàn)有技術(shù)的不足，提出一種部分均勻海雜波背景下的平滑廣義似然比檢測(cè)方法，該方法在不增加GLRT檢測(cè)器計(jì)算復(fù)雜度的前提下，提升了檢測(cè)器的性能。本發(fā)明解決其技術(shù)問(wèn)題所采取的技術(shù)方案是：一種部分均勻海雜波背景下的平滑廣義似然比檢測(cè)方法，該方法在不增加計(jì)算復(fù)雜度的前提下，能在實(shí)測(cè)海雜波數(shù)據(jù)實(shí)驗(yàn)中獲得更好的檢測(cè)性能。方法流程：步驟1：采用GLRT檢測(cè)器作為S-GLRT檢測(cè)器的數(shù)學(xué)原型；所述的GLRT檢測(cè)器數(shù)學(xué)表達(dá)式為：其中M表示雜波的協(xié)方差矩陣，p為多普勒導(dǎo)向矢量，z表示雷達(dá)接收待檢測(cè)單元的回波，H表示共軛轉(zhuǎn)置，β為尺度參數(shù)，ξ為判決門限。步驟2：以歸一化樣本協(xié)方差矩陣(normalizedsamplecovariancematrix,NSCM)估計(jì)或漸進(jìn)最大似然(approximatedmaximumlikelihood,AML)估計(jì)作為雜波協(xié)方差矩陣M的估計(jì)算法，以取中值估計(jì)作為平滑因子的估計(jì)算法；所述的雜波協(xié)方差矩陣M的歸一化樣本協(xié)方差矩陣(normalizedsamplecovariancematrix,NSCM)估計(jì)形式為：漸進(jìn)最大似然(approximatedmaximumlikelihood,AML)估計(jì)形式為：相應(yīng)的，取中值估計(jì)的平滑因子表達(dá)式為：步驟3，將GLRT檢測(cè)器數(shù)學(xué)表達(dá)式中的尺度參數(shù)β替換為得到GLRT檢測(cè)器的修正形式，即平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器的表達(dá)式；所述的平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器的表達(dá)式為：有益效果：本發(fā)明與現(xiàn)有技術(shù)比較具有以下優(yōu)點(diǎn)：(1)本發(fā)明提出的平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器與GLRT檢測(cè)器相比，在不增加計(jì)算復(fù)雜度的前提下，能在實(shí)測(cè)海雜波數(shù)據(jù)實(shí)驗(yàn)中獲得更好的檢測(cè)性能。(2)本發(fā)明提出的S-GLRT檢測(cè)器，其引入的平滑因子主要是為了削弱部分均勻海雜波對(duì)檢測(cè)器性能的影響。但不失通用性，對(duì)于均勻海雜波背景下的目標(biāo)檢測(cè)，S-GLRT檢測(cè)器仍具有與GLRT檢測(cè)器相近的檢測(cè)性能。符合實(shí)際的雜波環(huán)境要求。(3)本發(fā)明提出的S-GLRT對(duì)尺度參數(shù)具有恒虛警特性。(4)S-GLRT檢測(cè)器中的平滑因子采用取中值估計(jì)算法，在實(shí)際環(huán)境中具有較好的性能表現(xiàn)。附圖說(shuō)明圖1為本發(fā)明的方法流程圖。圖2為本發(fā)明提出的S-GLRT與GLRT在實(shí)測(cè)雜波情況下的性能比較圖。具體實(shí)施方式下面結(jié)合說(shuō)明書附圖對(duì)本發(fā)明創(chuàng)造作進(jìn)一步的詳細(xì)說(shuō)明。本發(fā)明在部分均勻海雜波背景下，提升了GLRT檢測(cè)器性能的方法，其中的主要技術(shù)問(wèn)題包括：(1)平滑因子估計(jì)算法的選擇。(2)S-GLRT檢測(cè)器數(shù)學(xué)表達(dá)式的導(dǎo)出。本發(fā)明所述的部分均勻海雜波中雷達(dá)目標(biāo)的平滑自適應(yīng)檢測(cè)算法包括以下技術(shù)措施：首先，給出GLRT檢測(cè)器的數(shù)學(xué)模型。然后，對(duì)雜波協(xié)方差矩陣M分別采用NSCM和AML估計(jì)算法，平滑因子采用取中值估計(jì)算法。最后，將GLRT檢測(cè)器數(shù)學(xué)表達(dá)式中的尺度參數(shù)β替換為得到平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器的數(shù)學(xué)模型。如圖1所示，本發(fā)明提供了一種部分均勻海雜波背景下的平滑廣義似然比檢測(cè)方法，該方法包括：步驟1：首先采用GLRT檢測(cè)器的數(shù)學(xué)表達(dá)式作為數(shù)學(xué)原型：公式(1)中M表示雜波的協(xié)方差矩陣，p為多普勒導(dǎo)向矢量，z表示雷達(dá)接收待檢測(cè)單元的回波，H表示共軛轉(zhuǎn)置，β為尺度參數(shù)，ξ為判決門限。步驟2：當(dāng)以歸一化樣本協(xié)方差矩陣(normalizedsamplecovariancematrix,NSCM)估計(jì)作為雜波協(xié)方差矩陣M的估計(jì)算法，以取中值估計(jì)算法作為的取值算法時(shí)，的表達(dá)式為公式(2)中的k表示的是雷達(dá)回波的樣本個(gè)數(shù)。取中值估計(jì)算法median表示的具體含義是，對(duì)花括號(hào)中的所有數(shù)值比較大小，取大小居中的那個(gè)數(shù)值作為取值。當(dāng)以漸進(jìn)最大似然(approximatedmaximumlikelihood,AML)估計(jì)作為雜波協(xié)方差矩陣M的估計(jì)算法，以取中值估計(jì)算法作為的取值算法時(shí)，的表達(dá)式為步驟3：針對(duì)公式(1)，將尺度參數(shù)β替換為相應(yīng)的GLRT檢測(cè)器修改為：公式(4)即為本發(fā)明提出的平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器。本發(fā)明提出的平滑GLRT(smoothGLRT,S-GLRT)檢測(cè)器的可以通過(guò)下面的實(shí)驗(yàn)進(jìn)一步驗(yàn)證。實(shí)驗(yàn)使用IPIX雷達(dá)采集的海雜波數(shù)據(jù)來(lái)分析S-GLRT的檢測(cè)性能，提供數(shù)據(jù)的網(wǎng)址：http://soma.mcmaster.ca/ipix.php，數(shù)據(jù)名為：19980223-170435(距離分辨率為15m)，HH極化，該數(shù)據(jù)共含有60000個(gè)時(shí)間脈沖，34個(gè)距離單元?？紤]到部分距離單元數(shù)據(jù)可能被污染，故而發(fā)明人選取了26個(gè)純海雜波單元的數(shù)據(jù)，目標(biāo)被加在第15個(gè)距離單元上。圖2是本發(fā)明提出的S-GLRT和傳統(tǒng)GLRT在不同協(xié)方差矩陣估計(jì)下的檢測(cè)性能比較。顯然，不論是NSCM估計(jì)器還是AML估計(jì)器，在實(shí)測(cè)雜波中，S-GLRT的檢測(cè)性能明顯優(yōu)于GLRT的檢測(cè)性能。