CN110412572B - P波段合成孔径雷达成像干扰抑制方法 - Google Patents
P波段合成孔径雷达成像干扰抑制方法 Download PDFInfo
- Publication number
- CN110412572B CN110412572B CN201910670721.6A CN201910670721A CN110412572B CN 110412572 B CN110412572 B CN 110412572B CN 201910670721 A CN201910670721 A CN 201910670721A CN 110412572 B CN110412572 B CN 110412572B
- Authority
- CN
- China
- Prior art keywords
- signal
- radar
- frequency
- rotation
- receiving
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种P波段合成孔径雷达成像干扰抑制方法,包括:设置P波段成像雷达发射机与接收机工作模式;根据雷达接收机在发送‑接收工作模式下每个雷达脉冲重复时间内收到的雷达回波构建旋转域信号;对所述旋转域信号进行滤波处理,得到滤波后的旋转域信号;在发送‑接收工作模式接收完毕后,从所述旋转域信号中减去滤波后的旋转域信号,得到剩余信号能量QS(ω);计算在一次静默接收工作模式下接收到的回波信号的信号能量QB(ω);基于所述剩余信号能量QS(ω)和所述回波信号的信号能量QB(ω)定义目标函数Jn,以目标函数Jn作为干扰抑制效果评估指标,不断调整所述滤波器的中心频率、截止频率进行迭代,直到||Jn+1‑Jn||≤ε。本发明解决了现有方法干扰信号参数估计不准确导致抑制程度不够或者过度的缺陷。
Description
技术领域
本发明涉及雷达成像技术领域,尤其涉及一种P波段SAR成像干扰抑制方法。
背景技术
合成孔径雷达(Syntheticaperture Radar,SAR)作为一种全天时、全天候的主动式对地观测系统,在地表环境监测、海洋观测、资源勘查等民用测绘方面具有重要作用。雷达波段(radar frequency band)指雷达发射电波的频率范围,其度量单位是赫兹(Hz)或周/秒(C/s),P波段的频率在0.23-1Ghz,波长范围为130-30cm,但P波段SAR由于工作频率低,需要提高雷达系统带宽才能实现高分辨率成像,这样导致P波段SAR工作带宽与电视信号、调频广播信号、移动通信信号等部分民用信号波段重叠,,这就造成了P波段SAR不可避免地遭到来自民用信号的干扰,严重影响P波段SAR成像质量,因此亟需有效的P波段SAR成像干扰抑制方法,图2为含有干扰信号的P波段SAR雷达回波频谱示意图。现有的P波段干扰抑制方法如要有以下几类:
1、陷波法。陷波法是利用雷达回波在距离频域的能量分布特性搜索可能的干扰频点,然后通过LMS时域自适应滤波器、频域置零、频域带阻滤波器、频域限波器、信号分解等方法实现干扰抑制。
该类方法的缺陷在于基于信号能量作为判别干扰的依据,容易受地物回波特性以及其它不确定性噪声的影响,导致无法准确判别。其次,该类方法关键点在于干扰检测门限的设置,当估计的干扰信号的频点不准确时会对雷达回波滤除,影响成像质量。另外,该类方法仅适用于射频窄带干扰的抑制,对于电视广播这种宽带信号,容易造成频谱断裂,将损失有用信号,并且导致旁瓣抬升,影响成像质量。
2、经验模态分解、小波分解、压缩感知等信号分离方法。但是当干扰信号频率相近时,或者干扰成分过多时,这类方法无法将雷达回波与干扰信号准确分离,此外,这类方法多维运算,计算量大,不利于雷达图像的快速处理,难以满足工程化需要。
3、设置专用的专用旁听通道用来实现干扰信号的参数估计,但这种方法将增加系统的复杂性,并且引入了通道间相互干扰。
发明内容
(一)要解决的技术问题
有鉴于此,本发明的主要目的在于提供一种P波段合成孔径雷达成像干扰抑制的方法,以抑制P波段合成孔径雷达成像干扰,解决现有方法干扰信号参数估计不准确导致抑制程度不够或者过度的问题。
(二)技术方案
为达到上述目的,本发明采用的技术方案如下:
一种P波段合成孔径雷达成像干扰抑制的方法,包括:
设置P波段合成孔径雷达发射机与接收机工作模式,发射机在连续进行“发送-接收”工作模式后,进行一次静默接收工作模式;
根据合成孔径雷达接收机在“发送-接收”工作模式下每个雷达脉冲重复时间内收到的雷达回波构建旋转域信号;
对所述旋转域信号进行滤波处理,得到滤波后的旋转域信号;
在“发送-接收”工作模式接收完毕后,从所述旋转域信号中减去滤波后的旋转域信号,得到剩余信号能量QS(ω),其中QS(ω)整体表示剩余信号能量,ω表示信号频率;
计算在一次静默接收工作模式下接收到的回波信号的信号能量QB(ω),其中QB(ω)整体表示接收到的回波信号的信号能量;
基于所述剩余信号能量QS(ω)和所述回波信号的信号能量QB(ω)定义目标函数Jn,以目标函数Jn作为干扰抑制效果评估指标,不断调整所述滤波器的中心频率、截止频率进行迭代,直到||Jn+1-Jn||≤ε,ε表示迭代终止参数,n表示迭代的次数。
(三)有益效果
(1)本发明提供的P波段合成孔径雷达成像干扰抑制的方法,发射机在连续进行n次“发送-接收”工作模式后,紧接着进行一次静默接收工作模式,雷达系统不需要设立专门的干扰信号接收天线,在准确提取干扰信号的同时能有效降低系统的复杂性,另外增强了信号相参性。
(2)本发明提供的P波段合成孔径雷达成像干扰抑制的方法,通过带通滤波器对所述旋转域信号进行处理,得到滤波后的旋转域信号/>此时/>中则主要包含线性调频信号分量,而混合信号中的调频广播和移动通信窄带干扰、电视信号宽带干扰被滤除,并且不会损失雷达回波中线性调频信号分量的信号能量,保证其频谱完整,解决了现有方法干扰信号抑制程度不够或者过度的缺陷。
(3)本发明提供的P波段合成孔径雷达成像干扰抑制的方法,无需进行复杂的矩阵运算而实现了P波段雷达回波中有用的线性调频信号与干扰信号的分离,计算量小,计算效率高,满足工程需求。
附图说明
图1是本发明实施例的P波段SAR成像干扰抑制方法的步骤示意图;
图2是本发明实施例的含有干扰信号的P波段合成孔径雷达回波频谱示意图;
图3是本发明实施例的P波段合成孔径雷达回波旋转域信号频谱示意图;
图4是本发明实施例的P波段合成孔径雷达发射机与接收机工作模式示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明公开了一种P波段合成孔径雷达成像干扰抑制的方法,图1为本发明实施例的步骤示意图,包括以下步骤:
S11,设置P波段SAR成像雷达发射机与接收机工作模式,发射机在连续进行“发送-接收”工作模式后,进行一次静默接收工作模式;
具体的,发射机在连续进行n次“发送-接收”工作模式后,紧接着进行一次静默接收工作模式,静默接收工作模式时间长度等于雷达脉冲重复时间,这样每隔n+1个雷达脉冲重复时间作为一次循环周期,以一次循环周期处理雷达回波,其中n≥7,可选的,n=9,如图4为本发明实施例的P波段合成孔径雷达发射机与接收机工作模式示意图。
S12根据合成孔径雷达接收机在“发送-接收”工作模式下每个雷达脉冲重复时间内收到的雷达回波构建旋转域信号;
具体的,根据P波段SAR成像雷达系统参数构建频率旋转算子其中j表示信号虚部,t表示时间,fc表示雷达信号载频,k表示雷达信号调频率。对雷达接收机在“发送-接收”工作模式下每个雷达脉冲重复时间内收到的雷达回波si(t)构建旋转域信号其计算公式如下所示:
式中i表示i个雷达脉冲重复时间内雷达接收机在“发送-接收”工作模式接收到的含有干扰的雷达回波混合信号,其中τ表示截取时间窗长度,ω表示信号频率,fc表示雷达发射信号载频,k表示雷达信号调频率,si(t)表示雷达回波,t表示时间,i表示i个雷达脉冲重复时间内雷达接收机在“发送-接收”工作模式接收到的含有干扰的雷达回波混合信号;频率旋转算子计算公式为图3表示雷达接收机在“发送-接收”工作模式下每个PRT时间内收到的雷达回波转域信号频谱示意图,将含有干扰的雷达回波混合信号在时频域旋转某一角度,达到改变混合信号中各个信号分量的时频分布的目的,在含有干扰的雷达回波混合信号中,由于线性调频信号分量与构建的频率旋转算子匹配,因此经过旋转后的时频分布的频谱能量最集中,并其能量分布平行于时间轴,并且带宽最窄,而对于混合信号中的调频广播和移动通信窄带干扰、电视信号宽带干扰,由于与频率旋转算子不匹配,因此,对应旋转后的时频分布将变分散,频谱能量集中度降低,这样达到了提高线性调频信号分量与干扰信号分量、噪声分量的区分度。
S13,对所述旋转域信号进行滤波处理,得到滤波后的旋转域信号;
具体的,设计带通滤波器第i个旋转域信号进行处理,经过带通滤波后得到第i个滤波后的旋转域信号/>此时/>中则主要包含线性调频信号分量,而混合信号中的调频广播和移动通信窄带干扰、电视信号宽带干扰被滤除,并且不会损失雷达回波中线性调频信号分量的信号能量,保证其频谱完整,然后对第i个滤波后的旋转域信号/>进行旋转反变换,从而复原该雷达回波中的线性调频信号分量;其中ω表示信号频率,fc表示雷达发射信号载频,k表示雷达信号调频率。
进一步的,带通滤波器选用Vold-Kalman滤波器,该滤波器对于步骤三中旋转反变换过程,具有很好的保相性,并且该滤波器的中心频率、截止频率可自适应调节,调节依据由目标函数迭代终止条件决定。
S14,在“发送-接收”工作模式接收完毕后,从所述旋转域信号中减去滤波后的旋转域信号,得到剩余信号能量QS(ω),其中ω表示信号频率;
其中,ω表示信号频率,fc表示雷达发射信号载频,k表示雷达信号调频率。
S15,计算在一次静默接收工作模式下接收到的回波信号的信号能量QB(ω),其中ω表示信号频率;
根据频率旋转算子对雷达接收机在静默接收工作模式下每个雷达脉冲重复时间内收到的雷达回波si(t)构建旋转域信号其中τ表示截取时间窗长度,ω表示信号频率,fc表示雷达发射信号载频,k表示雷达信号调频率,siB(τ)表示雷达回波,i表示i个雷达脉冲重复时间内雷达接收机在静默接收工作模式下接收到的含有干扰的雷达回波混合信号。
S16,基于所述剩余信号能量QS(ω)和所述回波信号的信号能量QB(ω)定义目标函数Jn,以目标函数Jn作为干扰抑制效果评估指标,不断调整所述滤波器的中心频率、截止频率进行迭代,直到时,ε表示迭代终止参数,n表示迭代的次数。
具体的,定义目标函数Jn作为干扰抑制效果评估指标,对干扰抑制效果进行客观、定量的评价,当||Jn+1-Jn||≤ε时,ε表示迭代终止参数,完成步骤三中带通滤波器参数达到最优,保持不变,否则调整所述Vold-Kalman滤波器的中心频率、截止频率,实现滤波器参数的自适应调节。
所述目标函数Jn的计算公式为:
在本实施例中发射机在连续进行9次“发送-接收”工作模式后,紧接着进行一次静默接收工作模式,采用这种策略,雷达系统不需要设立专门的干扰信号接收天线,在准确提取干扰信号的同时能有效降低系统的复杂性,另外增强了信号相参性。此外,通过带通滤波器对旋转域信号进行处理,得到滤波后的旋转域信号/>此时/>中则主要包含线性调频信号分量,而混合信号中的调频广播和移动通信窄带干扰、电视信号宽带干扰被滤除,并且不会损失雷达回波中线性调频信号分量的信号能量,保证其频谱完整,解决了现有方法干扰信号抑制程度不够或者过度的缺陷。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种P波段合成孔径雷达成像干扰抑制的方法,其特征在于,包括:
设置P波段合成孔径雷达发射机与接收机工作模式,发射机在连续进行“发送-接收”工作模式后,进行一次静默接收工作模式;
根据合成孔径雷达接收机在“发送-接收”工作模式下每个雷达脉冲重复时间内收到的雷达回波构建旋转域信号;
对所述旋转域信号进行滤波处理,得到滤波后的旋转域信号;
在“发送-接收”工作模式接收完毕后,从所述旋转域信号中减去滤波后的旋转域信号,得到剩余信号能量QS(w),其中QS(w)整体表示剩余信号能量,ω表示信号频率;
计算在一次静默接收工作模式下接收到的回波信号的信号能量QB(w),其中QB(w)整体表示接收到的回波信号的信号能量;
基于所述剩余信号能量QS(w)和所述回波信号的信号能量QB(w)定义目标函数Jn,以目标函数Jn作为干扰抑制效果评估指标,不断调整所述滤波器的中心频率、截止频率进行迭代,直到||Jn+1-Jn||≤ε,ε表示迭代终止参数,n表示迭代的次数。
2.根据权利要求1所述的方法,其特征在于,所述雷达发射机与接收机工作模式为:
发射机在连续进行n次“发送-接收”工作模式后,紧接着进行一次静默接收工作模式,静默接收工作模式时间长度等于雷达脉冲重复时间,这样每隔n+1个雷达脉冲重复时间作为一次循环周期,以一次循环周期处理雷达回波,其中n≥7。
3.根据权利要求2所述的方法,其特征在于,所述发射机在连续进行n次“发送-接收”工作模式,n=9。
6.根据权利要求5所述的方法,其特征在于,所述带通滤波器为Vold-Kalman滤波器。
7.根据权利要求1所述的方法,其特征在于,所述计算在一次静默接收工作模式下接收到的回波信号的信号能量QB(w),包括:
根据频率旋转算子对雷达接收机在静默接收工作模式下每个雷达脉冲重复时间内收到的雷达回波si(t)构建旋转域信号其中τ表示截取时间窗长度,ω表示信号频率,fc表示雷达发射信号载频,k表示雷达信号调频率,siB(τ)表示雷达回波,i表示i个雷达脉冲重复时间内雷达接收机在静默接收工作模式下接收到的含有干扰的雷达回波混合信号;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910670721.6A CN110412572B (zh) | 2019-07-23 | 2019-07-23 | P波段合成孔径雷达成像干扰抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910670721.6A CN110412572B (zh) | 2019-07-23 | 2019-07-23 | P波段合成孔径雷达成像干扰抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110412572A CN110412572A (zh) | 2019-11-05 |
CN110412572B true CN110412572B (zh) | 2023-03-24 |
Family
ID=68362787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910670721.6A Active CN110412572B (zh) | 2019-07-23 | 2019-07-23 | P波段合成孔径雷达成像干扰抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110412572B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111948653B (zh) * | 2020-07-31 | 2022-03-18 | 上海卫星工程研究所 | 基于p波段合成孔径雷达的林下目标检测方法及系统 |
CN113064122B (zh) * | 2021-02-25 | 2023-02-17 | 上海卫星工程研究所 | P波段sar干扰抑制算法的性能评价方法、系统及介质 |
CN113176543B (zh) * | 2021-03-26 | 2022-06-24 | 上海卫星工程研究所 | 多维度信息联合的射频干扰抑制方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014048193A1 (zh) * | 2012-09-28 | 2014-04-03 | 北京理工大学 | 一种用于舰艇编队情况下同型雷达同频干扰抑制方法 |
CN107132536A (zh) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | 一种消除目标微动对雷达成像干扰的方法 |
CN108776336A (zh) * | 2018-06-11 | 2018-11-09 | 电子科技大学 | 一种基于emd的自适应穿墙雷达静止人体目标定位方法 |
-
2019
- 2019-07-23 CN CN201910670721.6A patent/CN110412572B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014048193A1 (zh) * | 2012-09-28 | 2014-04-03 | 北京理工大学 | 一种用于舰艇编队情况下同型雷达同频干扰抑制方法 |
CN107132536A (zh) * | 2017-04-10 | 2017-09-05 | 中国科学院国家空间科学中心 | 一种消除目标微动对雷达成像干扰的方法 |
CN108776336A (zh) * | 2018-06-11 | 2018-11-09 | 电子科技大学 | 一种基于emd的自适应穿墙雷达静止人体目标定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110412572A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110412572B (zh) | P波段合成孔径雷达成像干扰抑制方法 | |
CN114050951B (zh) | 超短波跳频信号参数盲估计方法 | |
WO2015172622A1 (zh) | 一种高频地波雷达射频干扰抑制的方法 | |
CN110519003B (zh) | 一种基于信号特征差异的无人机上下行通信链路识别方法及装置 | |
CN109962714B (zh) | 一种电磁频谱伞罩数字域自干扰抑制方法及系统 | |
Ford et al. | RFI mitigation techniques in radio astronomy | |
Gurgel et al. | Radio frequency interference suppression techniques in FMCW modulated HF radars | |
Ads et al. | A comprehensive sounding of the ionospheric HF radio link from Antarctica to Spain | |
EP2662704B1 (en) | Method and device for non-uniform sampling of singularity point of multi-channel synthetic-aperture radar (SAR) system | |
CN106324598B (zh) | 一种多通道子带信号的均衡、同步方法及其系统 | |
CN106680790B (zh) | 一种伪码调相连续波雷达干扰检测方法 | |
FR2975193A1 (fr) | Procede et systeme de localisation d'interferences affectant un signal de radionavigation par satellite | |
JP2011117899A (ja) | レーダ装置 | |
EP2096768B1 (en) | Apparatus and method for filtering narrow band interference in received pulsed radio signal | |
Lamont-Smith et al. | Filtering approaches for interference suppression in low‐frequency SAR | |
Belgibaev et al. | Software-Defined Radio Ionosonde for Diagnostics of Wideband HF Channels with the Use of USRP Platform | |
CN103176184B (zh) | 一种结合干扰抑制方法的p波段sar成像处理方法 | |
CN103959087B (zh) | 探地雷达的自适应干扰抑制 | |
CN106154236B (zh) | 一种基于cemd抑制高频地波雷达射频干扰的方法 | |
US8737909B2 (en) | Digital processing for co-site interference mitigation | |
van der Merwe et al. | High-rate DFT-based data manipulator (HDDM) algorithm for effective interference mitigation | |
CN103248337A (zh) | 一种空时级联抑制斜向探测电离图外部干扰的方法 | |
Ruan et al. | Wide band noise interference suppression for SAR with dechirping and eigensubspace filtering | |
CN113238194B (zh) | 基于分数域-频域处理的宽带相控阵雷达抗诱骗干扰方法 | |
CN112630737B (zh) | 对雷达中频回波信号的预处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |