CN113820668B - 毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 - Google Patents
毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 Download PDFInfo
- Publication number
- CN113820668B CN113820668B CN202111168979.XA CN202111168979A CN113820668B CN 113820668 B CN113820668 B CN 113820668B CN 202111168979 A CN202111168979 A CN 202111168979A CN 113820668 B CN113820668 B CN 113820668B
- Authority
- CN
- China
- Prior art keywords
- interference
- signal
- echo
- unit
- envelope
- 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
- 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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种毫米波调频连续波雷达的强干扰抑制系统及其抑制方法,涉及新能源汽车及新一代信息技术领域,其技术方案在于:包括干扰包络检测单元及干扰抑制单元,所述干扰包络检测单元包括低通滤波单元、信号统计分析单元及干扰信号定位单元,所述低通滤波单元用于处理雷达回波,所述信号统计分析单元用于统计并判断雷达回波中是否存在强干扰,所述干扰信号定位单元用于对雷达回波中的干扰位置进行定位,所述干扰抑制信号单元用于形成干扰信号区域带有缺失数据的回波信号,并通过回波恢复算法恢复缺失数据。通过干扰包络检测单元判断当前雷达回波是否存在干扰,没有干扰时将不处理,减少了抑制方法的计算复杂度和提高了实时信号处理信息效率。
Description
技术领域
本发明涉及新能源汽车及新一代信息技术领域,尤其涉及一种毫米波调频连续波雷达的强干扰抑制系统及其抑制方法。
背景技术
因毫米波调频连续波雷达的工作带宽大,测距测速及测角度的精度高,体积小,成本低,尤其可在夜晚以及恶劣天气条件可靠稳定工作等等优点,该类型的雷达已广泛应用于车辆的主动防撞、辅助驾驶系统以及无人驾驶系统中。广泛的市场应用使得雷达在同时工作在同一区域的密度增高,雷达发射信号将对另一邻近雷达都是强干扰来源。虽然国内外已有部分研究,有的研究集中在非调频连续波雷达的方法,与当前主流调频连续波雷达芯片背道而驰,无法抑制毫米波调频连续波雷达间的干扰。另外已有对干扰区域简单置零的方法,虽然降低了干扰强度,但丢失了雷达目标信息,在雷达后续的距离、速度及角度压缩中产生了许多强的旁瓣。其他传统的信号处理方法,比如基于小波分析的方法也未能有效地抑制强干扰。
汽车毫米波雷达之间的干扰对汽车雷达对可靠性提出了颠覆性的挑战,目前欧美发达国家仍然没有解决汽车雷达间的干扰难题。
发明内容
本发明的目的是为了解决现有汽车雷达抗干扰能力不足的技术瓶颈问题。
为了实现上述目的,本发明采用了如下技术方案:
一种毫米波调频连续波雷达的强干扰抑制系统,包括干扰包络检测单元及干扰抑制单元,其中所述干扰包络检测单元包括低通滤波单元、信号统计分析单元及干扰信号定位单元,所述低通滤波单元用于处理雷达回波,所述信号统计分析单元用于统计并判断雷达回波中是否存在强干扰,所述干扰信号定位单元用于对雷达回波中的干扰位置进行定位,所述干扰抑制单元用于形成干扰信号区域带有缺失数据的回波信号,并通过回波恢复算法恢复缺失数据。
优选的,所述低通滤波单元用于对回波信号进行滤波得到回波包络,所述低通滤波单元的通带截止数字归一化频率为0.05-0.1,阻带的上边带开始归一化数字频率为0.1-0.25,所述低通滤波单元的阶数为10-30阶。
优选的,所述信号统计分析单元用于对干扰包络信号进行平均统计得到当前的包络均值,并统计得到包络信号的最大幅度峰值。
优选的,所述干扰信号定位单元用于对干扰位置进行定位,当包络信号的最大幅度峰值小于干扰包络信号的均值的三倍,则判定为当前回波不存在强干扰区域,反之,则判定为存在强干扰区域,且干扰发生的位置为高于包络信号均值的区域标记定位为干扰发生的位置。
本申请还提供了一种毫米波调频连续波雷达的强干扰抑制方法,使用上述所述的一种毫米波调频连续波雷达的强干扰抑制系统,其步骤在于:
S1:通过干扰包括检测单元对雷达信号进行取模,并对信号进行处理和分析,
S2:根据干扰包络的统计结果对相应的干扰区域形成缺失数据,并基于L1范数正则化稀疏方法,进而对干扰所在的回波数据进行恢复。
优选的,所述S1中具体步骤为:
然后,进行回波包络的均值统计为信号为u=E(ye),最大幅度峰值为am=max(ye)。
优选的,所述S1中,包络是否存在干扰的判断依据:如果am>3*u,则把满足|y(ni)|>u的部位全标记为干扰区域,形成干扰区域标识向量s,其中ni表示干扰所在的回波中的位置序号,即s(ni)=1,对于没有干扰的位置n,则s(n)=0;反之,如果am<3*u,设判定回波中不存在干扰。
优选的,所述S2中具体步骤为:
首先,对已标记为干扰区域的位置,即s(n)=1处,对应的回波信号删除,形成缺失信号ys为
ys=y*s (1)
假设不含干扰的回波信号在傅立叶变换域中的系数向量为c,则缺失数据信号可表示成
ys=Ac+n (2)
其中A为离散傅立叶变换矩阵,n为系统噪声向量。
因为回波在傅立叶变换的系数c具有良好的稀疏性,可运用L1范数正则化方法得到回波恢复公式为
上述所述的一种毫米波调频连续波雷达的强干扰抑制系统,能极大地抑制干扰并可靠地检测出干扰中的目标,可把信号与干扰加噪声比从-5.2d B,提升至了20.3dB,共提升了25.5dB。恢复目标的结果与没有发生干扰时的结果非常接近。并且通过干扰包络检测单元中利用统计分析可判断当前雷达回波是否存在干扰,没有干扰时将不进行干扰抑制处理,大大减少了抑制方法的计算复杂度和提高了实时信号处理信息效率。
附图说明
图1为本发明的实施例中一种毫米波调频连续波雷达的强干扰抑制系统的工作流程示意图;
图2为本发明的实施例中一种毫米波调频连续波雷达的强干扰抑制系统在雷达信号处理器中数据处理链路顺序原理框图;
图3为本发明的实施例中一种强干扰抑制效果雷达图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下结合具体实施例,对本发明作进一步地详细说明。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用不同于在此描述的其他方式来实施,因此,本发明并不限于下面公开说明书的具体实施例的限制。
一种毫米波调频连续波雷达的强干扰抑制系统,应用于雷达数字中频采样接收机,其包括干扰包络检测单元及干扰抑制单元,其中所述干扰包络检测单元用于对雷达数字中频采样接收机中接收到的雷达回波进行统计并分析,并在雷达回波中存有干扰时,通过干扰抑制单元进行干扰抑制处理。
请参阅图1,在一实施方式中,所述干扰包络检测单元包括低通滤波单元、信号统计分析单元及干扰信号定位单元,其中所述低通滤波单元用于对雷达数字中频采样接收机中接收到的雷达回波进行处理,在一实施方式中,所述低通滤波机的通带截止数字归一化频率在0.05-0.1范围内,阻带的上边带开始归一化数字频率范围为0.1-0.25,所述低通滤波单元的阶数为10阶至30阶。在一实施方式中,所述低通滤波机可以根据雷达的扫描频率来确定相关数值,所述低通滤波单元的相频响应为与阶数相关的线性相位表达,具体设计可以根据窗函数及最小均方误差方法进行设计得到低通滤波单元的系数。
具体的,在一实施方式中,以最小均方误差方法举例,进行设计得到20阶数字低通滤波单元的系数为[0.46,0.88,1.59,2.53,3.68,4.96,6.26,7.47,8.44,9.08,9.30,9.08,8.44,7.47,6.26,4.96,3.68,2.53,1.59,0.88,0.46]。
所述信号统计分析单元用于对干扰包络信号进行统计及分析,在一实施方式中,雷达回波进过低通滤波单元处理后输出为包络信号,通过信号统计分析单元对包络信号进行平均统计得到当前包络均值,并统计得到包络的最大幅度峰值。
所述干扰信号定位单元用于对干扰位置进行定位,在一实施方式中,当包络信号的最大幅度峰值小于干扰包络信号的均值的三倍,则判定为当前回波不存在强干扰区域,反之,则判定为存在强干扰区域,且干扰发生的位置为高于包络信号均值的区域标记定位为干扰发生的位置。
请参阅图2,在一实施方式中,所述干扰抑制信号单元用于形成干扰信号区域带有缺失数据的回波信号,并通过回波恢复算法恢复缺失数据。在一实施方式中,所述干扰抑制信号通过对干扰区域的雷达数据进行稀疏方法处理。具体的,在一实施方式中,根据雷达回波在离散傅里叶交换中具有稀疏的特征,将含有缺失数据的回波信号投影到傅里叶交换域中,通过L1范数约束进行回波恢复。
本申请还提供了一种毫米波调频连续波雷达的强干扰抑制方法,包含以下步骤:
S1:通过干扰包括检测单元对雷达信号进行取模:,并对信号进行处理和分析,具体步骤为:
然后,进行回波包络的均值统计为信号为u=E(ye),最大幅度峰值为am=max(ye)。
接下来,如果am>3*u,则把满足|y(ni)|>u的部位全标记为干扰区域,形成干扰区域标识向量s,其中ni表示干扰所在的回波中的位置序号,即s(ni)=1,对于没有干扰的位置n,则s(n)=0;反之,如果am<3*u,设判定回波中不存在干扰。
在一实施方式中,取模后通过低通滤波单元对雷达信号进行滤波从而得到干扰的包络,并根据包络的统计特性判断是否存在干扰并定位干扰所在位置,然后将干扰位置信息送达干扰抑制单元。
在一实施方式中,包络是否存在干扰的判断依据:当包络的最大幅度峰值小于包络均值的三倍时,则当前回波不存在强干扰,反之,则存在强干扰,并且幅度峰值高于包络均值的区域标记定位为干扰发生的位置。
S2:根据干扰包络的统计结果对相应的干扰区域形成缺失数据,并基于L1范数正则化稀疏方法,进而对干扰所在的回波数据进行恢复。
具体步骤为:
首先,对已标记为干扰区域的位置,即s(n)=1处,对应的回波信号删除,形成缺失信号ys为
ys=y*s (1)
假设不含干扰的回波信号在傅立叶变换域中的系数向量为c,则缺失数据信号可表示成
ys=Ac+n (2)
其中A为离散傅立叶变换矩阵,n为系统噪声向量。
因为回波在傅立叶变换的系数c具有良好的稀疏性,可运用L1范数正则化方法得到回波恢复公式为
请参阅图3,图3为本申请系统装置在强干扰环境下的一个干扰抑制实验结果图。其中雷达的起始工作频率为77GHz,带宽为600MHz,相干积累脉冲数为128,主雷达的信号扫描周期为51.2us,两个干扰雷达的带宽分别为300MHz和500MHz,调频率和扫描周期都跟主雷达不同,但工作频带存在重叠部分。
图3(b)为没有使用干扰抑制的幅度示意图,可以看出强干扰把其中一个目标淹没,并产生了许多强的虚假目标。图3(c)运用现有的小波抑制方法不能有效抑制干扰恢复目标。图3(d)使用了毫米波调频连续波雷达的强干扰抑制系统,通过附图可以看出,本申请所提供的一种毫米波调频连续波雷达的强干扰抑制系统,能极大地抑制干扰并可靠地检测出干扰中的目标,可把信号与干扰加噪声比从-5.2d B,提升至了20.3dB,共提升了25.5dB。恢复目标的结果与没有发生干扰时的结果非常接近。并且通过干扰包络检测单元中利用统计分析可判断当前雷达回波是否存在干扰,没有干扰时将不进行干扰抑制处理,大大减少了抑制方法的计算复杂度和提高了实时信号处理信息效率。
以上所述仅为本发明的实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围。
Claims (7)
1.一种毫米波调频连续波雷达的强干扰抑制系统,其特征在于:包括干扰包络检测单元及干扰抑制单元,其中所述干扰包络检测单元包括低通滤波单元、信号统计分析单元及干扰信号定位单元,所述低通滤波单元用于处理雷达回波,所述信号统计分析单元用于统计并判断雷达回波中是否存在强干扰,所述干扰信号定位单元用于对雷达回波中的干扰位置进行定位,当包络信号的最大幅度峰值小于干扰包络信号的均值的三倍,则判定为当前回波不存在强干扰区域,反之,则判定为存在强干扰区域,且干扰发生的位置为高于包络信号均值的区域标记定位为干扰发生的位置,所述干扰抑制单元用于形成干扰信号区域带有缺失数据的回波信号,并通过回波恢复算法恢复缺失数据。
2.根据权利要求1所述的毫米波调频连续波雷达的强干扰抑制系统,其特征在于:所述低通滤波单元用于对回波信号进行滤波得到回波包络,低通滤波单元的通带截止数字归一化频率为0.05-0.1,阻带的上边带开始归一化数字频率为0.1-0.25,所述低通滤波单元的阶数为10-30阶。
3.根据权利要求1所述的毫米波调频连续波雷达的强干扰抑制系统,其特征在于:所述信号统计分析单元用于对干扰包络信号进行平均统计得到当前的包络均值,并统计得到包络信号的最大幅度峰值。
4.一种毫米波调频连续波雷达的强干扰抑制方法,其特征在于:使用权利要求1-3所述的一种毫米波调频连续波雷达的强干扰抑制系统,其步骤在于:
S1:通过干扰包括检测单元对雷达信号进行取模,并对信号进行处理和分析,
S2:根据干扰包络的统计结果对相应的干扰区域形成缺失数据,并基于L1范数正则化稀疏方法,进而对干扰所在的回波数据进行恢复。
6.根据权利要求5所述的毫米波调频连续波雷达的强干扰抑制方法,其特征在于:所述S1中,包络是否存在干扰的判断依据:如果am>3*u,则把满足|y(ni)|>u的部位全标记为干扰区域,形成干扰区域标识向量s,其中ni表示干扰所在的回波中的位置序号,即s(ni)=1,对于没有干扰的位置n,则s(n)=0;反之,如果am<3*u,设判定回波中不存在干扰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111168979.XA CN113820668B (zh) | 2021-10-08 | 2021-10-08 | 毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111168979.XA CN113820668B (zh) | 2021-10-08 | 2021-10-08 | 毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113820668A CN113820668A (zh) | 2021-12-21 |
CN113820668B true CN113820668B (zh) | 2022-08-26 |
Family
ID=78916233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111168979.XA Active CN113820668B (zh) | 2021-10-08 | 2021-10-08 | 毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113820668B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114152915B (zh) * | 2021-12-30 | 2024-08-30 | 广州极飞科技股份有限公司 | 雷达信号处理设备、雷达设备和可移动平台 |
CN114859303B (zh) * | 2022-04-12 | 2024-09-27 | 南通大学 | 一种具有低运算复杂度的毫米波雷达干扰快速抑制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914411A (zh) * | 2015-05-21 | 2015-09-16 | 中国人民解放军63892部队 | 一种单脉冲雷达对交叉眼干扰的检测判别及抑制方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044686B (zh) * | 2015-08-03 | 2017-05-10 | 中国电子科技集团公司第二十八研究所 | 一种雷达密集假目标干扰抑制方法 |
CN106226744B (zh) * | 2016-07-06 | 2018-06-08 | 中国电子科技集团公司第二十八研究所 | 一种基于投影变换的干扰特征提取方法 |
US10408932B2 (en) * | 2016-12-16 | 2019-09-10 | Automotive Research & Testing Center | Environment recognition system using vehicular millimeter wave radar |
CN106959439B (zh) * | 2017-05-26 | 2020-02-07 | 南通大学 | 汽车调频连续波雷达的强干扰抑制方法及系统 |
CN111190147B (zh) * | 2020-02-12 | 2023-04-21 | 南京理工大学 | 一种针对毫米波汽车防撞雷达的干扰抑制方法 |
CN113433523B (zh) * | 2021-06-17 | 2022-09-16 | 同济大学 | 一种基于双框架循环rls的汽车雷达干扰抑制方法 |
-
2021
- 2021-10-08 CN CN202111168979.XA patent/CN113820668B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914411A (zh) * | 2015-05-21 | 2015-09-16 | 中国人民解放军63892部队 | 一种单脉冲雷达对交叉眼干扰的检测判别及抑制方法 |
Non-Patent Citations (1)
Title |
---|
舰载无源综合脉冲孔径雷达干扰抑制;张雅斌等;《电子与信息学报》;20080115(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113820668A (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113820668B (zh) | 毫米波调频连续波雷达的强干扰抑制系统及其抑制方法 | |
CN111399002B (zh) | 一种基于两级神经网络的gnss接收机组合干扰分类识别方法 | |
CN111624574A (zh) | 弱目标探测的目标检测方法、系统、存储介质和装置 | |
CN110133632B (zh) | 一种基于cwd时频分析的复合调制信号识别方法 | |
WO2008082973A1 (en) | System and method for reducing the effect of a radar interference signal | |
CN114217315B (zh) | 一种毫米波雷达动态目标快速检测方法 | |
CN109856600A (zh) | 一种毫米波防撞雷达抗干扰的方法 | |
CN112731307B (zh) | 基于距离-角度联合估计的ratm-cfar检测器及检测方法 | |
CN110398719B (zh) | 基于杂波图原理的雷达杂波信号抑制方法和雷达探测系统 | |
CN113176542B (zh) | 一种时域自适应抗主瓣自卫式干扰方法及系统 | |
CN114114184B (zh) | 一种汽车雷达之间的交叉干扰抑制系统及方法 | |
CN112327260B (zh) | 一种sar回波数据中脉冲式干扰信号的抑制方法和装置 | |
CN113866750A (zh) | 一种基于毫米波雷达的行人目标检测和跟踪方法 | |
CN112285656B (zh) | 一种基于航迹特征的抗时敏目标欺骗干扰方法 | |
CN111707994B (zh) | 一种毫米波汽车防撞雷达二维干扰抑制方法 | |
CN114265015A (zh) | 一种基于杂波特性感知的精细化雷达检测方法 | |
CN110488238B (zh) | 一种基于压缩感知的高分辨宽测绘带解距离模糊方法 | |
CN114137523A (zh) | 车载毫米波雷达原始点云的噪点排除方法和系统 | |
CN118244210B (zh) | 一种间歇采样转发干扰主动对抗方法、装置及系统 | |
Lopez-Risueno et al. | Two digital receivers based on time-frequency analysis for signal interception | |
CN116224280B (zh) | 雷达目标检测方法、装置、雷达设备及存储介质 | |
CN113640767B (zh) | 一种基于方差的低分辨雷达目标识别方法 | |
CN115494461B (zh) | 一种基于维特比算法的脉内转发式干扰时频域抑制方法 | |
CN118131165B (zh) | 一种基于自回归模型的干扰信号恢复方法 | |
KR102066579B1 (ko) | 프랙탈 차원을 이용한 표적 인식용 밀리미터파 탐색기 및 이를 이용한 표적 인식방법 |
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 |