CN114779199A - 端射阵机载雷达自适应互耦补偿与杂波抑制方法 - Google Patents
端射阵机载雷达自适应互耦补偿与杂波抑制方法 Download PDFInfo
- Publication number
- CN114779199A CN114779199A CN202210432661.6A CN202210432661A CN114779199A CN 114779199 A CN114779199 A CN 114779199A CN 202210432661 A CN202210432661 A CN 202210432661A CN 114779199 A CN114779199 A CN 114779199A
- Authority
- CN
- China
- Prior art keywords
- clutter
- representing
- airborne radar
- time
- mutual coupling
- 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
- 230000001629 suppression Effects 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000006880 cross-coupling reaction Methods 0.000 title claims abstract description 21
- 230000003044 adaptive effect Effects 0.000 title claims description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000010168 coupling process Methods 0.000 claims abstract description 36
- 238000005859 coupling reaction Methods 0.000 claims abstract description 36
- 230000008878 coupling Effects 0.000 claims abstract description 34
- 239000013598 vector Substances 0.000 claims abstract description 32
- 238000009826 distribution Methods 0.000 claims description 5
- 238000002592 echocardiography Methods 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000001808 coupling effect Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004148 unit process Methods 0.000 description 1
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/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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,包括,步骤1:构造互耦情况下的端射阵机载雷达杂波回波数据;步骤2:对步骤1构造的数据进行空时滑窗处理,并估计其峰值功率;步骤3:将步骤2估计的峰值功率与当前时刻的峰值功率进行比较处理;步骤4:利用步骤3的处理结果计算得到归一化阻抗矩阵;步骤5:利用步骤4的处理结果,对空域导向矢量失配进行补偿;步骤6:利用步骤4的处理结果,构造杂波协方差矩阵;步骤7:根据步骤6的处理结果形成空时自适应权值,并针对步骤5的结果进行杂波抑制处理。
Description
技术领域
本发明属于信号处理技术领域,更具体地,涉及一种端射阵机载雷达自适应互耦补偿与杂波抑制方法。
背景技术
端射阵天线是指最大辐射方向指向阵列排布轴向的一类天线,因其具有独特的低剖面特性和强定向辐射特性,所以在机载预警雷达远距离探测领域具有广阔的应用前景,近年来受到广泛关注。由于端射阵在最大辐射方向的方向系数不再与等效口径尺寸成正比,这便有效地解决了传统侧射阵天线由于大阵列口径尺寸而导致风阻较大的问题,因而特别适合应用于要求小风阻、低安装高度的平台,尤其是各种高速移动载体如飞机、战车等。例如美国研发的E-737“楔尾”预警机,便是采用端射阵天线进行前后向补盲。
在理想情况下我们通常假设天线各阵元为点源,而在实际情况下阵元有一定的物理尺寸,此时互耦效应不可避免。互耦效应是阵列天线固有的重要特性之一,它改变了孤立单元应该有的电流分布,不但对天线阵列的增益、方向图、波束宽度等参数产生一定的影响,而且还会影响到阵列天线发射信号和接收信号的幅度和相位。当阵列的结构和组成阵列的天线单元确定后,无论馈以何种信号序列,即无论该阵列工作在侧射阵状态还是端射阵状态,其阻抗矩阵都是相同的。唯一的区别是侧射激励下各阶耦合间存在一种近似相抵的关系,而端射激励下各阶耦合则是同相叠加的关系。因此端射阵互耦效应更加明显。
现有解决机载雷达互耦问题的方法主要有两种,一是假设归一化阻抗矩阵已知,在此基础上形成互耦补偿矩阵实现对互耦效应的有效补偿;二是基于实时测量得到的空域导向矢量形成空时自适应权值,实现存在互耦情况下的杂波抑制和目标有效检测。但是上述两种方法均假设互耦效应引起的归一化阻抗矩阵已知,而在实际工程中互耦效应会随时间的变化而变化,通常无法实时精确已知。此外,基于杂波协方差矩阵重构的互耦效应补偿方法没有考虑距离模糊的影响。
因此,急需发明一种有效的端射阵机载雷达自适应互耦补偿与杂波抑制方法,以提高端射阵机载雷达对运动目标的检测性能。
发明内容
为此,本发明提供一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,用以克服现有技术中存在的问题。
为实现上述目的,本发明提供一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,包括如下步骤,
步骤1:构造互耦情况下的端射阵机载雷达杂波回波数据;
步骤2:对步骤1构造的数据进行空时滑窗处理,并估计其峰值功率;
步骤3:将步骤2估计的峰值功率与当前时刻的峰值功率进行比较处理;
步骤4:利用步骤3的处理结果计算得到归一化阻抗矩阵;
步骤5:利用步骤4的处理结果,对空域导向矢量失配进行补偿;
步骤6:利用步骤4的处理结果,构造杂波协方差矩阵;
步骤7:根据步骤6的处理结果形成空时自适应权值,并针对步骤5的结果进行杂波抑制处理;
其中随时间变化的互耦误差矢量为
表示矩阵转置操作,j表示虚数单位,和分别表示幅度误差和相位误差,其中幅度误差服从高斯分布,相位误差服从均匀分布,即,;Z ij (i=1,2,...,M; j=1,2,...,M)表示第i个阵元和第j个阵元之间由互耦引起的互阻抗,Z ii 表示第i个天线阵元的自阻抗;表示相距n-1个列间距的两列的第i个阵元和第j个阵元之间的互阻抗,表示每个端口所接阻抗已知的固定负载;
则第q次搜索得到的存在互耦情况下的端射阵机载雷达杂波回波信号为
表示Hadamard乘积,和分别表示第l个距离单元第i个杂波块的时域导向矢量和空域导向矢量,N c表示杂波块的数目,N r表示模糊距离数,表示考虑互耦后的第q次搜索得到的第l个距离单元第i个杂波块回波幅度组成的列矢量。
进一步地,在所述步骤3中,判断的大小,其中表示当前时刻回波数据通过滑窗处理后估计得到的在空时平面上的峰值功率,表示步骤2搜索得到的峰值功率;若,则停止搜索,否则进行第q+1次搜索,直至满足要求为止,最终搜索到。
进一步地,在所述步骤5中,形成的互耦补偿矩阵为
其中补偿后的回波信号为
其中
进一步地,在所述步骤6中,构造的杂波协方差矩阵为
进一步地,在所述步骤7中,形成的空时自适应权值为
其中
经过杂波抑制处理后的输出信号为
(1)本发明通过提供一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,首先通过回波数据自适应迭代估计端射面阵的归一化阻抗矩阵,其次利用互耦补偿矩阵对互耦导致的空域导向矢量失配进行补偿,再次利用估计得到的归一化阻抗矩阵和系统参数构造杂波协方差矩阵,最后形成空时自适应权值,进行杂波抑制处理,提高了端射阵机载雷达在互耦误差情况下的运动目标检测性能;
(2)本发明以空时平面上的峰值功率值作为参考,通过自适应迭代方式估计互耦参数,峰值功率值同时包含了回波幅度和相位信息;
(3)本发明对互耦导致的空域导向矢量失配现象进行了补偿;
(4)本发明对于互耦导致的回波功率失配和端射阵特有的杂波非平稳分布现象,通过构造杂波协方差矩阵的方式进行了有效解决。
附图说明
图1为本发明所述端射阵机载雷达自适应互耦补偿与杂波抑制方法的结构框图;
图2为本发明所述端射阵机载雷达自适应互耦补偿与杂波抑制方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参阅图1-2所示,本发明通过提供一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,包括如下步骤:
步骤1:构造互耦情况下的端射阵机载雷达杂波回波数据;
其中随时间变化的互耦误差矢量为
表示矩阵转置操作,j表示虚数单位,和分别表示幅度误差和相位误差,其中幅度误差服从高斯分布,相位误差服从均匀分布,即,;Z ij (i=1,2,...,M; j=1,2,...,M)表示第i个阵元和第j个阵元之间由互耦引起的互阻抗,Z ii 表示第i个天线阵元的自阻抗;表示相距n-1个列间距的两列的第i个阵元和第j个阵元之间的互阻抗,表示每个端口所接阻抗已知的固定负载;
则第q次搜索得到的存在互耦情况下的端射阵机载雷达杂波回波信号为
表示Hadamard乘积,和分别表示第l个距离单元第i个杂波块的时域导向矢量和空域导向矢量,N c表示杂波块的数目,N r表示模糊距离数,表示考虑互耦后的第q次搜索得到的第l个距离单元第i个杂波块回波幅度组成的列矢量;
步骤2:对步骤1构造的数据进行空时滑窗处理,并估计其峰值功率;
步骤3:将步骤2估计的峰值功率与当前时刻的峰值功率进行比较处理;
峰值功率比较单元判断的大小,其中表示当前时刻回波数据通过滑窗处理后估计得到的在空时平面上的峰值功率,表示步骤2搜索得到的峰值功率;若,则停止搜索,否则进行第q+1次搜索,直至满足要求为止,最终搜索到。
步骤4:利用步骤3的处理结果计算得到归一化阻抗矩阵;
步骤5:利用步骤4的处理结果,对空域导向矢量失配进行补偿;
空时导向矢量失配补偿单元形成的互耦补偿矩阵为
其中补偿后的回波信号为
其中
步骤6:利用步骤4的处理结果,构造杂波协方差矩阵;
杂波协方差矩阵构造单元中构造的杂波协方差矩阵为
步骤7:根据步骤6的处理结果形成空时自适应权值,并针对步骤5的结果进行杂波抑制处理;
杂波抑制处理单元形成的空时自适应权值为
其中
经过杂波抑制处理后的输出信号为
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种端射阵机载雷达自适应互耦补偿与杂波抑制方法,其特征在于,包括如下步骤,
步骤1:构造互耦情况下的端射阵机载雷达杂波回波数据;
步骤2:对步骤1构造的数据进行空时滑窗处理,并估计其峰值功率;
步骤3:将步骤2估计的峰值功率与当前时刻的峰值功率进行比较处理;
步骤4:利用步骤3的处理结果计算得到归一化阻抗矩阵;
步骤5:利用步骤4的处理结果,对空域导向矢量失配进行补偿;
步骤6:利用步骤4的处理结果,构造杂波协方差矩阵;
步骤7:根据步骤6的处理结果形成空时自适应权值,并针对步骤5的结果进行杂波抑制处理;
其中随时间变化的互耦误差矢量为
表示矩阵转置操作,j表示虚数单位,和分别表示幅度误差和相位误差,其中幅度误差服从高斯分布,相位误差服从均匀分布,即,;Z ij (i=1,2,...,M; j=1,2,...,M)表示第i个阵元和第j个阵元之间由互耦引起的互阻抗,Z ii 表示第i个天线阵元的自阻抗;表示相距n-1个列间距的两列的第i个阵元和第j个阵元之间的互阻抗,表示每个端口所接阻抗已知的固定负载;
则第q次搜索得到的存在互耦情况下的端射阵机载雷达杂波回波信号为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432661.6A CN114779199A (zh) | 2022-04-24 | 2022-04-24 | 端射阵机载雷达自适应互耦补偿与杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432661.6A CN114779199A (zh) | 2022-04-24 | 2022-04-24 | 端射阵机载雷达自适应互耦补偿与杂波抑制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114779199A true CN114779199A (zh) | 2022-07-22 |
Family
ID=82433790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210432661.6A Pending CN114779199A (zh) | 2022-04-24 | 2022-04-24 | 端射阵机载雷达自适应互耦补偿与杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114779199A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116047462A (zh) * | 2023-03-31 | 2023-05-02 | 中国人民解放军空军预警学院 | 端射阵机载雷达最优阵元数目和阵元间距选取方法与装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116311A (zh) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 |
CN110374804A (zh) * | 2019-07-03 | 2019-10-25 | 西安交通大学 | 一种基于深度确定性策略梯度补偿的变桨距控制方法 |
CN113253223A (zh) * | 2021-03-30 | 2021-08-13 | 北京理工大学 | 基于步进频信号的非平稳杂波抑制的目标检测方法 |
-
2022
- 2022-04-24 CN CN202210432661.6A patent/CN114779199A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116311A (zh) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 |
CN110374804A (zh) * | 2019-07-03 | 2019-10-25 | 西安交通大学 | 一种基于深度确定性策略梯度补偿的变桨距控制方法 |
CN113253223A (zh) * | 2021-03-30 | 2021-08-13 | 北京理工大学 | 基于步进频信号的非平稳杂波抑制的目标检测方法 |
Non-Patent Citations (1)
Title |
---|
YUANYI XIONG 等: "Adaptive Mutual Coupling Compensation Method for Airborne STAP Radar With End-Fire Array", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116047462A (zh) * | 2023-03-31 | 2023-05-02 | 中国人民解放军空军预警学院 | 端射阵机载雷达最优阵元数目和阵元间距选取方法与装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111308437B (zh) | 一种毫米波mimo交通雷达求熵解速度模糊的方法 | |
CN111537966B (zh) | 一种适用于毫米波车载雷达领域的阵列天线误差校正方法 | |
CN109975807B (zh) | 一种适用于毫米波车载雷达的降维子空间测角方法 | |
CN110412559B (zh) | 分布式无人机mimo雷达的非相参融合目标检测方法 | |
CN107831480B (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN105929371A (zh) | 一种基于协方差矩阵估计的机载雷达杂波抑制方法 | |
CN101226236A (zh) | 基于旁瓣约束的天波超视距雷达自适应干扰抑制方法 | |
CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN109375179A (zh) | 一种基于稀疏表示的机载雷达近程强杂波抑制方法 | |
CN115469286B (zh) | 基于毫米波汽车雷达最小冗余mimo阵列的超分辨测角方法 | |
CN107490788A (zh) | 一种适于mimo机载雷达非平稳杂波抑制的空时自适应处理方法 | |
CN113238211B (zh) | 一种干扰条件下参数化自适应阵列信号检测方法与系统 | |
CN115372925B (zh) | 一种基于深度学习的阵列稳健自适应波束形成方法 | |
Shen et al. | Efficient adaptive approach for airborne radar short-range clutter suppression | |
CN114779199A (zh) | 端射阵机载雷达自适应互耦补偿与杂波抑制方法 | |
CN106125039B (zh) | 基于局域联合处理的改进空时自适应单脉冲测角方法 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN108896963B (zh) | 机载雷达空时自适应降维处理方法 | |
CN108196238B (zh) | 高斯背景下基于自适应匹配滤波的杂波图检测方法 | |
CN109633600B (zh) | 一种最小冗余线阵mimo-othr的doa估计方法 | |
CN118276026A (zh) | 一种基于能量轨迹提取的海面运动目标航向角估计方法 | |
CN109597034B (zh) | 一种基于欧几里得距离的空时自适应处理方法 | |
CN116106841A (zh) | 一种机载双基地mimo雷达高分宽幅成像方法 | |
CN113126076B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20220722 |