一种干扰协方差矩阵稳健重构方法技术

本发明专利技术属于阵列信号处理领域，主要涉及干扰噪声协方差矩阵重构算法的稳健性和基于该算法的最差性能最优化Capon自适应波束形成的稳健性。本发明专利技术首先采用最差性能最佳化准则建立期望信号和干扰信号导向矢量的误差模型，然后采用鲁棒Capon波束形成(Robust Capon Beamforming,RCB)来分别估计D-1个干扰信号的功率和导向矢量同时对样本协方差矩阵进行EVD来估计阵列接收高斯白噪声功率由此得到考虑干扰信号导向矢量误差的干扰噪声协方差矩阵重构可以有效的针对求和式干扰噪声协方差矩阵重构的固有不足，有效提高波束形成算法的稳健性。

【技术实现步骤摘要】
一种干扰协方差矩阵稳健重构方法
本专利技术属于阵列信号处理领域，主要涉及干扰噪声协方差矩阵重构算法的稳健性和基于该算法的最差性能最优化Capon自适应波束形成的稳健性。
技术介绍
Capon自适应波束形成算法可以在保证对期望信号无失真输出的条件下，使阵列输出功率最小，最大限度的提高输出信干噪比(Signal-to-Interference-plus-NoiseRatio，SINR)、最大限度的提高阵列增益，具有较好的方位分辨力和较强的干扰抑制能力。但是，Capon波束形成是建立在对期望信号导向矢量和干扰噪声协方差矩阵均精确已知的假想基础上的，对期望信号导向矢量和干扰噪声协方差矩阵的误差比较敏感，且在实际应用中，干扰噪声协方差矩阵一般是难以得到的，往往以阵列接收数据样本协方差矩阵来代替。在阵列接收数据快拍数有限的情况下，Capon自适应波束形成算法的性能会不可避免的有所下降，尤其是当阵列接收数据中包含有期望信号之时，性能下降的尤为严重。为解决这一问题，AmirLeshem在2012年提出了一种干扰协方差矩阵重构算法(Interference-plus-NoiseCovarianceMatrixReconstruction,INCMR)，即利用阵列接收数据和干扰噪声协方差矩阵的特性重新构造干扰噪声协方差矩阵，将其带入Capon波束形成优化问题，可大大提高算法的性能。假设阵列接收到D个来自远场信源的信号，各个信号的来波方向分别为θd,d＝1,…,D，不失一般性，假设第1个信号为期望信号，其余D-1个均为干扰信号.INCMR算法的核心思想是首先将信号来波方向的整个区间Θ＝[-90°，90°]划分为两个区间Θ1,Θ2，且两个区间满足其中期望信号来波方向在区间Θ1中，即θ1∈Θ1，而所有的干扰信号来波方向均在区间Θ2中，即θd∈Θ2,d＝2,…,D；然后在区间Θ2上按照进行干扰噪声协方差矩阵重构(a(θ)为来波方向角θ的导向矢量，为来波方向角θ的Capon空间谱估计，为阵列接收数据样本协方差矩阵)。但是，在区间Θ2中进行积分运算是非常复杂的，为有效降低计算量，将将区间Θ2进行区间离散化为包含L个元素的角度集合区间离散化的要求是包含各个干扰信号来波方向，即用离散区间求和来代替区间Θ2积分得到干扰噪声协方差矩阵的重构然而，INCMR算法存在其固有的不足，其离散区间必须包含所有干扰信号来波信号方向角度即中的求和项必须要包含所有干扰信号项这就需要精确已知阵列接收到的所有干扰信号的导向矢量a(θd),d＝2,…,D，否则其算法性能会大大降低；但是在实际应用中，干扰信号的导向矢量是未知的，需要采用类似于期望信号导向矢量估计的方法进行估计。因此，研究针对干扰信号导向矢量误差的干扰噪声协方差矩阵稳健重构算法是非常有必要的。
技术实现思路
本专利技术的目的在于提供一种针对干扰信号导向矢量误差的干扰噪声协方差矩阵稳健重构算法(aRobustAlgorithmforInterference-plus-NoiseCovarianceMatrixReconstructionagainstInterferenceVectorErrors,RA-INCMR)，类似于期望信号导向矢量误差的建模，构建干扰信号导向矢量误差的基于最差性能最佳化准则的误差模型，并在此基础上利用干扰噪声协方差矩阵的结构特性来重构考虑干扰信号导向矢量误差的干扰噪声协方差矩阵提高干扰噪声协方差矩阵重构算法对干扰信号导向矢量误差的稳健性。本专利技术的思路是：本专利技术基于理想干扰噪声协方差矩阵的结构特性是干扰信号d的导向矢量，干扰信号d的功率，σ2是阵列接收高斯白噪声功率，IN是N×N单位矩阵，首先采用最差性能最佳化准则建立阵列接收信号(包括期望信号和干扰信号)导向矢量误差的模型，用阵列接收数据来估计样本协方差矩阵在此基础上采用鲁棒Capon波束形成(RobustCaponBeamforming,RCB)来估计干扰信号d的功率和导向矢量同时对样本协方差矩阵进行特征值分解(EVD)得到阵列接收高斯白噪声的功率从而得到重构的干扰噪声协方差矩阵，然后用代替原始最差性能最佳化算法中的样本协方差矩阵按照WCPO准则使得最差情况下的波束输出OSINR最大化。本专利技术的目的通过如下步骤实现：S1、由M个阵元构成的均匀线阵接收到D个来自远场信源的信号，所述D个信号的来波方向分别为θd，设第1个信号为期望信号，其余D-1个均为干扰信号，且假设各个信号之间互不相关，且信号与噪声之间也互不相关，则第n个快拍下阵列接收数据记为其中，A＝[a(θ1),…,a(θD)]为阵列流型矩阵，s(n)为阵列接收到的信号源矢量，v(n)表示阵列接收到的噪声矢量，v(n)为零均值高斯白噪声，d＝1,…,D，则阵列接收到的N个快拍数据为由阵列接收数据矩阵X可以得到阵列接收数据的样本协方差矩阵S2、利用S1所述阵列接收数据的样本协方差矩阵来估计阵列接收的高斯白噪声功率S3、采用鲁棒Capon波束形成(RobustCaponBeamforming,RCB)来估计D-1个干扰信号的功率和导向矢量S31、干扰信号d的预估计导向矢量为真实的信号导向矢量a(θd)位于椭圆不确定集合中，利用S1所述阵列接收数据的样本协方差矩阵来构建干扰信号d的RCB波束形成优化问题：将所述干扰信号d的RCB波束形成优化问题进行一定整理之后转换为如下的半定规划问题：求解所述半定规划问题，可以得到干扰信号d的功率和导向矢量S32、分别取d＝2,…,D，重复步骤S31即可得到干扰噪声协方差矩阵中的干扰信号项S33、根据S2所述和S32所述干扰信号项得到考虑干扰信号导向矢量误差的干扰噪声协方差矩阵重构S4、期望信号的预估计导向矢量为其真实导向矢量a(θ1)位于椭圆不确定集合利用S3中估计的干扰噪声协方差矩阵来构建期望信号的最差性能最优化(WCPO)波束形成优化问题：将其进行一定整理之后转换为如下的二阶锥规划问题：采用已有的SeDuMi软件或CVX软件进行求解，得到其稳健的阵列加权wRA-INCMR。进一步地，S2所述估计阵列接收高斯白噪声功率具体为：对进行特征值分解(EVD)得到从大到小排列的特征值其中，D个大特征值对应于阵列接收到的D个信源信号部分，剩余的M-D个小特征值对应于阵列接收到的噪声部分，则噪声功率可用下式进行估计：本专利技术的有益效果是：首先采用最差性能最佳化准则建立期望信号和干扰信号导向矢量的误差模型，然后采用鲁棒Capon波束形成(RobustCaponBeamforming,RCB)来分别估计D-1个干扰信号的功率和导向矢量同时对样本协方差矩阵进行EVD来估计阵列接收高斯白噪声功率由此得到考虑干扰信号导向矢量误差的干扰噪声协方差矩阵重构可以有效的针对求和式干扰噪声协方差矩阵重构的固有不足，有效提高波束形成算法的稳健性。附图说明图1是本专利技术方法的流程图。图2是针对干扰信号导向矢量误差的干扰噪声协方差矩阵稳健重构算法的波束输出SINR随期望信号输入SNR的变化曲线图。图3是针对干扰信号导向矢量误差的干扰噪声协方差矩阵稳健重构算法的波束输出SINR随阵列接收数据快拍数的变化曲线图。具体实施方式下面结合实施例和附图，详细说明本专利技术的技术方案。图1是本专利技术本文档来自技高网...

