CN113325378B - 基于俯仰tda的距离模糊杂波抑制方法 - Google Patents

基于俯仰tda的距离模糊杂波抑制方法 Download PDF

Info

Publication number
CN113325378B
CN113325378B CN202110402564.8A CN202110402564A CN113325378B CN 113325378 B CN113325378 B CN 113325378B CN 202110402564 A CN202110402564 A CN 202110402564A CN 113325378 B CN113325378 B CN 113325378B
Authority
CN
China
Prior art keywords
distance
echo signal
clutter
tda
echo
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
Application number
CN202110402564.8A
Other languages
English (en)
Other versions
CN113325378A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202110402564.8A priority Critical patent/CN113325378B/zh
Publication of CN113325378A publication Critical patent/CN113325378A/zh
Application granted granted Critical
Publication of CN113325378B publication Critical patent/CN113325378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/41Details 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/414Discriminating targets with respect to background clutter
    • 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
    • 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/41Details 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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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

本发明公开了一种基于俯仰TDA的距离模糊杂波抑制方法,包括以下步骤:获取俯仰TDA的原始回波信号;对原始回波信号进行滤波处理,以对距离模糊进行抑制,得到期望回波信号;利用APC方法对期望回波信号进行解码,得到无模糊区域回波信号;对无模糊区域回波信号进行杂波补偿,得到目标回波信号。本发明提供的方法不需要正交波形即可有效解决距离模糊问题,且具有良好的杂波抑制性能,能够用于实现广域监视和目标跟踪。

Description

