李家強(qiáng)，姜慶剛，陳金立，葛俊祥

(1. 南京信息工程大學(xué) a. 氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心；　b. 電子與信息工程學(xué)院　南京 210044)

(2. 江蘇省氣象探測(cè)與信息處理重點(diǎn)實(shí)驗(yàn)室，　南京 210044)

?

船舶導(dǎo)航雷達(dá)的海雜波自適應(yīng)抑制算法

李家強(qiáng)1a,1b,2，姜慶剛1b，陳金立1a,1b,2，葛俊祥1a,1b,2

(1. 南京信息工程大學(xué) a. 氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心；b. 電子與信息工程學(xué)院南京 210044)

(2. 江蘇省氣象探測(cè)與信息處理重點(diǎn)實(shí)驗(yàn)室，南京 210044)

摘要：船舶導(dǎo)航雷達(dá)的海雜波抑制技術(shù)是提高海面目標(biāo)檢測(cè)性能的關(guān)鍵和難點(diǎn)技術(shù)，具有很高的軍事和民用價(jià)值。文中基于數(shù)理統(tǒng)計(jì)理論，提出了多模型選擇的海雜波自適應(yīng)抑制方法，通過(guò)對(duì)雷達(dá)多次掃描的回波數(shù)據(jù)建立海雜波模型庫(kù)，統(tǒng)計(jì)、估算各模型相應(yīng)參數(shù)選擇與雷達(dá)回波數(shù)據(jù)最佳匹配的模型，作抑制海雜波的對(duì)消處理。計(jì)算機(jī)仿真及實(shí)測(cè)數(shù)據(jù)處理，都驗(yàn)證該方法的有效性和可行性。

關(guān)鍵詞：船舶導(dǎo)航雷達(dá)；海雜波建模；海雜波抑制；目標(biāo)檢測(cè)

0引言

船舶導(dǎo)航雷達(dá)是船艇航海必備的導(dǎo)航設(shè)備之一，用于探測(cè)海面船只、礁石、冰山等多種物體,給船舶駕駛員提供近程和遠(yuǎn)程導(dǎo)航,以便及時(shí)發(fā)現(xiàn)、避讓海面各類障礙物,防止碰撞事故,提高船舶交通的安全性[1]。然而，船舶導(dǎo)航雷達(dá)接收到的回波信號(hào)中不僅包含目標(biāo)的回波信號(hào),也包含來(lái)自海面的各類動(dòng)態(tài)雜波，嚴(yán)重制約了海面目標(biāo)的可檢測(cè)性，尤其是海雜波的影響。因此，對(duì)海雜波的抑制是提高海面目標(biāo)檢測(cè)性能的關(guān)鍵和難點(diǎn)技術(shù)。

海雜波是雷達(dá)回波中最為復(fù)雜的一種雜波,不僅與天氣、風(fēng)速、海況等海面環(huán)境有關(guān)，還受雷達(dá)的入射角、發(fā)射頻率、分辨率、極化等工作狀態(tài)影響，國(guó)內(nèi)外很多學(xué)者對(duì)其進(jìn)行了深入研究。文獻(xiàn)[2-4]是基于混沌理論的海雜波分析方法，能夠在海雜波背景下有效地檢測(cè)到部分微弱信號(hào)，但目前這種理論還處于理論分析和實(shí)驗(yàn)仿真階段，理論比較復(fù)雜，且運(yùn)算量比較大，不易于工程應(yīng)用；文獻(xiàn)[5]是基于圖像處理角度的分析方式，運(yùn)用形態(tài)學(xué)濾波和圖像分割的原理來(lái)抑制海雜波,但這種算法適用范圍有限，只能是低分辨率船舶雷達(dá)，而且檢測(cè)的目標(biāo)通常是大型或中大型船只。針對(duì)上述問(wèn)題，本文提出一種基于數(shù)理統(tǒng)計(jì)理論的船舶導(dǎo)航雷達(dá)海雜波自適應(yīng)抑制方法，首先，建立海雜波模型庫(kù)；然后，根據(jù)雷達(dá)多次掃描回波數(shù)據(jù)統(tǒng)計(jì)模型庫(kù)中各模型的參數(shù)特征，得到各模型分布的概率密度分布曲線；其次，與雷達(dá)回波的分布曲線進(jìn)行對(duì)比，選擇與之最佳匹配的海雜波模型；最后，利用最佳匹配的海雜波模型對(duì)回波信號(hào)進(jìn)行對(duì)消，抑制海雜波的干擾。

