CN113422629A - 一种协方差矩阵重构自适应波束形成方法及系统 - Google Patents

一种协方差矩阵重构自适应波束形成方法及系统 Download PDF

Info

Publication number
CN113422629A
CN113422629A CN202110673467.2A CN202110673467A CN113422629A CN 113422629 A CN113422629 A CN 113422629A CN 202110673467 A CN202110673467 A CN 202110673467A CN 113422629 A CN113422629 A CN 113422629A
Authority
CN
China
Prior art keywords
covariance matrix
data
array
virtual
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110673467.2A
Other languages
English (en)
Other versions
CN113422629B (zh
Inventor
陈鹏
高婧洁
王威
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changan University
Original Assignee
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN202110673467.2A priority Critical patent/CN113422629B/zh
Publication of CN113422629A publication Critical patent/CN113422629A/zh
Application granted granted Critical
Publication of CN113422629B publication Critical patent/CN113422629B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B15/00Suppression or limitation of noise or interference
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开一种基于虚拟线列阵的协方差矩阵重构自适应波束形成方法及系统,首先利用均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并按照递推的方式依次对虚拟阵元位置出所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;然后,利用包含实际阵元和虚拟阵元的数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,最终添加入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;最后,通过同样的方式重构虚拟阵列对应的期望信号协方差矩阵,并提取出最大特征值对应的特征向量作为导向矢量,从而计算得到虚拟阵列对应的波束形成器加权向量和波束输出。

Description

一种协方差矩阵重构自适应波束形成方法及系统
技术领域
本发明涉及阵列信号处理领域,涉及一种协方差矩阵重构自适应波束形成方法及系统。
背景技术
在传感器阵列的自适应波束形成技术在过去的几十年里,得到了迅速的发展,相关成果已经广泛应用于各个领域。作为最著名的自适应波束形成器,Capon波束形成器能够在训练数据不含期望信号时,可以达到其理论性能的最优。然而,在实际的应用环境下,阵列的接收数据中往往存在期望信号,尤其是在信噪比较高时,会导致经典的自适应波束形成方法对期望信号导向矢量的失配更加敏感。
为了提高自适应波束形成器的稳健性,近年来国内外学者提出了许多稳健波束形成器,根据使用的协方差矩阵归纳起来有以下两类:
1)文献1“Robust Adaptive Beamforming,”IEEE Transactions on AcousticsSpeech and Signal Processing,vol.35,no.10,pp.1365-1376,Oct,1987.”和文献2“Robust Capon beamforming,”IEEE Signal Processing Letters,vol.10,no.6,pp.172-175,Jun,2003.”中公开的对角加载方法和不确定集约束方法能够较好的改善稳健性,但对角加载量和不确定集大小难以根据实际情况确定,难以实际使用。
2)文献3“Robust Adaptive Beamforming Based on Interference CovarianceMatrix Reconstruction and Steering Vector Estimation,IEEE Trans.SignalProcess.,vol.60,no.7,pp.3881-3885,Jul.2012.”中公开的重构干扰加噪声协方差矩阵(interference-plus-noise covariance matrix,INCM)算法能够有效地去除协方差矩阵中的期望信号分量,从而降低波束形成器对期望信号导向矢量误差的敏感性,但同时重构的INCM中可能存在干扰的导向矢量误差,从而降低了干扰抑制能力。
综上所述,第1)类方法直接利用接收数据协方差矩阵(Sample covariancematrix,SCM),为了保证稳健性,对干扰抑制能力有极大地牺牲;第2)类方法能够有效保证稳健性,但当存在模型失配时,其干扰抑制能力会下降。这两类算法都基于阵元接收数据,其输出信干噪比与输入信噪比的差值,不会超过阵增益。
发明内容
为了避免现有自适应波束形成的不足之处,本发明提出一种基于虚拟线列阵的协方差矩阵重构自适应波束形成方法,在所有入射信号方向未知的前提下,利用阵列接收数据和线性预测技术,拟合得到虚拟阵元的接收数据,从而提高自适应波束形成的自由度,最终提高自适应波束形成器对干扰和噪声的抑制能力。
为达到上述目的,本发明采用的技术方案为:一种协方差矩阵重构自适应波束形成方法,具体过程如下:
获取均匀线列阵接收时域数据,分别计算前向拟合系数和后向拟合系数,并采用递推的方式依次对虚拟阵元所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;
利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取较大特征值的方式剔除重构误差,最终引入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;
采用与计算改善干扰加噪声协方差矩阵的同样的方式,重构虚拟阵列对应的期望信号协方差矩阵,并提取出期望信号协方差矩阵最大特征值对应的特征向量,作为期望信号的导向矢量,根据所述导向矢量计算虚拟阵列对应的波束形成器加权向量,进而根据所述加权向量和虚拟线列阵的接收数据矩阵得到波束输出。
获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,具体过程如下:
均匀线列阵接收来自空间的信号和噪声的时域快拍数据;
计算均匀线列阵中除第1个振元以外其余阵元数据所对应的前向样本协方差矩阵,再计算前向拟合的拟合向量;
根据前向拟合的拟合向量计算得到前向拟合系数;
计算均匀线列阵中除最后一个振元以外其余振元数据所对应的前向样本协方差矩阵,再计算前向拟合的拟合向量;
根据所述拟合向量计算得到后向拟合系数。
按照递推的方式依次对虚拟阵元位置处所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据过程如下:
取虚拟振列中除最后一个振元以外其他振元数据,计算0号振元的前向拟合数据,依次计算到第n个前向拟合数据;
取虚拟振列中除第一个振元以外其他振元数据,计算后向编号最大振元的数据,依次分别计算到第n个后向拟合数据。
利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取较大特征值的方式剔除重构误差,最终添加入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵和重构虚拟阵列对应的期望信号协方差矩阵,具体过程如下:
将前向拟合数据、原始阵元接收数据以及后向拟合数据组合在一起,得到虚拟线列阵的接收数据矩阵;
根据所述接收数据矩阵计算虚拟线列阵的样协方差矩阵;
根据虚拟线列阵的阵元编号,得到对应的角度的导向矢量,
基于所述导向矢量计算各个角度对应的对角加载Capon空间谱;
根据所述Capon空间谱,划分期望信号扇区,然后划定除期望扇区以外的干扰加噪声扇区;
在所述干扰加噪声扇区内初步计算干扰加噪声协方差矩阵;
对所述干扰加噪声协方差矩阵进行特征分解;
根据特征分解结果得到修正的干扰加噪声协方差矩阵;
在期望信号扇区内重构期望信号协方差矩阵。
获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并按照递推的方式依次对虚拟阵元所接收到的数据进行拟合,预测各个虚拟阵元的快拍数据,具体包括以下步骤:
M个换能器组成的均匀线列阵接收来自空间的信号和噪声的时域快拍数据x(k)=[x1(k),x2(k),…,xM(k)]T,其中信号与噪声不相关,xm(k)表示第m个阵元接收到的时域数据,计算第2至M个阵元数据xF(k)=[x2(k),…,xM(k)]T所对应的前向样本协方差矩阵
Figure BDA0003119731180000041
Figure BDA0003119731180000042
计算前向拟合的拟合向量dF
Figure BDA0003119731180000043
最后得到前向拟合系数αF
Figure BDA0003119731180000044
计算第1至M-1个阵元数据XB(k)=[x1(k),…,xM-1(k)]T所对应的前向样本协方差矩阵,
Figure BDA0003119731180000045
再计算前向拟合的拟合向量dB
Figure BDA0003119731180000051
最后得到后向拟合系数β:
Figure BDA0003119731180000052
取第1至M-1个阵元数据xtemp(k)=[x1(k),…,xM-1(k)]T,计算0号阵元的前向拟合数据,
Figure BDA0003119731180000053
按照此规律,依次分别计算到第n个前向拟合数据:
Figure BDA0003119731180000054
取第2至M个阵元数据xtemp(k)=[x2(k),…,xM(k)]T,计算后向M+1号阵元的数据,
Figure BDA0003119731180000055
按照此规律,依次分别计算到第n个后向拟合数据:
Figure BDA0003119731180000056
将前向拟合数据、原始阵元接收数据以及后向拟合数据组合在一起,得到虚拟线列阵的接收数据矩阵,
xLP(k)=[x1-N(k),…,x1(k),…,xM(k),…,xM+N(k)]T
接着,计算采样协方差矩阵,
Figure BDA0003119731180000061
利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,最终结合理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;具体如下:
根据虚拟线列阵的阵元编号,得到对应的θ角度的导向矢量aLP(θ),计算各个角度对应的对角加载Capon空间谱p(θ),
Figure BDA0003119731180000062
其中I是单位对角矩阵,
根据Capon空间谱,划分期望信号扇区
Figure BDA0003119731180000063
其中B为设置的区域宽度,然后划定除期望扇区以外的干扰加噪声扇区
Figure BDA0003119731180000064
为Θ的补集,在干扰加噪声扇区内初步计算干扰加噪声协方差矩阵
Figure BDA0003119731180000065
Figure BDA0003119731180000066
其中θl
Figure BDA0003119731180000067
内的第l个离散角度,接着进行特征分解
Figure BDA0003119731180000068
Figure BDA0003119731180000069
其中:λm(m=1,…,M)为按照从大到小顺序排列的特征值,vm是对应的特征向量;Dint表示包含所有较大特征值的对角阵,Vint包含了较大特征值对应的特征向量,Re表示由较小特征值和特征向量组成的误差协方差矩阵,所述较大特征值是根据
Figure BDA00031197311800000610
估计干扰信号个数为Q,加上1个信号数为Q+1,Q+1是使得上式成立的最小正整数;由大至小排列的Q个干扰信号个数加上1个信号数Q+1组成较大特征值;此时,干扰加噪声协方差矩阵修正为:
Figure BDA0003119731180000071
其中
Figure BDA0003119731180000072
是估计的噪声功率,为最小特征值,
在期望信号扇区内,重构期望信号协方差矩阵
Figure BDA0003119731180000073
Figure BDA0003119731180000074
其中θl是Θ内的第l个离散角度,估计期望信号的导向矢量为,
Figure BDA0003119731180000075
其中,
Figure BDA0003119731180000076
为矩阵
Figure BDA0003119731180000077
的最大特征值对应的特征向量;
计算虚拟阵列对应的波束形成器加权向量wLP
Figure BDA0003119731180000078
根据该加权向量和虚拟线列阵的接收数据矩阵计算波束输出y(k):
Figure BDA0003119731180000079
一种协方差矩阵重构自适应波束形成系统,包括快拍数据获取模块、协方差矩阵计算模块、
快拍数据获取模块用于获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并采用递推的方式依次对虚拟阵元所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;
协方差矩阵计算模块用于根据包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,引入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;
输出波束计算模块用于采用和计算改善干扰加噪声协方差矩阵同样的方式重构虚拟阵列对应的期望信号协方差矩阵,并提取出最大特征值对应的特征向量作为期望信号的导向矢量,根据所述导向矢量计算虚拟阵列对应的波束形成器加权向量,进而根据所述加权向量和虚拟线列阵的接收数据矩阵得到波束输出。
一种计算机设备,包括处理器以及存储器,存储器用于存储计算机可执行程序,处理器从存储器中读取部分或全部所述计算机可执行程序并执行,处理器执行部分或全部计算可执行程序时能实现本发明所述协方差矩阵重构自适应波束形成方法。
一种计算机可读存储介质,计算机可读存储介质中存储有计算机程序,所述计算机程序被处理器执行时,能实现本发明所述的协方差矩阵重构自适应波束形成方法。
与现有技术相比,本发明至少具有以下有益效果:本发明将线性预测与协方差矩阵重构技术相结合,首先以线性预测拟合出前后向的虚拟阵元的接收数据,从而提升阵列的通道数。线性预测与协方差矩阵重构技术结合后,不仅能够提高自适应波束形成器的阵增益,还能提高自适应波束形成器的复杂度,一方面使得波束形成器在阵列误差较小时能够在干扰方向形成更深的凹陷,另一方面也能够在阵列误差较大时通过较低的旁瓣抑制干扰,最终提高自适应波束形成器对干扰和噪声的抑制能力。
附图说明
图1为所提虚拟阵列波束形成方法的输出SINR随输入SNR的变化情况。
图2为所提虚拟阵列波束形成方法的输出SINR随输入快拍数的变化情况。
图3为所提虚拟阵列波束形成方法的输出SINR随前后向拟合数N的变化情况。
具体实施方式
下面结合附图对本发明作进一步的描述。
本发明利用已有的均匀线列阵接收数据,通过线性预测的方式拟合出虚拟阵元,从而构成虚拟阵列,从而提升阵列的通道数,提高自适应波束形成的性能。
实施例
仿真1:水下声速为1500m/s,考虑一个10元线列阵,假设的阵元间距为0.5m,期望信号假设位于15°,信噪比从-10dB到20dB变化,干扰假设位于-25°和35°,干噪比为20dB,噪声为高斯均匀白噪声。期望信号和干扰信号均为中心频率为1500Hz的窄带不相干信号。前后向均拟合N=5个虚拟阵元,快拍数为K=30个。
为了阐述本方法的统计性能,共仿真200次,最终得到输出新干噪比随输入信噪比变化的统计结果。对于每次仿真,重复本方法的以下操作步骤以得到加权向量:
步骤1:10个水听器组成的均匀线列阵接收来自空间的信号和噪声的时域快拍数据x(k)=[x1(k),x2(k),…,x10(k)]T,其中信号与噪声不相关,xm(k)表示第m个阵元接收到的时域数据。虚拟振元和实际振元共同组成虚拟线列阵,计算第2至10个阵元数据xF(k)=[x2(k),…,x10(k)]T所对应的前向样本协方差矩阵,
Figure BDA0003119731180000091
再计算前向拟合的拟合向量,
Figure BDA0003119731180000092
最后得到前向拟合系数
Figure BDA0003119731180000093
接着,计算第1至9个阵元数据XB(k)=[x1(k),…,x9(k)]T所对应的前向样本协方差矩阵,
Figure BDA0003119731180000094
再计算前向拟合的拟合向量,
Figure BDA0003119731180000095
最后得到后向拟合系数
Figure BDA0003119731180000101
步骤2:取第1至9个阵元数据xtemp(k)=[x1(k),…,x9(k)]T,计算0号阵元的前向拟合数据,
Figure BDA0003119731180000102
按照此规律,依次分别计算到第5个前向拟合数据:
Figure BDA0003119731180000103
取第2至10个阵元数据xtemp(k)=[x2(k),…,x10(k)]T,计算后向11号阵元的数据,
Figure BDA0003119731180000104
按照此规律,依次分别计算到第5个后向拟合数据:
Figure BDA0003119731180000105
步骤3:已经前向和后向均拟合了5个数据,将前向拟合数据、原始阵元接收数据以及后向拟合数据组合在一起,得到虚拟线列阵的接收数据矩阵,
xLP(k)=[x-4(k),…,x1(k),…,x10(k),…,x15(k)]T
接着,计算采样协方差矩阵,
Figure BDA0003119731180000106
设定离散角度间隔为1°,那么[-90°,90°]内共181个离散角度。根据虚拟线列阵的阵元编号,得到对应的θ角度的导向矢量aLP(θ),,其包含20个元素,计算各个角度对应的对角加载Capon空间谱,
Figure BDA0003119731180000111
其中I是单位对角矩阵。
步骤4:根据Capon空间谱估计期望信号角度为
Figure BDA0003119731180000112
设置的区域宽度为B=10,那么划分期望信号扇区Θ=[10°,20°],然后划定除期望扇区以外的干扰加噪声扇区
Figure BDA0003119731180000113
共有L=181-11=170个离散角度。在干扰加噪声扇区内初步计算干扰加噪声协方差矩阵
Figure BDA0003119731180000114
其中θl
Figure BDA0003119731180000115
内的第l个离散角度。接着进行特征分解,
Figure BDA0003119731180000116
其中:λm(m=1,…,10)为按照从大到小顺序排列的特征值,vm是对应的特征向量;Dint表示包含所有较大特征值的对角阵,Vint包含了较大特征值对应的特征向量,Re表示由较小特征值和特征向量组成的误差协方差矩阵,所述较大特征值是根据
Figure BDA0003119731180000117
估计干扰信号个数为Q,加上1个信号数为Q+1,Q+1是使得上式成立的最小正整数;由大至小排列的Q个干扰信号个数加上1个信号数Q+1组成较大特征值;此时,干扰加噪声协方差矩阵可以修正为,
Figure BDA0003119731180000118
其中
Figure BDA0003119731180000119
是估计的噪声功率,通常可以认为是最小特征值;
步骤5:在期望信号扇区内共有11个离散角度,按照步骤4所述方式重构期望信号协方差矩阵,从而得到,
Figure BDA00031197311800001110
其中θl是Θ内的第l个离散角度,估计期望信号的导向矢量为,
Figure BDA0003119731180000121
所述
Figure BDA0003119731180000122
为矩阵
Figure BDA0003119731180000123
的最大特征值对应的特征向量;
步骤6:计算虚拟阵列对应的波束形成器加权向量,
Figure BDA0003119731180000124
根据该加权向量和虚拟线列阵的接收数据矩阵计算波束输出,
Figure BDA0003119731180000125
图1给出了当快拍数为30时,所提出的LP-REB与REB和RCB波束形成方法的对比,可以看出,RCB在高信噪比下性能较差,而REB由于重构误差的存在,失踪低于最优上限(optimal)。而LP-REB由于使用了虚拟阵列,增加了总的阵元数(10个真实阵元,10个虚拟阵元),能够有效地提高输出信干噪比,甚至优于10元阵列时的最优上限。
仿真2:输入信噪比在保持在10dB,让快拍数K从10到100变化,间隔为10.其余的仿真条件与仿真1相同,进行200次蒙特卡洛实验,并得到输出信干噪比的平均值随快拍数的变化图。步骤与仿真1对应步骤相同。
图2给出了各个波束形成器的输出信干噪比随快拍数变化的情况,可以看出,RCB算法由于样本协方差矩阵在低样本数时存在估计误差,所以性能逐渐增强,但依然远远低于20dB,REB和LP-REB由于采用了重构的协方差矩阵,所以对快拍数不敏感。其中,LP-REB由于采用了更多的虚拟阵元数据,其输出信干燥比远远高于REB。
仿真3:输入信噪比固定为10dB,快拍数固定为30,让前后向拟合数量N从0到9以1为间隔进行变化,其余的仿真条件与仿真1相同,进行200次蒙特卡洛实验,并得到输出信干噪比的平均值随前后向拟合数量N的变化图。
图3给出了LP-REB的SINR随前后向拟合数量的变化,可以看出,随着前后向拟合数的增大,虚拟阵元数也增大,从而有效地提高了虚拟阵列的增益,最终提高了SINR。
本发明提供一种协方差矩阵重构自适应波束形成系统,包括快拍数据获取模块、协方差矩阵计算模块、
快拍数据获取模块用于获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并采用递推的方式依次对虚拟阵元所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;
协方差矩阵计算模块用于根据包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,引入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;
输出波束计算模块用于采用和计算改善干扰加噪声协方差矩阵同样的方式重构虚拟阵列对应的期望信号协方差矩阵,并提取出最大特征值对应的特征向量作为期望信号的导向矢量,根据所述导向矢量计算虚拟阵列对应的波束形成器加权向量,进而根据所述加权向量和虚拟线列阵的接收数据矩阵得到波束输出。
本发明还可以提供一种计算机设备,包括处理器以及存储器,存储器用于存储计算机可执行程序,处理器从存储器中读取部分或全部所述计算机可执行程序并执行,处理器执行部分或全部计算可执行程序时能实现本发明所述协方差矩阵重构自适应波束形成方法。
另一方面,本发明提供一种计算机可读存储介质,计算机可读存储介质中存储有计算机程序,所述计算机程序被处理器执行时,能实现本发明所述的协方差矩阵重构自适应波束形成方法。
所述计算机设备可以采用车载计算机、笔记本电脑、平板电脑、桌面型计算机、手机或工作站。
处理器可以是中央处理器(CPU)、数字信号处理器(DSP)、专用集成电路(ASIC)或现成可编程门阵列(FPGA)。
对于本发明所述存储器,可以是车载计算机、笔记本电脑、平板电脑、桌面型计算机、手机或工作站的内部存储单元,如内存、硬盘;也可以采用外部存储单元,如移动硬盘、闪存卡。
计算机可读存储介质可以包括计算机存储介质和通信介质。计算机存储介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机可读存储介质可以包括:只读存储器(ROM,Read Only Memory)、随机存取记忆体(RAM,Random Access Memory)、固态硬盘(SSD,Solid State Drives)或光盘等。其中,随机存取记忆体可以包括电阻式随机存取记忆体(ReRAM,Resistance Random Access Memory)和动态随机存取存储器(DRAM,Dynamic Random Access Memory)。