基于俯仰TDA的距离模糊杂波抑制方法
技术领域
本发明属于信号处理技术领域,具体涉及一种基于俯仰TDA的距离模糊杂波抑制方法。
背景技术
在高脉冲重复频率(HPRF)或广域监视雷达领域,距离模糊是一个不可避免的问题,因为多普勒模糊和距离模糊之间的根本矛盾无法同时解决。距离模糊问题在系统参数设计和信号处理算法中起着重要的作用,特别是在高速平台系统中。通常,来自无模糊距离和模糊距离区域的杂波回波具有不同的后向散射特性。如果距离模糊和距离依赖(由前视几何引起)并存,杂波特性将变得极其复杂,这将进一步降低STAP(Space-Time AdaptiveProcessing,空时自适应处理)方法的性能。
为了解决距离模糊的问题,现有技术提供了一种方位相位编码波形技术,其利用多普勒带通滤波器能够可靠地分离每个发射波形对应的距离模糊回波;然后利用波形分集的优势,将多输入多输出(MIMO)系统与频率分集阵列(FDA)、俯仰维滤波等技术相结合来解决距离模糊问题。在地面运动目标检测中,MIMO可以形成距离-角度相关的波束模式,从而提高距离模糊杂波的抑制性能。
然而,GMTI雷达(地面动目标显示雷达)的大多数方法都假定多个发射机发射不相关的正交波形,并通过匹配滤波处理将这些发射波形在每个接收机上完美分离。在实际中,很难找到完全满足正交性要求的波形族。波形之间的互相关通常会产生较高的旁瓣电平,从而对杂波抑制产生不利影响。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于俯仰TDA的距离模糊杂波抑制方法。本发明要解决的技术问题通过以下技术方案实现:
一种基于俯仰TDA的距离模糊杂波抑制方法,包括以下步骤:
获取俯仰TDA的原始回波信号;
对所述原始回波信号进行滤波处理,以对距离模糊进行抑制,得到期望回波信号;
利用APC方法对所述期望回波信号进行解码,得到无模糊区域回波信号;
对所述无模糊区域回波信号进行杂波补偿,得目标回波信号。
在本发明的一个实施例中,所述原始回波信号表示为:
Figure GDA0004287839280000021
其中,ynk(t,tk)表示第n个阵元第k个脉冲接收到的信号,ρlq表示杂波散射系数,wa(tk)表示方位窗,χ(tk)表示方位相位编码,tk表示方位角时间,τT,m表示发射延迟时间,τR,m表示接受延迟时间,Rl表示斜距,ψlq表示天线圆锥角,βlq表示速度圆锥角,
Figure GDA0004287839280000031
表示接收空间频率,/>
Figure GDA0004287839280000032
表示归一化多普勒频率。
在本发明的一个实施例中,对所述原始回波信号进行滤波处理,以对距离模糊进行抑制,得到期望回波信号,包括:
对所述原始回波信号进行角度频域的匹配滤波;
对匹配滤波后的回波信号进行距离频域带通滤波,得到期望回波信号。
在本发明的一个实施例中,对所述原始回波信号进行角度频域的匹配滤波,包括:
构建匹配滤波函数如下:
Figure GDA0004287839280000033
其中,fr表示距离频域,B为信号带宽,μ表示调频信号调频率,pe(frl)表示俯仰TDA的空频图;
按照所述匹配滤波函数对回波信号进行匹配滤波,得到滤波后的回波信号,表示为:
Figure GDA0004287839280000034
其中,Ynk(fr,tk)表示ynk(t,tk)的距离域FFT。
在本发明的一个实施例中,对匹配滤波后的回波信号进行距离频域带通滤波,得到期望回波信号,包括:
构建带通滤波函数;
根据所述带通滤波函数对所述匹配滤波后的回波信号的每个距离单元进行距离模糊抑制和距离压缩,得到每个距离单元的压缩结果;
将所有距离单元的压缩结果合成,得到期望回波信号。
在本发明的一个实施例中,所述带通滤波函数的表达式为:
Figure GDA0004287839280000041
其中,
Figure GDA0004287839280000042
和/>
Figure GDA0004287839280000043
为pe(frl)的左右第一个零点。
在本发明的一个实施例中,所述每个距离单元的压缩结果表示为:
Figure GDA0004287839280000044
其中,Ar(·)为IFFT处理后的回波幅度,其表达式为:
Figure GDA0004287839280000045
在本发明的一个实施例中,利用APC方法对所述期望回波信号进行解码,以对距离模糊进行进一步抑制,得到无模糊区域回波信号,包括:
构建译码函数;
利用所述译码函数对所述期望回波信号进行解码,以实现模糊信号的多普勒频移;
对解码后的期望回波信号进行多普勒带通滤波,以滤除模糊信号的多普勒域,得到无模糊区域回波信号。
在本发明的一个实施例中,所述译码函数表示为:
Figure GDA0004287839280000051
其中,χ(tk)为方位相位编码,tk=(k-1)fPRF为方位角时间(慢时间),fPRF为脉冲重复频率;k(k=1,…,K)为发射脉冲数,K为累计相干脉冲总数,Mγ≥2为正整数,表示APC移位因子。
在本发明的一个实施例中,对所述无模糊区域回波信号进行杂波补偿,得目标回波信号,包括:
对所述无模糊区域回波信号进行多普勒频移补偿;
采用局域联合自适应处理方法对经过多普勒频移补偿的无模糊区域回波信号进行杂波抑制,得到目标回波信号。
本发明的有益效果:
1、本发明提供的基于俯仰TDA的距离模糊杂波抑制方法利用俯仰TDA雷达的杂波信号特性,在距离频域对距离模糊回波进行分离,并利用APC技术抑制残留的旁瓣模糊能量,最后再利用杂波补偿方法缓解距离依赖性,有效解决了距离模糊问题,具有良好的杂波抑制性能;
2、本发明提供的基于俯仰TDA的距离模糊杂波抑制方法不需要正交波形即可很容易的实现距离模糊杂波抑制,且杂波抑制性能较好,能够用于实现广域监视和目标跟踪。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的基于俯仰TDA的距离模糊杂波抑制方法流程示意图;
图2a-2b是传统相控阵雷达的空间频域和多普勒频域联合杂波谱图;
图3a-3b是本发明仿真实验中俯仰TDA雷达仅进行距离频率滤波的第一区域的杂波谱;
图4a-4b是本发明仿真中基于二维滤波的俯仰TDA雷达第一区域的杂波谱;
图5a-5d是本发明仿真中俯仰TDA雷达杂波补偿后的其他区域的杂波光谱;
图6a-6c是本发明仿真中SCNR损失随归一化多普勒频率变化图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
实施例一
本实施例首先对俯仰TDA阵列(Time diversity Arrays时间分集阵)的回波信号进行匹配滤波,再通过距离频域带通滤波将期望回波从距离模糊中分离出来,实现距离频域的距离压缩和模糊抑制;然后通过相位编码技术将模糊区域的多普勒频谱进行平移,通过多普勒带通滤波器进行残差距离模糊抑制;最后进行多普勒频移补偿,使主距离区域的杂波数据近似独立且同分布,利用JDL自适应处理方法对杂波进行抑制,实现距离模糊杂波抑制。
具体地,请参见图1,图1是本发明实施例提供的基于俯仰TDA的距离模糊杂波抑制方法流程示意图,包括以下步骤:
S1:获取俯仰TDA的原始回波信号。
在本实施例中,首先要对俯仰TDA的发射信号进行方位相位编码。
具体地,编码公式如下:
Figure GDA0004287839280000071
其中,tk=(k-1)fPRF是方位向时间,fPRF是脉冲重复频率;k(k=1,…,K)是发射脉冲数,K表示发射最大数量;Mγ≥2是一个正整数,表示方位相位编码技术中的位移因子。
然后,接收该经过编码的发射信号的回波信号。
具体地,本实施例构建TDA阵列第n个阵元第k个脉冲的接收回波信号如下:
Figure GDA0004287839280000072
其中,ynk(t,tk)表示第n个阵元第k个脉冲接收到的信号,ρlq表示杂波散射系数,wa(tk)表示方位窗,χ(tk)表示方位相位编码,τT,m表示发射延迟时间,τR,m表示接受延迟时间,Rl表示斜距,ψlq表示天线圆锥角,βlq表示速度圆锥角,
Figure GDA0004287839280000073
表示接收空间频率,/>
Figure GDA0004287839280000074
表示归一化多普勒频率。
S2:对原始回波信号进行滤波处理,以对距离模糊进行抑制,得到期望回波信号。具体包括:
S21:对原始回波信号进行角度频域的匹配滤波。
理论上,传统的距离频域匹配函数等于
Figure GDA0004287839280000081
然而,距离频谱和空间导向矢量之间的耦合使得传统的匹配滤波器不适用。因此,本实施例构建了一个新的滤波函数实现角度频域匹配滤波,以将脉冲压缩和等效发射波束形成结合到一个过程中。
具体地,构建匹配滤波函数如下:
Figure GDA0004287839280000082
其中,fr为距离频域,B为信号带宽,μ调频信号调频率,(·)*为共轭算子,pe(frl)表示俯仰TDA的空频图,其表达式为:
Figure GDA0004287839280000083
其中,M为TDA阵列的行数,fc为载频,Δt是两个相邻行的时间增量,间距都是d,φl表示仰角。
然后,按照上述匹配滤波函数对回波信号进行匹配滤波,得到滤波后的回波信号,表示为:
Figure GDA0004287839280000091
其中,Ynk(fr,tk)表示ynk(t,tk)的距离域FFT。
S22:对匹配滤波后的回波信号进行距离频域带通滤波,得到期望回波信号。
在本实施例中,角度频域匹配滤波后,对距离模糊进行抑制。为了通过距离频域带通滤波将期望回波从距离模糊中分离出来,需要从仰角φl得到期望回波的主瓣距离谱信息。根据上述分析,利用pe(frl)的负、正第一个零点的值可以获得主瓣谱的范围,因此,为了消除距离模糊,可以采用fr,peakl)和
Figure GDA0004287839280000092
分别作为通带滤波器的中心频率和截止频率。
具体地,步骤S22包括:
S22-1:构建带通滤波函数如下:
Figure GDA0004287839280000093
其中,
Figure GDA0004287839280000094
和/>
Figure GDA0004287839280000095
为pe(frl)的左右第一个零点。
S22-2:根据带通滤波函数对所述匹配滤波后的回波信号的每个距离单元进行距离模糊抑制和距离压缩,得到每个距离单元的压缩结果。
由于步骤S22-1构建的带通滤波其函数会随着感兴趣的距离单元的仰角而变化,这意味着带通滤波器是与角度相关的。换句话说,一个带通滤波器和一个快速傅里叶反变换(IFFT)操作不能得到不同距离单元的距离压缩结果。在这种情况下,必须在相应的距离频率通带上通过距离IFFT处理来独立完成每个距离单元的距离模糊抑制和距离压缩。则每个距离单元的压缩结果可以表示为:
Figure GDA0004287839280000101
其中,Ar(·)为IFFT处理后的回波幅度,其表达式为:
Figure GDA0004287839280000102
S22-3:将所有距离单元的压缩结果合成,得到期望回波信号。
通过上述处理,即可将从明确的距离区域提取出预期的回波,也即在距离频域实现了对主距离模糊回波的分离。
S3:利用APC(Azimuth phase coding,方位相位编码)方法对期望回波信号进行解码,得到无模糊区域回波信号。
虽然步骤S2完成了主距离模糊回波的分离,然而,来自其他距离区域的残余旁瓣能量仍然混杂在期望回波信号中,这将降低后续处理的杂波抑制性能。
为了减轻混合在当前期望回波中的残留模糊能量的影响,需要在多普勒域进一步抑制这些信号。对于残差模糊杂波抑制,以p阶距离模糊为例,其距离频率滤波后的信号表达式可以表示为:
Figure GDA0004287839280000103
其中,Rpl=Rl+p·Ru为模糊倾斜范围,τpl为时延;ρplq、ψplq和βplq分别为其对应的散射系数、天线锥角和速度锥角。
接下来,执行脉冲相位解码程序。
首先,构建译码函数。
在本实施例中,对于第k个脉冲的接收信号的译码函数可以构造为:
Figure GDA0004287839280000111
然后,利用译码函数对期望回波信号进行解码,以实现模糊信号的多普勒频移。
由于距离模糊,第p次模糊信号相对于无模糊信号有一个额外的往返延迟pfPRF。通过编解码,p次距离模糊的多普勒频移的形式为:
Figure GDA0004287839280000112
其中,
Figure GDA0004287839280000113
是模算子,它说明了离散傅里叶变换的周期性和Δfp在[-fPRF2,fPRF2]区间内的局限性。
具体地,将译码函数Hde应用到回波信号,可以得到回波信号为:
Figure GDA0004287839280000114
Figure GDA0004287839280000115
其中,
Figure GDA0004287839280000121
和/>
Figure GDA0004287839280000122
分别为回波信号中无距离模糊回波与距离模糊回波信号,将他们分别做慢时间FFT得到距离多普勒域信号
Figure GDA0004287839280000123
Figure GDA0004287839280000124
本实施例将译码函数应用于得到的期望回波信号,可以将p阶范围模糊的多普勒谱通过Δfp进行多普勒频移,使得无模糊信号和距离模糊信号在多普勒域,不再精确地叠加,且具有不同的支持区域。
最后,对解码后的期望回波信号进行多普勒带通滤波,以滤除模糊信号的多普勒域,得到无模糊区域回波信号。
具体地,本实施例设置滤波器的处理带宽为BD,考虑到前视几何、杂波的多普勒频谱只占总多普勒带宽的一半。设ζ表示方位角过采样因子,即ζ=PRF/BD。而通常,PRF(pulserepetition frequency,即脉冲重复频率)被设置为高于多普勒带宽,即ζ≥1。当ζ=1.5时,处理带宽为BD的滤波器可以部分去除处理带宽以外的剩余模糊能量,因此本实施例可通过APC技术获得较好的距离模糊抑制性能。
S4:对无模糊区域回波信号进行杂波补偿,得目标回波信号。
S41:对无模糊区域回波信号进行多普勒频移补偿。
在步骤2和步骤3中,利用距离频域和多普勒域的滤波运算,消除了来自其他距离区域的距离模糊杂波的影响,从而提取出包含杂波和目标信号的主距离区域混合回波。对于得到的无模糊区域回波,要抑制原距离区域的残留杂波,以便对目标进行检测。
对所有二维带通滤波器的输出进行叠加和重新排列,杂波片的时空快照向量可以表示为
Figure GDA0004287839280000131
其中,
Figure GDA0004287839280000132
为Kronecker积,/>
Figure GDA0004287839280000133
和/>
Figure GDA0004287839280000134
分别为接收矢量和时间矢量,表示为:
Figure GDA0004287839280000135
Figure GDA0004287839280000136
其中,上标(·)T表示转置操作。
地杂波可以看作是同一距离单元所有散射体回波的总和。因此,经过上述操作后,同一距离单元内的数据向量可表示为
Figure GDA0004287839280000137
其中,Nc表示在同一距离单元内独立的杂波散射体的数量。
同理,对于复振幅为ε0,径向速度为v0的运动目标,其返回的数据快照sl可以表示为
Figure GDA0004287839280000138
其中,fDt=(2vcos(β0)+2v0)/(λ·fPRF)为归一化多普勒频率,ψ0和β0分别为运动目标的天线和速度锥角。
因此,提取的由运动目标、杂波和高斯噪声组成的回波数据矢量可以重写为:
xl=cl+sl+nl
其中,nl为白高斯分布,信道和脉冲之间独立的系统噪声,即
Figure GDA0004287839280000141
INK和/>
Figure GDA0004287839280000142
分别表示NK×NK单位矩阵和噪声功率。
需要说明的是,由于前视几何,所提取的杂波仍然是与距离相关的。距离依赖性在近距离区域显著,在远距离区域逐渐减弱。为了解决这一问题,本实施例采用多普勒频移补偿法,即:
Figure GDA0004287839280000143
其中,ΔfD,l为第l距离单元的多普勒频移量。
S42:采用局域联合自适应处理方法对经过多普勒频移补偿的无模糊区域回波信号进行杂波抑制,得到目标回波信号。
由于经过多普勒频移补偿后,主距离区域的杂波数据近似独立且同分布,然后便可直接将常规的STAP算法应用到补偿后的数据中进行目标检测。
本实施例采用局域联合(JDL)自适应处理方法对杂波进行抑制。
考虑到主距离区域和圆锥角为ψ0的处理器,本实施采用在接收空间域形成Q1波束以进行降维自适应处理。
具体地,空间降维变换矩阵可表示为:
Figure GDA0004287839280000144
接下来,考虑Q2连续多普勒单元对应假定的速度v0,时间降维变换矩阵可以表示为:
Figure GDA0004287839280000145
由此,可以得到空时降维变换矩阵为:
Figure GDA0004287839280000151
进一步,可以得到基于JDL的自适应处理的快拍
然后,基于JDL的自适应处理可以写成以下优化问题的解决方案:
Figure GDA0004287839280000153
其中,
Figure GDA0004287839280000154
为降维杂波加噪声协方差矩阵,l0为目标所在的距离单元,
Figure GDA0004287839280000155
是降维之后导向向量。
本发明利用俯仰TDA雷达的杂波信号特性,在距离频域对距离模糊回波进行分离;利用APC技术抑制残留的旁瓣模糊能量;再对不同距离区域的回波进行提取后,可以采用传统的杂波补偿方法来缓解距离依赖性。该框架可用于实现广域监视和目标跟踪。与传统的STAP方法相比,该框架能够有效地解决距离模糊问题,具有良好的杂波抑制性能。此外,与基于FDA-MIMO的方法相比,该方法不需要正交波形,因此更容易实现。
实施例二
下面通过仿真实验对本发明的有益效果进行进一步说明。
1、实验条件
本实施例的仿真实验的硬件平台为:Intel(R)Core(TM)i5-8265U CPU@1.60GHz,频率为1.8GHz,Nvidia GeForce MX250。
本实施例的仿真实验的软件使用matlab2016b。
本实施例考虑了一种前视机载俯仰TDA雷达系统,主要仿真参数如表1所示,杂噪比为30dB,训练样本为300个。假设目标锥角为90度,径向速度为30m/s。
表1仿真参数一览表
Figure GDA0004287839280000161
2、仿真内容及结果分析
请参见图2a-2b,图2a-2b是传统相控阵雷达的空间频域和多普勒频域联合杂波谱图,其中,图2a为杂波补偿前,图2b为杂波补偿后。由图2a可见,不同距离区域的杂波光谱是无法区分的。此外,不同距离区域的杂波谱严重依赖于距离,这意味着IID条件不满足。需要注意的是,传统的杂波补偿程序只能在不存在距离模糊的情况下执行。在仿真中,高PRF导致了距离模糊问题,传统的杂波补偿技术不再有效。距离依赖性在近距离区域非常严重,在远距离区域则比较缓和,如图2a所示。由于各距离区域的杂波补偿函数差异较大,对第一个距离区域进行补偿后,距离依赖性减弱,而对其他距离模糊区域的杂波补偿减弱,如图2b所示。
请参见图3a-3b,图3a-3b是本发明仿真实验中俯仰TDA雷达仅进行距离频率滤波的第一区域的杂波谱,其中,图3a为杂波补偿前,图3b为杂波补偿后。使用本发明提出的距离频率滤波技术,大部分距离模糊杂波能量被消除,而一小部分旁瓣能量仍然存在,如图3a所示。通过距离相关补偿,对第一个距离区域的不同距离箱杂波能量进行聚焦。但从图3b可以看出,空时平面上存在散焦散斑,这是由于残留的旁瓣距离不明确信号造成的。
为了缓解这种现象,在俯仰TDA雷达中引入了APC技术。由于前视几何结构的存在,杂波谱在联合空间频域和多普勒频域只占据右半空间。利用APC技术,部分残留的模糊杂波可以转移到空时平面的空白区域。因此,可以采用多普勒带通滤波器来抑制剩余的模糊能量。请参见图4a-4b,图4a-4b是本发明仿真中基于二维滤波的俯仰TDA雷达其一区域的杂波谱;其中,图4a为杂波补偿前,图4b为杂波补偿后。从图4a可以看出,剩余的旁瓣杂波能量大部分已经被消除。经过杂波补偿后,聚焦杂波频谱如图4b所示。显然,与距离频率滤波方法相比,二维滤波方法可以获得更好的旁瓣杂波抑制性能。
同样,俯仰TDA雷达也可以提取第二和第三距离区域的距离模糊杂波信号。第二和第三模糊区域杂波的距离依赖性较弱,杂波补偿的效果可以忽略。请参见图5a-5d,图5a-5d是本发明仿真中俯仰TDA雷达杂波补偿后的其他区域的杂波光谱,其中,图5a为进行距离频率滤波的第二区域的杂波谱,图5b为二维滤波第二区域的杂波谱,图5c为进行距离频率滤波的第三区域的杂波谱,图5d为二维滤波第三区域的杂波谱。
由于PRF有限,第三量程区域的剩余模糊能量只占正半轴的一小部分。在相应的多普勒滤波过程中,我们缩小了滤波器的支持区域。即部分位于正半轴的杂波被消除,包括第三区域的原始杂波和其他区域的模糊杂波。
请参见图6a-6c,图6a-6c是本发明仿真中SCNR损失随归一化多普勒频率变化图。图6a-6c给出了三个区域不同方法的SCNR损耗曲线,包括相控阵传统的STAP方法、基于俯仰TDA的两种方法(带APC和不带APC)以及相应的杂波补偿方法,此外,还给出了理想的无干扰回波曲线作为基准。其中,图6a-6c中的SCNR损耗曲线对应于三个范围模糊区域。由于距离模糊,传统的相控阵分段方法性能严重下降。在这种情况下,不同距离区域的杂波补偿函数不同,会使其他距离区域的曲线切口变宽。相反相控阵,俯仰TDA分离距离模糊的杂波频域范围和传统杂波补偿程序可以应用于每一个独立的距离单元。因此,该方法的俯仰TDA雷达性能得到了很好的维护。然而,距离模糊区域的旁瓣能量混合在原距离区域会导致杂波抑制性能下降。利用APC技术可以将剩余的模糊能量转移到其他方向,并利用多普勒带通滤波器消除部分旁瓣能量。因此,可以在一定程度上提高杂波抑制性能。在第二和第三个距离区域,消除了杂波的距离依赖性,基本满足IID条件,对是否进行杂波补偿几乎没有影响。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (9)

1.一种基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,包括以下步骤:
获取俯仰时间分集阵的原始回波信号;
对所述原始回波信号进行滤波处理,以对距离模糊进行抑制,得到期望回波信号;包括:
对所述原始回波信号进行角度频域的匹配滤波;
对匹配滤波后的回波信号进行距离频域带通滤波,得到期望回波信号;
利用方位相位编码方法对所述期望回波信号进行解码,得到无模糊区域回波信号;
对所述无模糊区域回波信号进行杂波补偿,得到目标回波信号。
2.根据权利要求1所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,所述原始回波信号表示为:
Figure FDA0004106118860000011
其中,ynk(t,tk)表示第n个阵元第k个脉冲接收到的信号,ρlq表示杂波散射系数,wa(tk)表示方位窗,χ(tk)表示方位相位编码,tk表示方位角时间,τT,m表示发射延迟时间,τR,m表示接受延迟时间,Rl表示斜距,ψlq表示天线圆锥角,βlq表示速度圆锥角,
Figure FDA0004106118860000012
表示接收空间频率,/>
Figure FDA0004106118860000021
表示归一化多普勒频率;M表示TDA阵列的行数,fc表示载频,d表示两个相邻行的时间增量的间距,φl表示仰角,Δt表示两个相邻行的时间增量,fPRF表示脉冲重复频率。
3.根据权利要求2所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,对所述原始回波信号进行角度频域的匹配滤波,包括:
构建匹配滤波函数如下:
Figure FDA0004106118860000022
其中,fr表示距离频域,B为信号带宽,μ表示调频信号调频率,pe(frl)表示俯仰TDA的空频图;
按照所述匹配滤波函数对回波信号进行匹配滤波,得到滤波后的回波信号,表示为:
Figure FDA0004106118860000023
其中,Ynk(fr,tk)表示ynk(t,tk)的距离域FFT。
4.根据权利要求3所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,对匹配滤波后的回波信号进行距离频域带通滤波,得到期望回波信号,包括:
构建带通滤波函数;
根据所述带通滤波函数对所述匹配滤波后的回波信号的每个距离单元进行距离模糊抑制和距离压缩,得到每个距离单元的压缩结果;
将所有距离单元的压缩结果合成,得到期望回波信号。
5.根据权利要求4所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,所述带通滤波函数的表达式为:
Figure FDA0004106118860000031
其中,
Figure FDA0004106118860000032
和/>
Figure FDA0004106118860000033
为pe(frl)的左右第一个零点。
6.根据权利要求5所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,所述每个距离单元的压缩结果表示为:
Figure FDA0004106118860000034
其中,Ar(·)为IFFT处理后的回波幅度,其表达式为:
Figure FDA0004106118860000035
7.根据权利要求1所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,利用APC方法对所述期望回波信号进行解码,以对距离模糊进行进一步抑制,得到无模糊区域回波信号,包括:
构建译码函数;
利用所述译码函数对所述期望回波信号进行解码,以实现模糊信号的多普勒频移;
对解码后的期望回波信号进行多普勒带通滤波,以滤除模糊信号的多普勒域,得到无模糊区域回波信号。
8.根据权利要求7所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,所述译码函数表示为:
Figure FDA0004106118860000041
其中,χ(tk)为方位相位编码,tk=(k-1)/fPRF为方位角时间,fPRF为脉冲重复频率;k,k=1,…,K为发射脉冲数,K为累计相干脉冲总数,Mγ≥2为正整数,表示APC移位因子。
9.根据权利要求1所述的基于俯仰TDA的距离模糊杂波抑制方法,其特征在于,对所述无模糊区域回波信号进行杂波补偿,得目标回波信号,包括:
对所述无模糊区域回波信号进行多普勒频移补偿;
采用局域联合自适应处理方法对经过多普勒频移补偿的无模糊区域回波信号进行杂波抑制,得到目标回波信号。
CN202110402564.8A 2021-04-14 2021-04-14 基于俯仰tda的距离模糊杂波抑制方法 Active CN113325378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110402564.8A CN113325378B (zh) 2021-04-14 2021-04-14 基于俯仰tda的距离模糊杂波抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110402564.8A CN113325378B (zh) 2021-04-14 2021-04-14 基于俯仰tda的距离模糊杂波抑制方法

Publications (2)

Publication Number Publication Date
CN113325378A CN113325378A (zh) 2021-08-31
CN113325378B true CN113325378B (zh) 2023-07-14

Family

ID=77414987

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110402564.8A Active CN113325378B (zh) 2021-04-14 2021-04-14 基于俯仰tda的距离模糊杂波抑制方法

Country Status (1)

Country Link
CN (1) CN113325378B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977571A (zh) * 2015-06-25 2015-10-14 西安电子科技大学 基于俯仰频率分集stap的距离模糊杂波抑制方法
EP3040737A2 (en) * 2014-12-30 2016-07-06 Propagation Research Associates, Inc. Using orthogonal space projections to generate a constant false alarm rate control parameter
CN109738889A (zh) * 2018-12-19 2019-05-10 柳毅 一种认知型线调频脉冲航海雷达系统的探测方法
CN111693983A (zh) * 2020-06-12 2020-09-22 哈尔滨工业大学(威海) 互补波形构建方法及模块、认知雷达系统及波形发射方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107561535A (zh) * 2017-09-07 2018-01-09 中国科学院电子学研究所 一种合成孔径雷达距离模糊抑制方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3040737A2 (en) * 2014-12-30 2016-07-06 Propagation Research Associates, Inc. Using orthogonal space projections to generate a constant false alarm rate control parameter
CN104977571A (zh) * 2015-06-25 2015-10-14 西安电子科技大学 基于俯仰频率分集stap的距离模糊杂波抑制方法
CN109738889A (zh) * 2018-12-19 2019-05-10 柳毅 一种认知型线调频脉冲航海雷达系统的探测方法
CN111693983A (zh) * 2020-06-12 2020-09-22 哈尔滨工业大学(威海) 互补波形构建方法及模块、认知雷达系统及波形发射方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Dan P. Scholnik.Range-ambiguous clutter suppression with Pulse-diverse waveforms.《2011 IEEE RadarCon (RADAR)》.2011,336-341. *
Jingwei Xu et al..Space-Time Adaptive Processing With Vertical Frequency Diverse Array for Range-Ambiguous Clutter Suppression.《IEEE Transactions on Geoscience and Remote Sensing》.2016,第54卷(第9期),5352-5364. *
吴旭姿等.基于俯仰频率分集技术的波束形成方法.《电子与信息学报》.2016,第38卷(第12期),3070-3077. *
祝晓静等.基于改进方位相位编码的全极化SAR距离模糊抑制方法.《雷达学报》.2017,6(4),420-431. *
郭一帆等.基于张量稀疏恢复的MIMO雷达杂波抑制方法.《系统工程与电子技术》.2017,第39卷(第12期),2677-2682. *

Also Published As

Publication number Publication date
CN113325378A (zh) 2021-08-31

Similar Documents

Publication Publication Date Title
CN110412568B (zh) 基于扩展方位角相位编码的距离模糊杂波抑制方法
CN111880171B (zh) 一种消除雷达目标盲速的脉冲分段编码方法
CN102520395B (zh) 基于双基地多输入多输出雷达的杂波抑制方法
Wang et al. MIMO SAR using chirp diverse waveform for wide-swath remote sensing
DE102020111533A1 (de) Mimo-radar-vorrichtungen und mimo-radar-verfahren
CN110703209B (zh) 高重频机载前视阵雷达地面距离模糊杂波的抑制方法
CN110471033B (zh) 基于相干fda的二维空时编码干扰抑制方法
CN112596033A (zh) 阵元-脉冲编码mimo雷达欺骗式干扰抑制方法
CN110412570B (zh) 基于空间脉冲相位编码的hrws-sar成像方法
CN109814070B (zh) 基于辅助脉冲的距离模糊杂波抑制方法
US11693088B2 (en) Method and device for separating echo signals of space-time waveform encoding synthetic aperture radar in elevation
He et al. Range ambiguous clutter suppression approach with elevation time diverse array
CN106680797A (zh) 基于宽带模糊函数的目标参数估计新方法
Chen et al. Range-dependent ambiguous clutter suppression for airborne SSF-STAP radar
CN110940953A (zh) 地波雷达海杂波中目标的三维检测方法
CN113238212B (zh) 基于空时编码的频率分集阵列雷达距离分辨率增强方法
He et al. Near-range clutter suppression with elevation element multifrequency subpulse coding array radar
Wu et al. Interference suppression method with MR-FDA-MIMO radar
Qiu et al. Range-Ambiguous Clutter Suppression for Space-Based Early Warning Radar Using Vertical FDA and Horizontal EPC
CN113325378B (zh) 基于俯仰tda的距离模糊杂波抑制方法
He et al. An adaptive coding-angle-Doppler clutter suppression approach with extended azimuth phase coding array
CN109597034B (zh) 一种基于欧几里得距离的空时自适应处理方法
CN113109785A (zh) 一种mimo雷达多通道分离方法
CN112198510B (zh) 双接收通道混合全极化雷达系统抑制方位模糊的方法
Landi et al. Time-orthogonal-waveform-space-time adaptive processing for distributed aperture radars

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