[發明專利]一種基于包絡直方圖分布的雷達窄帶干擾識別方法

在審

• 申請號: 202210083162.0
• 申請日: 2022-01-25
• 公開（公告）號: CN114488031A
• 公開（公告）日: 2022-05-13
• 發明（設計）人: 王志剛;李穎;鄒成曉;耿利祥;楊志昆
• 申請（專利權）人: 中國船舶重工集團公司第七二四研究所
• 主分類號: G01S7/36
• 分類號: G01S7/36
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于包絡直方圖分布的雷達窄帶干擾識別方法，其特征在于：

S1：回波幅度包絡構建；

S2：等幅直方圖分布統計；以回波幅度包絡對數化結果最大值和最小值為邊界、設定等幅閾值為步長，將對數化結果映射至直方圖格子；

S3：臨近窗精細化補償；對每個直方圖格子進行后向補償，得到補償結果：

Z(i)＝S(i)+S(i+1),i∈[1,N_unit-1] (1)

其中，Z(i)表示直方圖臨近窗精細化補償結果；

S4：自適應分類識別；設定直方圖檢測門限T_jam：

T_jam＝mean(y_norm)+3×std(y_norm) (2)

其中，y_norm表示雷達未受干擾時回波求模對數化結果，mean(·)表示求均值，std(·)表示求方差；與檢測門限T_jam對應的直方圖格子索引N_Tj：

N_Tj＝floor((T_jam-y_min)/T_unit)+1 (3)

對直方圖索引i∈[N_Tj+1,N_unit-1]的補償后數據Z(i)進行分類識別：

S41：遍歷所有格子，如果存在Z(i)M/2，則為假目標干擾；否則跳轉S42；

S42：統計格子中Z(i)2的格子個數N_min2、所有格子中元素的和N_sum＝sum(Z(i))，其中，sum(·)表示求和；如果N_min2(N_unit-N_Tj)/2且N_sumM/2，則為窄脈沖干擾；否則跳轉S43；

S43：計算補償后數據Z(i)的峭度K_z，如果K_z小于設定閾值T_z，則為噪聲類干擾；否則為卷積類干擾。

2.根據權利要求1所述的一種基于包絡直方圖分布的雷達窄帶干擾識別方法，其特征在于：所述步驟S1中的回波幅度包絡構建包括：雷達原始零中頻回波數據x取模后再進行對數化處理：

y＝A×logB(|x|) (4)

其中，y表示回波幅度包絡對數化結果，A表示幅度系數，logB(·)表示以B為底數求對數，|·|表示取向量元素的模值。

3.根據權利要求1所述的一種基于包絡直方圖分布的雷達窄帶干擾識別方法，其特征在于：所述步驟S2的等幅直方圖分布統計，求解回波幅度包絡對數化結果y的最大值y_max和最小值y_min；設定等幅直方圖的等幅步長為T_unit，將直方圖劃分為N_unit個格子：

N_unit＝floor((y_max-y_min)/T_unit)+1 (5)

其中，floor(·)表示向負無窮大方向取整；

以y_min為直方圖格子起始、T_unit為直方圖格子步長，遍歷回波幅度包絡對數化結果y中每個元素，并映射至直方圖格子中：

n＝floor((y_m-y_min)/T_unit)+1,m∈[1,M] (6)

S(n)＝S(n)+1,n∈[1,N_unit] (7)

其中，M表示回波幅度包絡對數化結果y中元素總個數，n表示等幅直方圖格子索引，S(n)表示第n個直方圖格子中元素個數。