Claims (9)

1.一种协方差矩阵重构自适应波束形成方法,其特征在于,具体过程如下:
获取均匀线列阵接收时域数据,分别计算前向拟合系数和后向拟合系数,并采用递推的方式依次对虚拟阵元所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;
利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取较大特征值的方式剔除重构误差,最终引入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;
采用与计算改善干扰加噪声协方差矩阵的同样的方式,重构虚拟阵列对应的期望信号协方差矩阵,并提取出期望信号协方差矩阵最大特征值对应的特征向量,作为期望信号的导向矢量,根据所述导向矢量计算虚拟阵列对应的波束形成器加权向量,进而根据所述加权向量和虚拟线列阵的接收数据矩阵得到波束输出。
2.根据权利要求1所述的协方差矩阵重构自适应波束形成方法,其特征在于,获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,具体过程如下:
均匀线列阵接收来自空间的信号和噪声的时域快拍数据;
计算均匀线列阵中除第1个振元以外其余阵元数据所对应的前向样本协方差矩阵,再计算前向拟合的拟合向量;
根据前向拟合的拟合向量计算得到前向拟合系数;
计算均匀线列阵中除最后一个振元以外其余振元数据所对应的前向样本协方差矩阵,再计算前向拟合的拟合向量;
根据所述拟合向量计算得到后向拟合系数。
3.根据权利要求1所述的协方差矩阵重构自适应波束形成方法,其特征在于,按照递推的方式依次对虚拟阵元位置处所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据过程如下:
取虚拟振列中除最后一个振元以外其他振元数据,计算0号振元的前向拟合数据,依次计算到第n个前向拟合数据;
取虚拟振列中除第一个振元以外其他振元数据,计算后向编号最大振元的数据,依次分别计算到第n个后向拟合数据。
4.根据权利要求1所述的协方差矩阵重构自适应波束形成方法,其特征在于,利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取较大特征值的方式剔除重构误差,最终添加入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵和重构虚拟阵列对应的期望信号协方差矩阵,具体过程如下:
将前向拟合数据、原始阵元接收数据以及后向拟合数据组合在一起,得到虚拟线列阵的接收数据矩阵;
根据所述接收数据矩阵计算虚拟线列阵的样协方差矩阵;
根据虚拟线列阵的阵元编号,得到对应的角度的导向矢量,
基于所述导向矢量计算各个角度对应的对角加载Capon空间谱;
根据所述Capon空间谱,划分期望信号扇区,然后划定除期望扇区以外的干扰加噪声扇区;
在所述干扰加噪声扇区内初步计算干扰加噪声协方差矩阵;
对所述干扰加噪声协方差矩阵进行特征分解;
根据特征分解结果得到修正的干扰加噪声协方差矩阵;
在期望信号扇区内重构期望信号协方差矩阵。
5.根据权利要求1所述的协方差矩阵重构自适应波束形成方法,其特征在于,获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并按照递推的方式依次对虚拟阵元所接收到的数据进行拟合,预测各个虚拟阵元的快拍数据,具体包括以下步骤:
M个换能器组成的均匀线列阵接收来自空间的信号和噪声的时域快拍数据x(k)=[x1(k),x2(k),…,xM(k)]T,其中信号与噪声不相关,xm(k)表示第m个阵元接收到的时域数据,计算第2至M个阵元数据xF(k)=[x2(k),…,xM(k)]T所对应的前向样本协方差矩阵
Figure FDA0003119731170000038
Figure FDA0003119731170000031
计算前向拟合的拟合向量dF
Figure FDA0003119731170000032
最后得到前向拟合系数αF
Figure FDA0003119731170000033
计算第1至M-1个阵元数据XB(k)=[x1(k),…,xM-1(k)]T所对应的前向样本协方差矩阵,
Figure FDA0003119731170000034
再计算前向拟合的拟合向量dB
Figure FDA0003119731170000035
最后得到后向拟合系数β:
Figure FDA0003119731170000036
取第1至M-1个阵元数据xtemp(k)=[x1(k),…,xM-1(k)]T,计算0号阵元的前向拟合数据,
Figure FDA0003119731170000037
按照此规律,依次分别计算到第n个前向拟合数据:
Figure FDA0003119731170000041
取第2至M个阵元数据xtemp(k)=[x2(k),…,xM(k)]T,计算后向M+1号阵元的数据,
Figure FDA0003119731170000042
按照此规律,依次分别计算到第n个后向拟合数据:
Figure FDA0003119731170000043
将前向拟合数据、原始阵元接收数据以及后向拟合数据组合在一起,得到虚拟线列阵的接收数据矩阵,
xLP(k)=[x1-N(k),…,x1(k),…,xM(k),…,xM+N(k)]T
接着,计算采样协方差矩阵,
Figure FDA0003119731170000044
6.根据权利要求1所述的协方差矩阵重构自适应波束形成方法,其特征在于,利用包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,最终结合理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;具体如下:
根据虚拟线列阵的阵元编号,得到对应的θ角度的导向矢量aLP(θ),计算各个角度对应的对角加载Capon空间谱p(θ),
Figure FDA0003119731170000051
其中I是单位对角矩阵,
根据Capon空间谱,划分期望信号扇区
Figure FDA0003119731170000052
其中B为设置的区域宽度,然后划定除期望扇区以外的干扰加噪声扇区
Figure FDA0003119731170000053
为Θ的补集,在干扰加噪声扇区内初步计算干扰加噪声协方差矩阵
Figure FDA0003119731170000054
Figure FDA0003119731170000055
其中θl
Figure FDA0003119731170000056
内的第l个离散角度,接着进行特征分解
Figure FDA0003119731170000057
Figure FDA0003119731170000058
其中:λm(m=1,…,M)为按照从大到小顺序排列的特征值,vm是对应的特征向量;Dint表示包含所有较大特征值的对角阵,Vint包含了较大特征值对应的特征向量,Re表示由较小特征值和特征向量组成的误差协方差矩阵,所述较大特征值是根据
Figure FDA0003119731170000059
估计干扰信号个数为Q,加上1个信号数为Q+1,Q+1是使得上式成立的最小正整数;由大至小排列的Q个干扰信号个数加上1个信号数Q+1组成较大特征值;此时,干扰加噪声协方差矩阵修正为:
Figure FDA00031197311700000510
其中
Figure FDA00031197311700000511
是估计的噪声功率,为最小特征值,
在期望信号扇区内,重构期望信号协方差矩阵
Figure FDA00031197311700000512
Figure FDA00031197311700000513
其中θl是Θ内的第l个离散角度,估计期望信号的导向矢量为,
Figure FDA00031197311700000514
其中,
Figure FDA0003119731170000061
为矩阵
Figure FDA0003119731170000062
的最大特征值对应的特征向量;
计算虚拟阵列对应的波束形成器加权向量wLP
Figure FDA0003119731170000063
根据该加权向量和虚拟线列阵的接收数据矩阵计算波束输出y(k):
Figure FDA0003119731170000064
7.一种协方差矩阵重构自适应波束形成系统,其特征在于,包括快拍数据获取模块、协方差矩阵计算模块、
快拍数据获取模块用于获取均匀线列阵接收时域数据,分别计算前向和后向的线性预测系数,并采用递推的方式依次对虚拟阵元所接收到的数据进行拟合,从而预测出各个虚拟阵元的快拍数据;
协方差矩阵计算模块用于根据包含实际阵元和虚拟阵元的快拍数据计算更高维度的干扰加噪声协方差矩阵,并采用特征分解和提取大特征值的方式剔除重构误差,引入理想的噪声协方差矩阵,得到虚拟阵列对应的改善干扰加噪声协方差矩阵;
输出波束计算模块用于采用和计算改善干扰加噪声协方差矩阵同样的方式重构虚拟阵列对应的期望信号协方差矩阵,并提取出最大特征值对应的特征向量作为期望信号的导向矢量,根据所述导向矢量计算虚拟阵列对应的波束形成器加权向量,进而根据所述加权向量和虚拟线列阵的接收数据矩阵得到波束输出。
8.一种计算机设备,其特征在于,包括处理器以及存储器,存储器用于存储计算机可执行程序,处理器从存储器中读取部分或全部所述计算机可执行程序并执行,处理器执行部分或全部计算可执行程序时能实现权利要求1~6中任一项所述协方差矩阵重构自适应波束形成方法。
9.一种计算机可读存储介质,其特征在于,计算机可读存储介质中存储有计算机程序,所述计算机程序被处理器执行时,能实现如权利要求1~6中任一项所述的协方差矩阵重构自适应波束形成方法。
CN202110673467.2A 2021-06-17 2021-06-17 一种协方差矩阵重构自适应波束形成方法及系统 Active CN113422629B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110673467.2A CN113422629B (zh) 2021-06-17 2021-06-17 一种协方差矩阵重构自适应波束形成方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110673467.2A CN113422629B (zh) 2021-06-17 2021-06-17 一种协方差矩阵重构自适应波束形成方法及系统

