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

1.本发明涉及阵列信号处理技术领域，特别是一种基于子阵的数字模拟合成自适应抗干扰方法。


背景技术：

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


技术实现要素：

7.本发明的目的在于提供一种基于子阵的数字模拟合成自适应抗干扰方法。
8.实现本发明目的的技术解决方案为：一种基于子阵的数字模拟合成自适应抗干扰方法，包括以下步骤：
9.步骤一：建立n元线性阵列天线，采用不重叠非均匀划分将其划分为l个子阵，由阵元到子阵的降维矩阵设为t；线性阵列接收入射信号，入射信号由空间中p个互不相关的干扰信号和噪声组成；且n、l、p为正整数，；
10.步骤二：通过对子阵接收信号采样，在设定的快拍数下估计子阵级协方差矩阵，对子阵级协方差矩阵进行特征分解，得到分解后的干扰子空间和噪声子空间，并对该结果采用music算法进行空间谱估计，得到空间干扰的个数m和空间干扰导向矢量；
11.步骤三：对步骤二中得到的空间干扰的个数m进行判断，若，则执行步骤四；若，则执行步骤五；
12.步骤四：由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵c和约束响应向量f，采用lcmv波束形成算法，计算得到阵元级模拟权矢量，并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器，返回步骤二；
13.步骤五：在子阵级采用capon波束形成算法计算得到子阵级数字权矢量，并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器；
14.步骤六：将步骤五中各子阵通道的信号合成为一路波束作为输出。
15.进一步地，步骤一中使用不重叠子阵划分降维矩阵t：
16.分别为子阵级，阵元级在k时刻接收的样本，且子阵级协方差矩阵与阵元级协方差矩阵的关系为：
17.。
18.进一步地，步骤二中采用如下方法计算子阵级协方差矩阵并得到干扰子空间和噪声子空间：
19.子阵级在k时刻接收的样本为为正整数，为快拍数即采样个数，对子阵级协方差矩阵的估计为：
20.；
21.对进行特征值分解得到：
22.；
23.其中为的特征值，为特征值对应的特征向量，其中：
24.，且
25.为干扰子空间对应的特征值，
ꢀ
为噪声子空间对应的特征值, 即为估计的干扰子空间；即为估计的噪声子空间。
26.进一步地，步骤二中采用如下方法得到空间干扰导向矢量：
27.由music算法得到空间谱函数
28.，其中为信号的导向矢量，为估计的噪声子空间。
29.由的峰值对应的角度，得到空间干扰导向矢量。
30.进一步地，步骤四中由期望信号导向矢量和空间干扰导向矢量构建约束矩阵c，约束响应向量f为：
31.；
32.；
33.其中 。
34.进一步地，步骤三中由lcmv算法求得阵元级模拟权矢量方法为：
35.；
36.求得。
37.进一步地，步骤五中在子阵级采用capon波束形成算法得到子阵级数字权矢量的方法为：
38.将以下优化问题的解作为子阵级数字权矢量：
39.，其中为子阵级期望信号导向矢量，
40.求得。
41.进一步地，步骤六中得到的合成波束输出表达式如下：
42.。
43.本发明与现有方法相比，其显著优点是：（1）通过在子阵级对干扰信号个数进行循环检测与判断，并配置阵元级模拟权矢量，可以逐步减少输入到子阵级的干扰个数，最终可以抑制干扰的个数与全自适应阵列处理个数相当，从而相比常规的基于子阵的波束形成方法能够处理更多的干扰；（2）相比常规的全自适应阵列处理波束形成方法可以大大减少硬件资源的消耗，提高了计算收敛速度，适用于拥有成千上万个阵元天线的相控阵雷达。（3）该方法结合了在阵元级进行模拟波束形成和子阵级进行数字波束形成的优点，使得该系统相比普通的基于子阵的自适应波束形成方式抗干扰能力更强，测角精度更高。
附图说明
44.图1为本发明实施方式的信号处理流程图。
45.图2是在阵元级使用常规波束形成，在子阵级分别采用lcmv算法和capon算法的天线方向图仿真结果。
46.图3为全阵列使用lcmv算法的天线方向图仿真结果。
47.图4为在阵元级分别使用lcmv算法和常规波束形成，在子阵级采用capon算法的天线方向图仿真结果。
具体实施方式
48.下面结合附图和实施例对本发明的方法作进一步描述。
49.本发明的基本方案包括如下六个步骤，具体的信号处理流程如图1所示。
50.步骤一：建立n元线性阵列天线，采用不重叠非均匀划分将其划分为l个子阵，由阵元到子阵的降维矩阵设为t；线性阵列接收入射信号，入射信号由空间中p个互不相关的干扰信号和噪声组成；且n、l、p为正整数，；
51.步骤二：通过对子阵接收信号采样，在设定的快拍数下估计子阵级协方差矩阵，对子阵级协方差矩阵进行特征分解，得到分解后的干扰子空间和噪声子空间，并对该结果采用music算法进行空间谱估计，得到空间干扰的个数m和空间干扰导向矢量；
52.步骤三：对步骤二中得到的空间干扰的个数m进行判断，若，则执行步骤四；若，则执行步骤五；
53.步骤四：由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵c和约束响应向量f，采用lcmv波束形成算法，计算得到阵元级模拟权矢量，并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器，返回步骤二；
54.步骤五：在子阵级采用capon波束形成算法计算得到子阵级数字权矢量，并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器；
55.步骤六：将步骤五中各子阵通道的信号合成为一路波束作为输出。
56.基于上述基本方案，其中步骤一划分子阵的模型为：
57.先建立信号模型，设有一线性阵列具有阵元数为n，阵元之间的距离均为d，将该线性阵列划分为l个子阵，形成了l个接收通道。现假设存在p个互不相关的干扰信号，干扰信号的复包络为，干扰信号的导向矢量矩阵为，背景噪声为，则阵元接收信号。
58.子阵划分基于降维矩阵，其中为加权系数对角阵，为的子阵形成矩阵，在其第1列的所有元素中，与第l个子阵的阵元序号对应的元素值为1，其余均为0，其中l为1-l之间的所有正整数，则子阵接收信号为。
59.若在k时刻对信号进行采样，得到分别为子阵级，阵元级在k时刻接收的样本满足，且子阵级协方差矩阵与阵元级协方差矩阵的关系为：；子阵级导向矢量与阵元级导向矢量的关系为：
60.，
61.其中：
62.。
63.基于上述基本方案，其中步骤二计算接收阵列的子阵级协方差矩阵和特征分解的方法为：
64.求取协方差矩阵理论上需要准确地了解接收信号的统计特性，但实际上只能通过在一定的快拍数k下进行采样得到对协方差矩阵的估计。
65.对子阵接收信号进行采样，得到子阵在k个快拍数下的采样矩阵，则子阵级采样协方差矩阵为：。对进行特征值分解得到：
66.；
67.其中为子阵级采样协方差矩阵的特征值，为特征值对应的特征向量，其中
68.，且
69.为干扰子空间对应的特征值， 为噪声子空间对应的特征值，即为估计的干扰子空间；即为估计的噪声子空间。
70.基于上述基本方案，其中步骤二采用music算法得到空间干扰的个数m及其导向矢量的方法为：
71.理论上，信号空间与噪声空间是正交的，因此信号空间方向矩阵a的列向量与噪声空间的特征向量也是正交的，而a的列向量与信号的到达方向一一对应，根据此性质，便可以利用噪声空间的特征向量来求解信号的到达方向。首先构建噪声矩阵即步骤二中得到的噪声子空间，然后定义空间谱函数：
72.，其中为信号的导向矢量，为估计的噪声子空间。
73.由空间谱函数的定义可以看出，当和正交时，中的从遍历，于是会出现一个谱峰，谱峰对应的值便是信号到达方向的估计值，所以可以通过搜索空间谱函数的峰值来估计信号来波方向。
74.由的峰值对应的角度，便可以得到空间干扰的导向矢量。
75.基于上述基本方案，在进行完步骤二之后，得到了估计的空间干扰的导向矢量，在步骤四中为了得到阵元级模拟权矢量，需要构建约束矩阵和约束响应向量，使得算法能够自适应地抑制干扰信号并保持在波束指向的增益为常数，本实施方式构建lcmv算法流程如下：
76.由约束矩阵；
77.约束响应向量，
78.其中；
79.得到算法模型；
80.求得。
81.根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器；这样虽然上述步骤是在子阵级进行采样估计的，但是通过阵元级的配置可以带来不少的优点，如既能实现基于子阵方式的设备量较少，计算量少的优点又能够提高该系统的自由度，即可以提高由于子阵划分造成的对干扰抑制的自由度不足的问题。该方法经过一次循环能够抑制子阵个数级别的干扰，经过若干次循环后计算得到的阵元级模拟权矢量，最终能够抑制阵元个数级别的干扰，克服了其它基于子阵的自适应阵列信号处理方法可用自由度变小的缺点。
82.基于上述基本方案，其中步骤五中在子阵级采用capon波束形成算法得到子阵级数字权矢量的方法为：
83.将以下优化问题的解作为子阵级数字权矢量
[0084][0085]
求得。
[0086]
这样在空间干扰数目过多时，通过若干次循环而配置阵元级模拟权矢量，使得大部分的干扰得到了抑制，最后一次循环时子阵采样估计出的干扰已经很少，可以在子阵级进行干扰抑制，从而结合了在阵元级进行模拟波束形成和子阵级进行数字波束形成的优点，使得该系统相比普通的基于子阵的自适应波束形成方式抗干扰能力更强，测角精度更高。
[0087]
进行阵列方向图的仿真，仿真采用均匀线阵，仿真参数如表1所示。
[0088]
表1仿真参数设置
[0089]
[0090]
图2是在阵元级使用常规波束形成，在子阵级分别采用lcmv算法和capon算法的天线方向图仿真结果；
[0091]
图3为全阵列使用lcmv算法的天线方向图仿真结果；
[0092]
图4为在阵元级分别使用lcmv算法和常规波束形成，在子阵级采用capon算法的天线方向图仿真结果。
[0093]
由图2可以看到lcmv算法与capon算法有各自的优点，由仿真的结果图可以看出lcmv算法由于约束条件更精细，所以对干扰的抑制效果更好；而本发明设计的系统在阵元级使用lcmv算法得到阵元级模拟权矢量，在模拟部分就可以精确地将大部分干扰抑制掉，而在子阵级采取capon算法得到子阵级数字权矢量，相当于配合模拟部分抑制掉剩余的干扰，具有更强大的干扰抑制能力。
[0094]
由图3可以看到，虽然全阵列的lcmv波束形成方式具有很强的抗干扰能力，但是基于全阵列的实现方式，一般会造成设备量大，成本高等缺点。
[0095]
由图4可以看到相比于一般的基于子阵的自适应波束形成方法，如阵元cbf+子阵capon，对于20
°
方向的干扰没有抑制能力，而本系统的基于子阵的数字模拟合成自适应抗干扰方法明显具有更强的抗干扰能力。