CN113504509A - 一种波束域补偿的匀加速机载雷达杂波抑制方法 - Google Patents
一种波束域补偿的匀加速机载雷达杂波抑制方法 Download PDFInfo
- Publication number
- CN113504509A CN113504509A CN202110639482.5A CN202110639482A CN113504509A CN 113504509 A CN113504509 A CN 113504509A CN 202110639482 A CN202110639482 A CN 202110639482A CN 113504509 A CN113504509 A CN 113504509A
- Authority
- CN
- China
- Prior art keywords
- data
- pulse
- domain
- array element
- airborne 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.)
- Granted
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
-
- 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/414—Discriminating targets with respect to background clutter
-
- 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/415—Identification of targets based on measurements of movement associated with the target
-
- 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/418—Theoretical aspects
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information 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)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种波束域补偿的匀加速机载雷达杂波抑制方法,确定机载雷达,机载雷达检测范围内存在目标,机载雷达向其检测范围内发射信号并经目标反射后接收到的回波信号,记为雷达原始回波数据矩阵;采用阵元域FFT产生多个不同指向的波束,此时每个波束序号对应一个天线锥角余弦值,可以在方位上分离信号;按照天线锥角与方位角的关系,对不同波束序号的数据进行不同补偿;对波束域补偿后的数据进行降维STAP处理,得到波束域补偿的匀加速机载雷达杂波抑制处理结果。本发明抑制方法能够减小载机飞行加速度对地面散射点回波的影响,提高杂波抑制性能。
Description
技术领域
本发明属于雷达技术领域,涉及一种波束域补偿的匀加速机载雷达杂波抑制方法。
背景技术
在现代战争中,雷达因为其全天候、全天时等优势,逐渐成为军事上必不可少的核心电子装备。其中机载雷达更是因为其探测视野广、机动性强等特点,已成为现代战争中夺取制空权和制信息权,决定现代战争胜负的重要因素。
当机载雷达在下视工作时,往往会受到地杂波的影响。且由于载机和地面之间存在相对运动,不同方向的散射点与雷达的相对速度是不同的,杂波会在多普勒域上发生扩散,杂波抑制难度相较于地基固定雷达会明显增加。空时自适应处理(STAP)技术可以在空时二维域上对杂波和目标进行区分,实现杂波抑制和目标检测。
在杂波特性已知的假设下,全维STAP的性能是所有STAP处理方法的性能上限。但是在大天线系统情况下,全维STAP涉及较高的计算复杂度,且对独立同分布(i.i.d.)样本数的需求较大,难以在实际工程中应用。因此在实际应用中,通常会使用可以大大减少计算量和i.i.d.样本需求的降维STAP方法。降维STAP是以放弃部分系统自由度,损失杂波抑制性能为代价的,它只能达到杂波抑制的准最优性能。降维STAP算法分为固定结构降维算法和自适应降维算法,前者由于其实用性而更受关注。在固定结构的降维STAP方法中,由于时域滤波相较空域滤波能获得更低的副瓣,杂波的局域化效果更好,因此多普勒后降维STAP方法在实际应用中被更多使用。多普勒后降维STAP方法只需要单独考虑每个多普勒通道的杂波分量,分别形成空域或联合域自适应滤波器,能够在较少的i.i.d.样本下得到较好的杂波抑制性能。
对于在一个CPI内变速飞行的载机,其加速度往往会使雷达天线的采样点在空间上非均匀分布,导致目标或地面散射点回波的多普勒频率在一个CPI内发生变化,传统的针对匀速直线运动机载雷达的STAP技术无法直接应用,给杂波抑制和目标检测带来新的挑战。因此,对变速飞行机载雷达杂波抑制方法进行研究是十分必要的。
发明内容
本发明的目的是提供一种波束域补偿的匀加速机载雷达杂波抑制方法,能够减小载机飞行加速度对地面散射点回波的影响,提高杂波抑制性能。
本发明所采用的技术方案是,一种波束域补偿的匀加速机载雷达杂波抑制方法,具体按照以下步骤实施:
步骤1,确定机载雷达,且机载雷达检测范围内存在目标,设机载雷达向其检测范围内发射信号并经目标反射后接收到的回波信号,记为雷达原始回波数据矩阵X,雷达原始回波数据矩阵X为N×K×L维的阵元-脉冲-距离域数据,N表示机载雷达天线的阵元数,K表示机载雷达在一个相干积累时间内发射的脉冲数,L表示机载雷达接收到的回波数据包含的距离单元总个数;
步骤2,根据阵元加权FFT的切比雪夫加权值,确定阵元FFT加权系数矢量β;
步骤3,令l=1,l∈{1,2,…,L},l表示第l个距离单元;
步骤4,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl,根据第l个距离单元的斜距Rl,确定第l个距离单元上的散射点的俯仰角根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值,Q的大小和机载雷达天线的阵元数N相等;
步骤5,根据步骤4得到的归一化空间频率值和步骤2得到的阵元FFT加权系数矢量β,确定第l个距离单元的N×K维阵元-脉冲维数据的阵元FFT加权变换矩阵Fl;使用阵元FFT加权变换矩阵Fl对所述第l个距离单元的阵元-脉冲维数据Xl进行阵元域加权FFT变换,得到第l个距离单元波束-脉冲维回波数据Yl;
步骤6,令q=1,q∈{1,2,…,Q},q表示第q个波束;
步骤7,根据机载雷达飞行的加速度a及其它雷达参数,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql,并根据补偿矢量tql对第l个距离单元、第q个波束的脉冲维回波数据yql进行补偿,得到补偿后的脉冲维数据
本发明的特征还在于,
步骤2中,阵元FFT加权系数矢量β的表达式为:
β=[β0 β1 … βN-1]T (1)
式(1)中,βn,n=1,…,N表示阵元FFT加权系数的每一项,(·)T表示转置操作。
步骤4的具体过程为:
步骤4.1,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl;
机载雷达的距离分辨率的表达式为:
式(2)中,c表示光速,B表示机载雷达发射信号的带宽;
则第l个距离单元的斜距Rl的表达式为:
Rl=(l-1)ΔR (3)
第l个距离单元上的散射点的俯仰角的表达式为:
式(4)中,sin-1(·)表示反正弦函数,H表示载机飞行高度,Re=43re表示地球等效半径,re为地球真实半径;
步骤4.3,根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值,Q的大小和机载雷达天线的阵元数N相等;
假设多波束形成的角度范围为[π-θmax,θmax],其中θmax为形成多波束的波束中心指向的最大方位角;
对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值的表达式为:
步骤5的具体过程为:
步骤5.2,使用阵元FFT加权变换矩阵Fl对第l个距离单元的阵元-脉冲维数据Xl进行阵元域加权FFT变换,得到第l个距离单元波束-脉冲维回波数据Yl,表达式为:
式(7)中,(·)H表示共轭转置操作。
步骤7的具体过程为:
步骤7.1,根据机载雷达飞行的加速度a及其它雷达参数,得到散射点的时域导向矢量bq,通过散射点的时域导向矢量bq,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql;
其中,散射点的时域导向矢量bq的表达式为:
bq=bv,q⊙ba,q (8)
式(8)中,⊙表示Hadamard积;
式(9)-(11)中,v0表示载机飞行的初始速度,a表示载机飞行的加速度,λc表示机载雷达发射信号的波长,Tr表示机载雷达发射信号的脉冲重复周期,α表示机载雷达天线的安装角,θq表示第q个波束的波束中心指向的方位角;
步骤9的具体过程为:
步骤10的具体过程为:
其中,脉冲域FFT变换矩阵为:
步骤11的具体过程为:
步骤11.1,构造第m个多普勒通道的空时降维矩阵,即:
式(17)中,Bm=[Γm-1 Γm Γm+1]表示时域降维矩阵;IN表示N×N维的单位阵,是空域降维矩阵;Γm表示脉冲域FFT变换矩阵Γ的第m列;
用邻近第l个距离单元的P个训练单元的数据来估计第l个距离单元的杂波加噪声协方差矩阵,即:
步骤11.2,构造目标导向矢量,即:
步骤11.3,对于第l个距离单元,第m个多普勒通道的数据可以进行空时降维自适应处理,得到降维后的搜索空时导向矢量,即:
降维后的第l个距离单元,第m个多普勒通道的回波数据为:
降维后的杂波加噪声协方差矩阵为:
在保证目标信号增益不变的前提下,根据线性约束最小准则,计算第l个距离单元,第m个多普勒通道的空时自适应滤波器,表达式为:
本发明的有益效果是,本发明一种波束域补偿的匀加速机载雷达杂波抑制方法,能够对载机加速度导致的散射点杂波多普勒频率变化进行补偿,且对来自不同方向的杂波回波所补偿的量不同,有效降低杂波抑制后的剩余能量水平,提高输出信杂噪比。
附图说明
图1是本发明一种波束域补偿的匀加速机载雷达杂波抑制方法的流程图;
图2是雷达原始回波数据矩阵的距离多普勒谱图;
图3是雷达原始回波数据矩阵沿脉冲维进行加窗傅里叶变换后得到的结果示意图;
图4是本发明一种波束域补偿的匀加速机载雷达杂波抑制方法得到的自适应权值图;
图5是经过本发明抑制方法抑制后再进行脉冲压缩和距离多普勒处理后得到的结果示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明提供一种波束域补偿的匀加速机载雷达杂波抑制方法,如图1所示,具体按照以下步骤实施:
步骤1,确定机载雷达,且机载雷达检测范围内存在目标,设机载雷达向其检测范围内发射信号并经目标反射后接收到的回波信号,记为雷达原始回波数据矩阵X,雷达原始回波数据矩阵X为N×K×L维的阵元-脉冲-距离域数据;
其中,N表示机载雷达天线的阵元数,K表示机载雷达在一个相干积累时间内发射的脉冲数,L表示机载雷达接收到的回波数据包含的距离单元总个数;
步骤2,根据阵元加权FFT的切比雪夫加权值,确定阵元FFT加权系数矢量β,其表达式为:
β=[β0 β1 … βN-1]T (1)
式(1)中,βn,n=1,…,N表示阵元FFT加权系数的每一项,(·)T表示转置操作;
步骤3,令l=1,l∈{1,2,…,L},l表示第l个距离单元;
步骤4,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl,根据第l个距离单元的斜距Rl,确定第l个距离单元上的散射点的俯仰角根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值;
具体过程为:
步骤4.1,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl;
机载雷达的距离分辨率的表达式为:
式(2)中,c表示光速,B表示机载雷达发射信号的带宽;
则第l个距离单元的斜距Rl的表达式为:
Rl=(l-1)ΔR (3)
第l个距离单元上的散射点的俯仰角的表达式为:
式(4)中,sin-1(·)表示反正弦函数,H表示载机飞行高度,Re=43re表示地球等效半径,re为地球真实半径;
步骤4.3,根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值,Q的大小和机载雷达天线的阵元数N相等;
假设多波束形成的角度范围为[π-θmax,θmax],其中θmax为形成多波束的波束中心指向的最大方位角;
对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值的表达式为:
步骤5,根据步骤4得到的归一化空间频率值和步骤2得到的阵元FFT加权系数矢量β,确定第l个距离单元的N×K维阵元-脉冲维数据的阵元FFT加权变换矩阵Fl;使用阵元FFT加权变换矩阵Fl对所述第l个距离单元的阵元-脉冲维数据Xl进行阵元域加权FFT变换,得到第l个距离单元波束-脉冲维回波数据Yl;
具体过程为:
步骤5.2,使用阵元FFT加权变换矩阵Fl对第l个距离单元的阵元-脉冲维数据Xl进行阵元域加权FFT变换,得到第l个距离单元波束-脉冲维回波数据Yl,表达式为:
式(7)中,(·)H表示共轭转置操作;
步骤6,令q=1,q∈{1,2,…,Q},q表示第q个波束;
步骤7,根据机载雷达飞行的加速度a及其它雷达参数,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql,并根据补偿矢量tql对第l个距离单元、第q个波束的脉冲维回波数据yql进行补偿,得到补偿后的脉冲维数据
具体过程为:
步骤7.1,根据机载雷达飞行的加速度a及其它雷达参数,得到散射点的时域导向矢量bq,通过散射点的时域导向矢量bq,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql;
其中,散射点的时域导向矢量bq的表达式为:
bq=bv,q⊙ba,q (8)
式(8)中,⊙表示Hadamard积;
式(9)-(11)中,v0表示载机飞行的初始速度,a表示载机飞行的加速度,λc表示机载雷达发射信号的波长,Tr表示机载雷达发射信号的脉冲重复周期,α表示机载雷达天线的安装角,θq表示第q个波束的波束中心指向的方位角;
具体过程为:
具体过程为:
其中,脉冲域FFT变换矩阵为:
步骤11.1,构造第m个多普勒通道的空时降维矩阵,即:
式(17)中,Bm=[Γm-1 Γm Γm+1]表示时域降维矩阵;IN表示N×N维的单位阵,是空域降维矩阵;Γm表示脉冲域FFT变换矩阵Γ的第m列;
用邻近第l个距离单元的P个训练单元的数据来估计第l个距离单元的杂波加噪声协方差矩阵,即:
步骤11.2,构造目标导向矢量,即:
步骤11.3,对于第l个距离单元,第m个多普勒通道的数据可以进行空时降维自适应处理,得到降维后的搜索空时导向矢量,即:
降维后的第l个距离单元,第m个多普勒通道的回波数据为:
降维后的杂波加噪声协方差矩阵为:
在保证目标信号增益不变的前提下,根据线性约束最小准则,计算第l个距离单元,第m个多普勒通道的空时自适应滤波器,表达式为:
仿真实验:
(一)实验参数及实验条件
本次实验采用的参数如下:
1)机载雷达天线采用1行×64列的均匀线阵阵列,阵元间距为机载雷达发射波形半波长,则可以得到大小为N×K×L的雷达回波数据;雷达阵列正侧视阵安放。
2)载机1飞行高度为6公里,水平匀速飞行,载机速度为100m/s;载机2飞行高度为6公里,水平匀加速飞行,载机初始速度为100m/s,加速度为40m/s;两个载机雷达在同一个相干处理间隔内发射450个相干积累脉冲,脉冲重复频率是1.5kHz;距离采样频率为2MHz;主波束偏扫,与载机机头的夹角为80°;地球半径为6378公里。
(二)实验内容及结果分析
A.本次实验首先对匀速飞行的载机1的机载雷达原始回波数据和匀加速飞行的载机2的机载雷达原始回波数据进行脉冲压缩,并直接进行降维STAP处理,处理结果如图2和图3所示;其中,横坐标表示信号的多普勒通道数,纵坐标表示信号的距离门数,从图2和图3可以对比看出,对匀加速飞行的载机2的机载雷达接收的雷达回波直接进行降维STAP处理,相比对匀速机载雷达回波的降维STAP处理结果,输出结果的整体能量水平会变高。
B.对匀加速飞行的载机2的机载雷达原始回波数据进行脉冲压缩,并按照本发明的流程进行处理;图4和图5分别为对匀速机载雷达回波、匀加速机载雷达回波、波束域补偿数据降维STAP处理后的杂波剩余对比图以及改善因子对比图。可以看出本发明方法可以对匀加速机载雷达回波进行有效处理,可以有效降低杂波抑制后的剩余能量水平,提高输出信杂噪比。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。
Claims (8)
1.一种波束域补偿的匀加速机载雷达杂波抑制方法,其特征在于,具体按照以下步骤实施:
步骤1,确定机载雷达,且机载雷达检测范围内存在目标,设机载雷达向其检测范围内发射信号并经目标反射后接收到的回波信号,记为雷达原始回波数据矩阵X,雷达原始回波数据矩阵X为N×K×L维的阵元-脉冲-距离域数据,N表示机载雷达天线的阵元数,K表示机载雷达在一个相干积累时间内发射的脉冲数,L表示机载雷达接收到的回波数据包含的距离单元总个数;
步骤2,根据阵元加权FFT的切比雪夫加权值,确定阵元FFT加权系数矢量β;
步骤3,令l=1,l∈{1,2,…,L},l表示第l个距离单元;
步骤4,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl,根据第l个距离单元的斜距Rl,确定第l个距离单元上的散射点的俯仰角根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值,Q的大小和机载雷达天线的阵元数N相等;
步骤5,根据步骤4得到的归一化空间频率值和步骤2得到的阵元FFT加权系数矢量β,确定第l个距离单元的N×K维阵元-脉冲维数据的阵元FFT加权变换矩阵Fl;使用阵元FFT加权变换矩阵Fl对所述第l个距离单元的阵元-脉冲维数据Xl进行阵元域加权FFT变换,得到第l个距离单元波束-脉冲维回波数据Yl;
步骤6,令q=1,q∈{1,2,…,Q},q表示第q个波束;
步骤7,根据机载雷达飞行的加速度a及其它雷达参数,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql,并根据补偿矢量tql对第l个距离单元、第q个波束的脉冲维回波数据yql进行补偿,得到补偿后的脉冲维数据
2.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法,其特征在于,所述步骤2中,阵元FFT加权系数矢量β的表达式为:
β=[β0 β1 … βN-1]T (1)
式(1)中,βn,n=1,…,N表示阵元FFT加权系数的每一项,(·)T表示转置操作。
3.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法,其特征在于,步骤4的具体过程为:
步骤4.1,对第l个距离单元的N×K维阵元-脉冲维数据Xl,根据机载雷达的距离分辨率,确定第l个距离单元的斜距Rl;
机载雷达的距离分辨率的表达式为:
式(2)中,c表示光速,B表示机载雷达发射信号的带宽;
则第l个距离单元的斜距Rl的表达式为:
Rl=(l-1)ΔR (3)
第l个距离单元上的散射点的俯仰角的表达式为:
式(4)中,sin-1(·)表示反正弦函数,H表示载机飞行高度,Re=4/3re表示地球等效半径,re为地球真实半径;
步骤4.3,根据第l个距离单元上的散射点的俯仰角和多波束形成的角度范围,确定对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值,Q的大小和机载雷达天线的阵元数N相等;
假设多波束形成的角度范围为[π-θmax,θmax],其中θmax为形成多波束的波束中心指向的最大方位角;
对第l个距离单元的N×K维阵元-脉冲维数据Xl进行阵元加权FFT的Q个归一化空间频率值的表达式为:
5.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法,其特征在于,步骤7的具体过程为:
步骤7.1,根据机载雷达飞行的加速度a及其它雷达参数,得到散射点的时域导向矢量bq,通过散射点的时域导向矢量bq,确定对第l个距离单元、第q个波束的K×1维杂波数据yql的K×1维补偿矢量tql;
其中,散射点的时域导向矢量bq的表达式为:
bq=bv,q⊙ba,q (8)
式(8)中,⊙表示Hadamard积;
式(9)-(11)中,v0表示载机飞行的初始速度,a表示载机飞行的加速度,λc表示机载雷达发射信号的波长,Tr表示机载雷达发射信号的脉冲重复周期,α表示机载雷达天线的安装角,θq表示第q个波束的波束中心指向的方位角;
8.根据权利要求1所述的一种波束域补偿的匀加速机载雷达杂波抑制方法,其特征在于,步骤11的具体过程为:
步骤11.1,构造第m个多普勒通道的空时降维矩阵,即:
式(17)中,Bm=[Γm-1 Γm Γm+1]表示时域降维矩阵;IN表示N×N维的单位阵,是空域降维矩阵;Γm表示脉冲域FFT变换矩阵Γ的第m列;
用邻近第l个距离单元的P个训练单元的数据来估计第l个距离单元的杂波加噪声协方差矩阵,即:
步骤11.2,构造目标导向矢量,即:
步骤11.3,对于第l个距离单元,第m个多普勒通道的数据可以进行空时降维自适应处理,得到降维后的搜索空时导向矢量,即:
降维后的第l个距离单元,第m个多普勒通道的回波数据为:
降维后的杂波加噪声协方差矩阵为:
在保证目标信号增益不变的前提下,根据线性约束最小准则,计算第l个距离单元,第m个多普勒通道的空时自适应滤波器,表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110639482.5A CN113504509B (zh) | 2021-06-08 | 2021-06-08 | 一种波束域补偿的匀加速机载雷达杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110639482.5A CN113504509B (zh) | 2021-06-08 | 2021-06-08 | 一种波束域补偿的匀加速机载雷达杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113504509A true CN113504509A (zh) | 2021-10-15 |
CN113504509B CN113504509B (zh) | 2023-07-11 |
Family
ID=78009594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110639482.5A Active CN113504509B (zh) | 2021-06-08 | 2021-06-08 | 一种波束域补偿的匀加速机载雷达杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113504509B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114779198A (zh) * | 2022-04-24 | 2022-07-22 | 中国人民解放军空军预警学院 | 共形阵机载雷达空时杂波谱自适应补偿与杂波抑制方法 |
CN114966568A (zh) * | 2022-05-25 | 2022-08-30 | 西安电子科技大学 | 一种权矢量时变的匀加速飞行雷达空时自适应处理方法 |
CN116106852A (zh) * | 2023-04-12 | 2023-05-12 | 中国人民解放军63921部队 | 机载主杂波频道的确定方法、装置及电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107703490A (zh) * | 2017-09-29 | 2018-02-16 | 西安电子科技大学 | 基于fda‑mimo雷达的距离模糊杂波抑制方法 |
WO2018045566A1 (zh) * | 2016-09-09 | 2018-03-15 | 深圳大学 | 基于压缩感知的随机脉冲多普勒雷达角度-多普勒成像方法 |
-
2021
- 2021-06-08 CN CN202110639482.5A patent/CN113504509B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018045566A1 (zh) * | 2016-09-09 | 2018-03-15 | 深圳大学 | 基于压缩感知的随机脉冲多普勒雷达角度-多普勒成像方法 |
CN107703490A (zh) * | 2017-09-29 | 2018-02-16 | 西安电子科技大学 | 基于fda‑mimo雷达的距离模糊杂波抑制方法 |
Non-Patent Citations (1)
Title |
---|
刘锦辉;廖桂生;王泉;: "机载前视阵雷达阵元误差稳健的空时插值补偿方法", 系统工程与电子技术, no. 12 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114779198A (zh) * | 2022-04-24 | 2022-07-22 | 中国人民解放军空军预警学院 | 共形阵机载雷达空时杂波谱自适应补偿与杂波抑制方法 |
CN114966568A (zh) * | 2022-05-25 | 2022-08-30 | 西安电子科技大学 | 一种权矢量时变的匀加速飞行雷达空时自适应处理方法 |
CN116106852A (zh) * | 2023-04-12 | 2023-05-12 | 中国人民解放军63921部队 | 机载主杂波频道的确定方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN113504509B (zh) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113504509A (zh) | 一种波束域补偿的匀加速机载雷达杂波抑制方法 | |
CN108761419B (zh) | 基于组合空时主通道自适应处理的低空风切变风速估计方法 | |
CN104977571B (zh) | 基于俯仰频率分集stap的距离模糊杂波抑制方法 | |
CN110554391B (zh) | 基于ddd-gmb的低空风切变风速估计方法 | |
CN113253223B (zh) | 基于步进频信号的非平稳杂波抑制的目标检测方法 | |
CN111913157A (zh) | 一种基于雷达信号空-时去相关模型的海杂波抑制方法 | |
CN111220955A (zh) | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 | |
Zhu et al. | Deceptive jamming on space-borne SAR using frequency diverse array | |
CN110632573A (zh) | 一种机载宽带雷达空时二维keystone变换方法 | |
CN113917424B (zh) | 一种基于ddma与盲源分离的地基mimo雷达地杂波抑制方法 | |
CN111220986A (zh) | 回波功率筛选与dlcd辅助的低空风切变风速估计方法 | |
CN113534055B (zh) | 一种插值补偿的匀加速机载雷达杂波抑制方法 | |
CN113376603A (zh) | 宽带机载相控阵雷达的子带空时自适应处理方法 | |
CN110471040B (zh) | 一种基于fda天线的逆合成孔径雷达干扰方法 | |
CN115436943B (zh) | 基于脉内和脉间联合相位调制的可重构mimo-sar回波分离方法 | |
CN115144854A (zh) | 一种基于频率分集阵列的高分宽幅sar成像方法 | |
CN115407335A (zh) | 基于OFDM-chirp信号和DBF处理的MIMO-SAR解模糊方法及装置 | |
CN115113208A (zh) | 一种基于杂波特征精确认知的连续波雷达杂波对消方法 | |
CN114325595A (zh) | 联合发射阵列与脉冲编码的机载mimo-sar距离模糊分离方法 | |
CN113253222B (zh) | 机载fda-mimo双基地雷达距离模糊杂波抑制及降维搜索方法 | |
Zhu et al. | Multiple moving target resolution and imaging based on ISAR principle | |
CN117784078B (zh) | 机载雷达空时极化联合自适应处理杂波抑制方法与装置 | |
CN117269915A (zh) | 一种天基预警雷达波束通道联合降维方法 | |
CN115113161A (zh) | 一种目标补偿的匀加速飞行雷达空时自适应处理方法 | |
CN115343702B (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 |