Publications (2)

Publication Number Publication Date
CN113422629A true CN113422629A (zh) 2021-09-21
CN113422629B CN113422629B (zh) 2023-02-21

Family

ID=77788895

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110673467.2A Active CN113422629B (zh) 2021-06-17 2021-06-17 一种协方差矩阵重构自适应波束形成方法及系统

Country Status (1)

Country Link
CN (1) CN113422629B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114779176A (zh) * 2022-04-19 2022-07-22 四川大学 一种低复杂度的鲁棒自适应波束形成方法及装置

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050259006A1 (en) * 2004-05-24 2005-11-24 Samsung Electronics Co., Ltd. Beam forming apparatus and method using interference power estimation in an array antenna system
US20050271016A1 (en) * 2004-05-07 2005-12-08 Byoung-Yun Kim Beam forming apparatus and method for an array antenna system
US20160014393A1 (en) * 2014-07-09 2016-01-14 Massachusetts Institute Of Technology Methods and Apparatus for Virtual Sensor Array
CN105306117A (zh) * 2015-10-15 2016-02-03 哈尔滨工程大学 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法
US20160157828A1 (en) * 2014-06-05 2016-06-09 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
CN106405487A (zh) * 2016-10-17 2017-02-15 重庆邮电大学 一种通用基于扩展esprit技术的空间谱估计方法
CN106569181A (zh) * 2016-11-15 2017-04-19 大连大学 基于协方差矩阵重构稳健Capon波束形成的算法
CN109600152A (zh) * 2018-12-17 2019-04-09 西北工业大学 一种基于子空间基变换的自适应波束形成方法
US10574320B1 (en) * 2019-07-15 2020-02-25 Southwest Research Institute Matrix methods to speed processing for MVDR beamforming
CN111049556A (zh) * 2019-12-27 2020-04-21 中国科学技术大学 一种基于干扰协方差矩阵重构的互素阵稳健自适应波束形成方法
CN111123192A (zh) * 2019-11-29 2020-05-08 湖北工业大学 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN111651718A (zh) * 2020-05-25 2020-09-11 中国科学技术大学 针对高斯噪声中非高斯信号的稳健自适应波束形成方法
WO2020212569A1 (en) * 2019-04-17 2020-10-22 Université Du Luxembourg Method and device for beamforming in a mimo radar system

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050271016A1 (en) * 2004-05-07 2005-12-08 Byoung-Yun Kim Beam forming apparatus and method for an array antenna system
US20050259006A1 (en) * 2004-05-24 2005-11-24 Samsung Electronics Co., Ltd. Beam forming apparatus and method using interference power estimation in an array antenna system
US20160157828A1 (en) * 2014-06-05 2016-06-09 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
US20160014393A1 (en) * 2014-07-09 2016-01-14 Massachusetts Institute Of Technology Methods and Apparatus for Virtual Sensor Array
CN105306117A (zh) * 2015-10-15 2016-02-03 哈尔滨工程大学 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法
CN106405487A (zh) * 2016-10-17 2017-02-15 重庆邮电大学 一种通用基于扩展esprit技术的空间谱估计方法
CN106569181A (zh) * 2016-11-15 2017-04-19 大连大学 基于协方差矩阵重构稳健Capon波束形成的算法
CN109600152A (zh) * 2018-12-17 2019-04-09 西北工业大学 一种基于子空间基变换的自适应波束形成方法
WO2020212569A1 (en) * 2019-04-17 2020-10-22 Université Du Luxembourg Method and device for beamforming in a mimo radar system
US10574320B1 (en) * 2019-07-15 2020-02-25 Southwest Research Institute Matrix methods to speed processing for MVDR beamforming
CN111123192A (zh) * 2019-11-29 2020-05-08 湖北工业大学 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN111049556A (zh) * 2019-12-27 2020-04-21 中国科学技术大学 一种基于干扰协方差矩阵重构的互素阵稳健自适应波束形成方法
CN111651718A (zh) * 2020-05-25 2020-09-11 中国科学技术大学 针对高斯噪声中非高斯信号的稳健自适应波束形成方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
J. LI等: "Versatile robust Capon beamforming: theory and applications", 《PROCESSING WORKSHOP PROCEEDINGS, 2004 SENSOR ARRAY AND MULTICHANNEL SIGNAL》 *
鲁欢等: "基于IAA的协方差矩阵重构稳健波束形成方法", 《太赫兹科学与电子信息学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114779176A (zh) * 2022-04-19 2022-07-22 四川大学 一种低复杂度的鲁棒自适应波束形成方法及装置
CN114779176B (zh) * 2022-04-19 2023-05-05 四川大学 一种低复杂度的鲁棒自适应波束形成方法及装置