【技术保护点】
一种干扰协方差矩阵稳健重构方法，其特征在于，包括如下步骤：S1、由M个阵元构成的均匀线阵接收到D个来自远场信源的信号，所述D个信号的来波方向分别为θd，设第1个信号为期望信号，其余D‑1个均为干扰信号，且假设各个信号之间互不相关，且信号与噪声之间也互不相关，则第n个快拍下阵列接收数据记为x(n)=a(θ1)s1(n)+Σd=2Da(θd)sd(n)+v(n)=As(n)+v(n),]]>其中，A＝[a(θ1),…,a(θD)]为阵列流型矩阵，s(n)为阵列接收到的信号源矢量，v(n)表示阵列接收到的噪声矢量，v(n)为零均值高斯白噪声，d＝1,…,D，则阵列接收到的N个快拍数据为X=[x(1),...,x(N)]=AS+VS=[s(1),...,s(N)]V=[v(1),...,v(N)],]]>由阵列接收数据矩阵X可以得到阵列接收数据的样本协方差矩阵S2、利用S1所述阵列接收数据的样本协方差矩阵来估计阵列接收的高斯白噪声功率S3、采用鲁棒Capon波束形成(Robust Capon Beamforming,RCB)来估计D‑1个干扰信号的功率和导向矢量S31、干扰信号d的预估计导向矢量为真实的信号导向矢量a(θd)位于椭圆不确定集合中，利用S1所述阵列接收数据的样本协方差矩阵来构建干扰信号d的RCB波束形成优化问题：maxa(θd),σd2σd2s.t.R^x≥σd2a(θd)aH(θd),||a(θd)-a^(θd)||2≤ϵd]]>将所述干扰信号d的RCB波束形成优化问题进行一定整理之后转换为如下的半定规划问题：mina(θd),σd21/σd2]]>s.t.R^xa(θd)aH(θd)1/σd2≥0]]>ϵdIa(θd)-a^(θd)(a(θd)-a^(θd))H1≥0]]>求解所述半定规划问题，可以得到干扰信号d的功率和导向矢量S32、分别取d＝2,…,D，重复步骤S31即可得到干扰噪声协方差矩阵中的干扰信号项Σd=2Dσ~d2a~(θd)a~H(θd);]]>S33、根据S2所述和S32所述干扰信号项得到考虑干扰信号导向矢量误差的干扰噪声协方差矩阵重构R^i+n=Σd=2Dσ~d2a~(θd)a~H(θd)+σ~2IN;]]>S4、期望信号的预估计导向矢量为其真实导向矢量a(θ1)位于椭圆不确定集合S1=Δ{a(θ1)|a(θ1)=a^(θ1)+δ1,||δ1||2≤ϵ1},]]>利用S3中估计的干扰噪声协方差矩阵来构建期望信号的最差性能最优化(WCPO)波束形成优化问题：minwwHR^i+nw s.t.wHa~(θ1)≥ϵ||w||2+1]]>将其进行一定整理之后转换为如下的二阶锥规划问题：minw||Vw||s.t.ϵ||w||2≤wHa^(θ1)-1,R^i+n=VHV]]>采用已有的SeDuMi软件或CVX软件进行求解，得到其稳健的阵列加权wRA‑INCMR。...

