[發(fā)明專利]一種基于干擾規(guī)避的信噪比解算方法

在審

• 申請?zhí)枺?/td>: 202210020163.0
• 申請日: 2022-01-10
• 公開（公告）號: CN114374448A
• 公開（公告）日: 2022-04-19
• 發(fā)明（設計）人: 齊志剛;李強;慈國輝;張紅濤;申昱東;林興隆
• 申請（專利權）人: 中國電子科技集團公司第五十四研究所
• 主分類號: H04B17/336
• 分類號: H04B17/336
• 代理公司: 河北東尚律師事務所 13124
• 代理人: 王文慶;曲佳穎

1.一種基于干擾規(guī)避的信噪比解算方法，其特征在于，包括以下步驟：

步驟1：將模擬信號轉換成數字信號并進行下變頻，然后根據跟蹤信號最低信噪比要求選擇采樣數據，進行快速傅立葉變化形成頻譜數據；

步驟2：根據監(jiān)測帶寬選擇參加信噪比計算的數據，將頻譜數據進行拼接得到新的數據序列；

步驟3：在新的數據序列中找到最大功率值p_max和最小功率值p_min，并計算Δ_max＝10lgp_max-10lg p_min；

步驟4：如果Δ_max大于最大閾值，則梯度h_dB為1dB，如果Δ_max在最大閾值與最小閾值之間，則h_dB＝(10lg p_max-10lg p_min)/64，然后進入步驟6，如果Δ_max小于最小閾值,則直接設置信噪比為0；

步驟5：將功率值劃分為t個區(qū)間，利用梯度h_dB計算每個頻譜點所屬的區(qū)間j，并計算每個區(qū)間的平均功率P_j；其中，t為步驟4中所述的最大閾值；

步驟6：若跟蹤信號為單信標信號則噪聲譜點數為m_noise＝N*l，若跟蹤信號為寬帶信號則噪聲譜點數為m_noise＝(N-BW/Sample)*l，其中BW為跟蹤信號帶寬，Sample為數據頻率，N為新的數據序列中的數據個數，l為系數；

步驟7：按照直方圖統計，得到滿足并計算噪聲譜估計其中，q表示滿足要求的最小區(qū)間數；

步驟8：將所有頻譜與門限T＝K×noise進行比較，小于門限的譜線為噪聲譜線，選取一個噪聲譜線進行噪聲譜估計修正m＝m+1，其中，p_x為小于門限的頻譜點，K為門限因子；

步驟9：更新noise，選取另一個噪聲譜線重復步驟9，直到沒有譜線小于門限T或者達到允許的最大迭代次數；

步驟10：若跟蹤單信標信號，信噪比結果為若跟蹤寬帶信號，信噪比結果為其中∑p_i為寬帶信號內所有頻譜功率和，(N/2+BW/(2×Sample))＜i＜(N/2+BW/(2×Sample))，M為寬帶信號內所有頻譜個數，p_max-1和p_max+1分別為與最大功率值對應數據相鄰的前一個數據和后一個數據。

2.根據權利要求1所述的基于干擾規(guī)避的信噪比解算方法，其特征在于，步驟1中選擇的采樣數據個數C＝2^r，r為整數。

3.根據權利要求1所述的基于干擾規(guī)避的信噪比解算方法，其特征在于，步驟2中參加信噪比計算的數據個數其中，B為監(jiān)測帶寬，Sample為數據頻率，若N為偶數，則N的值加1。

4.根據權利要求1所述的基于干擾規(guī)避的信噪比解算方法，其特征在于，步驟2中，拼接的方法為：

1)按順序取頻譜數據L_i的后(N-1)/2個數據L_C-N/2+1,L_C-N/2+2…L_C；

2)按順序取頻譜數據L_i的前(N-1)/2+1個數據L₁,L₂…L_N/2+1；

3)將L_C-N/2+1,L_C-N/2+2…L_C,L₁,L₂…L_N/2+1作為參加信噪比計算的新的數據序列，總共頻譜數據為N個，每個頻譜數據的功率值用p_i，i＝1,2,3…N表示。

5.根據權利要求1所述的基于干擾規(guī)避的信噪比解算方法，其特征在于，步驟5具體包括以下步驟：

步驟5.1：將功率值劃分為t個區(qū)間，選取一個頻譜點p_i；其中，t為步驟4中所述的最大閾值；

步驟5.2：計算所選取頻譜點所屬的區(qū)間j：

并計算區(qū)間j的平均功率P_j：

其中，n_j為區(qū)間j內功率頻點數，j取整數，若j＜1，則取j為1，P_j-1為上一個頻點計算的平均功率；

步驟5.3：選取另一個頻譜點p_i，返回步驟5.2，直至所有頻點區(qū)間劃分完成。