Also Published As

Publication number Publication date
CN113422629B (zh) 2023-02-21

Similar Documents

Publication Publication Date Title
CN108809398B (zh) 一种基于信源个数约束的稳健自适应波束形成方法
CN108181507B (zh) 一种稳健自适应波束形成方法
CN111049556B (zh) 一种基于干扰协方差矩阵重构的互素阵稳健自适应波束形成方法
CN110045321B (zh) 基于稀疏和低秩恢复的稳健doa估计方法
CN108663668B (zh) 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法
CN109639333B (zh) 一种基于有效重构协方差矩阵的波束形成方法
CN105302936A (zh) 基于相关计算和协方差矩阵重构的自适应波束形成方法
CN113422630B (zh) 一种自适应聚焦宽带波束形成方法及系统
CN110988854A (zh) 基于交替方向乘子法的鲁棒自适应波束形成算法
CN108802705A (zh) 一种基于稀疏的空时自适应处理方法及系统
CN113422629B (zh) 一种协方差矩阵重构自适应波束形成方法及系统
CN115102597A (zh) 一种稳健的自适应波束形成方法
CN113673419B (zh) 适用于强干扰环境的波束域快速稀疏贝叶斯方位估计方法
CN112946564B (zh) 基于dnn的波束空间的doa估计方法、装置及计算机存储介质
CN111651718B (zh) 针对高斯噪声中非高斯信号的稳健自适应波束形成方法
CN111698011B (zh) 一种对称噪声中非对称信号的稳健自适应波束形成方法
CN111431574B (zh) 一种多重稳健自适应波束形成方法
CN114047481A (zh) 一种基于子空间正交性的稳健自适应波束形成方法
CN115808659A (zh) 基于低复杂度不确定集积分的稳健波束形成方法及系统
CN114779176B (zh) 一种低复杂度的鲁棒自适应波束形成方法及装置
CN114996653A (zh) 一种基于原子范数最小化的二维鲁棒自适应波束形成方法
CN114236471A (zh) 一种相关干扰源下的稳健自适应波束形成方法
CN114563764A (zh) 一种自适应波束生成方法及系统
CN114355293A (zh) 一种自适应波束生成方法及系统
CN107846241B (zh) 脉冲噪声环境下的波束形成方法、存储介质及波束形成器

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant