CN106019244A - 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 - Google Patents
基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 Download PDFInfo
- Publication number
- CN106019244A CN106019244A CN201610443068.6A CN201610443068A CN106019244A CN 106019244 A CN106019244 A CN 106019244A CN 201610443068 A CN201610443068 A CN 201610443068A CN 106019244 A CN106019244 A CN 106019244A
- Authority
- CN
- China
- Prior art keywords
- sequence
- shift
- negative slope
- target
- lfm waveform
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000010287 polarization Effects 0.000 title abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 238000010606 normalization Methods 0.000 claims abstract description 4
- 238000005388 cross polarization Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000000205 computational method Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims 1
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008485 antagonism Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012850 discrimination method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
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
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,包含以下步骤:(1)正交双极化雷达采用正负斜率LFM波形发射信号并接收回波序列,对序列进行匹配滤波并归一化,得到HH、HV、VH以及VV四个通道的输出序列,并将HV通道输出序列和VV通道输出序列进行时域上的反转;(2)对四个通道输出序列进行处理,得到序列一;(3)对序列一进行目标检测,得到待定目标位置,并从序列一中提取得到序列二;(4)根据待定目标位置,对四个通道输出序列进行处理,提取得到序列三;(5)计算序列三与序列二之比作为鉴别量;(6)通过鉴别量与门限值的对比,判断目标和移频干扰。本发明对匹配滤波后信号的处理不涉及复杂的运算,实时性较强。
Description
技术领域
本发明属于电子侦查领域,涉及一种基于正负斜率LFM波形正交极化雷达的抗移频干扰方法。
背景技术
移频干扰是目前对付线性调频(LFM)脉冲压缩雷达的一种有效手段。由于线性调频脉冲压缩雷达信号的模糊函数存在距离-多普勒频移耦合现象,这种距离与多普勒频移之间的相互影响会使得当目标与匹配滤波器存在多普勒失配时,对距离的估计也会产生一定的误差。对抗移频干扰的研究目前主要集中在采用时变的发射信号,如改变发射信号的调频斜率、初始相位等,但是这种方法会增加发射机的发射难度。除此之外还可以采用时频分析方法,但是对信噪比的要求比较高,并且有可能会出现信号混叠的情况。以上方法都需要相当大的计算量,无法保证目标检测的实时性。随着极化信息处理技术的发展,极化雷达的研制也越来越受到重视,其在目标监视和特征提取等领域的应用日益广泛。由于干扰同目标有着不同的极化特征,因此利用极化信息可以将目标同干扰区分出来。
发明内容
本发明要解决的技术问题是:针对当前电子侦查过程中所接收的线性调频脉冲压缩雷达信号中可能存在的移频干扰,利用正交极化雷达,根据干扰同真实目标回波中包含的不同的极化信息,对回波进行处理,从真实目标中鉴别出移频干扰。
本发明解决其技术问题所采用的技术方案是:基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,包括以下步骤:
步骤(1)、正交双极化雷达采用正负斜率LFM波形发射信号并接收回波序列,对序列进行匹配滤波并归一化,得到HH、HV、VH以及VV四个通道的输出序列,并将HV通道输出序列和VV通道输出序列进行时域上的反转,其中H和V分别表示水平和垂直极化;
步骤(2)、对四个通道输出序列进行处理,得到序列一;
步骤(3)、对序列一进行目标检测,得到待定目标位置,并从序列一中提取得到序列二;
步骤(4)、根据待定目标位置,对四个通道输出序列进行处理,提取得到序列三;
步骤(5)、计算序列三与序列二之比作为鉴别量;
步骤(6)、通过鉴别量与门限值的对比,判断目标和移频干扰。
所述步骤(2)中序列一提取方法如下所示:
其中,rHH(n)和rVH(n)分别为HH与VH通道的输出序列,为HV和VV通道的输出序列经时域反转后的结果;|·|表示对序列取绝对值,n为采样点。
所述步骤(3)中序列二提取方法如下所示:
v2(k)=v1(nk)
其中,序列二v2(k)为m个待定目标位置点组成的序列,即{nk|v1(nk)≥th1,k=0,1,...,m},th1为目标检测时所用门限值。
所述步骤(4)中序列三提取方法如下所示:
所述步骤(5)中鉴别量计算方法如下式:
ζ(k)=v3(k)/v2(k)
本发明与现有技术相比的优点在于:
(1)本发明使用的正交极化雷达采用了正负斜率LFM波形,其波形正交性满足仅发射一个脉冲即可分辨真假目标的要求;
(2)本发明对匹配滤波后信号的处理不涉及复杂的运算,实时性较强。
附图说明
图1为本发明基于正负斜率LFM波形正交极化雷达的抗移频干扰方法流程图;
图2为移频干扰叠加到目标回波上经匹配滤波处理后的仿真图;
图3为输出结果经步骤(2)处理后得到的序列一的仿真图。
具体实施方式
下面结合附图及具体实施方式详细介绍本发明。
干扰机的移频干扰通常可以通过数字射频存储器(Digital Radio Frequency Memory,DRFM)技术来实现。DRFM首先存储干扰机接收到的雷达发射信号,随后复制出原雷达信号,通过给上变频本振调制一个额外的频移,从而产生移频干扰信号。然而干扰信号的极化信息 同目标回波的极化信息有着很明显的不同,使用极化雷达发射并接收目标与干扰信号,经处理后可较易区分两种信号,本发明为基于同时极化测量体制的正交极化雷达对移频干扰的鉴别方法,其实施流程图如图1所示,具体包含如下4个步骤:
1、匹配滤波
正交极化雷达经由正交极化通道独立发射两路正交信号,然后同时接收目标回波的正交极化分量。两路正交发射信号为eH(t)、eV(t),在本发明中为一对正负斜率线性调频脉冲信号,不考虑调制幅度,可以表示为:
其中,tp为脉冲宽度,K为调频斜率,f0为载频,两极化通道的发射信号除调频斜率互为正负外其余参数完全一致。
接收回波信号后,按照抽样定理对回波时间序列进行采样,得到两个通道的N点时间序列,可以将其表示为rH(n)、rV(n)。对回波序列rH(n)、rV(n)进行如下滤波处理:
Rij(n)=IFFT[FFT(ri(n))·(FFT(ej(n)))*],(i,j=H,V) (2)
其中,RHH(n)、RHV(n)、RVH(n)和RVV(n)分别为HH、HV、VH和VV通道输出序列,经过归一化以及HV与VV通道输出序列的时域反转:
rij(n)=Rij(n)/Ej (3)
式(3)中Ej为发射信号ej(n)的能量,可由下式计算得到:
其中,N为序列长度。
此时我们所得到的四个匹配滤波输出峰值即为目标极化散射矩阵的测量结果,若rij(n)中存在多个峰值则表明有多个目标,其中可能包含移频干扰造成的假目标,需要对其进行处理。
2、计算序列一
图1中P1表示对输出序列进行如下操作:
通过上式即可得到序列一v1(n)。
3、目标检测
对序列一v1(n)进行目标检测,门限值由下式计算得到:
th1=mean[v1(n)]+a·std[v1(n)] (7)
其中,mean(·)表示计算均值,std(·)表示计算标准差,a为门限系数,其取值范围可取[0.5,2]。目标检测过程为:
若共有m个点检测到有目标,则提取得到待定目标位置点{nk|v1(nk)≥th1,k=0,1,...,m}。
4、计算序列二
图1中P2操作为取序列一中待定目标位置点处的值得到序列二,即:
v2(k)=v1(nk) (9)
5、计算序列三
图1中P3表示取rHH(n)、rVH(n)和四个序列在待定目标位置点nk处的值,对其进行如下处理:
得到序列三。
6、提取鉴别量
鉴别量ζ(k)提取方法如下所示:
7、目标判断
将鉴别量ζ(k)与门限值进行对比并判断真假目标
其中,th2为门限值。通过设置合理的门限值,一般而言,对于人造小目标,选择th2为0.2-0.4即可进行真实目标与移频干扰的有效判断。
下面以仿真实验来说明该鉴别方法的有效性。仿真条件设置如下:雷达发射信号脉冲宽度tp=10μs,带宽B=10MHz,载波频率f0=10GHz;目标的多普勒频移为fd=20kHz,目 标极化散射矩阵为信噪比SNR=10dB;干扰取两个固定移频干扰,频移量分别为ζ1=-0.1B、ζ2=0.2B,干扰极化状态Jones矢量为 干噪比分别为INR1=13dB,INR2=15dB;噪声符合零均值高斯分布,两极化通道接收信号各参数均相等。
图2为目标回波加入移频干扰后经过步骤2的输出,可以看到移频干扰形成了两个较为明显的假目标,影响了真实目标的判断。图3为序列一的计算结果,门限系数设置为2,由图中可以看出能够检测到3个目标,待定目标位置点nk=3001,4980,6000。由此计算序列二、序列三以及鉴别量如下表所示
表1序列二、序列三与鉴别量
由上表可知待定目标1和待定目标3的鉴别量远远小于待定目标2,此时门限值可取th2=0.3,对上面三个待定目标进行鉴别,得到鉴别结果为:待定目标2为真实目标,符合仿真条件。
本发明根据目标回波与移频干扰中蕴含的不同的极化信息,利用正交极化雷达发射正负斜率LFM波形信号并接收回波,经过处理回波来进行针对移频干扰的鉴别与对抗。根据仿真结果可以看出,在同时存在目标与干扰且雷达对目标和干扰信号接收功率相当的情况下,分别提取鉴别量,目标鉴别量与干扰鉴别量通常差距较大,据此可以有效区分真实目标与移频干扰。
Claims (5)
1.一种基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,其特征在于步骤如下:
步骤(1)、正交双极化雷达采用正负斜率LFM波形发射信号并接收回波序列,对序列进行匹配滤波并归一化,得到HH、HV、VH以及VV四个通道的输出序列,并将HV通道输出序列和VV通道输出序列进行时域上的反转,其中H和V分别表示水平和垂直极化;
步骤(2)、对四个通道输出序列进行处理,得到序列一;
步骤(3)、对序列一进行目标检测,得到待定目标位置,并从序列一中提取得到序列二;
步骤(4)、根据待定目标位置,对四个通道输出序列进行处理,提取得到序列三;
步骤(5)、计算序列三与序列二之比作为鉴别量;
步骤(6)、通过鉴别量与门限值的对比,判断目标和移频干扰。
2.根据权利要求1的基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,其特征在于:所述步骤(2)中得到序列一v1(n)的处理方法如下所示:
v1(n)=|rHH(n)|2+|r′HV(n)|2+|rVH(n)|2+|r′VV(n)|2
其中,rHH(n)和rVH(n)分别为HH与VH通道的输出序列,r′HV(n)和r′VV(n)分别为HV和VV通道的输出序列经时域反转后的结果;|·|表示对序列取绝对值,n为采样点。
3.根据权利要求2的基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,其特征在于:所述步骤(3)中序列二v2(k)提取方法如下所示:
v2(k)=v1(nk)
其中,序列二v2(k)为m个待定目标位置点组成的序列,即{nk|v1(nk)≥th1,k=1,…,m},th1为目标检测时所用门限值。
4.根据权利要求3的基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,其特征在于:所述步骤(4)中序列三v3(k)提取方法如下所示:
v3(k)=|rHH(nk)·r′VV(nk)-rVH(nk)·r′HV(nk)|。
5.根据权利要求4的基于正负斜率LFM波形正交极化雷达的抗移频干扰方法,其特征在于:所述步骤(5)中鉴别量ζ(k)的计算方法如下式:
ζ(k)=v3(k)/v2(k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443068.6A CN106019244B (zh) | 2016-06-20 | 2016-06-20 | 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610443068.6A CN106019244B (zh) | 2016-06-20 | 2016-06-20 | 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019244A true CN106019244A (zh) | 2016-10-12 |
CN106019244B CN106019244B (zh) | 2018-05-15 |
Family
ID=57088868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610443068.6A Expired - Fee Related CN106019244B (zh) | 2016-06-20 | 2016-06-20 | 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019244B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257835A (zh) * | 2020-02-17 | 2020-06-09 | 森思泰克河北科技有限公司 | 雷达的干扰抑制方法及终端设备 |
CN113009465A (zh) * | 2021-03-05 | 2021-06-22 | 中国人民解放军海军航空大学 | 一种基于两次相位补偿的稳健自适应脉冲压缩方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006177671A (ja) * | 2004-12-20 | 2006-07-06 | Mitsubishi Electric Corp | 妨害電波キャンセラ装置 |
CN104281854A (zh) * | 2014-09-28 | 2015-01-14 | 中国人民解放军海军航空工程学院 | 高分辨率极化sar舰船目标与干扰目标鉴别方法 |
CN104765020A (zh) * | 2015-04-12 | 2015-07-08 | 西安电子科技大学 | 有源假目标干扰的极化鉴别方法 |
-
2016
- 2016-06-20 CN CN201610443068.6A patent/CN106019244B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006177671A (ja) * | 2004-12-20 | 2006-07-06 | Mitsubishi Electric Corp | 妨害電波キャンセラ装置 |
CN104281854A (zh) * | 2014-09-28 | 2015-01-14 | 中国人民解放军海军航空工程学院 | 高分辨率极化sar舰船目标与干扰目标鉴别方法 |
CN104765020A (zh) * | 2015-04-12 | 2015-07-08 | 西安电子科技大学 | 有源假目标干扰的极化鉴别方法 |
Non-Patent Citations (3)
Title |
---|
WENGUANG WANG,YU JI,XIAOXIA LIN: "A Novel Fusion-Based Ship Detection Method from Pol-SAR Images", 《SENSOR》 * |
周利华: "雷达干扰信号分析技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李永祯 等: "有源转发式干扰的全极化单脉冲雷达抑制方法研究", 《电子与信息学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257835A (zh) * | 2020-02-17 | 2020-06-09 | 森思泰克河北科技有限公司 | 雷达的干扰抑制方法及终端设备 |
CN111257835B (zh) * | 2020-02-17 | 2022-02-18 | 森思泰克河北科技有限公司 | 雷达的干扰抑制方法及终端设备 |
CN113009465A (zh) * | 2021-03-05 | 2021-06-22 | 中国人民解放军海军航空大学 | 一种基于两次相位补偿的稳健自适应脉冲压缩方法 |
CN113009465B (zh) * | 2021-03-05 | 2022-05-27 | 中国人民解放军海军航空大学 | 一种基于两次相位补偿的稳健自适应脉冲压缩方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106019244B (zh) | 2018-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104765020B (zh) | 有源假目标干扰的极化鉴别方法 | |
Li et al. | Fast coherent integration for maneuvering target with high-order range migration via TRT-SKT-LVD | |
CN106908771B (zh) | 微波暗室中脉冲体制雷达目标探测模拟方法 | |
CN106353748B (zh) | 用于fmcw雷达测距系统的信号处理装置及方法 | |
CN105403875B (zh) | 双极化接收雷达的目标检测方法 | |
CN106125053A (zh) | 脉冲多普勒雷达极化抗干扰方法 | |
CN105974376B (zh) | 一种sar射频干扰抑制方法 | |
CN105807267A (zh) | 一种mimo雷达扩展目标的检测方法 | |
CN106405517B (zh) | 针对脉冲多普勒雷达的微动假目标生成方法 | |
CN106249212A (zh) | 主瓣压制干扰背景下有源假目标的极化鉴别方法 | |
CN106338715B (zh) | 一种全极化有源假目标的鉴别方法 | |
CN110376559B (zh) | 单通道雷达主瓣多源干扰分离方法、装置和设备 | |
CN106019243A (zh) | 一种基于三次初相和frft的drfm干扰检测和对抗方法 | |
CN106168661B (zh) | 随机捷变极化有源假目标的鉴别方法 | |
CN114578296B (zh) | 基于相位编码信号的间歇采样干扰抑制方法 | |
CN106019244A (zh) | 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法 | |
CN106249209A (zh) | 一种抗速度欺骗干扰的自适应迭代估计方法 | |
CN105954729B (zh) | 雷达抗距离-速度联合欺骗干扰的自适应迭代滤波方法 | |
Chen et al. | Airborne SAR suppression of blanket jamming based on second order blind identification and fractional order Fourier transform | |
Kim et al. | Extraction of micro‐Doppler characteristics of drones using high‐resolution time‐frequency transforms | |
CN104237864B (zh) | 基于匹配模糊函数的距离扩展目标检测方法 | |
CN106338714B (zh) | Pri间恒定极化有源假目标鉴别方法 | |
Ganveer et al. | SAR implementation using LFM signal | |
CN108055218B (zh) | 基于频率选择性编码多载波调制的泄漏电缆入侵探测方法 | |
CN108983189B (zh) | 一种振动目标的二维微动轨迹估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180515 |