基于子阵的数字模拟合成自适应抗干扰方法技术

本发明专利技术公开了一种基于子阵的数字模拟合成自适应抗干扰方法，包括：1）将N元线性阵不重叠非均匀划分为L个子阵；2）对接收信号采样，进行空间谱估计，得到空间干扰的个数及其导向矢量；3）对空间干扰的个数进行检测与判断；4）若干扰个数大于子阵个数的一半时，采用LCMV算法，得到阵元级模拟权矢量，进而配置阵元级移相器和衰减器；5）若干扰个数小于或等于子阵个数一半时，利用Capon波束形成算法，得到子阵级数字权矢量，进而配置子阵级数字移相器和衰减器；6）合成波束输出。本发明专利技术提出的方法不仅测角精度高，抗干扰能力强，还可以减少硬件资源的消耗，提高计算的收敛速度。提高计算的收敛速度。提高计算的收敛速度。

【技术实现步骤摘要】
基于子阵的数字模拟合成自适应抗干扰方法


[0001]本专利技术涉及阵列信号处理
，特别是一种基于子阵的数字模拟合成自适应抗干扰方法。

技术介绍

[0002]阵列天线常被称为“相控阵”，即相位控制阵列的简称。顾名思义，相控阵天线是由许多辐射单元排列而成，各个单元的馈电相位由计算机灵活控制的阵列。阵列天线通常是由多个按一定规则排列的全向天线排列而成，通过改变阵列中各个天线单元的幅度和相位关系，可以实现对阵列天线方向图的波束指向的控制。
[0003]在电子战环境中，雷达系统不可避免地受到各种无源和有源的电磁干扰，这些干扰信号通过接收天线方向图的旁瓣或主瓣进入接收系统，使接收机的信号检测能力大为降低。雷达的抗干扰性能已成为衡量其性能的重要指标，阵列天线技术是提高雷达抗干扰能力的重要手段。阵列天线的自适应数字波束形成技术（ADBF）使雷达根据干扰特性，自动采取相应的抗干扰方案，并能随其变化动态地调整工作参数，达到某种最佳性能。目前阵列天线的ADBF技术已成为雷达信号处理研究的一个重要内容。
[0004]近些年提出的各种波束形成方法，其应用大多是基于阵元级的，这对于只有几个天线阵元的小型阵列来说是可行的。但现在的阵列天线都朝着大型的稀疏阵列的方向发展，特别在相控阵雷达中，几百乃至上千个天线阵元已是屡见不鲜。如果在这种大型阵列中，在阵元级上应用波束形成方法，就得对每一个天线阵元的接收信号进行单独处理，即每一个阵元组成一个接收通道，每一个接收通道都必须包含几次放大、混频，最后是视频处理或模数（A/D）转换。可想而知，硬件成本将会成倍增加。因此必须减少信号处理器的维数，同时尽可能保持对阵列响应的良好控制。
[0005]减少信号处理器的维数的一般方法是采用基于子阵的自适应阵列信号处理。
[0006]对于一个N元自适应阵列，有N个自由度可以利用，其中一部分自由度用来满足特定的约束，以利用某些可以得到的先验信息，称为约束自由度；剩下的自由度用来自适应地抑制干扰和噪声，称为自适应自由度。自适应自由度越多，自适应算法的计算量越大，收敛速度越慢；反之，收敛速度加快。基于子阵的自适应阵列信号处理正是减少了系统的自适应自由度来加快收敛速度的，因此缺点就是可以自适应抑制干扰的能力下降了。

技术实现思路

