基于干扰加噪声协方差矩阵重构的自适应波束形成方法技术

本发明专利技术属于阵列信号处理技术领域，公开了一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，包括：建立均匀直线阵列，用于检测包含1个待检测期望信号和Q‑1个干扰信号的远场窄带信号，得到接收信号；对接收信号采样，得到K个接收信号样本，组成样本矩阵，计算样本矩阵的干扰加噪声采样协方差矩阵；构造加权矩阵，对干扰加噪声采样协方差矩阵加权，得到加权后的协方差矩阵；构造采样矩阵，进而构造第一干扰加噪声协方差矩阵；根据加权矩阵对第一干扰加噪声协方差矩阵再次加权，得到第二干扰加噪声协方差矩阵；根据第二干扰加噪声协方差矩阵计算自适应波束形成器的自适应权矢量；能够在保证期望信号的检测效果的基础上减少运算量。

【技术实现步骤摘要】
基于干扰加噪声协方差矩阵重构的自适应波束形成方法
本专利技术属于阵列信号处理
，尤其涉及一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，适用于阵列天线信号处理系统中自适应波束形成器的设计。
技术介绍
阵列信号处理在通信、雷达、声纳、电子对抗、医疗成像、射电天文等领域具有广泛的应用，波束形成技术是阵列信号处理的一个重要分支。随着相控阵天线在雷达以及通信系统中的普及，波束形成技术与方法也有着飞速的发展与改进。原始的波束形成技术通过对不同的阵列传感器施加不同的馈电相位使得阵列天线接收的复信号输出相位为同一个方向，使得阵列输出增益在该方向达到最大，这也是相控阵天线的基本工作原理。传统的波束形成方法虽然可以在需要的方向上形成高的增益，但是不具有抑制干扰信号的功能；在20世纪60年代，Capon提出的最小方差无失真响应(MVDR)波束形成器推导出通过理论上的干扰加噪声协方差矩阵计算阵列自适应权值的方法，从而在保证期望方向上的信号的增益的情况下在干扰方向上形成响应零陷，理论上可以有效抑制干扰信号；之后出现的使用接收信号采样协方差矩阵代替理论上的干扰加噪声协方差矩阵的采样矩阵求逆(SMI)波束形成器能够根据阵列天线的接收信号自适应地抑制干扰，是经典的自适应波束形成器。当前的波束形成方法的研究主要集中于接收信号采样协方差矩阵含有期望信号分量的情况，在考虑实际工程中的阵列幅相误差、信号观测误差以及阵列校准误差时，采样协方差矩阵中的期望信号分量的存在使得求解的自适应权值会对期望信号产生相消响应，从而导致自适应波束形成器性能的显著下降；针对该问题，2003年Vorobyov提出了基于椭圆不确定集优化方程的最差性能最优化方法，在一定程度上削减了采样协方差矩阵中期望信号分量对于自适应波束形成器性能的影响，然而该方法的性能受先验参数的影响，并且在高信噪比的信号环境下性能较差。2012年，Gu等人提出了基于曲线积分的干扰加噪声协方差矩阵重构的方法，该方法虽然能够较好地消除采样协方差矩阵中的期望信号分量，但是对于干扰信号的角度估计不够准确，并且大区域的积分运算造成较大的时间计算复杂度，使得该方法在实际工程中的应用范围有限。
技术实现思路
针对上述现有技术的不足，本专利技术提供一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，能够在保证期望信号的检测效果的基础上减少运算量。为达到上述目的，本专利技术采用如下技术方案予以实现。一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，所述方法包括如下步骤：步骤1，建立具有M个阵元的均匀直线阵列，所述均匀直线阵列用于接收Q个远场窄带信号，所述均匀直线阵列接收的所述Q个远场窄带信号包含1个待检测期望信号和Q-1个干扰信号，所述1个待检测期望信号的入射方向为θ0，所述Q-1个干扰信号的入射方向分别为θ1,θ2,…,θQ-1；根据所述具有M个阵元的均匀直线阵列、所述Q个远场窄带信号，得到所述均匀直线阵列的接收信号t表示时间变量；步骤2，对所述均匀直线阵列的接收信号进行采样，得到K个接收信号样本，分别记为所述K个接收信号样本组成M×K维样本矩阵XK，根据所述M×K维样本矩阵计算所述M×K维样本矩阵的M×M维干扰加噪声采样协方差矩阵t1,t2,…,tK为K个不同的采样时刻；步骤3，构造加权矩阵，根据所述加权矩阵对所述干扰加噪声采样协方差矩阵加权，得到加权后的协方差矩阵R；步骤4，构造采样矩阵，根据所述采样矩阵和所述加权后的协方差矩阵R，构造第一干扰加噪声协方差矩阵步骤5，根据所述加权矩阵对所述第一干扰加噪声协方差矩阵再次加权，得到第二干扰加噪声协方差矩阵Ri+n；步骤6，根据所述第二干扰加噪声协方差矩阵Ri+n计算自适应波束形成器的M×1维的自适应权矢量，从而根据所述自适应权矢量进行自适应波束形成。本专利技术技术方案通过空间功率谱采样来重构干扰加噪声协方差矩阵，用重构的干扰加噪声协方差矩阵计算自适应波束形成器的自适应权矢量，不仅避免了空间功率谱估计的步骤，减少了运算量，而且得到的输出信号能够保持高的输出信杂噪比，保证了期望信号的检测效果。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例提供的一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法的流程示意图；图2为样本数为60的情况下，分别使用传统的采样矩阵求逆方法、最差性能最优化方法和本专利技术方法设计的自适应波束形成器示意图；其中，横坐标为期望信号的输入信噪比，纵坐标为自适应波束形成器的输出信干噪比的大小，单位为分贝(dB)；图3为输入信噪比为10dB的情况下，分别使用传统的采样矩阵求逆方法、最差性能最优化方法和本专利技术方法设计的自适应波束形成器示意图；其中，横坐标为采样样本数，纵坐标为自适应波束形成器的输出信干噪比的大小，单位为分贝(dB)。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术实施例提供种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，如图1所示，所述方法包括如下步骤：步骤1，建立具有M个阵元的均匀直线阵列，所述均匀直线阵列用于接收Q个远场窄带信号，所述均匀直线阵列接收的所述Q个远场窄带信号包含1个待检测期望信号和Q-1个干扰信号，所述1个待检测期望信号的入射方向为θ0，所述Q-1个干扰信号的入射方向分别为θ1,θ2,…,θQ-1；根据所述具有M个阵元的均匀直线阵列、所述Q个远场窄带信号，得到所述均匀直线阵列的接收信号t表示时间变量。所述步骤1具体为：(1a)所述均匀直线阵列的接收信号其中，表示M×1维待检测期望信号分量，表示M×1维干扰信号分量，表示M×1维高斯白噪声分量，t表示时间变量；(1b)将入射方向为θ0的待检测期望信号记为X(θ0)，维数为M×1维，其表达式为：X(θ0)＝a(θ0)S(θ0)，从而得到待检测期望信号分量S(θ0)表示入射方向为θ0的待检测期望信号的幅度的复包络；a(θ0)表示入射方向为θ0的待检测期望信号的M×1维导向矢量，其表达为：其中，t表示时间变量，M表示均匀直线阵列包含的阵元个数，d表示均匀直线阵列的M个阵元的阵元间距，λ表示信号的波长，上标T表示转置；(1c)将入射方向为θq的干扰信号记为X(θq)，其表达式为：X(θq)＝a(θq)S(θq)，θq表示第q个干扰信号的入射方向，q∈{1,2,…,Q-1}，Q-1为均匀直线阵列接收不同干扰方向信号的个数，S(θq)表示入射方向为θq的干扰信号的幅度的复包络，a(θq)表示入射方向为θq的干扰信号的M×1维导向矢量；然后令q分别取1至Q-1，进而得到均匀直线阵列接收的Q-1个干扰信号分量其表达式为：(1d)从而所述均匀直线阵列的接收信号表示为：q'∈{0,1,…,Q-1}，A表示均匀直线阵列流本文档来自技高网...

