CN111948598B - 空域干扰信号检测方法与装置 - Google Patents
空域干扰信号检测方法与装置 Download PDFInfo
- Publication number
- CN111948598B CN111948598B CN202010747227.8A CN202010747227A CN111948598B CN 111948598 B CN111948598 B CN 111948598B CN 202010747227 A CN202010747227 A CN 202010747227A CN 111948598 B CN111948598 B CN 111948598B
- Authority
- CN
- China
- Prior art keywords
- interference
- value
- agc
- signal
- data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000001514 detection method Methods 0.000 claims abstract description 88
- 238000001228 spectrum Methods 0.000 claims abstract description 83
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000012545 processing Methods 0.000 claims abstract description 42
- 238000004364 calculation method Methods 0.000 claims abstract description 38
- 239000013598 vector Substances 0.000 claims abstract description 37
- 230000003595 spectral effect Effects 0.000 claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 23
- 238000007781 pre-processing Methods 0.000 claims abstract description 14
- 238000010276 construction Methods 0.000 claims description 17
- 239000011295 pitch Substances 0.000 claims description 16
- 238000009825 accumulation Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 238000003672 processing method Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims 1
- 230000001427 coherent effect Effects 0.000 claims 1
- 230000017105 transposition Effects 0.000 claims 1
- 238000013461 design Methods 0.000 description 13
- 230000000875 corresponding effect Effects 0.000 description 11
- 238000005259 measurement Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开的一种空域干扰信号检测方法与装置涉及阵列信号处理领域,本发明通过下述技术方案实现:采用FPGA芯片并行硬件处理架构,利用阵列采样信号数据形成数据协方差矩阵,基于MVDR算法获得空间谱值形成二维空间谱图。然后利用AGC值辅助空间谱搜索,实现干扰的来向检测。同时还通过数字AGC反检法实现对快速变化的干扰功率进行检测,得到干扰功率值,从而丰富了干扰信息的检测结果。本发明通过预处理方式获得导向矢量以及将谱峰计算转化为谱谷计算,进一步减少实现的复杂度,最后还利用空间谱噪底门限和局部最值联合判断来代替采样数据协方差矩阵的特征值求解来进行干扰个数的判定,获得干扰个数的检测,使得工程应用高效可行。
Description
技术领域
本发明涉及阵列信号处理领域,更具体的,涉及阵列信号处理领域中空域干扰信号 检测方法及其装置。
背景技术
通信设备在复杂电磁环境下,很容易受到包括单音、窄带、宽带等多种形式的有意或无意信号的干扰。如何有效的抑制干扰、改善信号接收质量,一直是雷达、通信、导航以及声纳等应用领域中特别引入关注的问题。当雷达处于强的有源干扰环境时,干扰信号从天 线旁瓣进入接收机的概率非常高,有时可淹没目标信号,造成雷达不能正常工作。为降低或 消除干扰信号对通信设备的影响,一种有效的方法是先对干扰信号进行检测获取干扰信息, 如干扰的来向、个数等,然后采取措施在物理上避免接收干扰方向上的信号或对接收到的信 号进行干扰抑制处理。
干扰信号的特征主要包含干扰的波达方向、干扰个数及干扰信号功率等。目前,工程上常采用阵列天线结合阵列信号处理算法来进行干扰信号检测。阵列信号处理主要的研究 内容就是:如何利用信号的空域特性来增强接收的有用信号,以及如何有效地提取包括信号 空域信息在内的信号的其它信息。其主要的研究领域可以分为波束形成技术、零点技术和空 间谱估计技术等几个方面,它们都是基于对信号进行空间采样的数据进行处理,因此,这些 技术是相互渗透和相互关联的。由于处理的目的不同,其着眼点有所差别,因此导致有不同 的算法。众所周知,空间谱估计测向的最终实现,是由测向处理软件完成的,而软件***的 核心,则是空间谱估计算法。这些算法基本上可以归为2大类,即基于经典估计理论的极大 似然估计和归结于信号参量估计的参数谱方法。将极大似然参数估计方法应用于超分辨估计, 由于这个估计似然函数是非线性的,求解其最优解需进行多维搜索,运算量巨大。
经过多年的发展,已经产生了大量性能优异的测向算法可资利用,典型的有MUSIC、ESPRIT、子空间拟合、多维MUSIC等。其中,MUSIC算法最为经典最 具代表性,由于具有超瑞利限特性,所以又被称作超分辨测向方法或高分辨测向算法,但要 求干扰个数准确已知,且在硬件实现上复杂度高。基于波束扫描的波达方向估计是最常用的 一种方法,该方法受瑞利限的限制,对干扰信号的空域分辨率低,检测精度差。另外,上述 算法均无法准确估计干扰信号个数以及干扰信号的功率,也限制了上述方法在干扰信号检测 中的应用。空间谱估计测向方法需经过复杂的计算才能得到待测信号的到达方向。
现有技术空域干扰检测方法干扰来向估计分辨率差、干扰个数估计不准确,尤其是 基于DSP芯片的空间谱构建、谱峰搜索的检测处理架构,无法满足大规模阵列空间谱构建 的复杂度要求及谱峰方位和俯仰二维域搜索实时性要求。
发明内容
本发明所要解决的技术问题是针对现有技术空域干扰检测方法干扰来向估计分辨率 差、干扰个数估计不准确,尤其是基于DSP芯片的空间谱构建、谱峰搜索的检测处理架构, 无法满足大规模阵列空间谱构建的复杂度要求及谱峰方位和俯仰二维域搜索实时性要求。本 发明提供解决上述问题的空域干扰信号检测方法。
本发明通过下述技术方案实现:一种空域干扰信号检测方法,其特征在于包括如下 步骤:
S1:干扰功率检测:采用自适应门限结合局部最值搜索的方法,对干扰的波达方向、干扰个 数进行检测,利用时域处理方法获得数字自动增益控制AGC值,通过数字AGC反检法对 干扰功率的检测,根据M为阵元个数,得到所有阵元输出干扰功率值power(1)、power(2)…power(M),然后对M个阵元的干扰功率值进行比较,选最大的那个干扰功率值作为干扰功率的输出值;即时测量干扰信号功率,实时对干扰功率进行精确估计;
S2:空间谱构建:将M个阵元的预处理后信号转化成空域信号,得到第M个天线阵元的信 号形成的空域信号,利用阵列采样信号数据形成数据协方差矩阵,进行最小方差无失真响应 MVDR算法空间谱值的计算;空间谱噪底门限和局部最值联合判断来代替采样数据协方差 矩阵的特征值求解来进行干扰个数的判定。
S3:干扰来向,个数检测:采用FPGA芯片并行硬件处理架构,对快速变化的干扰来向检测,采用预处理的方式计算导向矢量,通过方位值θ和俯仰值来计算RAM的地址,作为导向矢量RAM的地址进行查找,得到对应θ和的导向矢量基于步骤S1获得的信号数字AGC信息辅助,首先进行干扰有无判定,若无干扰 则不进行谱峰搜索,若有干扰进行谱峰搜索,以及利用空间谱噪底门限和局部最值联合判断 来代替采样数据协方差矩阵的特征值求解来进行干扰个数的判定,从而获得干扰来向及个数的检测。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明方法基于以上设计,解决了现有空域干扰检测方法干扰来向估计分辨率差、干 扰个数估计不准确、工程实现难及工程实时性等问题,干扰来向检测结果分辨率高,实时性 高,具有工程应用价值。设计中采用自适应门限结合局部最值搜索的方法实现对干扰的波达 方向、干扰个数的检测,避免传统方法对采样数据协方差矩阵特征值的求解,降低了计算复 杂度,可工程化应用;其次通过数字AGC反检法实时对干扰功率进行精确估计,实现对干 扰信号功率即时测量;最后采用FPGA芯片的硬件处理架构,利用芯片并行处理能力实现对 快速变化的干扰来向检测、空间谱构建及谱峰搜索。该发明具有复杂度低、稳定性高、测量 精度高等优点,满足工程适应性要求。通过以上算法优化及设计,使得该方法具有快速实现 对干扰信息的检测能力,整体优化了程序设计的复杂度,使得工程应用高效合理。
本发明采用自适应门限结合局部最值搜索的方法对干扰的波达方向、干扰个数进行检 测,利用时域处理方法获得数字自动增益控制AGC值,即时测量干扰信号功率,通过数字 AGC反检法对干扰功率的检测,得到所有阵元输出干扰功率值power(1)、 power(2)…power(M),并对M个阵元的干扰功率值进行比较,选择最大的那个干扰功率值作 为干扰功率的输出值,对干扰功率进行实时精确估计;利用AGC数值信息是否进行空间谱 搜索启动,避免了无干扰时错误判定,提升了工程的稳定性和可靠性。有效减小了运算量, 并且可根据精度设计需要调节预计算值存储量,减小了硬件实现的复杂度,工程实现简单。 本发明采用基于空间谱噪底门限和局部最值联合判决方法,相对比基于采样数据协方差矩阵 的特征值求解进行干扰个数的判定的方法,突破了工程实现的困难,简单可靠,突破了工程实现中矩阵特征值求解复杂度高的瓶颈。
其次通过数字AGC反检法实时对干扰功率进行精确估计,实现对干扰信号功率即时测量; 干扰功率检测通过数字AGC系数查表得到,克服直接对信号进行实时FFT分析时波动较大 的问题,并且较频域FFT变换来说其复杂度更低、稳定性和精度均更高,对硬件资源需求 更小。
本发明将M个阵元的预处理后信号转化成空域信号,得到第M个天线阵元的信号形成的空域信号,利用阵列采样信号数据形成数据协方差矩阵,基于最小方差无失真响应MVDR算法进行空间谱值的计算。将谱峰搜索转化为谱谷搜索,避免除法处理硬件实现, 在不影响算法性能的同时简化了计算复杂度。
本发明采用FPGA芯片并行硬件处理架构,对快速变化的干扰来向检测,基于步骤S1 获得的信号数字AGC信息辅助,进行干扰有无判定,若无干扰则不进行谱谷搜索,若有干扰进行谱谷搜索,将谱峰搜索转化为谱谷搜索,在不影响算法性能的同时简化了计算复杂度。 采用局部最小值代替全局极小值方法对空间谱搜索,加快了空间谱搜索处理速度,提高工程 实时性。避免了二维谱谷极值搜索进行主对角线和副对角线的比较,严重影响搜索速度的缺 陷。采用自适应门限结合局部最值搜索的方法实现对干扰的波达方向、干扰个数的检测,避 免传统方法对采样数据协方差矩阵特征值的求解,降低了计算复杂度,可工程化应用;
本发明采用基于空间谱噪底门限和局部最值联合判决方法,相对比基于采样数据协 方差矩阵的特征值求解进行干扰个数的判定的方法,突破了工程实现的困难,简单可靠,突 破了工程实现中矩阵特征值求解复杂度高的瓶颈。最后采用FPGA芯片的硬件处理架构,利 用芯片并行处理能力实现对快速变化的干扰来向检测、空间谱构建及谱峰搜索。该发明具有 复杂度低、稳定性高、测量精度高等优点,满足工程适应性要求。
本发明采用自适应门限结合局部最值搜索的方法实现对干扰的波达方向、干扰个数 的检测,解决了空域干扰信号检测算法干扰来向估计分辨率差、干扰个数估计不准等难点。 避免了传统方法对采样数据协方差矩阵特征值的求解,降低了实现的复杂度。设计中还采用 数字AGC信息辅助方法值来估计是否有干扰,避免无干扰时进行空间谱搜索提升干扰检测 结果的稳定性。以及利用AGC信息反检获得干扰信号功率,避免频域处理,减少复杂度。 本发明采用FPGA芯片并行硬件处理架构,对快速变化的干扰来向检测,基于步骤S1获得 的信号数字AGC信息辅助,利用FPGA芯片硬件处理架构的并行处理能力实现对快速变化的干扰来向检测、空间谱构建及谱峰搜索。利用FPGA芯片的特点,硬件实现中导向矢量的计算通过预处理的方式进行预计算,并存储于RAM中,通过方位值θ和俯仰值来计算导 向矢量RAM的地址,然后进行查找得到对应导向矢量,有效减小了运算量,并且可根据设 计需要来调节所需要的精度进行预计算值存储,减小了硬件实现的复杂度。使得该方法在基 于FPGA芯片硬件处理架构上具有具有复杂度低、稳定性高、测量精度高等优点,使得工程 应用高效合理。利用AGC数值信息是否进行空间谱搜索启动,避免了无干扰时错误判定, 提升了工程的稳定性和可靠性。
发明方法基于以上设计,解决了现有空域干扰检测方法干扰来向估计分辨率差、干扰 个数估计不准确、工程实现难及工程实时性等问题,干扰来向检测结果分辨率高,实时性高, 具有工程应用价值。设计中采用自适应门限结合局部最值搜索的方法实现对干扰的波达方向、 干扰个数的检测,避免传统方法对采样数据协方差矩阵特征值的求解,降低了计算复杂度, 可工程化应用;其次通过数字AGC反检法实时对干扰功率进行精确估计,实现对干扰信号 功率即时测量;最后采用FPGA芯片的硬件处理架构,利用芯片并行处理能力实现对快速变 化的干扰来向检测、空间谱构建及谱峰搜索。该发明具有复杂度低、稳定性高、测量精度高 等优点,满足工程适应性要求。通过以上算法优化及设计,使得该方法具有快速实现对干扰 信息的检测能力,整体优化了程序设计的复杂度,使得工程应用高效合理。
本发明解决了空域干扰信号检测算法干扰来向估计分辨率差、干扰个数估计不准等 难点。采用自适应门限结合局部最值搜索的方法实现对干扰的波达方向、干扰个数的检测, 避免传统方法对采样数据协方差矩阵特征值的求解,降低了实现的复杂度。设计中还采用数 字AGC信息辅助方法值来估计是否有干扰,避免无干扰时进行空间谱搜索提升干扰检测结 果的稳定性。以及利用AGC信息反检获得干扰信号功率,避免频域处理,减少复杂度。
本发明采用FPGA芯片的硬件处理架构,利用芯片并行处理能力实现对快速变化的干扰来向检测、空间谱构建及谱峰搜索。利用FPGA芯片的特点,硬件实现中导向矢量的计算通过预处理的方式进行预计算,并存储于RAM中,通过方位值和俯仰值来计算导向矢量RAM的地址,然后进行查找得到对应导向矢量,有效减小了运算量,并且可根据设计需要 来调节所需要的精度进行预计算值存储,减小了硬件实现的复杂度。使得该方法在基于 FPGA芯片硬件处理架构上具有具有复杂度低、稳定性高、测量精度高等优点,使得工程应 用高效合理。
本发明基于硬件实现的一种AGC信息辅助的空域干扰信号检测装置,解决了空域干 扰信号检测算法干扰来向估计分辨率差、干扰个数估计不准等难点。采用自适应门限结合局 部最值搜索的方法实现对干扰的波达方向、干扰个数的检测,避免传统方法对采样数据协方 差矩阵特征值的求解,降低了实现的复杂度。设计中还采用数字AGC信息辅助方法值来估 计是否有干扰,避免无干扰时进行空间谱搜索提升干扰检测结果的稳定性。以及利用AGC 信息反检获得干扰信号功率,避免频域处理,减少复杂度。
本发明还采用FPGA芯片的硬件处理架构,利用芯片并行处理能力实现对快速变化的干扰来向检测、空间谱构建及谱峰搜索。利用FPGA芯片的特点,硬件实现中导向矢量的计算通过预处理的方式进行预计算,并存储于RAM中,通过方位值和俯仰值来计算导向矢量RAM的地址,然后进行查找得到对应导向矢量,有效减小了运算量,并且可根据设计需 要来调节所需要的精度进行预计算值存储,减小了硬件实现的复杂度。使得该方法在基于FPGA芯片硬件处理架构上具有具有复杂度低、稳定性高、测量精度高等优点,使得工程应用高效合理。
本发明可以应用于同类或者其他产品上,不能规避此技术手段。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分, 并不构成对本发明实施例的限定。在附图中:
图1是本发明干扰功率检测流程图。
图2为本发明预处理流程图。
图3为本发明空间谱构建的实现框图。
图4为本发明中步骤S3空间谱搜索的详细步骤流程图。
图5本发明支持空域干扰信号检测方法的空域干扰信号检测装置电路原理示意图。 为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进 一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明 的限定。
具体实施方式
参阅图1。根据本发明,所述的空域干扰信号检测方法,包括:
干扰功率检测:采用自适应门限结合局部最值搜索的方法,对干扰的波达方向、干扰个数进 行检测,利用时域处理方法获得数字自动增益控制AGC值,通过数字AGC反检法对干扰 功率的检测,得到所有阵元输出干扰功率值power(1)、power(2)…power(M),然后对M个阵 元的干扰功率值进行比较,选最大的那个干扰功率值作为干扰功率的输出值;即时测量干扰 信号功率,实时对干扰功率进行精确估计;
S2:空间谱构建:将M个阵元的预处理后信号转化成空域信号,得到第M个天线阵元的信 号形成的空域信号,利用阵列采样信号数据形成数据协方差矩阵,进行最小方差无失真响应 MVDR算法空间谱值的计算;空间谱噪底门限和局部最值联合判断来代替采样数据协方差 矩阵的特征值求解来进行干扰个数的判定。
S3:干扰来向,个数检测:采用FPGA芯片并行硬件处理架构,对快速变化的干扰 来向检测,采用预处理的方式计算导向矢量,通过方位值θ和俯仰值来计算RAM的地址,作为导向矢量RAM的地址进行查找,得到对应θ和的导向矢量基于步骤S1获得的信号数字AGC信息辅助,首先进行干扰有无判定,若无干扰 则不进行谱峰搜索,若有干扰进行谱峰搜索,以及利用空间谱噪底门限和局部最值联合判断 来代替采样数据协方差矩阵的特征值求解来进行干扰个数的判定,从而获得干扰来向及个数的检测。
具体实施时,图1为基于硬件实现的高性能抗干扰方法的流程图,该图对整个处理流程进行了阐述,本发明首先通过AD对接收阵元的中频信号进行采样,得到每个阵元的数字AGC值和数据预处理后IQ信号;其次根据得到数字AGC值进行干扰功率检测并输出, 同时进行采样数据协方差矩阵RXX的计算。然后根据RXX的值,利用MVDR算法计算每个空 间来向上的谱值,构建出空间谱搜索的二维空间谱图。最后根据数字AGC值判定是否进行 二维空间谱图的谱搜索,无干扰则退出搜索,否则按照谱搜索的方法进行空间谱搜索输出干 扰来向及个数。
进一步地,参阅图2。在干扰功率检测步骤中,首先对接收到的阵元信号进行数字AGC计算,其次根据AGC值进行反检,得到干扰功率值,最后输出所有阵元中功率最强的 值作为干扰功率的检测结果。1-M个阵元接收信号经模数AD采样,得到时域中频采样数据data,将采样数据data和初始AGC系数值agc_value相乘,得到相乘后的信号值,相乘后的信号值取平方运算进行累加得到累加后的值sum_data,这里的累加点为N,然后将sum_data值除以累加点数N,即对数据求和sumdata/N,再通过数字AGC得到1-M个阵元的AGC值, 将得到的值mean_data与目标值tg_data进行比较,如果值比tg_data大则AGC系数偏大,需 要往小调整,否则往大调整,直到信号相干累加后的mean_data值在稳定的目标值tgdata上下 小范围内波动。经过AGC系数值agc_value和信号data相乘后功率保持一个相对稳定值,可 以利用AGC值进行查表得到相应的干扰功率值,不断重复上述过程,得到所有阵元输出干 扰功率值power(1)、power(2)…power(M),M为阵元个数,然后对M个阵元的干扰功率值 进行比较,选最大的那个作为干扰功率检测的干扰功率输出值。
进一步地,参阅图3。在干扰来向检测中:将模数AD采样的数据送入数据预处理模块,预处理后得到带有幅度、相位信息的IQ两路信号,IQ两路信号送入正交变换滤波器, 对IQ两路信号进行幅度和相位数据校正,得到校正后的IQ两路信号,然后,将M个阵元 的预处理后信号转化成空域信号矩阵,得到第M个天线阵元的信号形成的空域信号 X(n)=[x1(n)x2(n) ... xM(n)]H,其中,xM(n)表示第M个天线在时刻n的采样 数据。
利用空域信号x(n)形成采样数据协方差矩阵Rxx,采样数据协方差矩阵Rxx的计算公为:
在导向矢量计算中,方位角θ和俯仰角按照1度步进进行空间遍历,预先计算所有 方位和俯仰组合得到的导向矢量设空间谱搜索的范围为θ∈[0,90],共 有32760=91*360个值,将计算结果按a(0,0)、a(1,0)...a(359,0)、a(0,1)....a(359,90)顺序存入 FPGA的RAM中,然后根据作为导向矢量RAM的地址进行查找,得 到对应θ和的导向矢量构建空间谱计算公式为:
为了减小计算的复杂度,空间谱函数计算公式简化为:
pij为公式(4)计算得到的值,然后按列,即按照 p1,1,...,pS1,1,p1,2,…,pS1,2,...,p1,S2,…,pS1,S2的顺序写入RAM,也就是说最先写入RAM的是 PMVDR矩阵的第1列,即其中,S1=91,S2=360。因此RAM的地 址addr_ram与空间谱值pij是一一对应的,且RAM的地址addr_ram与pij的下标i,j的对应关 系为:
j=mod(addr_ram,S1)+1(8)
其中表示取整,mod取余数,而pij的下标i,j的值就是来波信号的方位角和俯仰角, 因此通过搜索最值pij对应的RAM地址addr_ram就可以获得来波信号的方位角sig_cokf和俯仰 角sig_cokt。
进一步地,在二维谱谷搜索过程中,搜索二维谱矩阵PMVDR中的最值,从空间谱函数计算公式所示的矩阵PMVDR中寻找二维极小值,极小值对应的列 和行,即为来波信号的方位角sig_cokf和俯仰角sig_cokt,如果pij为一个二维极小值点,则 pij应同时小于等于pi-1,j-1、pi-1,j、pi-1,j+1、pi,j-1、pi,j+1、pi+1,j-1、pi+1,j和pi+1,j+1。如用数学表达式表达谱谷的位置,即{(i,j)|pij≤pi1,j1,i1∈[i-1,i+1],j1∈[j-1,j+1],(i1,j1)≠(i,j)}
(6)也就是说,二维谱谷搜索,除了行和列的比较之外还要进行主对角线和副对角线 的比较。如果完全按照式(6)进行空间谱搜索,严重影响搜索速度。
进一步地,参阅图5,对空间谱搜索采用局部最小值代替极小值,具体步骤包括:
步骤A:将当前时刻输入AGC值赋值给d_amp,及只有噪声输入AGC值赋值给 d_amp0,比较d_amp值和d_amp0,确定当前是否有信号输入,如果没有信号则直接结 束搜索,否则转步骤B;
步骤B:取二维谱矩阵PMVDR的列求均值得到噪底水平;然后取矩阵第j~j+9列中元素的 最小值,若该值在第j+1~j+8列,且小于噪底水平4倍以上,则记下该值所在的RAM的地址addr_ram和值的大小,其中j=1,9,17,25,…353,当j=353时,取第353~360,1,2列。遍历完所有j的取值,转步骤C;
步骤C:若最后满足要求的局部最小值大于M-1个,M为阵元个数,则选择所有最值中最小的M-1个的所对应RAM的地址作为最后结果。
参阅图5。另一方面,本发明还提供了一种AGC信息辅助的空域干扰信号检测装置,该装置支持所述的一种AGC信息辅助的空域干扰信号检测方法,该装置包括:连接阵列信号接收单元的干扰来向检测单元和AGC信息辅助的干扰功率检测单元,其中:所述干扰来向检测单元包括顺次串联的预处理模块、协方差矩阵生成模块、空间谱峰构建模块和空间谱 峰搜索模块,干扰功率检测单元包括连接空间谱峰搜索模块及干扰信息输出模块的AGC模 块和干扰功率计算模块。
阵列信号接收单元将接收到的阵列天线信号送入模数转换器进行模数AD采样处理, 得到数字中频信号;干扰功率检测单元采用时域数字AGC反检方法,对干扰功率的检测, 将AD采样数字中频信号送入AGC模块,计算输入数字中频信号的AGC值,干扰功率计 算模块根据AGC模块获得的AGC值,进行干扰功率的检测、空间谱搜索模块根据空间谱 构建模块获得空间谱值及及所述AGC模块给出的AGC值进行联合处理,检测出干扰的来 向和个数,干扰信息输出模块根据所述干扰功率检测单元和空间谱峰搜索模块给出的干扰功 率、干扰来向及个数,输出干扰信息检测结果。
预处理模块将中频信号进行正交处理得到IQ信号,并进行幅度和相位校正,得到校 正后的IQ信号,协方差矩阵生成模块利用所有阵元的采样数据进行自相关运算,获得采样 数据协方差矩阵,输出协方差矩阵值;空间谱峰构建模块采用MVDR空间谱估计方法,根据所述协方差矩阵生成模块获得采样数据协方差矩阵Rxx及导向矢量来计算空间谱值, 其计算公式为式中为的共轭转置,为Rxx的逆。
进一步的谱峰构建模块还包括导向矢量计算模块和协方差矩阵的求逆模块,导向矢量 计算模块是为了获得来波信号所在方向的约束矢量,协方差矩阵求逆是空间谱计算所需的关 键值。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详 细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的 保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含 在本发明的保护范围之内。
Claims (7)
1.一种空域干扰信号检测方法,其特征在于包括如下步骤:
S1:干扰功率检测:采用自适应门限结合局部最值搜索的方法,对干扰的波达方向、干扰个数进行检测,利用时域处理方法获得数字自动增益控制AGC值,通过数字AGC反检法对干扰功率进行检测,首先对接收到的阵元信号进行数字AGC计算,根据1-M个阵元接收信号经模数AD采样,得到时域中频采样数据data,将采样数据data和初始AGC系数值agc_value相乘,得到相乘后的信号值,相乘后的信号值取平方运算进行累加得到累加后的值sum_data,然后以N为累加点,将sum_data值除以累加点数N,对数据求和sumdata/N,再通过数字AGC得到1到M个阵元的AGC值,将得到的值mean_data与目标值tg_data进行比较,如果值比tg_data大则AGC系数偏大,需要往小调整,否则往大调整,直到信号相干累加后的mean_data值在稳定的目标值tgdata上下小范围内波动;经过AGC系数值agc_value和信号data相乘后功率保持一个相对稳定值,利用AGC值进行查表,即根据AGC值进行反检,得到相应的干扰功率值,不断重复上述过程,得到所有阵元输出干扰功率值power(1)、power(2)···power(M),M为阵元个数,然后对M个阵元的干扰功率值进行比较,选最大的那个作为干扰功率检测的干扰功率输出值;即时测量干扰信号功率,实时对干扰功率进行精确估计,输出所有阵元中功率最强的值作为干扰功率的检测结果;
S2:空间谱构建:将M个阵元预处理后的信号转化成空域信号,得到空域信号矩阵,利用阵列采样信号数据形成数据协方差矩阵,进行最小方差无失真响应MVDR算法空间谱值的计算;再利用空间谱噪底门限和局部最值联合判断来代替采样数据协方差矩阵的特征值求解过程,进行干扰个数判定;
2.如权利要求1所述的空域干扰信号检测方法,其特征在于:在干扰来向检测中,模数AD采样的数据送入数据预处理模块,预处理后得到带有幅度、相位信息的IQ两路信号,IQ两路信号送入正交变换滤波器,对IQ两路信号进行幅度和相位数据校正,得到校正后的IQ两路信号校正参数,这里的校正参数是硬件通道间的相位和幅度差,然后,将M个阵元的预处理后信号转化成空域信号矩阵,得到第M个天线阵元的信号形成的空域信号X(n)=[x1(n)x2(n) ... xM(n)]H,其中,xM(n)表示第M个天线在时刻n的采样数据,H表示转置。
7.一种使用权利要求1所述空域干扰信号检测方法的装置,包括:连接阵列信号接收单元的干扰来向检测单元和AGC信息辅助的干扰功率检测单元,其中:所述干扰来向检测单元包括顺次串联的预处理模块、协方差矩阵生成模块、空间谱峰构建模块和空间谱峰搜索模块,干扰功率检测单元包括连接空间谱峰搜索模块及干扰信息输出模块的AGC模块和干扰功率计算模块,其特征在于:阵列信号接收单元将接收到的阵列天线信号送入模数转换器进行模数AD采样处理,得到数字中频信号;干扰功率检测单元采用时域数字AGC反检方法,对干扰功率的检测,将AD采样数字中频信号送入AGC模块,计算输入数字中频信号的AGC值,干扰功率计算模块根据AGC模块获得的AGC值,进行干扰功率的检测、空间谱搜索模块根据空间谱构建模块获得空间谱值及所述AGC模块给出的AGC值进行联合处理,检测出干扰的来向和个数,干扰信息输出模块根据所述干扰功率检测单元和空间谱峰搜索模块给出的干扰功率、干扰来向及个数,输出干扰信息检测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010747227.8A CN111948598B (zh) | 2020-07-30 | 2020-07-30 | 空域干扰信号检测方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010747227.8A CN111948598B (zh) | 2020-07-30 | 2020-07-30 | 空域干扰信号检测方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111948598A CN111948598A (zh) | 2020-11-17 |
CN111948598B true CN111948598B (zh) | 2022-11-29 |
Family
ID=73339782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010747227.8A Active CN111948598B (zh) | 2020-07-30 | 2020-07-30 | 空域干扰信号检测方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111948598B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112769469B (zh) * | 2021-01-23 | 2023-02-24 | 成都振芯科技股份有限公司 | 一种基于波束形成控制运算阵元数的方法及装置 |
CN113156450B (zh) * | 2021-03-17 | 2024-02-27 | 北京理工大学 | 无人机上的主动旋转激光雷达***及其控制方法 |
CN113156392B (zh) * | 2021-05-26 | 2022-04-08 | 哈尔滨工业大学 | 一种基于俯仰域自适应处理的杂波抑制方法 |
CN113534198B (zh) * | 2021-06-16 | 2023-05-23 | 北京遥感设备研究所 | 基于协方差矩阵重构的卫星导航动态抗干扰方法及其*** |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403791A (zh) * | 2008-11-06 | 2009-04-08 | 哈尔滨工程大学 | 快速实时空间谱估计超分辨测向装置及方法 |
CN101990216A (zh) * | 2009-08-04 | 2011-03-23 | 鼎桥通信技术有限公司 | 一种对td-scdma***的小区外干扰源检测方法 |
CN103427871A (zh) * | 2013-07-22 | 2013-12-04 | 西安电子科技大学 | 一种抑制SmartAGC抗干扰能力的转发式干扰方法 |
CN104536014A (zh) * | 2015-01-23 | 2015-04-22 | 哈尔滨工业大学 | 一种基于Relax算法的抗转发式欺骗干扰方法 |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
CN106970374A (zh) * | 2017-05-27 | 2017-07-21 | 中国人民解放军63892部队 | 干扰条件下单脉冲雷达目标和干扰源角度的测量方法 |
CN106970361A (zh) * | 2017-05-27 | 2017-07-21 | 中国人民解放军63892部队 | 干扰条件下双极化单脉冲雷达目标角度的估计方法 |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN109799516A (zh) * | 2018-12-22 | 2019-05-24 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种用于gnss干扰检测定位的装置 |
CN109856647A (zh) * | 2019-02-15 | 2019-06-07 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种空地协同的gnss干扰检测定位***及方法 |
CN110879380A (zh) * | 2019-08-19 | 2020-03-13 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于手机的gnss干扰源定位方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2495128B (en) * | 2011-09-30 | 2018-04-04 | Skype | Processing signals |
GB2495129B (en) * | 2011-09-30 | 2017-07-19 | Skype | Processing signals |
US20130332156A1 (en) * | 2012-06-11 | 2013-12-12 | Apple Inc. | Sensor Fusion to Improve Speech/Audio Processing in a Mobile Device |
US10845463B2 (en) * | 2015-07-17 | 2020-11-24 | Origin Wireless, Inc. | Method, apparatus, and system for wireless object scanning |
-
2020
- 2020-07-30 CN CN202010747227.8A patent/CN111948598B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403791A (zh) * | 2008-11-06 | 2009-04-08 | 哈尔滨工程大学 | 快速实时空间谱估计超分辨测向装置及方法 |
CN101990216A (zh) * | 2009-08-04 | 2011-03-23 | 鼎桥通信技术有限公司 | 一种对td-scdma***的小区外干扰源检测方法 |
CN103427871A (zh) * | 2013-07-22 | 2013-12-04 | 西安电子科技大学 | 一种抑制SmartAGC抗干扰能力的转发式干扰方法 |
CN104536014A (zh) * | 2015-01-23 | 2015-04-22 | 哈尔滨工业大学 | 一种基于Relax算法的抗转发式欺骗干扰方法 |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
CN106970374A (zh) * | 2017-05-27 | 2017-07-21 | 中国人民解放军63892部队 | 干扰条件下单脉冲雷达目标和干扰源角度的测量方法 |
CN106970361A (zh) * | 2017-05-27 | 2017-07-21 | 中国人民解放军63892部队 | 干扰条件下双极化单脉冲雷达目标角度的估计方法 |
CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
CN109799516A (zh) * | 2018-12-22 | 2019-05-24 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种用于gnss干扰检测定位的装置 |
CN109856647A (zh) * | 2019-02-15 | 2019-06-07 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种空地协同的gnss干扰检测定位***及方法 |
CN110879380A (zh) * | 2019-08-19 | 2020-03-13 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于手机的gnss干扰源定位方法 |
Non-Patent Citations (1)
Title |
---|
"GNSS干扰检测与识别技术研究";朱鹏程;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20190415;第I136-422页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111948598A (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111948598B (zh) | 空域干扰信号检测方法与装置 | |
CN109459744B (zh) | 一种实现多干扰抑制的稳健自适应波束形成方法 | |
RU2382379C2 (ru) | Способ многосигнальной пеленгации источников радиоизлучения на одной частоте | |
KR102001394B1 (ko) | 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법과 이를 위한 장치 | |
CN109991519B (zh) | 基于神经网络和无线传感阵列的局部放电测向方法及*** | |
CN112612010A (zh) | 一种基于波瓣***预处理的米波雷达低仰角测高方法 | |
CN111628790A (zh) | 一种基于干扰带宽检测的高精度抗干扰方法及装置 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN115436896A (zh) | 快速的雷达单快拍music测角方法 | |
CN111896913A (zh) | 高频雷达单极子/交叉环天线通道增益校准方法及装置 | |
CN114137476A (zh) | 一种基于极化修正的短波测向方法 | |
CN111665469B (zh) | 一种基于空间时频分布的水下多径信号参数估计方法 | |
Okane et al. | Resolution improvement of wideband direction-of-arrival estimation" Squared-TOPS" | |
CN103901421A (zh) | 基于对角减载的水声阵列smi-mvdr空间谱估计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN111722178B (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
CN115248413A (zh) | 一种适用于非均匀线阵的离格信号波达方向估计方法 | |
CN109861770B (zh) | 一种基于波束形成输出功率组合的宽带信号检测方法 | |
CN113589223A (zh) | 基于互耦情况下嵌套阵列的测向方法 | |
JP4119719B2 (ja) | 移動局方向推定方法及び装置 | |
CN114690114A (zh) | 一种信号测向方法、装置、***及存储介质 | |
CN111323750A (zh) | 一种基于声矢量阵列网络的直接定位方法 | |
CN112649787B (zh) | 一种基于低频圆环阵的目标方位估计方法 | |
JP4391771B2 (ja) | 測角装置 |
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 |