1多參數(shù)選擇的海雜波模型庫(kù)

海雜波的幅度分布是建模中重要的一個(gè)內(nèi)容，經(jīng)國(guó)內(nèi)外學(xué)者深入研究與實(shí)驗(yàn)，先后建立多種海雜波幅度分布模型，典型的模型主要包括以下四種：瑞利分布(Rayleigh)、對(duì)數(shù)正態(tài)分布(LogNormal)、韋布爾分布(Weibull)和K分布[6]，現(xiàn)對(duì)四種模型作簡(jiǎn)要分析。

1.1瑞利分布

當(dāng)散射源的數(shù)目足夠多而且彼此相互獨(dú)立時(shí)，海雜波的包絡(luò)振幅服從瑞利分布。這種雜波模型通常描述低分辨率雷達(dá)(天線波束寬度大于2°，脈沖寬度大于1 μs)接收到的海雜波，其概率密度函數(shù)為

(1)

式中：x為海雜波幅度；σ為瑞利分布的參數(shù)；最大似然估計(jì)為[7]

(2)

式中：N為數(shù)據(jù)點(diǎn)的個(gè)數(shù)；xi為回波的包絡(luò)振幅。不同σ值下的瑞利分布曲線如圖1所示。

1.2對(duì)數(shù)正態(tài)分布

對(duì)數(shù)正態(tài)分布常用于表示平坦地區(qū)雷達(dá)分辨率相對(duì)較高、入射角相對(duì)較小(脈寬Φ<0.5 μs、角φ<5°)情況下的海雜波，也適用于描述海面形成“海尖峰”現(xiàn)象時(shí)的海雜波幅度。其概率密度函數(shù)為

x>0,σ>0，μ>0

(3)

式中：x為海雜波幅度；u為尺度參數(shù)；σ為形狀參數(shù)；最大似然估計(jì)為[7]

(4)

式中：N為數(shù)據(jù)點(diǎn)的個(gè)數(shù)；xi為回波的包絡(luò)振幅。不同σ、μ值下的對(duì)數(shù)正態(tài)分布曲線如圖2所示。

圖2　不同σ、μ值下的對(duì)數(shù)正態(tài)分布曲線

1.3韋布爾分布

韋布爾分布是一種比較常用的雜波模型，與瑞利分布、對(duì)數(shù)正態(tài)分布相比，有著更寬廣的動(dòng)態(tài)范圍和更高的精確度,常用于描述高分辨率雷達(dá)、低入射角和起伏較為均勻情況下的海雜波[8]。其概率密度函數(shù)

x≥0，p>0，q>0

(5)

式中：x為海雜波幅度；q為尺度參數(shù)；p為形狀參數(shù)；Menon[9]提出了一種參數(shù)有效估計(jì)方法

(6)

(7)

式中：N為數(shù)據(jù)點(diǎn)的個(gè)數(shù)；xi為回波的包絡(luò)振幅。不同p、q值下的韋布爾分布曲線如圖3所示。

圖3　不同p、q值下的韋布爾分布曲線

1.4

K分布

作為一種復(fù)合分布模型，K分布可用兩個(gè)相互獨(dú)立的斑點(diǎn)分量和基本幅度調(diào)制分量相乘描述[10]，其中，斑點(diǎn)分量又稱為快變化分量，服從瑞利分布，基本幅度調(diào)制分量又稱為慢變化分量，符合Gamma分布。K分布通常描述復(fù)雜海況和高分辨率雷達(dá)情況下非均勻的海雜波[11]，其概率密度函數(shù)為

x>0，a>0，v>0

(8)

