[發明專利]基于高精度頻率估計的噪聲調幅干擾抑制方法

有效

• 申請號: 201510523710.7
• 申請日: 2015-08-24
• 公開（公告）號: CN105158740B
• 公開（公告）日: 2017-07-28
• 發明（設計）人: 金艷;閔柏成;姬紅兵
• 申請（專利權）人: 西安電子科技大學
• 主分類號: G01S7/35
• 分類號: G01S7/35;G01S7/36
• 代理公司: 陜西電子工業專利中心61205
• 代理人: 田文英,王品華

技術領域

本發明屬于通信技術領域，更進一步涉及雷達信號抗干擾技術領域中的一種基于高精度頻率估計的噪聲調幅干擾抑制方法。本發明通過比較頻率偏差值與門限值的大小估計最佳的頻率偏差，從而得到噪聲調幅干擾中心頻率的高精度估計，利用信號的解調以及頻域對消，實現了強干擾環境下噪聲調幅干擾的抑制。

背景技術

噪聲調幅干擾是雷達對抗中一種重要的干擾方式，具有信號產生簡單，帶寬可變，壓制效果明顯等優點，已成為瞄準式及復合式干擾的重要組成部分。隨著干擾技術的發展，進入雷達接收機的干擾信號能量可超過雷達回波信號達數十分貝以上，雷達回波信號完全淹沒于干擾信號中。在強干擾背景下，通過估計干擾的特征參數，采用對消方法抑制干擾，提高信號處理增益是現代雷達的一個重要抗干擾手段。因此，有效提取干擾信號的特征參數是實現干擾對消的重要前提。

南京郵電大學申請的專利“基于DFT正弦信號的頻率估計方法”(申請號201410151762.1，申請日2014.04.15，申請公開號CN 103941089 A，申請公開日2014.07.23)中提出一種基于DFT的正弦信號頻率估計方法。該方法通過分析Candon算法和2N點DFT算法的性能，在對原始信號進行離散化預處理后，粗估計階段利用Candon算法估計頻率偏差，然后對修正后的原始信號進行2N點DFT算法精估計，增加了對原始信號的頻率修正步驟，可實現在相對頻偏為任意值時，頻率估計的均方根誤差均接近克拉美羅界下限。該方法存在的不足之處是，噪聲調幅干擾中調制噪聲對譜線的幅值影響較大，使Candon算法和2N點DFT算法的頻率估計精度下降，導致利用該專利不能準確估計噪聲調幅干擾的中心頻率，無法實現干擾對消。

中國人民解放軍后勤工程學院申請的專利“一種正弦信號頻率估計的相位匹配方法”(申請號201410342480.X，申請日2014.07.13，申請公開號CN 104142425 A，申請公開日2014.11.12)中提出了一種基于初相匹配的正弦信號頻率估計方法。該方法首先利用正弦信號的線性預測性質、自相關計算同頻信號的自相關序列，消除初相位不同對正弦信號信息融合的影響；然后，對自相關序列進行平域平均處理，提高信號的信噪比；最后，根據自相關序列的特點生成零初相信號序列，通過二者的相關關系構造誤差函數，使誤差函數最小獲得頻率估計值。該方法存在的不足之處是，噪聲調幅干擾中調制噪聲的影響存在相位模糊問題，使得頻率估計誤差較大，導致利用該方法無法實現頻率的高精度估計，不能有效抑制噪聲調幅干擾。

綜上所述，針對噪聲調幅干擾中心頻率的估計問題，已有分析方法僅僅實現了頻率的粗估計，噪聲調幅干擾的中心頻率估計精度低，嚴重影響信號的解調及頻域對消效果，無法有效地抑制干擾。

發明內容

本發明目的在于克服上述已有噪聲調幅干擾中心頻率估計方法的不足，提出一種基于高精度頻率估計的噪聲調幅干擾抑制方法。本發明充分考慮噪聲調幅干擾中調制噪聲的影響，估計最佳的頻率偏差，從而獲得高精度的頻率估計并能有效抑制噪聲調幅干擾。

實現本發明目的的具體思路是：首先對頻率偏差值作線性加權處理，消除干擾中的調制噪聲對干擾頻率估計精度的影響；然后比較頻率偏差值與門限值的大小，估計最佳的頻率偏差，從而得到噪聲調幅干擾的中心頻率，通過解調和頻域對消完成噪聲調幅干擾的抑制。

實現本發明目的的具體步驟如下：

(1)采集信號：

信號采集系統通過連續波雷達的接收機設備，采集連續波雷達天線中任意一段含有線性調頻回波和噪聲調幅干擾的接收信號，將接收信號作為采集信號；

(2)搜索信號頻譜信息：

(2a)對采集信號做快速傅里葉變換，搜索快速傅里葉變換后頻譜峰值的位置；

(2b)搜索快速傅里葉變換后與峰值相鄰譜線處的幅值；

(2c)計算快速傅里葉變換后的量化頻率；

(3)估計頻率偏差值：

(3a)采用雙線幅度公式，估計量化頻率與真實頻率間的第一個頻率偏差值α；

(3b)采用三譜線頻偏校正公式，估計量化頻率與真實頻率間的第二個頻率偏差值β；

(4)估計最佳頻率偏差值：

(4a)比較頻率偏差值α、β與門限值A、B的大小，其中，α表示量化頻率與真實頻率的第一個頻率偏差值，β表示量化頻率與真實頻率間的第二個頻率偏差值，A＝0.06，B＝0.16，當|α|＞B時，執行步驟(4b)，當|α|≤B且|β|＜A時，執行步驟(4c)，當|α|≤B且A＜|β|＜B時，執行步驟(4d)；