【技术特征摘要】
1.一种干扰协方差矩阵稳健重构方法，其特征在于，包括如下步骤：S1、由M个阵元构成的均匀线阵接收到D个来自远场信源的信号，所述D个信号的来波方向分别为θd，设第1个信号为期望信号，其余D-1个均为干扰信号，且假设各个信号之间互不相关，且信号与噪声之间也互不相关，则第n个快拍下阵列接收数据记为其中，A＝[a(θ1),…,a(θD)]为阵列流型矩阵，s(n)为阵列接收到的信号源矢量，v(n)表示阵列接收到的噪声矢量，v(n)为零均值高斯白噪声，d＝1,…,D，s1(n)为期望信号的时域窄带信号波形，sd(n)为第d个干扰信号的时域窄带信号波形，则阵列接收到的N个快拍数据为由阵列接收数据矩阵X可以得到阵列接收数据的样本协方差矩阵S2、利用S1所述阵列接收数据的样本协方差矩阵来估计阵列接收的高斯白噪声功率S3、采用鲁棒Capon波束形成(RobustCaponBeamforming,RCB)来估计D-1个干扰信号的功率和导向矢量S31、第d个干扰信号的预估计导向矢量为真实的信号导向矢量a(θd)位于椭圆不确定集合中，利用S1所述阵列接收数据的样本协方差矩阵来构建第d个干扰信号的RCB波束形成优化问题，其中，δd表示第d个干扰信号的真实导向矢量a(θd)与其预估的导向矢量为之间的估计误差，εd表示该导向矢量估计误差δd的范数约束上限值，表示第d个干扰信号的信号功率，“I”表示维度为M×M的单位矩阵：将所述第d个干扰信号的RCB波束形成优化问题进行一定整理之后转换为如下的半定规划问题：

【专利技术属性】
技术研发人员：袁晓垒，甘露，杜继萍，廖红舒，张花国，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51