式中：x為海雜波幅度；a為尺度參數(shù)；v為形狀參數(shù)；Γ(·)為Gamma函數(shù)；kv(·)為第二類修正Bessel函數(shù)。根據(jù)二階/四階矩估計(jì)法參數(shù)估算為

(9)

圖4　不同a、v值下的K分布曲線

2多模型選擇海雜波抑制算法原理

本文提出的一種基于數(shù)理統(tǒng)計(jì)理論的船舶導(dǎo)航雷達(dá)海雜波自適應(yīng)抑制方法，可以根據(jù)不同的海情，自適應(yīng)地選擇與之匹配最佳的雜波分布模型，通過(guò)對(duì)消處理達(dá)到抑制海雜波的目的，其關(guān)鍵在于分布曲線的擬合、雜波模型的選取。算法的流程圖如圖5所示。主要包括如下步驟。

圖5　海雜波抑制的流程圖

步驟1建立海雜波模型庫(kù)，本文選擇典型的四種分布模型(Rayleigh分布、LogNormal分布、Weibull分布和K分布)建立海雜波模型庫(kù)，并根據(jù)需要可做新的雜波模型添加；

步驟2對(duì)雷達(dá)回波中前N幀數(shù)據(jù)f1的幅度特性進(jìn)行分析，統(tǒng)計(jì)出模型庫(kù)中各模型的參數(shù)特征，從而得到各模型分布的概率密度分布曲線；

步驟3將得到的模型的分布曲線與雷達(dá)回波幅度分布直方圖進(jìn)行對(duì)比，選擇與雷達(dá)回波信號(hào)直方圖相差最小、匹配效果最佳的海雜波模型，作為雜波對(duì)消處理的背景；

步驟4將雷達(dá)回波數(shù)據(jù)f2與模型作對(duì)消處理，達(dá)到抑制海雜波的效果。

3仿真與實(shí)測(cè)處理結(jié)果與分析

3.1仿真實(shí)驗(yàn)

具體實(shí)施如下：

(1)通過(guò)仿真得到一組雷達(dá)回波數(shù)據(jù)，包含兩目標(biāo)，其多普勒頻移分別對(duì)應(yīng)100 Hz和150 Hz。回波信號(hào)的時(shí)域和頻域的波形如圖6a)、6b)所示，從圖中可以看出目標(biāo)信號(hào)淹沒(méi)在雜波中，很難辨別。

圖6　歸一化后的回波信號(hào)

(2)對(duì)雷達(dá)回波中前N幀數(shù)據(jù)f1的幅度特性進(jìn)行分析，分別統(tǒng)計(jì)出模型庫(kù)中各模型的參數(shù)特征，從而得到瑞利分布模型、對(duì)數(shù)正態(tài)分布模型、韋布爾分布模型、K分布模型的概率密度分布曲線。然后將得到的四種模型的分布曲線分別與雷達(dá)回波幅度分布直方圖進(jìn)行對(duì)比，選擇與回波信號(hào)直方圖相差最小、匹配效果最佳的海雜波模型。通過(guò)實(shí)驗(yàn)對(duì)比發(fā)現(xiàn)，在尺度參數(shù)a=0.6，形狀參數(shù)v=3時(shí)，K分布模型的概率密度分布曲線與雷達(dá)回波擬合效果最佳，結(jié)果如圖7所示，因此選擇模型庫(kù)中在尺度參數(shù)a=0.6，形狀參數(shù)v=3情況下的K分布模型作為算法對(duì)消處理的背景。

圖7　K分布概率密度分布曲線與回波信號(hào)的擬合圖

(3)將雷達(dá)回波中數(shù)據(jù)f2與K分布模型作對(duì)消處理，能夠有效地抑制海雜波。結(jié)果如圖8所示，能夠清楚地辨別目標(biāo)信號(hào)，即在頻率100 Hz和150 Hz處的兩個(gè)譜峰。

圖8　海雜波抑制結(jié)果

3.2實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)