[0007]本专利技术的目的在于提供一种基于子阵的数字模拟合成自适应抗干扰方法。
[0008]实现本专利技术目的的技术解决方案为：一种基于子阵的数字模拟合成自适应抗干扰方法，包括以下步骤：
[0009]步骤一：建立N元线性阵列天线，采用不重叠非均匀划分将其划分为L个子阵，由阵元到子阵的降维矩阵设为T；线性阵列接收入射信号，入射信号由空间中P个互不相关的干扰信号和噪声组成；且N、L、P为正整数，；
[0010]步骤二：通过对子阵接收信号采样，在设定的快拍数下估计子阵级协方差矩阵，对子阵级协方差矩阵进行特征分解，得到分解后的干扰子空间和噪声子空间，并对该结果采用MUSIC算法进行空间谱估计，得到空间干扰的个数M和空间干扰导向矢量；
[0011]步骤三：对步骤二中得到的空间干扰的个数M进行判断，若，则执行步骤四；若，则执行步骤五；
[0012]步骤四：由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵C和约束响应向量f，采用LCMV波束形成算法，计算得到阵元级模拟权矢量，并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器，返回步骤二；
[0013]步骤五：在子阵级采用Capon波束形成算法计算得到子阵级数字权矢量，并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器；
[0014]步骤六：将步骤五中各子阵通道的信号合成为一路波束作为输出。
[0015]进一步地，步骤一中使用不重叠子阵划分降维矩阵T：
[0016]分别为子阵级，阵元级在k时刻接收的样本，且子阵级协方差矩阵与阵元级协方差矩阵的关系为：
[0017]。
[0018]进一步地，步骤二中采用如下方法计算子阵级协方差矩阵并得到干扰子空间和噪声子空间：
[0019]子阵级在k时刻接收的样本为为正整数，为快拍数即采样个数，对子阵级协方差矩阵的估计为：
[0020]；
[0021]对进行特征值分解得到：
[0022]；
[0023]其中为的特征值，为特征值对应的特征向量，其中：
[0024]，且
[0025]为干扰子空间对应的特征值，
ꢀ
为噪声子空间对应的特征值, 即为估计的干扰子空间；即为估计的噪声子空间。
[0026]进一步地，步骤二中采用如下方法得到空间干扰导向矢量：
[0027]由MUSIC算法得到空间谱函数
[0028]，其中为信号的导向矢量，为估计的噪声子空间。
[0029]由的峰值对应的角度，得到空间干扰导向矢量。
[0030]进一步地，步骤四中由期望信号导向矢量和空间干扰导向矢量构建约束矩阵C，约束响应向量f为：
[0031]；
[0032]；
[0033]其中 。
[0034]进一步地，步骤三中由LCMV算法求得阵元级模拟权矢量方法为：
[0035]；
[0036]求得。
[0037]进一步地，步骤五中在子阵级采用Capon波束形成算法得到子阵级数字权矢量的方法为：
[0038]将以下优化问题的解作为子阵级数字权矢量：
[0039]，其中为子阵级期望信号导向矢量，
[0040]求得。
[0041]进一步地，步骤六中得到的合成波束输出表达式如下：
[0042]。
[0043]本专利技术与现有方法相比，其显著优点是：（1）通过在子阵级对干扰信号个数进行循环检测与判断，并配置阵元级模拟权矢量，可以逐步减少输入到子阵级的干扰个数，最终可以抑制干扰的个数与全自适应阵列处理个数相当，从而相比常规的基于子阵的波束形成方法能够处理更多的干扰；（2）相比常规的全自适应阵列处理波束形成方法可以大大减少硬件资源的消耗，提高了计算收敛速度，适用于拥有成千上万个阵元天线的相控阵雷达。（3）该方法结合了在阵元级进行模拟波束形成和子阵级进行数字波束形成的优点，使得该系统相比普通的基于子阵的自适应波束形成方式抗干扰能力更强，测角精度更高。
附图说明
[0044]图1为本专利技术实施方式的信号处理流程图。
[0045]图2是在阵元级使用常规波束形成，在子阵级分别采用LCMV算法和Capon算法的天线方向图仿真结果。
[0046]图3为全阵列使用LCMV算法的天线方向图仿真结果。
[0047]图4为在阵元级分别使用LCMV算法和常规波束形成，在子阵级采用Capon算法的天线方向图仿真结果。
具体实施方式
[0048]下面结合附图和实施例对本专利技术的方法作进一步描述。
[0049]本专利技术的基本方案包括如下六个步骤，具体的信号处理流程如图1所示。
[本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于子阵的数字模拟合成自适应抗干扰方法，其特征在于，包括如下步骤：步骤一：建立N元线性阵列天线，采用不重叠非均匀划分将其划分为L个子阵，由阵元到子阵的降维矩阵设为T；线性阵列接收入射信号，入射信号由空间中P个互不相关的干扰信号和噪声组成；且N、L、P为正整数，；步骤二：通过对子阵接收信号采样，在设定的快拍数下估计子阵级协方差矩阵，对子阵级协方差矩阵进行特征分解，得到分解后的干扰子空间和噪声子空间，并对该结果采用MUSIC算法进行空间谱估计，得到空间干扰的个数M和空间干扰导向矢量；步骤三：对步骤二中得到的空间干扰的个数M进行判断，若，则执行步骤四；若，则执行步骤五；步骤四：由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵C和约束响应向量f，采用LCMV波束形成算法，计算得到阵元级模拟权矢量，并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器，返回步骤二；步骤五：在子阵级采用Capon波束形成算法计算得到子阵级数字权矢量，并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器；步骤六：将步骤五中各子阵通道的信号合成为一路波束作为输出。2.根据权利要求1所述的基于子阵的数字模拟合成自适应抗干扰方法，其特征在于，所述步骤一中的降维矩阵T满足：，分别为子阵级、阵元级在k时刻接收的样本，且子阵级协方差矩阵与阵元级协方差矩阵的关系为：，子阵级导向矢量与阵元级导向矢量的关系为：。3.根据权利要求1所述的基于子阵的数字模拟合成自适应抗干扰方法，其特征在于，所述步骤二中采用如下方法估计子阵级协方差矩阵并得到干扰...

【专利技术属性】
技术研发人员：蔡凌萍，邹阳，余其旺，李洪涛，田巳睿，邱林康，钱浩楠，邢灵尔，黄雪琴，狄儒霄，初瑞雪，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：