CN107102307A - 基于均匀相位分集体制的极化雷达目标检测方法 - Google Patents
基于均匀相位分集体制的极化雷达目标检测方法 Download PDFInfo
- Publication number
- CN107102307A CN107102307A CN201710443094.3A CN201710443094A CN107102307A CN 107102307 A CN107102307 A CN 107102307A CN 201710443094 A CN201710443094 A CN 201710443094A CN 107102307 A CN107102307 A CN 107102307A
- Authority
- CN
- China
- Prior art keywords
- polarization
- radar
- target
- spectrum
- interference
- 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
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/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
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
技术领域
本发明属于极化雷达探测技术领域,主要涉及雷达天线主瓣内同时存在多个射频干扰或主瓣内存在变极化射频干扰情况下的极化雷达目标检测问题。
背景技术
极化雷达一般指的是具有双极化或多极化天线,并能够快速改变发射、接收极化状态的雷达。由于极化雷达能够通过使其接收极化匹配于雷达目标回波极化或正交于干扰信号极化,从而提升信干噪比(信号与干扰、噪声的功率比),因此,极化雷达具有更好的目标检测能力。
极化雷达的目标检测性能优势在单个主瓣干扰情况下表现的尤为突出,其基本原理是通过极化对消抑制主瓣干扰从而显著提升检测前的信干噪比。然而,当主瓣同时存在多个干扰时其检测性能迅速下降。
随着民用射频环境的日益复杂和雷达工作频段的不断延伸,雷达受到同频段民用射频干扰的威胁越来越大,其中与雷达目标方向临近的射频干扰会在雷达天线主瓣内形成多个极化不同、相关性较弱的干扰,使得合成干扰信号的极化度迅速降低,导致极化雷达干扰抑制和目标检测能力严重下降。特别是当主瓣内存在多个空间紧邻的干扰或主瓣内存在极化状态随机变化的干扰(称为变极化干扰)时,现有方法很难有效检测目标。
主瓣同时存在多个空间紧邻干扰或变极化干扰背景下的雷达目标检测问题,是雷达界当前面临的难点问题,还未有公开文献报道过相关技术。
发明内容
本发明提出了一种极化雷达目标检测方法,这种方法采用均匀相位分集的雷达发射极化分集体制,并基于极化状态分布特征差异,能够在主瓣内同时存在多个空间紧邻干扰或变极化干扰背景下有效检测雷达目标。
本发明的技术方案是,雷达采用正交双极化阵列接收天线和均匀相位分集的发射极化分集体制,对每一个PRI(Pulse Repeated Interval,脉冲重复周期)的接收信号生成极化-空间联合谱;然后将所有极化-空间联合谱叠加形成合成谱。对合成谱取模后,采用四维单元平均恒虚警检测(Cell-Average Constant False Alarm Rate,CA-CFAR)和点迹合并方法得到疑似目标点。对合成谱在方位-俯仰空间中进行方位-俯仰网格划分,计算每个方位-俯仰网格内所有疑似目标点的极化状态距离标准差特征量,并与检测门限进行比较,判定疑似目标是否为雷达目标。
本发明的技术效果:
本发明采用均匀相位分集的方式实现雷达发射极化分集,进而利用目标极化状态随发射极化有规律变化而干扰极化状态不服从该规律的差异,通过计算合成谱,并且利用合成谱的极化分布特征,有效地判别目标的“存在”,因此本发明能够显著提升主瓣同时存在多个紧邻干扰或主瓣存在变极化干扰情况下的雷达目标检测能力。
附图说明
图1为本发明所提方法的信号处理流程示意图;
图2为利用本发明所提方法计算得到的目标、主瓣多个紧邻干扰、主瓣变极化干扰的“极化状态距离标准差”检测特征量与检测门限对比图。
具体实施方式
以下结合附图对本发明的实施方式作进一步说明。
本发明的技术方案是,一种基于均匀相位分集体制的极化雷达目标检测方法,极化雷达采用正交双极化阵列接收天线,不失一般性,设为H(Horizontal,水平)极化、V(Vertical,垂直)极化,设H极化阵列和V极化阵列的阵元数均为K,阵元间距不大于极化雷达发射信号半波长;极化雷达采用均匀相位分集的发射极化分集体制,每个PRI发射极化的相位描述子为(γn,ηn),其中 n=1,…,N,N为发射极化分集数,其特征在于,包括以下三个步骤:
第一步,形成合成谱。
首先,对第n个PRI对应的双极化阵列接收信号进行极化-空间联合谱估计,可以利用现有的任意一种方法进行极化-空间联合谱估计。优先使用极化-空间MUSIC谱估计方法,参见文献:徐有根,刘志文,龚晓峰,极化敏感阵列信号处理[M],北京:北京理工大学出版社,2013。设得到的极化-空间联合谱Pn(θ,φ,γ,η),其中θ为俯仰角,φ为方位角,极化角γ表示相对幅度信息,极化角η表示相对相位信息,极化角γ和极化角η一起构成极化矢量相位描述子。
将N个PRI对应的极化-空间联合谱叠加,得到合成谱P(θ,φ,γ,η):
第二步,提取疑似目标。
本步骤包括两个主要环节,一是进行恒虚警检测,二是进行点迹合并。具体处理过程如下:
对合成谱的幅度谱|P(θ,φ,γ,η)|(“| |”表示取模)利用四维CA-CFAR进行检测,保护单元长度和检测门限值根据具体情况确定。四维CA-CFAR检测方法是一维CA-CFAR检测方法的简单拓展,即在(θ,φ,γ,η)构成的四维空间中,在待检测点周围选择与其在四个坐标轴上均相距保护单元长度的点。在四维空间中,存在16个点,与待检测点的距离在每个坐标轴上均相距保护单元长度。将16个点的幅度值进行平均作为背景参考电平。
将四维CA-CFAR检测的过门限点在四维空间进一步进行点迹合并处理,得到疑似目标点。四维空间的点迹合并处理是一维空间点迹合并处理的简单拓展,一维空间点迹合并处理主要方法是雷达从一维坐标轴上对连续多个过门限点进行选择并只输出其中一个过门限点,一般选择的是连续过门限点中的最大值点,本发明中,四维空间点迹合并处理是对方位角φ、俯仰角θ、极化角γ、极化角η构成四维空间上的连续过门限点选择最大值点进行输出。
设经过上述点迹合并处理后得到的疑似目标共有I个,设第i个疑似目标在合成谱P(θ,φ,γ,η)上的坐标位置(θi,φi,γi,ηi),i=1,…,I。
第三步,确认检测。
确认检测通过计算“极化状态距离标准差”检测特征量并与检测门限进行比较,判定疑似目标是否为极化雷达目标。具体处理过程如下:
在(θ,φ,γ,η)四维空间的俯仰角θ轴和方位角φ轴上,以一定的角度间隔划分方位-俯仰网格(该角度间隔根据角估计精度、系统可忍受计算量共同决定,如果角估计精度较高、可忍受计算量较大,则角度间隔可以取得更小)。
计算落入每个方位-俯仰网格内的疑似目标点数目,对数目大于或等于N的网格进一步进行如下的确认检测处理:
设落入方位-俯仰网格(θj,φl)中的疑似目标点集合为Cj,l,计算集合Cj,l中所有疑似目标点在极化角η轴上的极化状态距离标准差STDj,l,求取公式如下:
其中Δηi=ηi+1-ηi(i=1,…,Ij,l-1),ηi为疑似目标点极化角η的值,为集合Cj,l中所有Δηi的均值,Ij,l为集合Cj,l中疑似目标点的数目。
STDj,l表征了在网格(θj,φl)中疑似目标点极化状态距离的起伏情况,即当接收信号来自于雷达目标回波时,极化角η具有接近于等间距变化的特点,对STDj,l而言,则表现为取值较小,而当接收信号来自于多个紧邻干扰或变极化干扰信号时,极化角η不具有等间距变化的特点,对STDj,l而言,则表现为取值较大。
因此,将STDj,l作为检测统计量:当STDj,l小于检测门限Lstd时,判为目标;当STDj,l超过检测门限Lstd时,判为干扰。
其中,检测门限Lstd的设置主要由极化雷达所关心目标的极化散射矩阵特性、分集数目决定,随目标极化散射矩阵的两个共极化分量之间相对差异增大而增大、随交叉极化分量增大而增大。检测门限Lstd通过数值仿真得到:设目标散射矩阵为根据关心目标类型,设定目标散射矩阵两个共极化分量幅值比和交叉极化分量相对值的取值范围,在给定发射极化分集数目N的情况下,在和的取值范围中通过数值仿真计算极化状态距离标准差STDj,l,找到STDj,l的最大值作为检测门限Lstd。
采用本发明进行仿真实验,参数设置如下:极化雷达发射极化分集数目N=4,相对相位分别取值为对应的发射极化相位描述子分别为(45°,0°),(45°,30°),(45°,60°),(45°,90°);目标散射矩阵的取值范围为根据数值仿真可得检测门限为Lstd=8.955;方位-俯仰网格的角度间隔均设为1度。
仿真中,共设置一个目标、三个紧邻干扰(干扰1、干扰2、干扰3)和一个变极化干扰(干扰4),目标和所有干扰的俯仰角均为0度,方位角和极化状态分别为:目标方位角40度;干扰1方位角为37.5度,极化状态为(90°,30°);干扰2方位角为38度,极化状态为(45°,60°);干扰3方位角为38.5度,极化状态为(45°,30°);干扰4方位角为42度,极化状态为随机变极化。所有干扰信号波形为相互独立的噪声压制干扰,每个干扰信号的干信比(即单个噪声压制干扰源的功率与目标回波功率的比值)均为30dB。
图2为利用本发明所提方法计算得到的目标、主瓣多个紧邻干扰、主瓣变极化干扰的“极化状态距离标准差”检测特征量与检测门限对比图。图中,横轴为蒙特卡洛仿真序号(总次数为1000),纵轴为“极化状态距离标准差”检测特征量,图中,“△”表示目标、“+”表示紧邻干扰、“*”表示变极化干扰、“-”表示判决门限。由图2可以看出,目标明显位于判决门限(图中所示纵坐标为8.955的黑实线)以下,而紧邻干扰、变极化干扰则位于判决门限以上,这表明该方法能够较好地区分目标与主瓣内的紧邻干扰、变极化干扰,且检测门限设置较为合理。
Claims (1)
1.一种极化雷达目标检测方法,其特征在于,雷达采用正交双极化阵列接收天线和均匀相位分集的发射极化分集体制,对每一个脉冲重复周期的接收信号生成极化-空间联合谱;然后将所有极化-空间联合谱叠加形成合成谱;对合成谱取模后,采用四维单元平均恒虚警检测和点迹合并方法得到疑似目标点;对合成谱在方位-俯仰空间中进行方位-俯仰网格划分,计算每个方位-俯仰网格内所有疑似目标点的极化状态距离标准差特征量,并与检测门限进行比较,判定疑似目标是否为雷达目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443094.3A CN107102307B (zh) | 2017-06-13 | 2017-06-13 | 基于均匀相位分集体制的极化雷达目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443094.3A CN107102307B (zh) | 2017-06-13 | 2017-06-13 | 基于均匀相位分集体制的极化雷达目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107102307A true CN107102307A (zh) | 2017-08-29 |
CN107102307B CN107102307B (zh) | 2019-06-07 |
Family
ID=59660701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710443094.3A Active CN107102307B (zh) | 2017-06-13 | 2017-06-13 | 基于均匀相位分集体制的极化雷达目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107102307B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375761A (zh) * | 2018-02-08 | 2018-08-07 | 电子科技大学 | 针对多发多收雷达系统的单目标异步信号检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941239A (zh) * | 2014-05-08 | 2014-07-23 | 西安电子科技大学 | 网络化雷达极化对消抑制主瓣压制式干扰方法 |
CN104376330A (zh) * | 2014-11-19 | 2015-02-25 | 西安电子科技大学 | 基于超像素散射机制的极化sar图像舰船目标检测方法 |
CN105403875A (zh) * | 2015-12-31 | 2016-03-16 | 中国人民解放军国防科学技术大学 | 双极化接收雷达的目标检测方法 |
-
2017
- 2017-06-13 CN CN201710443094.3A patent/CN107102307B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941239A (zh) * | 2014-05-08 | 2014-07-23 | 西安电子科技大学 | 网络化雷达极化对消抑制主瓣压制式干扰方法 |
CN104376330A (zh) * | 2014-11-19 | 2015-02-25 | 西安电子科技大学 | 基于超像素散射机制的极化sar图像舰船目标检测方法 |
CN105403875A (zh) * | 2015-12-31 | 2016-03-16 | 中国人民解放军国防科学技术大学 | 双极化接收雷达的目标检测方法 |
Non-Patent Citations (4)
Title |
---|
BO REN等: "Statistical properties of the polarisation ratio for dual-polarisation radar operating at simultaneous reception mode", 《IET RADAR, SONAR & NAVIGATION》 * |
XU CHENG等: "Novel polarimetric detector for target detection in heterogeneous clutter", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 * |
刘勇等: "有源压制干扰背景下的雷达目标空域极化检测", 《电波科学学报》 * |
张莹等: "干扰条件下极化雷达目标检测性能分析", 《电子信息对抗技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375761A (zh) * | 2018-02-08 | 2018-08-07 | 电子科技大学 | 针对多发多收雷达系统的单目标异步信号检测方法 |
CN108375761B (zh) * | 2018-02-08 | 2020-04-07 | 电子科技大学 | 针对多发多收雷达系统的单目标异步信号检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107102307B (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11714199B2 (en) | System and method for detecting spoofing of GNSS signals | |
CN101533091A (zh) | 空时二维阻塞窄带干扰方法 | |
CN103728597A (zh) | 基于辅助阵列的网络化雷达抑制压制式主瓣干扰方法 | |
CN109490845B (zh) | 多站雷达抑制主瓣压制式干扰的方法 | |
CN102298140A (zh) | 基于探测概率几何分布曲线的雷达网抗干扰效能评估方法 | |
CN104330787A (zh) | 水下运动阵列多目标检测和方位估计一体化方法 | |
CN103278818B (zh) | 基于相位统计参数加权的穿墙雷达成像方法 | |
Li et al. | Combining sum-difference and auxiliary beams for adaptive monopulse in jamming | |
CN104977585A (zh) | 一种鲁棒的运动声纳目标检测方法 | |
CN109143265B (zh) | 一种基于空间相关性识别的gnss抗欺骗方法 | |
CN103728596B (zh) | 双基地mimo雷达抑制分布式干扰的方法 | |
CN110967673A (zh) | 一种多域联合抗主瓣干扰方法 | |
CN103558584A (zh) | 一种波达方向的检测方法及装置 | |
CN114089288B (zh) | 一种相控阵雷达抗干扰方法、装置及存储介质 | |
Xianrong et al. | Adaptive ionospheric clutter suppression based on subarrays in monostatic HF surface wave radar | |
CN108919205B (zh) | 基于辅助阵元的频率分集阵列欺骗干扰抑制方法 | |
CN107102307A (zh) | 基于均匀相位分集体制的极化雷达目标检测方法 | |
Son et al. | Passive polarimetrie multistatic radar for ground moving target | |
CN104536014A (zh) | 一种基于Relax算法的抗转发式欺骗干扰方法 | |
Wu et al. | Performance analysis of polarization-space-time three-domain joint processing for clutter suppression in airborne radar | |
Wei et al. | An empirical exploration of a simple deghosting method for multistatic radars | |
Cai et al. | A method for countering dense false target jamming based on correlation sample selection | |
Cattenoz et al. | Adaptive processing methods for MIMO radar experimental signals | |
Yang et al. | A GPS Spoofing Detection Method Based on Compressed Sensing | |
Yin et al. | Range-angle two-dimensional anti jamming against dense false targets |
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 |