本文采用的實(shí)測(cè)數(shù)據(jù)是來(lái)自本單位自主研發(fā)的X波段的船舶導(dǎo)航雷達(dá)實(shí)測(cè)回波數(shù)據(jù)，本實(shí)驗(yàn)將實(shí)測(cè)雷達(dá)數(shù)據(jù)通過(guò)MATLAB軟件處理和顯示。圖9為雷達(dá)回波的原始數(shù)據(jù)顯示，目標(biāo)淹沒(méi)在雜波中，圖10為采用本文方法抑制結(jié)果的結(jié)果圖，為了便于比較，給出了國(guó)外某型號(hào)船用導(dǎo)航雷達(dá)顯示結(jié)果，如圖11所示。通過(guò)圖9、圖10和圖11對(duì)比可見，本文提出的海雜波抑制方法能夠有效地抑制海雜波。

圖9　雷達(dá)回波的原始數(shù)據(jù)

圖10　海雜波抑制結(jié)果

圖11　國(guó)外某船舶導(dǎo)航雷達(dá)目標(biāo)回波圖

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

本文依據(jù)雷達(dá)雜波特性進(jìn)行統(tǒng)計(jì)分析，提出一種基于數(shù)理統(tǒng)計(jì)理論的多模型選擇海雜波自適應(yīng)抑制方法。通過(guò)對(duì)雷達(dá)回波數(shù)據(jù)的分析，從海雜波模型庫(kù)中選擇相適應(yīng)的雜波模型作為雜波對(duì)消背景。通過(guò)仿真分析，表明該方法能夠根據(jù)雷達(dá)回波信息，自適應(yīng)地選擇與之最佳匹配的雜波模型，對(duì)海雜波進(jìn)行有效抑制。最后通過(guò)實(shí)測(cè)雷達(dá)數(shù)據(jù)的處理與對(duì)比，也說(shuō)明了本文方法抑制雜波可行且效果明顯，從而證明此方法的有效性，對(duì)工程實(shí)踐具有一定的理論指導(dǎo)和參考價(jià)值。

參 考 文 獻(xiàn)

[1]姚元飛, 何奎, 錢延軍, 等. 降低船用連續(xù)波雷達(dá)虛警概率的海雜波抑制方法[J] . 電訊技術(shù), 2013, 53(5): 592-596.

YAO Yuanfei, HE Kui, QIAN Yanjun, et al. A sea clutter suppression method for reducing marine CW radar false alarm probability[J]. Telecommunication Engineering, 2013, 53(5) : 592-596.

[2]謝紅梅, 俞卞章. 基于神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)器的混沌海雜波弱信號(hào)檢測(cè)[J] . 現(xiàn)代雷達(dá), 2004, 26(9): 50-52.

XIE Hongmei, YU Bianzhang. Detecting weak signal in chaotic sea clutter using RBF neural network predictor[J]. Modern Radar, 2004, 26(9): 50-52.

[3]行鴻彥, 朱清清, 徐偉. 一種混沌海雜波背景下的微弱信號(hào)檢測(cè)方法[J]. 物理學(xué)報(bào), 2014, 63(10) : 1-7.

XING Hongyan, ZHU Qingqing, XU Wei. A method of weak target detection based on the sea clutter[J]. Acta Physica Sinica, 2014, 63(10) : 1-7.

[4]李正周, 陳靜, 沈美容, 等. 基于混沌神經(jīng)網(wǎng)絡(luò)的海上目標(biāo)圖像的海雜波抑制方法[J]. 光電子·激光, 2014, 25(3) : 588-594.

LI Zhengzhou, CHEN Jing, SHEN Meirong, et al. Sea clutter suppression approach for target images at sea based on chaotic neural network[J]. Journal of Optoeletronics·Laster, 2014, 25(3) : 588-594.

[5]於建生, 馬紅星, 朱征宇. 基于雷達(dá)回波相關(guān)特征的海雜波抑制方法[J]. 艦船電子工程, 2009, 29(6) : 121-123.

YU Jiansheng, MA Hongxing, ZHU Zhengyu. Correlation character of radar echo-based sea clutter suppression algorithm[J]. Ship Electronic Engineering, 2009, 29(6): 121-123.

