一种用于相控阵卫通天线的波达角零陷展宽方法及系统

摘要

本发明提供了一种用于相控阵卫通天线的波达角零陷展宽方法及系统，属于卫星通信抗干扰技术领域，方法包括：运用相控阵卫通天线接收辐射信号，估计干扰信号的波达方向；确定干扰信号的零陷展宽区间；利用牛顿科茨公式求积分对接收阵元的导向矢量相关矩阵进行简化；对接收阵元的导向矢量相关矩阵特征分解，构建特征值子空间在零陷展宽区间内做线性约束抑制；采用投影变换矩阵对辐射信号的协方差矩阵重构，并加入修正因子对角加载；计算自适应零限展宽最优加权权值；根据最优加权权值生成干扰抑制波形。本发明通过零陷展宽和加深方法，克服了高速运动平台阵列的自适应波束的零陷角度和干扰信号的入射角度失配而造成阵列抗干扰性能迅速下降的缺点。

说明书

技术领域

本发明属于卫星通信抗干扰技术领域，更具体地，涉及一种用于相控阵卫通天线的波达角零陷展宽方法及系统。

背景技术

相控阵卫通天线通常采用多频段集成一体的卫通天线接收多个频点的下行信号，具有低剖面、高精度控制的特点，在多阵元合成高增益主波束的同时，具有较高的副瓣电平，因而同频段空间电磁干扰很有可能通过副瓣区域混入接收天线，影响卫通天线的接收信噪比。由于相控阵天线的集成化射频前端设计，导致传统相控阵天线在信号对准过程中，缺乏面对高动态场景下稳健的波达角信号处理方法。鉴于此，有必要提供一种用于相控阵卫通天线的波达角处理方法。

发明内容

针对现有技术的缺陷，本发明的目的在于提供一种用于相控阵卫通天线的波达角零陷展宽方法及系统，旨在解决由于相控阵天线的集成化射频前端设计，导致传统相控阵天线在信号对准过程中，缺乏面对高动态场景下稳健的波达角信号处理方法的问题。

为实现上述目的，一方面，本发明提供了一种用于相控阵卫通天线的波达角零陷展宽方法，包括以下步骤：

运用相控阵卫通天线接收辐射信号，并预估计辐射信号中干扰信号的波达方向；

根据预估计干扰信号的波达方向，确定干扰信号的零陷展宽区间；

将接收阵元的导向矢量相关矩阵表示为接收阵元的导向矢量及其共轭转置的乘积在干扰信号的零陷展宽区间内的积分，利用牛顿科茨公式求积分对接收阵元的导向矢量相关矩阵进行简化；

对简化后的接收阵元的导向矢量相关矩阵进行特征分解，将特征值从大到小排列，筛选部分特征值对应的特征向量构建特征值子空间；

基于特征值子空间与接收阵元的权重系数向量之间的乘积为0，结合接收阵元的导向矢量，在零陷展宽区间内做线性约束抑制；

采用特征值子空间中的特征向量，构建投影变换矩阵；并采用投影变换矩阵对辐射信号的协方差矩阵进行重构；

对重构后的辐射信号的协方差矩阵加入修正因子进行对角加载；

利用修正后的辐射信号的协方差矩阵、特征值子空间以及期望信号的导向矢量，获取自适应零陷展宽最优加权权值；

根据最优加权权值生成干扰抑制波形；

其中，接收阵元等间距排布，阵元间距

进一步优选地，接收阵元的导向矢量相关矩阵为：

简化后的接收阵元的导向矢量相关矩阵为：

其中，

进一步优选地，筛选部分特征值的方法为：

根据

其中，

进一步优选地，线性约束抑制为

进一步优选地，辐射信号的协方差矩阵重构方式为

其中，

进一步优选地，自适应零陷展宽最优加权权值为：

其中，

另一方面，本发明提供了一种用于相控阵卫通天线的波达角零陷展宽系统，包括：

相控阵卫通天线，用于接收辐射信号；其中，辐射信号包括干扰信号和期望信号；

波达方向预估模块，用于预估计辐射信号中干扰信号的波达方向；

导向矢量相关矩阵的简化模块，用于将接收阵元的导向矢量相关矩阵表示为接收阵元的导向矢量及其共轭转置的乘积在干扰信号的零陷展宽区间内的积分，利用牛顿科茨公式求积分对接收阵元的导向矢量相关矩阵进行简化；

特征值子空间的构建模块，用于对简化后的接收阵元的导向矢量相关矩阵进行特征分解，将特征值从大到小排列，筛选部分特征值对应的特征向量构建特征值子空间；

线性约束抑制的构建模块，用于基于特征值子空间与接收阵元的权重向量之间的乘积为0，结合接收阵元的导向矢量，在零陷展宽区间内做线性约束抑制；

投影变换矩阵的构建模块，用于采用特征子子空间中的特征向量，构建投影变换矩阵；

协方差矩阵的重构模块，用于采用投影变换矩阵对辐射信号的协方差矩阵进行重构；

协方差矩阵的修正模块，用于对重构后辐射信号的协方差矩阵加入修正因子进行对角加载；

自适应零限展宽最优加权权值的计算模块，用于利用修正后辐射信号的协方差矩阵、特征子空间以及期望信号的导向矢量，获取自适应零限展宽最优加权权值；

干扰抑制波形的生成模块，用于根据最优加权权值生成干扰抑制波形；

其中，接收阵元等间距排布，阵元间距

进一步，接收阵元的导向矢量相关矩阵为：

简化后的接收阵元的导向矢量相关矩阵为：

其中，

进一步优选地，筛选部分特征值的方法为：

根据

其中，

进一步优选地，辐射信号的协方差矩阵重构方式为

其中，

自适应零陷展宽最优加权权值为：

其中，

总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：

本发明提供了自适应零陷展宽最优加权权值

本发明通过投影变换与线性约束抑制结合的零陷展宽方法，可实时针对干扰信号的形成鲁棒的波束方向图，具有快拍影响小、零陷深、零陷宽度可控的性能优势。

附图说明

图1是本发明实施例提供的相控阵卫通天线的拓扑示意图；

图2是本发明实施例提供的相控阵卫通天线的波达角零陷展宽方法流程图；

图3是本发明实施例提供的几种不同算法获取的归一化波束图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

图1为相控阵卫通天

步骤一：运用相控阵卫通天线接收辐射信号，并设置初始时刻各接收阵元的权值系数为1；波达角零陷展宽方法适用于相控阵形式卫通接收阵列，；接收阵元等间距排布，阵元间距

步骤二：预估计干扰信号的波达方向；

波达方向预估计利用波达角估计算法获取，可以为典型MUSIC算法或其他先验估计算法实现；

步骤三：根据预估计干扰信号的波达方向，确定干扰信号的零陷展宽区间，计算接收阵元的导向矢量相关矩阵

为了简化算法复杂度，接收阵元的导向矢量相关矩阵可简化为

步骤四：对

特征分解为

步骤五：在零陷展宽区间内做线性约束抑制；

线性约束抑制为

步骤六：采用特征值子空间

对辐射信号的协方差矩阵重构方式为

步骤七：利用修正后的协方差矩阵

权值系数更新准则一般选用LCMV算法或MVDR算法；

步骤八：根据最优加权权值生成干扰抑制波形。

另一方面，本发明提供了一种用于相控阵卫通天线的波达角零陷展宽系统，包括：

相控阵卫通天线，用于接收辐射信号；其中，辐射信号包括干扰信号和期望信号；

波达方向预估模块，用于预估计辐射信号中干扰信号的波达方向；

导向矢量相关矩阵的简化模块，用于将接收阵元的导向矢量相关矩阵表示为接收阵元的导向矢量及其共轭转置的乘积在干扰信号的零陷展宽区间内的积分，利用牛顿科茨公式求积分对接收阵元的导向矢量相关矩阵进行简化；

特征值子空间的构建模块，用于对简化后的接收阵元的导向矢量相关矩阵进行特征分解，将特征值从大到小排列，筛选部分特征值对应的特征向量构建特征值子空间；

线性约束抑制的构建模块，用于基于特征值子空间与接收阵元的权重向量之间的乘积为0，结合接收阵元的导向矢量，在零陷展宽区间内做线性约束抑制；

投影变换矩阵的构建模块，用于采用特征子子空间中的特征向量，构建投影变换矩阵；

协方差矩阵的重构模块，用于采用投影变换矩阵对辐射信号的协方差矩阵进行重构；

协方差矩阵的修正模块，用于对重构后辐射信号的协方差矩阵加入修正因子进行对角加载；

自适应零限展宽最优加权权值的计算模块，用于利用修正后辐射信号的协方差矩阵、特征子空间以及期望信号的导向矢量，获取自适应零限展宽最优加权权值；

干扰抑制波形的生成模块，用于根据最优加权权值生成干扰抑制波形；

其中，接收阵元等间距排布，阵元间距

进一步，接收阵元的导向矢量相关矩阵为：

简化后的接收阵元的导向矢量相关矩阵为：

其中，

进一步优选地，筛选部分特征值的方法为：

根据

其中，

进一步优选地，辐射信号的协方差矩阵重构方式为

其中，

自适应零陷展宽最优加权权值为：

其中，

为验证算法的有效性，假设有2个独立的干扰信号从俯仰角[-30°，40°]，方位角为90°入射平面阵列；平面阵为10×10的平面阵列。其阵元间距为半波长（

设置需要展宽的零陷区间为Δ

综上所述，本发明与现有技术相比，存在以下优势：

（1）本发明提供了自适应零陷展宽最优加权权值

（2）本发明通过投影变换与线性约束抑制结合的零陷展宽方法，可实时针对干扰信号的形成鲁棒的波束方向图，具有快拍影响小、零陷深、零陷宽度可控的性能优势。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。