CN112731304A - 一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 - Google Patents
一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 Download PDFInfo
- Publication number
- CN112731304A CN112731304A CN202011484658.6A CN202011484658A CN112731304A CN 112731304 A CN112731304 A CN 112731304A CN 202011484658 A CN202011484658 A CN 202011484658A CN 112731304 A CN112731304 A CN 112731304A
- Authority
- CN
- China
- Prior art keywords
- azimuth
- theta
- angle domain
- azimuth angle
- arc array
- 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.)
- Granted
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000001629 suppression Effects 0.000 title claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 16
- 230000006835 compression Effects 0.000 claims abstract description 11
- 238000007906 compression Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000002592 echocardiography Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
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
技术领域
本发明属于雷达信号处理领域,特别涉及了一种弧形阵列雷达杂波抑制方法。
背景技术
弧形阵列天线通常采用频率调制连续波收发分置的等效MIMO机制快速形成弧形孔径,针对机载或车载平台附近是待探测区域进行运动目标探测。在弧形阵列天线雷达***中,在平行于电磁波传播的方向,发射信号决定了距离向分辨率;沿弧型角度设置若干天线阵元形成弧形孔径,回波信号进行相干积累以决定方位向分辨率。该***依靠微波高速开关实现信号的等效多发多收,不间断取得目标信息,同时弧形阵列天线受平台振动影响较小,且不依靠平台移动以实现合成孔径,适用条件较为广泛。该雷达***具有受平台振动影响小、无需平台移动的优点。
综上所述,弧形阵列天线雷达***能够获得相较一般线阵天线更大的观测范围,在视场观测、目标检测等方面都有优势。但是弧形阵列天线在使用过程中会受到杂波的干扰,影响阵列天线的性能,而传统的滤波方法对杂波的抑制效果并不理想。
发明内容
为了解决上述背景技术提到的技术问题,本发明提出了一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法。
为了实现上述技术目的,本发明的技术方案为:
一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法,包括以下步骤:
(1)对弧形阵列雷达的回波信号进行距离向处理;
(2)方位角度域滤波;
(3)对运动目标进行方位向压缩,完成杂波抑制。
进一步地,在步骤(1)中,弧形阵列雷达采用调频连续波FMCW信号,其发射信号的形式如下:
上式中,ST(t)为发射信号,fc表示发射信号载频,K表示调频率,t表示采样时间,j为虚数单位;
对回波信号进行解调频Dechirp处理、补偿残余视频相位、距离频率替换和距离向IFFT后,得回波信号的形式如下:
上式中,SIFFT(t,θn)为回波信号,θn为弧形阵列波束指示角,σ为目标的雷达截面积,pr为距离向的点扩散函数,wa为方位角度域的矩形窗函数,f为信号中心频率,c为电磁波传播速度,η为方位向时间,且θ为方位角度域,ω为波束扫描角速度,Rmv(η)表示运动目标与弧形阵列雷达中心间的距离。
进一步地,在步骤(2)中,设弧形阵列天线波束的3dB宽度为θa,波束扫描范围是[-θb,θb],天线波束指向角为θx∈[-θb,θb],且位于天线波束的3dB宽度θa以内的目标才有雷达回波信息;
滤除回波信号在方位角度域中内角度域θinner部分的回波,保留外角度域θouter中包含运动目标信息的回波信号,方位角度域滤波过程如下:
其中,θclutter为静止目标相对弧形阵列雷达中心的角度。
进一步地,在步骤(3)中,对运动目标进行方位向压缩,使得运动目标落入斜距为R0、方位角为θclutter的成像单元内,该成像单元相对于每个波束扫描角对应的天线阵元的距离如下:
其中,Rarc为阵列弧半径;
采用上述技术方案带来的有益效果:
弧形阵列雷达具有特殊的结构,建立方位向角频率与方位向扫描角度之间的对应关系,其方位向角频率带宽决定了弧形阵列雷达的方位向分辨能力。本发明借助方位向角频率和方位向扫描角度的对应关系,通过构造相应的方位角度域滤波器,对回波进行方位角度域滤波处理,实现杂波抑制。仿真结果验证了本发明的可行性和有效性。
附图说明
图1是本发明的方法流程图;
图2是弧形阵列天线波束扫描示意图;
图3是实施例中接收回波信号仿真图;
图4是实施例中未经杂波抑制处理的目标情况图;
图5是实施例中未经杂波抑制处理的目标情况(方位向FFT)图;
图6是实施例中频域滤波结果图;
图7是实施例中频域滤波结果(方位向FFT)图;
图8是实施例中方位角度域滤波结果图;
图9是实施例中方位角度域滤波结果(方位向FFT)图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明设计了一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法,如图1所示,步骤如下:
步骤1:对弧形阵列雷达的回波信号进行距离向处理;
步骤2:方位角度域滤波;
步骤3:对运动目标进行方位向压缩,完成杂波抑制。
每个步骤的具体实现过程如下:
1、回波信号距离向处理
弧形阵列雷达采用调频连续波FMCW信号,发射信号为:
其中,fc表示发射信号载频,K表示调频率,t表示采样时间。
经过解调频Dechirp处理、补偿残余视频相位、距离频率替换、距离向IFFT后,回波信号的形式如下:
其中,θn为弧形阵列波束指示角,σ为目标的雷达截面积RCS,pr为距离向的点扩散函数,wa为方位角度域的矩形窗函数,f为信号中心频率,c为电磁波传播速度,η为方位向时间,且θ为方位角度域,ω为波束扫描角速度,Rmv(η)表示运动目标与弧形阵列雷达中心间的距离。
假设某静止目标(可产生杂波)成像单元与雷达弧形阵列中心间的距离由Rclutter(η)表示,有:
其中,Rarc为阵列弧半径,R0和θclutter分别为静止目标相对雷达弧形阵列中心的距离和角度。
2、方位角度域滤波
杂波的相位历程可以表示为:
杂波方位角度对应的方位角频率ρclutter可以表示为:
如图2所示,弧形阵列天线波束的3dB宽度为θa,波束扫描范围是在[-θb,θb]之间,设天线波束指向角为θx∈[-θb,θb],且位于天线波束的3dB宽度θa以内的目标才有雷达回波信息。
θall=θinner+θouter
若要避免杂波谱混叠、抑制静止目标产生的杂波,应滤除回波信号在方位角度域中内角度域θinner部分的回波,保留外角度域θouter中包含运动目标信息的回波信号,角度域滤波过程如下:
3、方位向压缩
对运动目标进行方位向压缩,可以使得运动目标落入斜距为R0、方位角为θclutter的成像单元内。该成像单元相对于每个波束扫描角对应的天线阵元的距离,可以计算为:
方位向压缩可表示为将方位角度域滤波后的回波信号与滤波函数的相位相乘,匹配过程如下:
完成以上步骤,可获得较好的杂波抑制效果。
以下参照上述步骤,通过建立MATLAB软件仿真模型以验证发明的有效性,具体的检测仿真***参数如表1。
表1检测仿真***参数
为展现滤波算法的杂波抑制效果,布置6个运动点目标和1个静止点目标用于回波仿真,目标参数如表2。
表2点目标仿真参数
注:(1)4号点目标为预设静止目标,产生杂波信号;(2)在弧形阵列机载雷达***中,运动目标切向速度vq对聚焦效果影响微小,同时也基本不会引起运动目标散焦,为便于仿真,设置目标切向速度vq均为0m/s。
首先进行步骤1,回波信号距离向处理,回波仿真如图3所示。将信号进行匹配滤波处理,可观察全角度域θall中、没有经过任何杂波抑制处理的目标情况,如图4、图5所示。
(1)常规频域滤波法
将频域滤波后的信号进行傅里叶逆变换,恢复为方位向时域信号。将信号进行方位向压缩处理,可观测该方法的杂波抑制效果,如图6、图7所示。
将图6与图4对比可知,频域滤波后,杂波(4号点目标)信号能量衰减显著,分析信号数据可得,杂波能量衰减约为38.27dB。
(2)本发明方位向角度域滤波法
进行步骤2,滤除回波信号在方位角度域中内角度域θinner部分的回波,保留外角度域θouter中包含运动目标信息的回波信号。
之后,进行步骤3,将信号进行方位向压缩处理,观测该方法的杂波抑制效果。以上两个步骤的计算过程如下:
观测方位向角度域滤波后的杂波抑制效果如图8、图9所示。
对比图8与图4分析可知,该算法对于杂波(4号点目标)的抑制效果较好,分析信号数据可得,杂波能量衰减约为35.20dB,本发明的有效性和可行性得到验证。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (4)
1.一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法,其特征在于,包括以下步骤:
(1)对弧形阵列雷达的回波信号进行距离向处理;
(2)方位角度域滤波;
(3)对运动目标进行方位向压缩,完成杂波抑制。
2.根据权利要求1所述基于方位角度域滤波的弧形阵列雷达杂波抑制方法,其特征在于,在步骤(1)中,弧形阵列雷达采用调频连续波FMCW信号,其发射信号的形式如下:
上式中,ST(t)为发射信号,fc表示发射信号载频,K表示调频率,t表示采样时间,j为虚数单位;
对回波信号进行解调频Dechirp处理、补偿残余视频相位、距离频率替换和距离向IFFT后,得回波信号的形式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011484658.6A CN112731304B (zh) | 2020-12-16 | 2020-12-16 | 一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011484658.6A CN112731304B (zh) | 2020-12-16 | 2020-12-16 | 一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112731304A true CN112731304A (zh) | 2021-04-30 |
CN112731304B CN112731304B (zh) | 2024-04-02 |
Family
ID=75602733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011484658.6A Active CN112731304B (zh) | 2020-12-16 | 2020-12-16 | 一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112731304B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113820673A (zh) * | 2021-10-19 | 2021-12-21 | 中安锐达(北京)电子科技有限公司 | 一种基于伺服转速的机扫雷达方位波束中心的修正方法 |
CN114627205A (zh) * | 2022-03-31 | 2022-06-14 | 昆明理工大学 | 一种煤层平行电磁波束地质探测的多尺度反投影滤波及可视化解释图像重建算法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1580573A1 (en) * | 2002-03-13 | 2005-09-28 | Raytheon Canada Limited | System and method for noise suppression in pre-processed radar data |
CN104237860A (zh) * | 2014-09-10 | 2014-12-24 | 电子科技大学 | 一种双基地mimo雷达跟踪定位及欺骗干扰识别方法 |
US20160061944A1 (en) * | 2014-08-29 | 2016-03-03 | Furuno Electric Co., Ltd. | Clutter suppressing device and radar apparatus provided with the same |
CN106872954A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 一种高超声速平台杂波抑制和运动目标成像方法 |
CN109633598A (zh) * | 2019-01-25 | 2019-04-16 | 西安电子科技大学 | 基于特征分析的阵列雷达目标检测方法 |
CN109946665A (zh) * | 2019-03-07 | 2019-06-28 | 西安电子科技大学 | 基于阵列雷达的获取真实目标的方法 |
-
2020
- 2020-12-16 CN CN202011484658.6A patent/CN112731304B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1580573A1 (en) * | 2002-03-13 | 2005-09-28 | Raytheon Canada Limited | System and method for noise suppression in pre-processed radar data |
US20160061944A1 (en) * | 2014-08-29 | 2016-03-03 | Furuno Electric Co., Ltd. | Clutter suppressing device and radar apparatus provided with the same |
CN104237860A (zh) * | 2014-09-10 | 2014-12-24 | 电子科技大学 | 一种双基地mimo雷达跟踪定位及欺骗干扰识别方法 |
CN106872954A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 一种高超声速平台杂波抑制和运动目标成像方法 |
CN109633598A (zh) * | 2019-01-25 | 2019-04-16 | 西安电子科技大学 | 基于特征分析的阵列雷达目标检测方法 |
CN109946665A (zh) * | 2019-03-07 | 2019-06-28 | 西安电子科技大学 | 基于阵列雷达的获取真实目标的方法 |
Non-Patent Citations (3)
Title |
---|
ZHENG MJ 等: "Research on Strong Clutter Suppression for Gaofen-3 Dual-Channel SAR/GMTI", 《SENSORS》, vol. 18, no. 4, pages 1 - 15 * |
ZHIQI GAO 等: "Arc Array Radar IAA-STAP Algorithm Based on Sparse Constraint", 《2020 INTERNATIONAL CONFERENCE ON INFORMATION SCIENCE, PARALLEL AND DISTRIBUTED SYSTEMS (ISPDS)》, pages 277 - 282 * |
项建弘 等: "基于虚拟天线的自适应波束形成零陷改善方法", 《无线电工程》, vol. 48, no. 11, pages 953 - 959 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113820673A (zh) * | 2021-10-19 | 2021-12-21 | 中安锐达(北京)电子科技有限公司 | 一种基于伺服转速的机扫雷达方位波束中心的修正方法 |
CN114627205A (zh) * | 2022-03-31 | 2022-06-14 | 昆明理工大学 | 一种煤层平行电磁波束地质探测的多尺度反投影滤波及可视化解释图像重建算法 |
Also Published As
Publication number | Publication date |
---|---|
CN112731304B (zh) | 2024-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107121670B (zh) | 一种基于合成孔径雷达的反无人机防御方法 | |
Zhao et al. | Multipath clutter rejection for digital radio mondiale-based HF passive bistatic radar with OFDM waveform | |
CN109932716B (zh) | 一种低空目标微多普勒特征提取方法 | |
CN112014806B (zh) | 一种复杂干扰场景下的机载雷达无意干扰抑制方法 | |
CN112881982B (zh) | 一种捷变频雷达抑制密集转发干扰方法 | |
CN108196241B (zh) | 一种基于Hough变换的高速动目标速度估计方法 | |
CN111551925B (zh) | 基于fft的脉间频率捷变雷达的目标速度估计方法 | |
CN112731304A (zh) | 一种基于方位角度域滤波的弧形阵列雷达杂波抑制方法 | |
Zuo et al. | UAV detection via long-time coherent integration for passive bistatic radar | |
Zhao et al. | Reduced complexity multipath clutter rejection approach for DRM-based HF passive bistatic radar | |
CN113884989A (zh) | 一种基于ofdm信号的外辐射源雷达时变杂波抑制方法 | |
CN110632573A (zh) | 一种机载宽带雷达空时二维keystone变换方法 | |
CN105005032A (zh) | 一种基于级数反演成像算法的sar移频干扰方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN113625273A (zh) | 混叠数字信号合成孔径定位方法 | |
CN113885002B (zh) | 基于尺度变换的外辐射源雷达捷变频信号相参积累方法 | |
Root | HF radar ship detection through clutter cancellation | |
CN104391288B (zh) | 拟合干涉相位的虚警剔除方法 | |
CN113359131A (zh) | 一种sar低截获射频隐身***及其设计方法 | |
CN109143235A (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN116068498A (zh) | 基于主副瓣空域极化特性的多干扰源抑制方法 | |
CN112710999B (zh) | 基于径向速度搜索的弧形阵列雷达运动目标聚焦检测方法 | |
CN112630744B (zh) | 一种多相参积累方法融合的海上小目标检测方法及*** | |
CN114397649A (zh) | 一种基于子载波的目标匹配方法 | |
Chen et al. | Robust mismatched filtering algorithm for passive bistatic radar using worst-case performance optimization |
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 |