加权子空间自适应天线方向图旁瓣形状精确控制方法技术

本发明专利技术公开了一种加权子空间自适应天线方向图旁瓣形状精确控制方法，包括以下步骤：根据波束指向和期望的旁瓣电平分布确定旁瓣加权函数；使用旁瓣加权函数对旁瓣区导向性矢量进行加权，采用加权后的旁瓣区导向性矢量计算旁瓣协方差矩阵，取其主特征向量构建旁瓣子空间矩阵；将自适应权重向旁瓣子空间矩阵投影，通过不等式约束其模值；将该不等式约束和MVDR优化问题联立，形成改进的MVDR代价函数；采用内点法求解自适应权重。本发明专利技术在保证主瓣准确指向期望信号，自适应抑制旁瓣区连续干扰的同时，能够精确地控制自适应方向图峰值旁瓣电平以抑制突发脉冲干扰，并具备复杂的旁瓣区电平分布控制能力，以适用于特殊的空域干扰和杂波抑制场景。

An Accurate Sidelobe Shape Control Method for Weighted Subspace Adaptive Antenna Patterns

The invention discloses an accurate sidelobe shape control method for weighted subspace adaptive antenna pattern, which includes the following steps: determining the sidelobe weighting function according to the beam direction and expected sidelobe level distribution; weighting the sidelobe steering vector by using the sidelobe weighting function, and calculating the sidelobe covariance matrix by using the weighted sidelobe steering vector, and taking its main characteristic direction. The adaptive weight is projected to the sidelobe subspace matrix, and its modulus is constrained by inequality. The inequality constraint and MVDR optimization problem are combined to form an improved MVDR cost function. The interior point method is used to solve the adaptive weight. While ensuring that the main lobe accurately points to the desired signal and adaptively suppresses the continuous interference in the sidelobe region, the method can accurately control the peak sidelobe level of the adaptive pattern to suppress the burst pulse interference, and has the complex level distribution control capability of the sidelobe region, so as to be suitable for special spatial interference and clutter suppression scenarios.

【技术实现步骤摘要】
加权子空间自适应天线方向图旁瓣形状精确控制方法
本专利技术属于阵列天线空域抗干扰领域，具体涉及一种加权子空间自适应天线方向图旁瓣形状精确控制方法。
技术介绍
自适应数字波束形成技术的显著优势在于它能够在保证天线方向图主瓣指向期望信号的同时自适应地在干扰方向产生零陷。最著名的自适应数字波束形成器是最小方差无失真响应(MinimumVarianceDistortionlessResponse,MVDR)波束形成器和由之推广而来的线性约束最小方差(LinearConstraintMinimumVariance，LCMV)波束形成器，它们都是在满足给定线性约束的条件下最小化波束形成器的输出功率从而自适应地抑制干扰。MVDR波束形成器可以表示为：w＝argminwHRxws.t.wHa(θ0)＝1。但是对一些时间快变干扰，尤其是突发脉冲形式的干扰，自适应波束形成器通常来不及或者无法及时更新权重系数来产生自适应零陷，造成波束形成器的输出性能下降。低旁瓣控制技术，在不更新权重系数的情况下，对此类干扰也具有一定的抑制作用。另一方面，高精度的目标角度估计和跟踪一般采用和差单脉冲测角技术。而单脉冲测角要求测角波束具有稳定的主瓣形状和指向，若主瓣变形或者峰值偏移，不仅影响测角精度，而且会导致输出信干噪比下降。传统的MVDR和LCMV波束形成器往往不具备稳健的方向图主瓣控制能力，当存在接收数据模型误差，如低快拍，主瓣信号，导向性矢量误差或者阵列误差情况下，方向图主瓣特性会严重恶化。因此，研究自适应天线方向图控制尤为重要，同时为了应对多样化的干扰和杂波，复杂的旁瓣形状控制也至关重要。对角加载是一种经典的旁瓣控制方法，通过在协方差矩阵中人为注入噪声提高协方差矩阵的稳健性，从而避免方向图大幅度抖动，但是对角加载量通常难以确定；罚函数方法通过让自适应方向图或者自适应权重逼近预先优化好的静态方向图或者静态权重来达到控制旁瓣电平的目的，然而此类方法一般无法精确控制自适应方向图的峰值旁瓣电平，更不具有复杂的旁瓣形状控制能力。文献1(R.Wu,Z.Bao,Y.L.Ma,“Controlofpeaksidelobelevelinadaptivearrays,”IEEETransactionsonAntennas&Propagation,vol.44,no.10,1996,pp.1341-1347.)通过对角加载方法求解其提出的罚函数模型，并推导出对角加载量和期望峰值旁瓣电平之间准确的数值关系，利用这一关系可以实现峰值旁瓣电平精确控制，但是其精确性依赖于合适的静态权重；文献2(J.Liu,A.B.Gershman,Z.Q.Luo,etal.,“Adaptivebeamformingwithsidelobecontrol:asecond-orderconeprogrammingapproach,”IEEESignalProcessingletters,vol.10,no.11,2003,pp.331-334.)采用多个二次不等式约束直接控制旁瓣电平，该方法保证优化出的峰值旁瓣电平位于期望值以下，但是计算量大，且方向图主瓣指向不稳定。
技术实现思路
本专利技术的目的在于提供一种基于加权子空间的自适应天线方向图旁瓣形状精确控制方法，通过精确地控制方向图峰值旁瓣电平以抑制旁瓣区突发脉冲干扰，同时通过自适应零陷抑制旁瓣区连续干扰，并且保证主瓣稳定准确指向期望信号。实现本专利技术目的的技术方案为：一种加权子空间自适应天线方向图旁瓣形状精确控制方法，包括以下步骤：步骤1，根据波束指向和期望的旁瓣电平分布确定旁瓣加权函数；步骤2，使用旁瓣加权函数对旁瓣区导向性矢量进行加权，采用加权后的旁瓣区导向性矢量计算旁瓣协方差矩阵，取其主特征向量构建旁瓣子空间矩阵；步骤3，将自适应权重向旁瓣子空间矩阵投影，通过不等式约束其模值；将该不等式约束和MVDR优化问题联立，形成改进的MVDR代价函数；步骤4，采用内点法求解自适应权重。本专利技术与现有技术相比，其显著优点为：(1)本专利技术能够精确控制自适应方向图峰值旁瓣电平，实现复杂旁瓣形状控制，具备稳健的主瓣形状和指向控制能力；(2)本专利技术通过单个子空间约束对整个旁瓣进行控制，增加了权重求解的数值稳定性；在构造旁瓣子空间时，引入由期望旁瓣电平分布确定的加权函数，达到拟合期望旁瓣电平分布的目的；(3)本专利技术通过子空间约束能显著减少约束维数，降低运算复杂度；算法可以转化为SOCP问题，能通过内点法高效求解；(4)该方法计算复杂度低，可广泛应用于雷达、通信、声呐、射电天文，语音信号处理等系统中的自适应阵列天线。附图说明图1是本专利技术的算法实现流程图。图2(a)、图2(b)、图2(c)分别是实施例中50阵元均匀直线阵在不同期望旁瓣电平分布下的自适应方向图。图3是实施例中存在一个主瓣信号和两个旁瓣干扰情况下输出信干噪比随采样快拍数的变化曲线图。具体实施方式结合图1，一种基于加权子空间的自适应天线方向图旁瓣形状精确控制方法，包括以下步骤：步骤1，根据波束指向和期望的旁瓣电平分布确定旁瓣加权函数。步骤2，使用旁瓣加权函数对旁瓣区导向性矢量进行加权，采用加权后的旁瓣区导向性矢量计算旁瓣协方差矩阵，取其主特征向量构建旁瓣子空间矩阵；步骤3，将自适应权重向旁瓣子空间矩阵投影，通过不等式约束其模值；将该不等式约束和MVDR优化问题联立，形成改进的MVDR代价函数；步骤4，采用内点法求解自适应权重。进一步的，步骤1具体为：步骤1-1，构建旁瓣加权子函数h1(θ)；假设，方向图期望波束指向为θ0，其旁瓣区范围为Θ，期望的旁瓣电平分布函数为DSL(θ)，θ∈Θ，单位为dB。那么，峰值旁瓣电平可以表示为DPSLA＝10(max(DSL(θ))/20)，θ∈Θ。取该加权子函数用于控制自适应方向图旁瓣电平均匀分布；步骤1-2，构建旁瓣加权子函数h2(θ)；取h2(θ)＝10(-DSL(θ)/20)，θ∈Θ，该加权子函数实现期望旁瓣形状控制；步骤1-3，构建最终旁瓣加权函数h(θ)＝h1(θ)×h2(θ)，θ∈Θ。进一步的，步骤2具体为：步骤2-1，构建加权旁瓣协方差矩阵RΘ；在旁瓣区Θ内均匀选取J个角度θj，j＝1,2,...,J，根据公式计算加权旁瓣协方差矩阵RΘ，a(θj)为阵列导向性矢量，J>>N，N为阵元数，保证RΘ为满秩矩阵。指数p根据期望的旁瓣电平分布函数DSL(θ)设置，若仅需要控制自适应方向图峰值旁瓣电平，对旁瓣形状无要求，则p＝0，即不进行加权；否则p＝1，保证加权函数h(θ)有效。步骤2-2，构建旁瓣区约束矩阵VΘ；对RΘ进行特征值分解，将特征值按降序排列，λn为RΘ的第n个特征值，vn为对应的归一化特征向量；取前M个特征向量构成旁瓣子空间VΘ＝[v1,v2,...,vM]，以VΘ作为旁瓣区约束矩阵。进一步的，步骤3具体为：确定最优波束形成器优化模型；优化问题的代价函数为上式中，Rx为阵列接收信号协方差估计矩阵，可以表示为其中K为采样快拍数，x(k)为阵列接收信号向量。wHa(θ0)＝1为期望信号单位增益约束，为旁瓣子空间不等式约束，DPSLA＝10(max(DSL(θ))/20)为期望的峰值旁瓣电平幅度。根据上述描述，总结本专利技术的本文档来自技高网...

【技术保护点】
1.一种基于加权子空间的自适应天线方向图旁瓣形状精确控制方法，其特征在于，包括以下步骤：步骤1，根据波束指向和期望的旁瓣电平分布确定旁瓣加权函数；步骤2，使用旁瓣加权函数对旁瓣区导向性矢量进行加权，采用加权后的旁瓣区导向性矢量计算旁瓣协方差矩阵，取其主特征向量构建旁瓣子空间矩阵；步骤3，将自适应权重向旁瓣子空间矩阵投影，通过不等式约束其模值；将该不等式约束和MVDR优化问题联立，形成改进的MVDR代价函数；步骤4，采用内点法求解自适应权重。

【技术特征摘要】
1.一种基于加权子空间的自适应天线方向图旁瓣形状精确控制方法，其特征在于，包括以下步骤：步骤1，根据波束指向和期望的旁瓣电平分布确定旁瓣加权函数；步骤2，使用旁瓣加权函数对旁瓣区导向性矢量进行加权，采用加权后的旁瓣区导向性矢量计算旁瓣协方差矩阵，取其主特征向量构建旁瓣子空间矩阵；步骤3，将自适应权重向旁瓣子空间矩阵投影，通过不等式约束其模值；将该不等式约束和MVDR优化问题联立，形成改进的MVDR代价函数；步骤4，采用内点法求解自适应权重。2.根据权利要求1所述的基于加权子空间的自适应天线方向图旁瓣形状精确控制方法，其特征在于，步骤1具体为：步骤1-1，构建旁瓣加权子函数h1(θ)；假设方向图期望波束指向为θ0，其旁瓣区范围为Θ，期望的旁瓣电平分布函数为DSL(θ)，θ∈Θ，单位为dB；那么，峰值旁瓣电平可表示为DPSLA＝10(max(DSL(θ))/20)，θ∈Θ；取该加权子函数用于控制自适应方向图旁瓣电平均匀分布；步骤1-2，构建旁瓣加权子函数h2(θ)；取h2(θ)＝10(-DSL(θ)/20)，θ∈Θ，该加权子函数实现期望旁瓣形状控制；步骤1-3，构建最终旁瓣加权函数h(θ)＝h1(θ)×h2(θ)，θ∈Θ。3.根据权利要求1所述的基于加...

【专利技术属性】
技术研发人员：马晓峰，周淼，盛卫星，韩玉兵，张仁李，崔杰，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32