CN106168661A - 随机捷变极化有源假目标的鉴别方法 - Google Patents

随机捷变极化有源假目标的鉴别方法 Download PDF

Info

Publication number
CN106168661A
CN106168661A CN201610707828.XA CN201610707828A CN106168661A CN 106168661 A CN106168661 A CN 106168661A CN 201610707828 A CN201610707828 A CN 201610707828A CN 106168661 A CN106168661 A CN 106168661A
Authority
CN
China
Prior art keywords
polarization
signal
subpulse
target
channel
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
Application number
CN201610707828.XA
Other languages
English (en)
Other versions
CN106168661B (zh
Inventor
施龙飞
宗志伟
李永祯
王雪松
肖顺平
杨勇
马佳智
毛楚乔
胥文泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201610707828.XA priority Critical patent/CN106168661B/zh
Publication of CN106168661A publication Critical patent/CN106168661A/zh
Application granted granted Critical
Publication of CN106168661B publication Critical patent/CN106168661B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

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

本发明提供一种随机捷变极化有源假目标的鉴别方法。技术方案是:雷达包括H极化发射通道和V极化发射通道,两个通道的发射信号均由M个无缝连接的子脉冲构成;H通道、V通道接收信号经过同步接收后,首先按照子脉冲波形进行匹配滤波,对其输出信号采用M点滑窗方法实现非相参积累,进而采用恒虚警检测方法实现目标检测,对检测出的采样位置提取目标极化特征后,利用该特征进行真、假目标属性的鉴别。本发明的技术效果是有源假目标鉴别能力强,并且功率利用率高。

Description

随机捷变极化有源假目标的鉴别方法
技术领域
本发明属于雷达抗干扰技术领域,主要涉及雷达对有源假目标干扰的鉴别方法。
背景技术
有源假目标干扰是一种重要的雷达欺骗干扰,通过模拟目标特征并发射假目标信号,使雷达出现虚假目标以扰乱情报雷达对空情态势的掌握,或破坏雷达目标跟踪以扰乱跟踪制导雷达对目标的跟踪。
先进的有源假目标干扰能够做到与目标回波信号在波形调制、重复周期、多普勒频率、航迹以及RCS(Radar Cross Section,雷达散射截面积)起伏特性等方面几乎完全一致,雷达难以利用时域、频域、空域、调制域的特征差异来鉴别并剔除有源假目标干扰。由于有源假目标与雷达目标在极化域普遍存在着物理差异,基于极化域特征差异鉴别有源假目标干扰已成为一条重要的技术途径。
近年来,随着干扰技术的发展,已出现各种变极化假目标干扰,随机捷变极化假目标干扰就是其中十分重要的一种,其在PRI(Pulse Repeated Interval,脉冲重复周期)间或PRI内的极化状态随机捷变而在脉冲持续时间内极化状态恒定。
对于随机捷变极化假目标干扰,通常要求雷达具有同时发射多极化的能力,即雷达通过同时发射多个极化状态的矢量信号,将雷达目标回波在脉冲持续时间内的极化多样性与假目标信号在脉冲持续时间内的极化单一性进行区分,实现假目标的鉴别。现有报道中对于随机捷变极化假目标干扰的鉴别方法,都是基于严格意义上的同时发射多极化体制,即雷达在两个正交极化通道同时发射多个波形相互正交的信号,通过正交化接收处理实现对不同极化回波的隔离。这种体制在工程应用中存在造价高、发射功率利用率低的问题。本发明基于相位分集准同时变极化体制,通过子脉冲间极化相位的捷变实现有限的发射变极化能力,既能够实现假目标鉴别,又能够以饱和功率发射而具有功率利用率高的特点。基于该体制的有源假目标鉴别方法,还未见有文献报道。
发明内容
本发明针对随机捷变极化假目标,提出了一种基于相位分集准同时变极化体制的假目标鉴别方法。
本发明的技术方案是,一种基于相位分集准同时变极化体制的随机捷变极化假目标鉴别方法,其特征在于:
(1)雷达采用正交双极化天线,不失一般性,设为H(Horizontal,水平)极化、V(Vertical,垂直)极化,雷达的H极化通道、V极化通道同时发射、同时接收信号;
(2)H极化通道、V极化通道发射信号均由M个紧密连接的子脉冲构成,子脉冲的波形均为s(t)、宽度均为Δt,M为子脉冲个数,一般大于2;
(3)H极化通道、V极化通道发射信号在每个子脉冲上的相对相位差为 m为子脉冲序号,m=1,…,M。
综上,H极化通道、V极化通道发射信号的合成矢量的表达式为:
其中分别表示H极化和V极化发射矢量。
本发明的信号处理过程包括以下五个步骤:
第一步,信号同步接收。
雷达接收的中频信号经过同步采样,得到两个极化通道的复信号数据序列xH(n)、xV(n),n=1,…,N,N为一个脉冲周期内的采样点数。
第二步,子脉冲匹配滤波。
依据下式分别计算H极化通道、V极化通道匹配滤波输出信号yH(n)、yV(n):
y H ( n ) = I F F T [ F F T [ x H ( n ) ] · U ( n ) ] y V ( n ) = I F F T [ F F T [ x V ( n ) ] · U ( n ) ] - - - ( 1 )
上式中,FFT[]表示信号的傅里叶变换,IFFT[]表示信号的傅里叶反变换;U(n)是对子脉冲信号s(t)的基准波形进行N点傅里叶变换得到的频谱。
第三步,滑窗积累与检测。
对yH(n)、yV(n)进行滑窗积累,得到信号z(n)的过程如下:
z ( n ) = Σ m = 0 M - 1 { | y H ( n + m · Δ N ) | + | y V , m ( n + m · Δ N ) | } , n = 1 , ... , N - Δ N · ( M - 1 )
其中,ΔN为子脉冲对应的采样点数。
对信号z(n)进行目标检测,由CFAR(Constant false alarm rate,恒虚警检测)完成,具体采用CA-CFAR(Cell averaging-Constant false alarm rate,简称单元平均恒虚警检测)。
设对z(n)中所有点进行CA-CFAR处理后,共得到K个目标,第k个目标对应的位置为Tk,k=1,…,K
第四步,极化特征提取。
利用下式计算第k个目标第m个子脉冲回波的极化比ρm(Tk):
ρ m ( T k ) = y V ( T k + ( m - 1 ) * Δ N ) y H ( T k + ( m - 1 ) * Δ N ) , R m ( T k ) = r e a l [ ρ m ( T k ) ] , I m ( T k ) = i m a g [ ρ m ( T k ) ] , k = 1 , ... , K , m = 1 , ... , M - - - ( 2 )
其中,real[·]为取实部,imag[·]为取虚部。
对第k个目标,计算均值
R ‾ k = Σ m = 1 M R m ( T k ) , I ‾ k = Σ m = 1 M I m ( T k ) , m = 1 , ... , M
利用下式计算第k个目标对应的鉴别特征量
d R - k ( m ) = | R m ( T k ) - R ‾ k | , d I - k ( m ) = | I m ( T k ) - I ‾ k | , m = 1 , ... , M
第五步,极化特征鉴别。
首先确定鉴别门限。鉴别门限D由下式给出:
D = η · Mσ V 2 Σ m = 1 M | y H ( T k + ( m - 1 ) * Δ N ) | 2
其中,在假目标干扰情况下是{Rm(Tk)}、{Im(Tk)}的正态分布方差,σV 2为V通道接收噪声功率,通过对匹配滤波后纯噪声数据的方差估计得到;η一般取3,此时有源假目标正确鉴别概率不小于99.7%(由统计学上关于正态分布的“3σ准则”给出)。
其次,根据门限D对第k个目标进行鉴别(k=1,…,K):
1).当鉴别特征量{dR-k(m)}和{dI-k(m)}(m=1,…,M)均小于门限D时,判定第k个目标为有源假目标;
2).当鉴别特征量{dR-k(m)}或{dI-k(m)}(m=1,…,M)中存在大于门限D的情况,则判定第k个目标为真实雷达目标。
本发明的技术效果:
一、有源假目标鉴别能力强。本发明基于相位分集准同时变极化体制,通过子脉冲间极化相位的捷变,实现在PRI内的发射变极化,能够有效鉴别随机捷变极化假目标,而对于PRI间恒定极化假目标亦能够提供鉴别能力,因此,本文提出的方法适用性较广,能够应对当前面临的大多数有源假目标干扰。
二、功率利用率高。本发明所采用的相位分集准同时变极化体制,相对于同时变极化体制,具有更为简单的信号波形,在多极化情况下H极化、V极化通道发射信号的幅度包络仍然可以保持恒定,因而能够以饱和功率发射,具有更高的功率利用率。
附图说明
图1为本发明所依赖的相位分集准同时发射变极化体制的发射信号示意图。
图2为本发明所提方法的信号处理流程示意图;
图3为采用本发明所提方法后雷达接收信号鉴别特征量随发射极化的变化。
图4为利用本发明所提方法对实体雷达目标和随机捷变极化假目标进行鉴别处理的效果。
具体实施方式
以下结合附图对本发明的实施方式作进一步说明。
图1为本发明所依赖的相位分集准同时变极化体制的发射信号示意图。图中,H通道、V通道分别表示雷达的H极化发射通道和V极化发射通道,两个通道的发射信号均由M个无缝连接的子脉冲构成,子脉冲波形均为s(t)、宽度均为Δt,H通道、V通道发射信号在每个子脉冲上的相对相位差为m=1,…,M。
图2为本发明所提方法的信号处理流程示意图。如图所示,H通道、V通道接收信号经过同步接收后,首先按照子脉冲波形进行匹配滤波,对其输出信号采用M点滑窗方法实现非相参积累,进而采用恒虚警检测(CA-CFAR)方法实现目标检测,对检测出的采样位置提取目标极化特征后,利用该特征进行真、假目标属性的鉴别。
图3为采用本发明所提方法后雷达接收信号鉴别特征量随发射极化的变化。仿真中参数设置如下:
a).雷达采用相位分集准同时发射变极化体制,相位分集数、子脉冲个数M=4,相对相位分别取值为0、π、
b).雷达发射信号脉冲宽度为100us,带宽为1MHz,载频为10GHz;
c).干扰机发射天线极化状态的Jones矢量为
d).实体目标的极化散射矩阵设置为
图中,横轴有8个点,表示8种发射极化,纵轴为鉴别特征量{dR-k(m)}和{dI-k(m)}的数值;“△”表示{dR-k(m)},“○”表示{dI-k(m)};左图为实体目标情况,右图为假目标干扰情况。如图所示,实体目标鉴别特征量随发射极化变化而改变,表现为{dR-k(m)}、{dI-k(m)}均出现起伏,而假目标信号鉴别特征量不随发射极化变化而改变,表现为{dR-k(m)}、{dI-k(m)}均恒定,这种差异为极化鉴别提供了基础。
图4为利用本发明所提方法对实体雷达目标和假目标进行鉴别处理的效果,其仿真参数设置同图3。图中,横轴为“信噪比”(即子脉冲匹配滤波及滑窗积累后,目标回波功率或假目标干扰功率与噪声功率的比值),单位为dB,纵轴为被判定为“假目标”的概率;左图是实体目标情况,图中五条曲线分别表示五种目标极化散射矩阵情况下,被正确判定为“目标”的概率,由图可见,当“信噪比”大于10dB时,目标正确鉴别率都已高于99%;右图是假目标情况,“△”表示鉴别正确率仿真统计值,由图可见,在仿真设置的“信噪比”范围内,假目标的正确鉴别概率较高,均超过97%。

Claims (1)

1.一种基于相位分集准同时变极化体制的随机捷变极化假目标鉴别方法,其特征在于:
雷达采用正交双极化天线,设为H极化、V极化,雷达的H极化通道、V极化通道同时发射、同时接收信号;
H极化通道、V极化通道发射信号均由M个紧密连接的子脉冲构成,子脉冲的波形均为s(t)、宽度均为Δt,M为子脉冲个数且大于2;
H极化通道、V极化通道发射信号在每个子脉冲上的相对相位差为 m为子脉冲序号,m=1,…,M;
综上,H极化通道、V极化通道发射信号的合成矢量的表达式为:
其中分别表示H极化和V极化发射矢量;
本发明的信号处理过程包括以下五个步骤:
第一步,信号同步接收:
雷达接收的中频信号经过同步采样,得到两个极化通道的复信号数据序列xH(n)、xV(n),n=1,…,N,N为一个脉冲周期内的采样点数;
第二步,子脉冲匹配滤波:
依据下式分别计算H极化通道、V极化通道匹配滤波输出信号yH(n)、yV(n):
yH(n)=IFFT[FFT[xH(n)]·U(n)]
yV(n)=IFFT[FFT[xV(n)]·U(n)]
上式中,FFT[]表示信号的傅里叶变换,IFFT[]表示信号的傅里叶反变换;U(n)是对子脉冲信号s(t)的基准波形进行N点傅里叶变换得到的频谱;
第三步,滑窗积累与检测:
对yH(n)、yV(n)进行滑窗积累,得到信号z(n)的过程如下:
z ( n ) = Σ m = 0 M - 1 { | y H ( n + m · Δ N ) | + | y V , m ( n + m · Δ N ) | } , n = 1 , ... , N - Δ N · ( M - 1 )
其中,ΔN为子脉冲对应的采样点数;
利用恒虚警检测对信号z(n)进行目标检测,设共得到K个目标,第k个目标对应的位置为Tk,k=1,…,K;
第四步,极化特征提取:
利用下式计算第k个目标第m个子脉冲回波的极化比ρm(Tk):
Rm(Tk)=real[ρm(Tk)],Im(Tk)=imag[ρm(Tk)],
其中,real[·]为取实部,imag[·]为取虚部;
对第k个目标,计算均值
R ‾ k = Σ m = 1 M R m ( T k ) , I ‾ k = Σ m = 1 M I m ( T k )
利用下式计算第k个目标对应的鉴别特征量
d R - k ( m ) = | R m ( T k ) - R ‾ k | , d I - k ( m ) = | I m ( T k ) - I ‾ k |
第五步,极化特征鉴别:
由下式给出鉴别门限D:
D = η · Mσ V 2 Σ m = 1 M | y H ( T k + ( m - 1 ) * Δ N ) | 2
其中,σV 2为V通道接收噪声功率,通过对匹配滤波后纯噪声数据的方差估计得到;η根据实际要求确定;
根据门限D对第k个目标进行鉴别:
当鉴别特征量{dR-k(m)}和{dI-k(m)}均小于门限D时,判定第k个目标为有源假目标;
当鉴别特征量{dR-k(m)}或{dI-k(m)}中存在大于门限D的情况,判定第k个目标为真实雷达目标。
CN201610707828.XA 2016-08-23 2016-08-23 随机捷变极化有源假目标的鉴别方法 Active CN106168661B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610707828.XA CN106168661B (zh) 2016-08-23 2016-08-23 随机捷变极化有源假目标的鉴别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610707828.XA CN106168661B (zh) 2016-08-23 2016-08-23 随机捷变极化有源假目标的鉴别方法

Publications (2)

Publication Number Publication Date
CN106168661A true CN106168661A (zh) 2016-11-30
CN106168661B CN106168661B (zh) 2018-07-24

Family

ID=57376598

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610707828.XA Active CN106168661B (zh) 2016-08-23 2016-08-23 随机捷变极化有源假目标的鉴别方法

Country Status (1)

Country Link
CN (1) CN106168661B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107290727A (zh) * 2017-06-14 2017-10-24 中国人民解放军国防科学技术大学 一种位置可控式isar假目标图像合成方法
CN108919206A (zh) * 2018-07-18 2018-11-30 武汉大学 一种基于子载波处理的外辐射源雷达极化滤波方法
CN110794374A (zh) * 2019-11-08 2020-02-14 哈尔滨工业大学 一种间歇采样转发干扰的参数辨识方法
CN112799028A (zh) * 2020-12-14 2021-05-14 中电科仪器仪表有限公司 一种基于rcs起伏统计特性差异的虚假目标识别方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007163318A (ja) * 2005-12-14 2007-06-28 Toshiba Corp レーダ信号処理装置とその一定誤警報確率検出方法
CN103018722A (zh) * 2012-10-15 2013-04-03 西安电子科技大学 利用组网雷达系统对抗欺骗式假目标的方法
CN104765020A (zh) * 2015-04-12 2015-07-08 西安电子科技大学 有源假目标干扰的极化鉴别方法
CN105403875A (zh) * 2015-12-31 2016-03-16 中国人民解放军国防科学技术大学 双极化接收雷达的目标检测方法
CN105652273A (zh) * 2016-03-17 2016-06-08 哈尔滨工程大学 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法
CN105676203A (zh) * 2016-01-23 2016-06-15 中国人民解放军国防科学技术大学 一种全极化雷达目标极化相干特征的可视化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007163318A (ja) * 2005-12-14 2007-06-28 Toshiba Corp レーダ信号処理装置とその一定誤警報確率検出方法
CN103018722A (zh) * 2012-10-15 2013-04-03 西安电子科技大学 利用组网雷达系统对抗欺骗式假目标的方法
CN104765020A (zh) * 2015-04-12 2015-07-08 西安电子科技大学 有源假目标干扰的极化鉴别方法
CN105403875A (zh) * 2015-12-31 2016-03-16 中国人民解放军国防科学技术大学 双极化接收雷达的目标检测方法
CN105676203A (zh) * 2016-01-23 2016-06-15 中国人民解放军国防科学技术大学 一种全极化雷达目标极化相干特征的可视化方法
CN105652273A (zh) * 2016-03-17 2016-06-08 哈尔滨工程大学 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107290727A (zh) * 2017-06-14 2017-10-24 中国人民解放军国防科学技术大学 一种位置可控式isar假目标图像合成方法
CN107290727B (zh) * 2017-06-14 2019-10-11 中国人民解放军国防科学技术大学 一种位置可控式isar假目标图像合成方法
CN108919206A (zh) * 2018-07-18 2018-11-30 武汉大学 一种基于子载波处理的外辐射源雷达极化滤波方法
CN108919206B (zh) * 2018-07-18 2022-03-04 武汉大学 一种基于子载波处理的外辐射源雷达极化滤波方法
CN110794374A (zh) * 2019-11-08 2020-02-14 哈尔滨工业大学 一种间歇采样转发干扰的参数辨识方法
CN112799028A (zh) * 2020-12-14 2021-05-14 中电科仪器仪表有限公司 一种基于rcs起伏统计特性差异的虚假目标识别方法

Also Published As

Publication number Publication date
CN106168661B (zh) 2018-07-24

Similar Documents

Publication Publication Date Title
CN104765020B (zh) 有源假目标干扰的极化鉴别方法
CN106249212A (zh) 主瓣压制干扰背景下有源假目标的极化鉴别方法
CN106125053A (zh) 脉冲多普勒雷达极化抗干扰方法
CN102866391B (zh) 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法
CN104267379B (zh) 一种基于波形设计的主被动雷达协同抗干扰方法
CN103728598B (zh) 用异地配置的主被动雷达网抑制航迹欺骗干扰的方法
CN103852756B (zh) 一种连续波雷达目标检测跟踪方法
CN106168661A (zh) 随机捷变极化有源假目标的鉴别方法
CN103760530B (zh) 基于信号锥的雷达压制式干扰和欺骗式干扰识别方法
CN106338715B (zh) 一种全极化有源假目标的鉴别方法
CN105259541B (zh) 一种多站雷达抗有源欺骗式干扰的方法
CN105137396B (zh) 一种smsp干扰和c&i干扰的检测方法
CN105403875B (zh) 双极化接收雷达的目标检测方法
CN105807267A (zh) 一种mimo雷达扩展目标的检测方法
CN103954943A (zh) 网络化雷达系统抗欺骗式干扰方法
CN106249208A (zh) 基于分数阶傅里叶变换的噪声调幅干扰下信号检测方法
CN106019243A (zh) 一种基于三次初相和frft的drfm干扰检测和对抗方法
CN105044686A (zh) 一种雷达密集假目标干扰抑制方法
CN105527614A (zh) 基于动目标检测的lfm脉冲雷达抗drfm干扰方法
CN106443660A (zh) 一种线性调频连续波雷达多运动目标匹配方法
CN105022041A (zh) 基于经验模态分解特征库的拖引干扰子类型识别方法
CN106338714B (zh) Pri间恒定极化有源假目标鉴别方法
CN105093199A (zh) 基于雷达时域回波的目标识别特征提取方法
CN106019250A (zh) 基于角闪烁转发式假目标鉴别方法
CN106209715A (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