CN109212489B - 一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 - Google Patents
一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 Download PDFInfo
- Publication number
- CN109212489B CN109212489B CN201811240433.9A CN201811240433A CN109212489B CN 109212489 B CN109212489 B CN 109212489B CN 201811240433 A CN201811240433 A CN 201811240433A CN 109212489 B CN109212489 B CN 109212489B
- Authority
- CN
- China
- Prior art keywords
- clutter
- distance
- representing
- auxiliary pulse
- auxiliary
- 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
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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于辅助多波形脉冲FDA‑MIMO雷达距离模糊杂波抑制方法。步骤是:1.选取与主发射波形具有超低互相关特性的波形;2.根据系统参数以及设计的OFDM信号,在相干积累时间内发射辅助脉冲;3.利用辅助脉冲得到主观测区间的模糊杂波协方差矩阵;4.基于最小方差准则抑制主观测区间的距离模糊杂波;5.对剩余杂波按照传统的非正侧视方式进行杂波抑制。本发明利用辅助脉冲与主观测区间距离模糊杂波的一致性,通过辅助脉冲设计匹配滤波器克服了传统FDA‑MIMO方法固定加权的非自适应性缺陷,提高了对阵列误差的稳定度。
Description
技术领域
本发明属于雷达技术领域,是一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,可用于星载/机载非正侧视大场景下的复杂杂波抑制。
背景技术
星载雷达下视观测时,受到观测场景大幅宽以及高脉冲重复频率的影响,存在非常严重的距离模糊,空时信号处理的杂波抑制性能会受到距离模糊杂波的影响,因此需要对模糊杂波进行抑制。
一类是基于脉冲重复频率(Pulse Repetition Frequency,PRF)优化的系统方法,能够降低空时自适应处理STAP(Space-time adaptive processing)的检测盲区,但是这类方法本质上并不能抑制距离模糊杂波,同时需要多组不同PRF的CPI(Coherent ProcessingInterval)回波实现目标检测。另一类则是基于平面阵列雷达的俯仰维自由度去抑制距离模糊杂波,利用近程杂波与远程杂波在俯仰维波束域的区别。这类方法需要比较宽的俯仰维波束角度,从而区分来自远程和近程不同俯仰角度的回波。但是对于远程杂波而言,多个距离模糊单元的角度区分度非常低,很难利用俯仰维的自由度实现区分。另外,这类方法能够形成俯仰维方向图,对于模糊的杂波存在一定的抑制,但是,对于距离模糊的情况,杂波独立同分布样本难以满足。
FDA(frequency diverse array)雷达不同于传统的相控阵雷达,FDA雷达沿着阵元分布方向采用了一个存在比较小的频率增量的载频发射信号,能够使回波存在距离-角度-时间的三维依赖。FDA-MIMO(Mutiple Input Multiple Output)雷达能够将距离角度的波束扫描移至接收端,在接收端得到距离和角度依赖的可控自由度。利用这个特点,已经有几种FDA-MIMO体制下的杂波抑制方法。其中基于发射空域的模糊杂波抑制方法受到模糊杂波空域自由度的影响,抑制能力有限。俯仰维-FDA方法能够抑制距离模糊杂波,但是该方法设计的滤波器是独立于数据的,对于阵列误差比较敏感,当存在阵列误差时,杂波抑制性能存在比较大的损失。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,解决在非正侧视环境中,因距离模糊杂波引入的杂波分布严重非均匀问题。
本发明的技术方案是:一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,步骤如下:
步骤1:选取辅助脉冲波形;假定第m个阵元发射的主脉冲波形对应的信号形式为其中t代表时间,p代表主观测区间;在OFDM体制下,辅助脉冲对应的发射信号形式为其中s代表辅助脉冲,表示第个辅助脉冲相对于参考载频的偏移量,表示OFDM信号的频率步进量;
步骤4:根据最小方差准则,利用拉格朗日乘子法得到主观测区的模糊杂波抑制自适应权,利用自适应权得到第p0个主观测区间的无杂波模糊回波数;
步骤6:重复步骤3-5,直到处理完所有感兴趣的观测区间。
所述步骤1中辅助脉冲波形的选取原则为:选取与主波形脉冲具有互相关低于-70db以上的辅助脉冲波形。
所述步骤2中在雷达系统的相干处理时间间隔期间发射辅助脉冲时,每一种正交波形在辅助脉冲发射阶段只能使用一次。
表示Kronecker积运算,diag()是将一个列矢量的各个元素排列在一个方阵中的对角线上,方阵中的其余元素全为0;其中rl表示主观测无模糊区间第l个距离单元对应的斜距,Δf表示雷达发射信号的频率步进量,c代表光速,M表示发射阵元总数目,1N是N列的单位列向量。
所述步骤4中的具体过程为:利用拉格朗日乘子法得到主观测区的模糊杂波抑制自适应权
其中
p0表示当前感兴趣的主观测区间,ψ0代表目标对应的空间锥角;ru=cTp/2表示系统的最大无模糊距离;λ0是参考载频f0对应的波长;
fR(ψ0)=dR cos(ψ0)/λ0。
本发明与现有技术相比有益效果为:
本发明相比于传统的非正侧视空时自适应处理方法,利用FDA体质杂波数据在距离和空域的耦合特性,能够自适应的设计出抑制非主观测区间的距离模糊杂波滤波器,改善非均匀分布杂波对主观测区间空时滤波器的影响;同时,利用辅助脉冲与主观测区间距离模糊杂波的一致性,通过辅助脉冲设计匹配滤波器克服了已有的FDA-MIMO方法固定加权的非自适应性缺陷,提高了对阵列误差的稳定度。
附图说明
图1为基于辅助脉冲的FDA-MIMO雷达距离模糊杂波抑制算法流程图;
图2为OFDM信号的功率谱分布图;
图3为不同观测区间的发射-接收空间平面杂波分布谱;图3(a)、图3(b)、图3(c)、图3(d)分别给出了剔除p=1、p=2、p=3和p=4观测区间的剩余距离模糊杂波两维空域分布。
图4为不同辅助脉冲数目下对第4个主观测区间进行模糊杂波抑制处理的信杂噪比损失。
具体实施方式
下面对本发明实施及效果作进一步的详细描述。
本发明的使用场景为:本发明可应用于非正侧视情况下的距离模糊杂波抑制。基于辅助多波形脉冲的FDA-MIMO雷达距离模糊杂波抑制方法,克服了现有方法对距离模糊杂波抑制能力的不足,能够显著提升距离模糊杂波的协方差矩阵估计精度,改善目标的输出信杂噪比,并且对阵列误差具有稳健性。如图1所示,其实现步骤如下:
步骤1:选取与主波形脉冲具有互相关低于-70db以上的辅助脉冲波形。假定第m个阵元发射的主脉冲波形对应的信号形式为其中p(primary)代表主观测区间。在OFDM体制下,辅助脉冲对应的发射信号形式为其中s(secondary)代表辅助脉冲,表示第个辅助脉冲相对于参考载频的偏移量,t代表时间。其中表示OFDM信号的频率步进量。
步骤3:对接收到的辅助脉冲串回波数据,利用匹配滤波器对数据进行滤波预处理,如图2所示,从而得到属于第个距离单元、第个辅助脉冲的回波数据其中Tp表示脉冲重复周期,fm=f0+(m-1)Δf表示第m个阵元的发射信号载频,其中f0是参考载频,Δf是各阵元发射载频的频率步进量。联合P个协方差矩阵,得到主观测区间距离模糊杂波的协方差矩阵
表示Knonecker积运算,diag()是将一个列矢量的各个元素排列在一个方阵中的对角线上,方阵中的其余元素全为0;其中rl表示主观测无模糊区间第l个距离单元对应的斜距,Δf表示雷达发射信号的频率步进量,c代表光速,M表示发射阵元总数目,1N是N列的单位列向量。
步骤4根据最小方差准则,利用拉格朗日乘子法可以得到主观测区的模糊杂波抑制自适应权为
其中
p0表示当前感兴趣的主观测区间,ψ0代表目标对应的空间锥角。ru=cTp/2是系统的最大无模糊距离。λ0表示参考载频对应的波长。
fR(ψ0)=dR cos(ψ0)/λ0。
dT和dR分别表示发射和接收阵元的阵元间距。利用自适应权可以得到第p0个主观测区间的无杂波模糊回波数其中是第k个(k=1,2,...,K)主发射脉冲,第l个(l=1,2,...,L)主观测区间距离单元对应的回波数据。
步骤6重复步骤3-5,直到处理完所有感兴趣的观测区间。
下面通过仿真实验进一步证明本发明的效果。
图3为星载雷达下视情况下,对海面目标观测得到的剔除主观测区间之后的杂波发射-接收两维平面分布图,图3(a)-图3(d)分别给出了剔除p=1、p=2、p=3和p=4区间的剩余距离模糊杂波分布。可见,采用所提方法能够很好的剔除主观测区间的杂波,从而实现模糊杂波区的高精度协方差矩阵估计。
图4展示了采用不同辅助脉冲数目进行杂波抑制的信杂噪比损失,从图中可以发现采用一个辅助脉冲也能够达到接近多个辅助脉冲的性能,因此在实际FDA-MIMO雷达系统中,考虑到有限的时间资源,使用2个辅助脉冲即可达到模糊杂波抑制的要求。
仿真结论:仿真结果表明,相比传统的非正侧视阵空时自适应处理方法,所提方法能够估计出主观测区间的模糊杂波协方差矩阵,剔除其非均匀特性对主观测区间杂波抑制性能的影响。而与已发表的FDA-MIMO雷达相比,能够自适应的生成模糊杂波区间的滤波器,对阵列误差不敏感,具有更佳的适应性和稳定性。
Claims (4)
1.一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,其特征在于步骤如下:
步骤1:选取辅助脉冲波形;假定第m个阵元发射的主脉冲波形对应的信号形式为其中t代表时间,p代表主观测区间;在OFDM体制下,辅助脉冲对应的发射信号形式为其中s代表辅助脉冲,表示第个辅助脉冲相对于参考载频的偏移量,表示OFDM信号的频率步进量;
表示Kronecker积运算,diag()是将一个列矢量的各个元素排列在一个方阵中的对角线上,方阵中的其余元素全为0;其中rl表示主观测无模糊区间第l个距离单元对应的斜距,Δf表示雷达发射信号的频率步进量,c代表光速,M表示发射阵元总数目,1N是N列的单位列向量;
步骤4:根据最小方差准则,利用拉格朗日乘子法得到主观测区的模糊杂波抑制自适应权,利用自适应权得到第p0个主观测区间的无杂波模糊回波数;
所述步骤4中的具体过程为:利用拉格朗日乘子法得到主观测区的模糊杂波抑制自适应权
其中
p0表示当前感兴趣的主观测区间,ψ0代表目标对应的空间锥角;ru=cTp/2表示系统的最大无模糊距离;λ0是参考载频f0对应的波长;
fR(ψ0)=dRcos(ψ0)/λ0。
步骤6:重复步骤3-5,直到处理完所有感兴趣的观测区间;
2.根据权利要求1所述的一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,其特征在于:所述步骤1中辅助脉冲波形的选取原则为:选取与主波形脉冲具有互相关低于-70db以上的辅助脉冲波形。
3.根据权利要求1所述的一种基于辅助脉冲的FDA-MIMO雷达模糊杂波抑制方法,其特征在于:所述步骤2中在雷达系统的相干处理时间间隔期间发射辅助脉冲时,每一种正交波形在辅助脉冲发射阶段只能使用一次。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811240433.9A CN109212489B (zh) | 2018-10-24 | 2018-10-24 | 一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811240433.9A CN109212489B (zh) | 2018-10-24 | 2018-10-24 | 一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109212489A CN109212489A (zh) | 2019-01-15 |
CN109212489B true CN109212489B (zh) | 2020-08-14 |
Family
ID=64979802
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811240433.9A Active CN109212489B (zh) | 2018-10-24 | 2018-10-24 | 一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109212489B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110082744A (zh) * | 2019-04-24 | 2019-08-02 | 西安电子科技大学 | 多普勒分级复用的mimo机载双基雷达杂波抑制方法 |
CN110208763B (zh) * | 2019-04-28 | 2022-12-27 | 西安电子科技大学 | Fda-mimo雷达频率偏差误差的估计方法 |
CN110221283A (zh) * | 2019-05-14 | 2019-09-10 | 电子科技大学 | 一种基于fda-mimo的运动目标检测方法 |
CN111352077B (zh) * | 2019-12-20 | 2020-11-10 | 湖北工业大学 | 基于多比例分式规划的低截获频控阵mimo雷达系统的设计方法 |
CN113189554B (zh) * | 2021-04-12 | 2022-06-07 | 中国人民解放军海军航空大学航空作战勤务学院 | 雷达实测回波数据的处理方法、电子设备及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140266868A1 (en) * | 2013-03-15 | 2014-09-18 | Src, Inc. | Methods And Systems For Multiple Input Multiple Output Synthetic Aperture Radar Ground Moving Target Indicator |
CN104155633B (zh) * | 2014-08-12 | 2017-01-25 | 西安电子科技大学 | 一种非正侧视双基地mimo雷达的杂波抑制方法 |
CN105738879B (zh) * | 2016-02-29 | 2018-06-05 | 西安电子科技大学 | 基于稀疏恢复的雷达杂波空时自适应预滤波方法 |
CN107703490A (zh) * | 2017-09-29 | 2018-02-16 | 西安电子科技大学 | 基于fda‑mimo雷达的距离模糊杂波抑制方法 |
-
2018
- 2018-10-24 CN CN201811240433.9A patent/CN109212489B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109212489A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212489B (zh) | 一种基于辅助脉冲的fda-mimo雷达模糊杂波抑制方法 | |
US11454702B2 (en) | Synthetic aperture radar method and synthetic aperture radar device | |
CN102520395B (zh) | 基于双基地多输入多输出雷达的杂波抑制方法 | |
CN109581362B (zh) | 合成孔径雷达高度计在可变脉冲簇模式下的信号处理方法 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN109581352B (zh) | 一种基于毫米波雷达的超分辨测角系统 | |
CN112444811A (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
CN111007503B (zh) | 基于频率谱精确定位的运动目标聚焦和定位方法及系统 | |
CN110412568A (zh) | 基于扩展方位角相位编码的距离模糊杂波抑制方法 | |
CN104714231B (zh) | 一种基于完全互补序列和相位补偿的mimo sar成像方法 | |
CN111220955B (zh) | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 | |
CN108387884B (zh) | 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 | |
CN105223560A (zh) | 基于杂波俯仰方位谱稀疏恢复的机载雷达目标检测方法 | |
CN104635214B (zh) | 机载前视频率分集阵列雷达距离模糊杂波抑制方法 | |
CN108020817A (zh) | 基于配准的机载前视阵雷达杂波抑制方法 | |
CN109655819B (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN110412533A (zh) | 基于三维角度多普勒补偿的杂波抑制方法 | |
CN110488283A (zh) | 一种用于多通道hrws-sar通道的误差校正方法 | |
CN114895261A (zh) | 基于多频子脉冲编码阵列的杂波抑制方法 | |
CN106802408B (zh) | 基于稀疏恢复的机载非正侧阵近程杂波距离模糊抑制方法 | |
CN114325700A (zh) | 一种星载多通道sar动目标成像方法 | |
CN108761417B (zh) | 基于知识辅助最大似然的机载雷达杂波抑制方法 | |
CN114428228B (zh) | 高重频和差天线雷达导引头的杂波抑制方法 | |
CN116381639A (zh) | 一种基于频率分集阵列的sar-mmti杂波抑制方法 | |
CN103576157A (zh) | 一种基于多维波形编码的合成孔径声纳成像方法及系统 |
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 |