CN103744067A - 一种非自适应的机载非正侧视雷达近程杂波抑制方法 - Google Patents
一种非自适应的机载非正侧视雷达近程杂波抑制方法 Download PDFInfo
- Publication number
- CN103744067A CN103744067A CN201410018068.2A CN201410018068A CN103744067A CN 103744067 A CN103744067 A CN 103744067A CN 201410018068 A CN201410018068 A CN 201410018068A CN 103744067 A CN103744067 A CN 103744067A
- Authority
- CN
- China
- Prior art keywords
- clutter
- space
- matrix
- time
- radar
- 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.)
- Pending
Links
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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—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
技术领域
本发明属于雷达信号处理技术领域,具体说是对机载非正侧视雷达,用一种非自适应的由空-时滤波和空-时匹配两部分级联构成的近程杂波抑制方法(Short-range ClutterSuppression Approach,SCSA),来抑制非正侧视雷达的近程杂波,实现对目标的检测。
背景技术
第二次世界大战之后,现代雷达技术已经历了80多年的蓬勃发展,一些重大技术的突破使雷达的性能更加完善。传统的地基雷达由于是静止放置的,杂波回波不存在空时耦合特性,因而地杂波主要分布在中心零频附近,广泛使用的动目标显示(MTI)技术和动目标检测(MTD)技术就可以实现对目标的检测。但是,对于机载预警雷达,由于其杂波呈现强烈的空时耦合特性,这些传统的成熟技术已不再适用。为了能够更好地抑制有空时耦合特性的杂波,人们提出了空时自适应处理(STAP)技术。STAP技术能够补偿载机的运动效应,通过空时二维滤波很好地抑制地杂波,提高机载预警雷达对动目标的检测性能,并能对通道误差、近程散射以及干扰进行自适应的抑制。Brennan和Reed率先将阵列自适应处理技术推广到机载雷达,提出了STAP的概念,并设计了最优STAP处理器。最优STAP处理器在通常情况下可以得到最优的杂波抑制效果,并获得最好的目标检测性能。但是,在实际应用中,最优STAP处理器的计算复杂度很高,并且估计杂波的协方差矩阵需要大量独立同分布(IID)的训练样本,因而不能满足实时处理的要求。为了降低STAP算法的运算量和对样本的要求,人们提出了降维STAP算法,其中最有代表性的是因子化(FA)算法、扩展因子化(EFA)算法和局域化(JDL)算法。FA算法是在一个多普勒通道内进行的空域自适应算法,它是先将训练样本投影到多普勒域,然后选取目标所在的多普勒通道进行空域自适应处理。EFA算法是用相邻的多普勒通道进行联合空域自适应处理,它除了要选取目标所在的多普勒通道之外,还要额外地选取相邻的多个多普勒通道进行联合自适应处理。JDL方法是在波束-多普勒域进行自适应处理,它是先将训练样本投影到多普勒-波束域,然后进行多普勒-波束域的联合自适应处理。但是,这些典型的自适应算法都不可避免地依赖于训练样本来估计杂波和噪声协方差矩阵。Reed,Brennan和Mallett提出的RMB准则表明,随着独立同分布样本数的增加,自适应算法会逐渐地收敛于最优的性能。因此,大量满足IID条件的样本对于自适应算法非常重要。然而,实际中绝大多数情况下,杂波样本都不满足IID条件,即杂波是非均匀的。只有在最简单的正侧视线阵的情况下,杂波样本才基本满足IID条件。因此,除正侧视雷达之外,对于一般的非正侧视雷达,其杂波都是非均匀的。如果将非均匀杂波的回波数据用作训练样本,由于用来估计杂波协方差矩阵的杂波样本与当前检测距离单元的杂波样本在角度-多普勒域的分布是不同的,这就会导致对杂波的协方差矩阵估计的误差,并最终会导致如FA算法、EFA算法和JDL算法等经典自适应算法性能的严重下降。特别地,由于杂波相对于载机的俯仰角在近程是快变的,近程杂波的非均匀性因而会更加严重,经典的自适应算法有时甚至不能正常使用。
发明内容
本发明的目的是:针对现有的STAP算法应用于机载雷达杂波抑制时存在的诸多不足,如计算复杂度高,估计杂波的协方差矩阵需要大量独立同分布的训练样本以及无法有效抑制近程杂波等缺陷,本发明提出了一种非自适应的机载雷达近程杂波抑制方法(SCSA)。本发明方法是由空-时滤波和空-时匹配两部分级联构成的非自适应方法,不会受到非均匀样本的影响,因而比自适应方法有更好的检测性能,对于近程杂波也有显著的抑制效果。此外,本发明方法的计算复杂度低,可以方便快速地实现对杂波的抑制。
本发明的技术方案概括为:首先将雷达各通道接收到的信号排列成一个列向量,然后对接收信号进行空-时滤波,实现对杂波能量的抑制,最后进行空-时匹配,实现对目标的检测。具体实现过程如下:
(1)将N个雷达通道在第k个雷达脉冲时刻分别接收到的杂波信号排列成一个长为N的列向量其中,k表示雷达脉冲序号,k=0,1,…,K-1,K表示在一个相干处理间隔内进行相干累积的脉冲个数,为修正空域矩阵,F为增益矩阵,Wk为多普勒矩阵,向量由U个杂波散射单元的回波幅度构成,U为杂波散射单元个数,上标T表示对矩阵或向量求转置;
(2)空-时滤波,以最小化剩余杂波能量:将第k个和第k+1个脉冲的杂波相减,得到剩余杂波能量为: 其中,Q∈CN×N是空-时滤波的系数矩阵,符号||||F表示矩阵或向量的F-范数,符号CA×B表示所有A×B维的复矩阵的集合;通过最小化剩余杂波能量得到然后由Q构造最终的空-时滤波的系数矩阵其中代表Moore-Penrose右伪逆操作,是由U个杂波散射单元的多普勒项构成的对角阵;
(3)空-时匹配,运用步骤(2)得到的空-时滤波系数矩阵求得空-时匹配的权矢量其中,st∈CK×1和ss∈CN×1分别为时域和空域的导向矢量,bt∈CK×1和bs∈CN×1分别是时域和空域的静态权矢量,符号⊙和分别表示Hadamard积和Kronecker积;用求得的权矢量进行空-时二维匹配,实现对目标的检测。
本发明与现有技术相比具有以下特点:
1、传统的自适应算法,如最优(Clairvoyant)算法、FA算法、EFA算法JDL算法,计算量巨大,工程实现困难,难以满足实时处理要求,表现在:***处理自由度(包括空域和时域两维的自由度)巨大,特别是在大规模阵列中,自由度更是十分庞大,这将会直接导致在杂波协方差矩阵的求逆运算中,运算量十分巨大,所需的设备量也十分惊人,满足不了实时处理的要求。而本发明方法中,如果方位角俯仰角θ以及载机飞行速度v均能被精确地测量,而且雷达参数能被提前获得,那么空-时滤波的系数矩阵Q就可以提前离线计算出来。如果载机速度和雷达参数发生变化,系数矩阵Q也会发生变化。为此,可以制成一个针对不同的情况的系数矩阵表格供使用时查询。因此它的计算量比传统的自适应算法要小很多,更有利于实时处理。
2、传统的自适应算法如最优(Clairvoyant)算法、FA算法、EFA算法JDL算法都不可避免地依赖于训练样本来估计杂波和噪声的协方差矩阵,要获得较好的估计性能就需要大量满足IID条件的杂波样本,即要求杂波均匀。如果将非均匀杂波的回波数据用作训练样本,由于用来估计杂波协方差矩阵的杂波样本与当前检测距离单元的杂波样本在角度-多普勒域的分布不同,因而会导致对杂波的协方差矩阵估计的误差,并最终会导致这些经典自适应算法的性能严重下降。特别地,由于杂波相对于载机的俯仰角在近程是快变的,近程杂波的非均匀性会更加严重,经典的自适应算法有时甚至不能正常使用。而本发明SCSA方法中,空-时滤波矩阵可以被提前计算出来,是一种非自适应方法,空-时匹配也是一种非自适应方法,因此SCSA方法是一种非自适应方法,不会受到非均匀样本的影响,因而比自适应方法有更好的检测性能。尤其是在抑制近程杂波时,由于近程杂波强烈的非均匀性导致经典自适应方法性能下降严重,而本发明的非自适应SCSA方法却可以获得一个较好的处理性能。如图9的(a)和图9的(b)所示,由于近程杂波的非均匀性,使得自适应方法的性能有所下降,而SCSA方法的性能相较于自适应方法有一定提升。在主瓣杂波区和旁瓣区,SCSA的信干噪比损失(SINR loss)曲线比FA、EFA和JDL算法均有一定的改善。
附图说明
图1是非正侧视雷达的杂波模型图;
图2是本发明的近程杂波抑制方法(SCSA)流程图;
图3是三种雷达杂波在角度-多普勒平面的分布轨迹图:其中,图3的(a)是正侧视雷达(α=0°)的杂波在角度-多普勒平面的分布轨迹图;图3的(b)是非正侧视雷达(α=30°)的杂波在角度-多普勒平面的分布轨迹图;图3的(c)是前视雷达(α=90°)的杂波在角度-多普勒平面的分布轨迹图;
图4是仿真试验一的信干噪比损失变化曲线图:图4的(a)是α=30°时,近程杂波的最优信干噪比损失(Clairvoyant SINR loss)随距离单元数的变化曲线;图4的(b)是α=30°时,样本估计下的信干噪比损失(Asymptotic SINR loss)随距离单元数的变化曲线;
图5是仿真试验一的功率、特征谱图:图5的(a)是α=30°时,原杂波最小方差无畸变响应(MVDR)功率谱图;图5的(b)是α=30°时,空-时滤波后MVDR功率谱图;图5的(c)是α=30°时,原杂波与空-时滤波后杂波的特征谱对比图;
图6是仿真试验一中本发明SCSA方法与较传统的自适应方法对比图:图6的(a)是α=30°时,最优(Clairvoyant)算法、FA、EFA、JDL算法和本发明SCSA方法在距离单元50的信干噪比损失(SINR loss)曲线;图6的(b)是α=30°时,Clairvoyant算法、FA、EFA、JDL算法和SCSA方法在距离单元95的SINR loss曲线;
图7是仿真试验二的信干噪比损失变化曲线图:图7的(a)是α=90°时,clairvoyantSINR loss曲线在前视雷达前512个距离单元随距离单元数变化的曲线;图7的(b)是α=90°时,asymptotic SINR loss曲线在前视雷达前512个距离单元随距离单元数变化的曲线;图8是仿真试验二的功率、特征谱图:图8的(a)是α=90°时,原杂波MVDR功率谱图;图8的(b)是α=90°时,空-时滤波后MVDR功率谱图;图8的(c)是α=90°时,原杂波与空-时滤波后杂波的特征谱对比图;
图9是仿真试验二中本发明SCSA方法与较传统的自适应方法对比图:图9的(a)是α=90°时,Clairvoyant算法、FA、EFA、JDL算法和SCSA方法在距离单元50的SINRloss曲线;图9的(b)是α=90°时,Clairvoyant算法、FA、EFA、JDL算法和SCSA方法在距离单元95的SINR loss曲线。
具体实施方式
下面参照附图说明本发明的方法实施过程。
为了更好地理解本发明,先介绍非正侧视雷达的杂波模型。图1为非正侧视雷达的杂波模型图。M×N的平面阵安装在载机上,天线是收发共置,并且M×N个天线在行和列方向上的间隔均为d。若将每列的M个阵元均进行微波合成,则可以得到N个合成通道。这N个合成通道可以看成是N个均匀线阵,它平行于地面(图1中的xoy平面)。设雷达的工作波长为λ,脉冲重复频率为fr。在一个相干处理间隔内,对K个脉冲进行相干累积。假设载机平行于地面飞行,速度为v。载机的速度与合成的线阵之间的夹角称为偏航角α。当偏航角α=0°时,阵列为正侧视阵(SLAR);当α≠0°时,阵列为非正侧视阵(non-SLAR)。其中,非正侧视阵在实际中更为常见,特别的,当偏航角α=90°时,阵列称为前视阵(FLAR)。
假设相比较于雷达与被检测距离单元间的距离,天线的尺寸很小,即雷达满足远场条件,那么,由方位角为俯仰角为θ的杂波散射单元反射回的由第n个通道在第k个脉冲接收到的杂波信号,可以写成如下形式:
其中,n=0,1,…,N-1表示接收天线序号,k=0,1,…,K-1表示脉冲序号;R(θ)表示雷达与被检测距离单元间的斜距;表示发射幅度增益,由发射方向图产生;gn(θ)表示接收幅度增益,由列合成的接收方向图产生,且gn(θ)只是俯仰角θ的函数。表示杂波散射单元的随机幅度,与散射单元的雷达反射截面(RCS)有关;是空间频率项;是多普勒时间频率项。
假定无距离模糊时,一次雷达回波是一个距离分辨单元内的所有杂波散射单元的回波之和。杂波散射单元共有俯仰角θ,且θ位于0°到180°之间。假设沿着方位角将一个距离分辨单元划分为U个杂波散射单元,那么,第k个脉冲内的杂波回波信号的离散形式可以写成如下形式:
其中,表示第i个杂波散射单元的方位角。由上式可见,对于最简单的正侧视线阵(α=0°),空间频率项与多普勒时间频率项成正比,因而杂波在角度-多普勒平面的分布轨迹不受俯仰角θ的影响,分布轨迹成一条直线。如图3的(a)所示,正侧视雷达的杂波在角度-多普勒平面的分布轨迹是一条直线,并且彼此重合。也就是说,无论是近程杂波还是远程杂波距离单元的杂波被用作训练样本时,样本都表现出相对的均匀一致性。但是,对于非正侧视雷达(α≠0°)而言,杂波在角度-多普勒平面的分布轨迹对于不同的俯仰角θ是不同的,分布轨迹会受到θ的影响,即杂波是非均匀的。此外,由于θ在近程是快变的,而在远程相对慢变,因此,近程杂波的非均匀性会更强烈一些。如图3的(b)和图3的(c),杂波的分布轨迹随着不同距离分辨单元(俯仰角θ)的变化是不同的,杂波是非均匀的。特别地,近程杂波的非均匀性更强烈一些。此外,由图3的(b)和图3的(c)还可以看出,前视雷达(α=90°)的杂波非均匀性比非正侧视雷达(α=30°)更强烈。
根据图2,本发明的非自适应近程杂波抑制方法如下:
1.将N个雷达通道在第k个雷达脉冲时刻分别接收到的杂波信号排列成一个长为N的列向量ck,形式如下:
其中,g(θ)=diag[g0(θ),g1(θ),…,gN-1(θ)]是由gn(θ)(n=0,1,…,N-1)构成的对角阵;表示阵列空间导向矢量,表示修正的阵列空间导向矢量;是与载机运动有关的多普勒相位项;由于在一个距离分辨单元内,所有的杂波散射点的距离R(θ)都是相同的,因此R(θ)是常数,是杂波散射单元的雷达反射截面积与常数R-2(θ)的乘积。
令
则式(3)可以写成如下矩阵形式:
其中,为修正空域矩阵,它包含了U个杂波散射单元的导向矢量;F∈CU×U为增益矩阵,它包含了U个杂波散射单元的发射幅度增益;Wk∈CU×U为多普勒矩阵,它是由U个杂波散射单元的多普勒项构成的对角矩阵;是由U个杂波散射单元的回波幅度构成的列向量。
第k个和第k+1个脉冲的杂波相减后的剩余杂波能量为:
其中,表示剩余杂波的能量;Q∈CN×N是空-时滤波的系数矩阵。
最小化剩余杂波能量,有:
令J(Q)对Q的共轭导数等于0,有:
其中,上标*和上标H分别表示对矩阵或向量取共轭和取复共轭转置。
求解(11)式可以得到:
将一个相干处理间隔内接收到的的信号排列成一个长为NK的空-时向量x,为:
x=[x1,1,x2,1,…,xN,K]T∈CNK×1,其中,xn,k=cn,k+wn,k为接收到的杂波和噪声信号,wn,k为高斯白噪声,那么,空-时滤波后的数据可以表示为:
I∈CN×N表示N阶单位阵。
3.利用空-时滤波的系数矩阵进行空-时二维匹配,实现对目标的检测。假设目标信号的空-时导向矢量为其中,符号表示Kronecker积,st∈CK×1为时域的导向矢量,ss∈CN×1为空域的导向矢量。空-时滤波后的目标导向矢量为:
相应的空-时匹配的权矢量为:
其中,符号⊙表示Hadamard积,即将向量的对应元素相乘;bt∈CK×1为时域的静态权矢量,bs∈CN×1为空域的静态权矢量,二者可以用来降低滤波器旁瓣。
仿真试验对比:
为了进一步说明本发明非自适应的近程杂波抑制方法(SCSA)较传统自适应方法(如FA算法、EFA算法JDL算法等)的优越性,做如下两个仿真试验。
***模型:载机上的相控阵为16×16的平面阵,即M=N=16,阵元在行方向和列方向上的阵元间距均为d=0.1m。假设雷达主波束指向方位角俯仰角θ0=0°,且有-35dB的Chebyshev加权。雷达工作波长λ=0.2m,脉冲重复频率fr=2000Hz。在一个相干处理间隔内,有K=16个脉冲进行相干累积。载机以速度v=100m/s平行于地面飞行,飞行高度H=8000m。每个距离单元沿方位角从0°到180°,被等间隔的划分为300个杂波散射单元,相邻散射单元的间隔为0.6°。每个杂波散射单元的幅度为高斯随机幅度,且被发射方向图和距离的-2次方调制。设噪声为高斯白噪声,功率为1,杂噪比CNR=60dB。假设第一个距离分辨单元从10km处开始。两个仿真试验中将给出信干噪比损失(SINR loss)随距离单元的变化曲线来验证近程杂波的非均匀性。由于SINR loss通常用来检测一种算法的性能,所以试验中将给出近程杂波的最优信干噪比损失曲线(clairvoyant SINR loss)和样本估计下的信干噪比损失曲线(asymptotic SINR loss)。其中,clairvoyant SINR loss曲线是最优算法的性能输出,而asymptotic SINR loss曲线是是采样协方差求逆(SMI)方法的性能输出。通过比较clairvoyant SINR loss曲线和asymptotic SINR loss曲线,就可以看出因杂波的非均匀性而导致的估计杂波协方差的误差,并最终导致的自适应算法性能的下降。
试验一:一般非正侧视雷达的情况(α=30°)
图4是仿真试验一的信干噪比损失变化曲线图:图4的(a)给出的是在零方向使用理想杂波协方差计算的Clairvoyant SINR loss随距离单元数变化的曲线。图中曲线取值较低的区域对应的是杂波的主杂波区域,从图中可以看出,前200个距离单元的主杂波区位置随距离单元的变化更强烈,而后面距离单元的主杂波区位置随距离单元的变化比较缓慢,这从一个方面说明了近场杂波的非均匀性。图4的(b)中给出的是在零方向使用估计的杂波协方差计算的asymptotic SINRloss随距离单元数变化的曲线,试验中训练样本数L=512(从距离单元1到距离单元512)。从图中可以看出,asymptotic SINR loss曲线在主杂波区明显增宽,且信干噪比损失明显增加。这是由于算法的性能输出响应对应的是平均后的训练样本,而非被检测距离单元本身。asymptotic SINR loss曲线性能的下降,说明了杂波协方差估计带来的误差,也说明了近程杂波的非均匀性。
图5是仿真试验一的功率、特征谱图:图5的(a)和图5的(b)中,待检测距离单元的主瓣杂波区旁***了一个归一化多普勒频率为2fd/fr=-0.1,方位角的目标。图5的(a)给出了原杂波的MVDR功率谱,从图中可以看出,杂波的分布轨迹为半椭圆(这是由于天线背板的遮挡效应,否则杂波轨迹为整个椭圆)。与强杂波相比,杂波旁的目标很微弱。图5的(b)中给出的是经过空-时滤波后剩余杂波的MVDR功率谱,从图中可以看出,经过空-时滤波后,杂波沿着其分布轨迹被极大地抑制,只在主瓣区有一些剩余杂波,此时目标可以很容易地被检测到。图5的(c)给出的是原杂波和空时滤波后杂波的特征谱对比图,图中的大特征值部分对应杂波的功率,小特征值部分对应噪声的功率。从图中可以看出,经过空-时滤波后,大的特征值得到明显抑制,尽管噪声功率有一点提高,但是整体上杂波功率得到了抑制。
图6是仿真试验一中本发明SCSA方法与较传统的自适应方法对比图:图6的(a)和图6的(b)分别是最优(Clairvoyant)算法、FA、EFA、JDL和SCSA方法在距离单元50和距离单元95的信干噪比损失(SINR loss)的零方向曲线,其中,最优算法的曲线给出的是检测性能的上界。从这两个图可以看出,SCSA方法的性能比自适应算法有一定的提升。在主瓣杂波区2fd/fr=-0.3处,SCSA方法比FA、EFA和JDL算法在图6的(a)中分别有大约45.03dB、29.55dB和6.96dB的改善,在图6的(b)中分别有36.15dB、4.56dB和2.85dB的改善。在旁瓣杂波区2fd/fr=0.6处,SCSA方法比FA、EFA和JDL算法在图6的(a)中分别有大约2.52dB、2.02dB和1.80dB的改善,在图6的(b)中分别有1.24dB、0.96dB和0.35dB的改善。
试验二:前视雷达的情况(α=90°)
图7是仿真试验二的信干噪比损失变化曲线图:图7的(a)和图7的(b)分别给出了clairvoyant SINR loss和asymptotic SINR loss曲线在前视雷达前512个距离单元随距离单元数变化的零方向曲线。从图7的(a)中可以看出,clairvoyant SINR loss曲线在近程是随距离单元变化的,且变化强度比试验一中的斜视雷达(α=30°)剧烈。从图7的(b)中可以看出,由于样本的非均匀性,导致asymptotic SINR loss曲线变宽,性能下降。
图8中主瓣杂波旁***了一个归一化多普勒频率2fd/fr=-0.7,方位角的目标。从图8的(a)、(b)、(c)三个图中可以看出,经过空-时滤波后,杂波的能量得到了极大的抑制。因此,本发明的SCSA方法更容易检测到目标。
图9是仿真试验二中本发明SCSA方法与较传统的自适应方法对比图:图9的(a)和图9的(b)给出了前视雷达的性能对比曲线。从图中可以看出,由于近程杂波的非均匀性,使得自适应方法的性能有所下降,而SCSA方法相较于自适应方法有一定的提升。在图9的(a)中,主瓣杂波区2fd/fr=-0.8处,SCSA方法的信干躁比损失比FA、EFA和JDL算法分别有45.34dB、29.36dB和15.00dB的改善。在旁瓣区2fd/fr=0.2处,SCSA方法的信干躁比损失比FA、EFA和JDL算法分别有2.15dB、1.36dB和1.25dB的改善。在图9的(b)中,主瓣杂波区2fd/fr=-0.86处,SCSA方法比FA、EFA和JDL算法分别有大约38.81dB、20.46dB和6.52dB的改善。旁瓣区2fd/fr=0.2处,SCSA方法比FA、EFA和JDL算法分别有1.22dB、0.96dB和0.83dB的改善。
Claims (4)
1.一种非自适应的机载非正侧视雷达近程杂波抑制方法,其特征是:用一种非自适应的由空-时滤波和空-时匹配两部分级联构成的近程杂波抑制方法来抑制近程杂波;首先将雷达各通道的接收信号排成一个列向量,然后对接收信号进行空-时滤波,实现对杂波能量的抑制,最后进行空-时二维匹配,实现对目标的检测;具体实现过程如下:
1)将N个雷达通道在第k个雷达脉冲时刻分别接收到的杂波信号排列成一个长为N的列向量其中,k表示雷达脉冲序号,k=0,1,…,K-1,K表示在一个相干处理间隔内进行相干累积的脉冲个数,为修正空域矩阵,F为增益矩阵,Wk为多普勒矩阵,向量由U个杂波散射单元的回波幅度构成,U为杂波散射单元的个数,上标T表示对矩阵或向量求转置;
2)空-时滤波,以最小化剩余杂波能量:将第k个和第k+1个脉冲的杂波信号相减,得到剩余杂波能量为: 其中,Q∈CN×N是空-时滤波的系数矩阵,符号||||F表示矩阵或向量的F-范数,符号CA×B表示所有A×B维的复矩阵的集合;通过最小化剩余杂波能量得到然后由Q构造最终的空-时滤波的系数矩阵其中代表Moore-Penrose右伪逆操作,是由U个杂波散射单元的多普勒项构成的对角阵;
2.根据权利要求1所述的一种非自适应的机载非正侧视雷达近程杂波抑制方法,其特征是:对雷达各通道在各脉冲时刻接收到的杂波信号排成的列向量进行空-时滤波,将某个脉冲时刻接收到的杂波信号排成一个列向量的具体过程如下:
{1}将N个雷达通道在第k个雷达脉冲时刻分别接收到的杂波信号排列成一个长为N的列向量ck,形式如下:
其中,θ表示杂波散射单元的俯仰角,为第i个杂波散射单元的方位角,i=1,…,U;R(θ)表示雷达与被检测距离单元间的斜距,g(θ)=diag[g0(θ),g1(θ),…,gN-1(θ)]是由gn(θ)构成的对角阵,n=0,1,…,N-1,gn(θ)表示接收幅度增益,由列合成的接收方向图产生;表示阵列空间导向矢量,是空间频率项,是修正的阵列空间导向矢量,为发射幅度增益,是与载机运动有关的多普勒相位项,是多普勒时间频率项;由于在一个距离分辨单元内,所有的杂波散射点的距离R(θ)都是相同的,因此R(θ)是常数,是杂波散射单元的雷达反射截面积与常数R-2(θ)的乘积;
{2}令
则将式(1)简化为如下矩阵形式:
3.根据权利要求2所述的一种非自适应的机载非正侧视雷达近程杂波抑制方法,其特征是:对接收信号进行空时-滤波以抑制杂波:将第k个和第k+1个脉冲的杂波相减,构造空-时滤波的系数矩阵以最小化剩余杂波能量,具体过程如下:
[1]由式(3)知,第k+1个脉冲接收的杂波有如下形式:
[2]第k个和第k+1个脉冲的杂波相减后的剩余杂波能量为:
[3]最小化剩余杂波能量,有:
[4]令J(Q)对Q的共轭导数等于0,有:
其中,上标*表示对矩阵或向量取共轭,上标H表示矩阵或向量的复共轭转置;
[5]求解式(9)得到:
[6]将一个相干处理间隔内接收到的信号排列成一个长为NK的空-时向量x,为:x=[x0,0,x1,0,…,xN-1,K-1]T∈CNK×1,其中,xn,k=cn,k+wn,k,n=0,…,N-1,k=0,1,…,K-1,cn,k为接收到的杂波信号,wn,k为高斯白噪声信号,那么,空-时滤波后的数据为:
I∈CN×N表示N阶单位矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410018068.2A CN103744067A (zh) | 2014-01-15 | 2014-01-15 | 一种非自适应的机载非正侧视雷达近程杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410018068.2A CN103744067A (zh) | 2014-01-15 | 2014-01-15 | 一种非自适应的机载非正侧视雷达近程杂波抑制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103744067A true CN103744067A (zh) | 2014-04-23 |
Family
ID=50501102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410018068.2A Pending CN103744067A (zh) | 2014-01-15 | 2014-01-15 | 一种非自适应的机载非正侧视雷达近程杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103744067A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155633A (zh) * | 2014-08-12 | 2014-11-19 | 西安电子科技大学 | 一种非正侧视双基地mimo雷达的杂波抑制方法 |
CN104408276A (zh) * | 2014-09-15 | 2015-03-11 | 电子科技大学 | 一种采样远场方向图诊断阵列天线失效阵元的方法 |
CN106093870A (zh) * | 2016-05-30 | 2016-11-09 | 西安电子科技大学 | 高超声速飞行器下降段的sar‑gmti杂波抑制方法 |
CN106353732A (zh) * | 2016-08-15 | 2017-01-25 | 西安电子科技大学 | 一种基于认知的机载雷达异构杂波抑制方法 |
CN109100690A (zh) * | 2018-07-10 | 2018-12-28 | 哈尔滨工业大学 | 一种基于Hilbert-Huang变换的超宽带时域RCS测量信号处理方法 |
CN109597034A (zh) * | 2018-12-12 | 2019-04-09 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
CN110531326A (zh) * | 2018-05-24 | 2019-12-03 | 南京锐达思普电子科技有限公司 | 低慢小雷达对地杂波抑制的发射波束控制算法 |
WO2020132972A1 (zh) * | 2018-12-26 | 2020-07-02 | 深圳市大疆创新科技有限公司 | 目标检测方法、***及计算机可读存储介质 |
CN112904289A (zh) * | 2021-01-18 | 2021-06-04 | 西安电子科技大学 | 基于对角加载的机载非正侧视阵雷达稳健杂波抑制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414002A (zh) * | 2008-12-01 | 2009-04-22 | 西安电子科技大学 | 机载雷达非自适应杂波对消方法 |
CN101907702A (zh) * | 2010-06-24 | 2010-12-08 | 西安电子科技大学 | 用于mimo雷达的两维多脉冲对消器 |
CN103176168A (zh) * | 2013-02-05 | 2013-06-26 | 西安电子科技大学 | 一种机载非正侧视阵雷达近程杂波对消方法 |
-
2014
- 2014-01-15 CN CN201410018068.2A patent/CN103744067A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414002A (zh) * | 2008-12-01 | 2009-04-22 | 西安电子科技大学 | 机载雷达非自适应杂波对消方法 |
CN101907702A (zh) * | 2010-06-24 | 2010-12-08 | 西安电子科技大学 | 用于mimo雷达的两维多脉冲对消器 |
CN103176168A (zh) * | 2013-02-05 | 2013-06-26 | 西安电子科技大学 | 一种机载非正侧视阵雷达近程杂波对消方法 |
Non-Patent Citations (1)
Title |
---|
曹杨 等: "机载MIMO雷达空时自适应杂波对消器", 《航空学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155633A (zh) * | 2014-08-12 | 2014-11-19 | 西安电子科技大学 | 一种非正侧视双基地mimo雷达的杂波抑制方法 |
CN104408276A (zh) * | 2014-09-15 | 2015-03-11 | 电子科技大学 | 一种采样远场方向图诊断阵列天线失效阵元的方法 |
CN106093870A (zh) * | 2016-05-30 | 2016-11-09 | 西安电子科技大学 | 高超声速飞行器下降段的sar‑gmti杂波抑制方法 |
CN106353732A (zh) * | 2016-08-15 | 2017-01-25 | 西安电子科技大学 | 一种基于认知的机载雷达异构杂波抑制方法 |
CN106353732B (zh) * | 2016-08-15 | 2019-02-19 | 西安电子科技大学 | 一种基于认知的机载雷达异构杂波抑制方法 |
CN110531326A (zh) * | 2018-05-24 | 2019-12-03 | 南京锐达思普电子科技有限公司 | 低慢小雷达对地杂波抑制的发射波束控制算法 |
CN109100690A (zh) * | 2018-07-10 | 2018-12-28 | 哈尔滨工业大学 | 一种基于Hilbert-Huang变换的超宽带时域RCS测量信号处理方法 |
CN109597034A (zh) * | 2018-12-12 | 2019-04-09 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
CN109597034B (zh) * | 2018-12-12 | 2021-08-31 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
WO2020132972A1 (zh) * | 2018-12-26 | 2020-07-02 | 深圳市大疆创新科技有限公司 | 目标检测方法、***及计算机可读存储介质 |
CN111727380A (zh) * | 2018-12-26 | 2020-09-29 | 深圳市大疆创新科技有限公司 | 目标检测方法、***及计算机可读存储介质 |
CN112904289A (zh) * | 2021-01-18 | 2021-06-04 | 西安电子科技大学 | 基于对角加载的机载非正侧视阵雷达稳健杂波抑制方法 |
CN112904289B (zh) * | 2021-01-18 | 2022-12-02 | 西安电子科技大学 | 基于对角加载的机载非正侧视阵雷达稳健杂波抑制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103744067A (zh) | 一种非自适应的机载非正侧视雷达近程杂波抑制方法 | |
CN102721947B (zh) | 机载数字阵列雷达高效三维空时自适应杂波抑制方法 | |
CN101561497B (zh) | 机载雷达杂波抑制方法 | |
CN101556328B (zh) | 基于杂波协方差矩阵的机载雷达空时二维滤波器构建方法 | |
CN106546965B (zh) | 基于雷达幅度和多普勒频率估计的空时自适应处理方法 | |
CN105929371B (zh) | 一种基于协方差矩阵估计的机载雷达杂波抑制方法 | |
CN104360325B (zh) | 机载前视阵雷达的空时自适应处理方法 | |
CN107703490A (zh) | 基于fda‑mimo雷达的距离模糊杂波抑制方法 | |
CN104977571B (zh) | 基于俯仰频率分集stap的距离模糊杂波抑制方法 | |
CN103728606A (zh) | 机载mimo雷达的多普勒通道关联两级降维方法 | |
CN110703209B (zh) | 高重频机载前视阵雷达地面距离模糊杂波的抑制方法 | |
CN103954941B (zh) | 一种机载相控阵雷达二维多脉冲的认知杂波抑制方法 | |
CN101907703B (zh) | 天波超视距雷达空时联合自适应抗干扰方法 | |
CN102288948B (zh) | 基于stap的高速平台高速空中动目标检测方法 | |
CN103760529B (zh) | 一种无源探测高效级联空时自适应处理方法 | |
CN103412290B (zh) | 知识辅助的apr非均匀样本检测方法 | |
CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN105403864B (zh) | 基于改进斜投影的二维船载高频地波雷达海杂波抑制方法 | |
CN103018727A (zh) | 一种基于样本训练的机载雷达非平稳杂波抑制方法 | |
CN109765529B (zh) | 一种基于数字波束形成的毫米波雷达抗干扰方法及*** | |
CN104345301A (zh) | 机载mimo雷达非自适应杂波预滤波空时二维对消方法 | |
CN110082744A (zh) | 多普勒分级复用的mimo机载双基雷达杂波抑制方法 | |
CN104635219B (zh) | 基于阵元‑脉冲域补偿的匀加速平台空时自适应处理方法 | |
CN101907702A (zh) | 用于mimo雷达的两维多脉冲对消器 | |
CN103267964A (zh) | 一种基于低秩矩阵恢复的弹载导引头雷达σδ-stap方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140423 |