【技术保护点】
一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，其特征在于，所述方法包括如下步骤：步骤1，建立具有M个阵元的均匀直线阵列，所述均匀直线阵列用于接收Q个远场窄带信号，所述均匀直线阵列接收的所述Q个远场窄带信号包含1个待检测期望信号和Q‑1个干扰信号，所述1个待检测期望信号的入射方向为θ

【技术特征摘要】
1.一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，其特征在于，所述方法包括如下步骤：步骤1，建立具有M个阵元的均匀直线阵列，所述均匀直线阵列用于接收Q个远场窄带信号，所述均匀直线阵列接收的所述Q个远场窄带信号包含1个待检测期望信号和Q-1个干扰信号，所述1个待检测期望信号的入射方向为θ0，所述Q-1个干扰信号的入射方向分别为θ1,θ2,…,θQ-1；根据所述具有M个阵元的均匀直线阵列、所述Q个远场窄带信号，得到所述均匀直线阵列的接收信号t表示时间变量；步骤2，对所述均匀直线阵列的接收信号进行采样，得到K个接收信号样本，分别记为所述K个接收信号样本组成M×K维样本矩阵XK，根据所述M×K维样本矩阵计算所述M×K维样本矩阵的M×M维干扰加噪声采样协方差矩阵t1,t2,…,tK为K个不同的采样时刻；步骤3，构造加权矩阵，根据所述加权矩阵对所述干扰加噪声采样协方差矩阵加权，得到加权后的协方差矩阵R；步骤4，构造采样矩阵，根据所述采样矩阵和所述加权后的协方差矩阵R，构造第一干扰加噪声协方差矩阵步骤5，根据所述加权矩阵对所述第一干扰加噪声协方差矩阵再次加权，得到第二干扰加噪声协方差矩阵Ri+n；步骤6，根据所述第二干扰加噪声协方差矩阵Ri+n计算自适应波束形成器的M×1维的自适应权矢量，从而根据所述自适应权矢量进行自适应波束形成。2.根据权利要求1所述的一种基于干扰加噪声协方差矩阵重构的自适应波束形成方法，其特征在于，所述步骤1具体为：(1a)所述均匀直线阵列的接收信号其中，表示M×1维待检测期望信号分量，表示M×1维干扰信号分量，表示M×1维高斯白噪声分量，t表示时间变量；(1b)将入射方向为θ0的待检测期望信号记为X(θ0)，维数为M×1维，其表达式为：X(θ0)＝a(θ0)S(θ0)，从而得到待检测期望信号分量S(θ0)表示入射方向为θ0的待检测期望信号的幅度的复包络；a(θ0)表示入射方向为θ0的待检测期望信号的M×1维导向矢量，其表达为：其中，t表示时间变量，M表示均匀直线阵列包含的阵元个数，d表示均匀直线阵列的M个阵元的阵元间距，λ表示信号的波长，上标T表示转置；(1c)将入射方向为θq的干扰信号记为X(θq)，其表达式为：X(θq)＝a(θq)S(θq)，θq表示第q个干扰信号的入射方向，q∈{1,2,…,Q-1}，Q-1为均匀直线阵列接收不同干扰方向信号的个数，S(θq)表示入射方向为θq的干扰信号的幅度的复包络，a(θq)表示入射方向为θq的干扰信号的M×1维导向矢量；然后令q分别取1至Q-1，进而得到均匀直线阵列接收的Q-1个干扰信号分量其表达式为：(1d)从而所述均匀直线阵列的接收信号表示为；q'∈{0,1,…,Q-1}，A表示均匀直线阵列流型矩阵，是由Q个不同方向信号的导向矢量组成的矩阵，A＝[a(θ0)a(θ1)a(θ2...

【专利技术属性】
技术研发人员：王彤，蔡启程，李博文，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61