一种两维数字阵列雷达快速干扰测向方法技术

本发明专利技术涉及一种两维数字阵列雷达快速干扰测向方法，通过增加雷达发射静默周期，侦听周围干扰，同时通过对接收的阵元数据进行分析，设置合理的门限阈值，判定是否有干扰进入，同时通过FFT的方法快速形成方向图，进行空域滤波搜索，从而快速获得方位及俯仰维度某一纬度的指向角度，之后通过对干扰入射方向的数据进行互相关估计，完成两个维度的干扰角度配对处理，获得配对后的干扰入射角的方位、俯仰两维角度。维角度。维角度。

【技术实现步骤摘要】
一种两维数字阵列雷达快速干扰测向方法


[0001]本专利技术属于两维数字阵列雷达
，具体涉及的技术包括方向图形成技术，波达方向估计技术，角度匹配技术。该专利技术可以快速测量来自空间的干扰方向，通过数字波束形成技术，形成指向干扰的波束，并测量干扰方向，为后续的干扰处理提供干扰指向作为先验知识。

技术介绍

[0002]随着技术的发展，军用雷达工作的电磁环境日益复杂，不但有来自空间的敌方干扰，还有来自本方其他电子设备的电磁辐射，在检测目标的同时，需要高效的对抗空间干扰。于此同时，大规模数字处理器件的快速发展，更大规模的数字阵列雷达，越来越多的进入到了工程研制阶段，为干扰测向提供了更多的可能性。
[0003]干扰入射方向的估计，属于波达方向估计的范畴。传统的雷达，由于系统的自由度不够，能够采用的波达方向估计技术并不多。在数字阵列中对于波达方向的估计，如MUISC，ESPRIT等算法，这些算法虽然测角精度高，但是考虑到雷达工作的波束驻留时间很短，工程应用中受限于运算时间的要求难以实现，而对于两维数字阵列，其复杂度进一步提高。
[0004]对于干扰的处理，首先要获得干扰的方向，对于一维数字阵列，估计获得的干扰指向角θ0，即可认为是干扰角度；对于两维数字阵列，如果有一个干扰时，获得干扰的方位指向角θ0，干扰的俯仰指向角那么即为干扰指向。但是对于大于两个干扰的时候，情况就复杂了很多，假设两个干扰方位指向角θ0，θ1；两个干扰的俯仰指向角那么对于最终的干扰角度就有了4种排列组合，对于最终的干扰角度就有了4种排列组合，从4个可能的指向中匹配出正确的指向，即角度配对问题成了要解决的实际问题。而对于更多干扰的情况，复杂度进一步提升。

技术实现思路

[0005]要解决的技术问题
[0006]为了解决上述两维数字阵列雷达中干扰角度测量及角度配对在工程中的应用问题，从而避免复杂的计算，减少运算时间。本专利技术提出一种两维数字阵列雷达快速干扰测向方法。
[0007]技术方案
[0008]一种两维数字阵列雷达快速干扰测向方法，其特征在于步骤如下：
[0009]步骤1：搭建数字阵列采样处理平台，获得每一个阵元完成采样及下变频后的基带信号，其全阵面接收的第n个采样快拍的数据矩阵为：
[0010][0011]其中：
[0012][0013]p，q为阵元坐标号，S(n)，N(n)，分别为采样信号与噪声，n＝0,1,2.....为采样快拍号，d
x
，d
y
为阵元间距；
[0014]步骤2：对两维子阵平面中，选取一个“十字架”，即中间的一行及一列，进行干扰测向；对于单个阵元，选取某一列的所有行，即：
[0015][0016]则在方位维可近似认为全向波束，俯仰维为一维线阵，可以进行俯仰维干扰角度测向；
[0017]同样的，选取某一列的所有行，即
[0018]X
p
(n)＝(X
0,p
(n) X
1,p
(n)...X
Q-1,p
(n))
ꢀꢀ
(3)
[0019]则在俯仰维可近似认为全向波束，方位维为一维线阵，可以进行方位维干扰角度测向；
[0020]步骤3：设计雷达系统工作时序，增加静默侦查区，所述的静默侦查区不辐射电磁信号，仅接收信号；
[0021]步骤4：在确定周围环境无干扰时，对行或列采样的快拍进行FFT，记录雷达正常工作时，雷达静默区的接收噪声，并设定适当门限为k；正常工作时，如雷达静默区侦听到的信号高于门限，判定为系统受到干扰，假设干扰的个数为M个；
[0022]步骤5：利用FFT对式2进行方向图形成，“列”单次采样快拍进行FFT运算，可以获得俯仰维指向干扰角度的FFT峰值点，记录其幅度是否大于接收系统噪声门限k；大于则记录干扰角度为：
[0023][0024]FFT后第m个谱峰对应的空间采样点为k
m
，n_FFT为FFT的点数；记录FFT后第m个谱峰的值
[0025]步骤6：利用FFT对式3进行方向图形成，对“行”单次采样快拍进行FFT运算，可以获得方位指向干扰方向角度的FFT峰值点，记录其幅度是否大于接收系统噪声门限k；大于则记录干扰角度为：
[0026]其中m＝0,1...M
[0027]FFT后第m个谱峰对应的空间采样点为k
m
，n_FFT为FFT的点数。记录FFT后第m个谱峰的值A
θ(m)
；
[0028]步骤7：假设共采样N个快拍，X
q
(n),n＝0,1....N，为天线阵面“列”采样点，X
p
(n),n＝0,1....N，为天线阵面“行”采样点；重复过程步骤5，计算获得俯仰维M个干扰方向的N次采样快拍的记录值，即M个干扰的入射角度为：
[0029][0030]此外则测得的FFT后的峰值位置也应一致；N个快拍M个干扰入射对应的谱峰位置的信号为：
[0031][0032]重复步骤6，计算获得方位维M个干扰方向的N次采样快拍的记录值：
[0033]θ
M
＝[θ(0)θ(1) ... θ(M)]ꢀꢀ
(7)对应的谱峰位置的信号为：
[0034][0035]步骤8：对俯仰维谱峰信号的矩阵及方位维谱峰信号的矩阵进行自相关运算；获得相关矩阵：
[0036][0037]步骤9：构造与相关矩阵对应的角度匹配矩阵，该矩阵并不计算具体数值，而是通过对矩阵θ1×
M
，进行与R
M
×
M
矩阵相同的矩阵运算过程，生成角度查找表；如下式：
[0038][0039]即S(i,j)该点的方位角为θ(i)，俯仰角为
[0040]相关矩阵R
M
×
M
幅度越大则表明相关性越大，即就是俯仰与方位角度配对程度越高；针对M个干扰通过对相关矩阵R
MXM
幅度排序后，获得M个极值；通过M个极值在R
M
×
M
矩阵中的获得该极值数据的下标(i,j)，对比相关矩阵R
M
×
M
与查找表S
M
×
M
的对应关系，从矩阵S
M
×
M
中，根据坐标(i,j)，并找到M个配对的干扰角度完成干扰角度配对。
[0041]有益效果
[0042]本专利技术提出的一种两维数字阵列雷达快速干扰测向方法，通过增加雷达发射静默周期，侦听周围干扰，同时通过对接收的阵元数据进行分析，设置合理的门限阈值，判定是否有干扰进入，同时通过FFT的方法快速形成方向图，进行空域滤波搜索，从而快速获得方位及俯仰维度某一纬度的指向角度，之后通过对干扰入射方向的数据进行互相关估计，完成两个维度的干扰角度配对处理，获得配对后的干扰入射角的方位、俯仰两维角本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种两维数字阵列雷达快速干扰测向方法，其特征在于步骤如下：步骤1：搭建数字阵列采样处理平台，获得每一个阵元完成采样及下变频后的基带信号，其全阵面接收的第n个采样快拍的数据矩阵为：其中：p，q为阵元坐标号，S(n)，N(n)，分别为采样信号与噪声，n＝0,1,2.....为采样快拍号，d
x
，d
y
为阵元间距；步骤2：对两维子阵平面中，选取一个“十字架”，即中间的一行及一列，进行干扰测向；对于单个阵元，选取某一列的所有行，即：则在方位维可近似认为全向波束，俯仰维为一维线阵，可以进行俯仰维干扰角度测向；同样的，选取某一列的所有行，即X
p
(n)＝(X
0,p
(n) X
1,p
(n) ... X
Q-1,p
(n))
ꢀꢀ
(3)则在俯仰维可近似认为全向波束，方位维为一维线阵，可以进行方位维干扰角度测向；步骤3：设计雷达系统工作时序，增加静默侦查区，所述的静默侦查区不辐射电磁信号，仅接收信号；步骤4：在确定周围环境无干扰时，对行或列采样的快拍进行FFT，记录雷达正常工作时，雷达静默区的接收噪声，并设定适当门限为k；正常工作时，如雷达静默区侦听到的信号高于门限，判定为系统受到干扰，假设干扰的个数为M个；步骤5：利用FFT对式2进行方向图形成，“列”单次采样快拍进行FFT运算，可以获得俯仰维指向干扰角度的FFT峰值点，记录其幅度是否大于接收系统噪声门限k；大于则记录干扰角度为：FFT后第m个谱峰对应的空间采样点为k
m
，n_FFT为FFT的点数；记录FFT后第m个谱峰的值步骤6：利用FFT对式3进行方向图形成，对“行”单次采样快拍进行FFT运算，可以获得方位指向干扰方向角度的FFT峰值点，记录其幅度是否大于接收系统噪声门限k；大于则记录干扰角度为：其中m＝0,1...MFFT后第m个谱峰对应的空间采样点...

【专利技术属性】
技术研发人员：崔向阳，付学斌，田欢，陈亮，田剑峰，高珊，
申请(专利权)人：西安电子工程研究所，
类型：发明
国别省市：