[6]鄭琦, 陳棟, 王其沖. 艦載雷達(dá)復(fù)雜環(huán)境下的CFAR選擇[J]. 現(xiàn)代雷達(dá), 2009, 31(4): 15-18.

ZHENG Qi, CHEN Dong, WANG Qichong. CFAR selection for ship-borne radarin complex environment[J]. Modern Radar, 2009, 31(4): 15-18.

[7]鄧澤懷. 基于實(shí)測(cè)數(shù)據(jù)的海雜波統(tǒng)計(jì)建模[D] . 西安: 西安電子科技大學(xué), 2014.

DENG Zehuai. Statistical modeling of sea clutter based on measured data[D]. Xi′an: Xidian University, 2014.

[8]姚云萍, 段昶, 方慶. 一種船用雷達(dá)的海雜波抑制算法[J]. 火控雷達(dá)技術(shù), 2011, 40(3) : 32-36.

YAO Yunping, DUAN Chang, FANG Qing. A sea clutter suppression algorithm for marine radar[J]. Fire Control Radar Technology, 2011, 40(3) : 32-36.

[9]MENON M V. Estimation of the shape and scale parameters of the Weibull distribution[J]. Technometrics, 1963, 5(2): 175-182.

[10]WARD K D. Compound representation of high resolution sea clutter[J]. Electronics Letters, 1981, 17(22) : 561-563.

[11]夏棟, 察豪, 秦華, 等. 基于實(shí)采數(shù)據(jù)的海雜波與氣象雜波模擬[J]. 現(xiàn)代雷達(dá), 2012, 34(5): 77-80.

XIA Dong, CHA Hao, QIN Hua, et al. Simulation of sea clutter and weather clutter based on real data sampled[J]. Modern Radar, 2012, 34(5): 77-80.

李家強(qiáng)男，1976年生，博士，副教授，碩士生導(dǎo)師。研究方向?yàn)槔走_(dá)信號(hào)/數(shù)據(jù)處理。

姜慶剛男，1990年生，碩士研究生。研究方向?yàn)槔走_(dá)信號(hào)處理。

陳金立男，1982年生，博士。研究方向?yàn)槔走_(dá)信號(hào)/數(shù)據(jù)處理。

葛俊祥男，1960年生，教授，博士生導(dǎo)師。研究方向?yàn)槔走_(dá)系統(tǒng)、微波毫米波與天線技術(shù)、電磁散射與繞射理論。

A Sea Clutter Adaptive Suppression Algorithm forMarine Navigation Radar

LI Jiaqiang1a,1b,2，JIANG Qinggang1b，CHEN Jinli1a,1b,2，GE Junxiang1a,1b,2

(1a. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters;b. College of Electronic and Information Engineering,Nanjing University of Information Science and Technology,Nanjing 210044, China)(2. Jiangsu Key Laboratory of Meteorological Observation and Information Processing,Nanjing 210044, China)

Abstract：Sea clutter suppression technique based on marine navigation radar is the key and challenge for improvement of targets detection above the sea surface. It is of great value for military and civilian applications. Based on the theory of mathematical statistics, a sea clutter adaptive suppression algorithm via multi-model selection is proposed in this paper. Multiple scans of echo signals are utilized for establishment of the sea clutter model library. With statistics and estimation of corresponding parameters, the optimal matched model for echo data is applied to sea clutter suppression via cancellation. Results of simulations and experiments have demonstrated the validity and feasibility of the proposed algorithm.

Key words：the marine navigation radar; sea clutter modeling; sea clutter suppression; target detection

中圖分類號(hào)：TN957.51

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

文章編號(hào)：1004-7859(2016)02-0023-04

收稿日期：2015-10-20

修訂日期：2015-12-31

通信作者：李家強(qiáng)Email：lijiaqiang@sina.com

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目(61372066，61302188)；江蘇省自然科學(xué)基金資助項(xiàng)目(BK20131005)；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目(PAPD-II)

DOI：·信號(hào)處理· 10.16592/ j.cnki.1004-7859.2016.02.006