1.基于协方差矩阵估计和子空间构造的波束形成方法,其特征在于,包括下列步骤:步骤1,将阵列划分为两个阵元数相等的子阵列,两个子阵列所对应的阵列接收数据矩阵分别定义为X1和X2,即阵列的阵列接收数据矩阵基于子阵列的阵列接收数据矩阵计算两个子阵列的采样自协方差矩阵 以及计算两个子阵列之间的采样互协方差矩阵基于采样自协方差矩阵 和采样互协方差矩阵 构成阵列的阵列接收数据矩阵X的采样协方差矩阵步骤2,对 和 分别进行协方差矩阵估计处理,得到估计的自协方差矩阵 和估计的互协方差矩阵步骤3,对 和 进行特征值分解,得到 的特征空间E11和特征值矩阵D11,得到 的特征空间E21,并构造矩阵G:H
步骤4,对GG进行特征值分解,并将特征值和其特征向量按升序方式排列,得到M个升序排列的特征值,记为λ1~λM,其对应的特征向量依次记为:v1,…,vM;其中,M表示阵元数;
步骤5,对采样协方差矩阵 进行特征值分解,并将得到的M个特征值降序排列,记为记为γ1~γM;
基于设置的门限η,从大到小依次将P个特征值γ1~γP相加并与M个特征值γ1~γM的和作比值,若该比值超过门限η,则基于P的值确定信源数;其中P的初始值为1,每与门限η比对一次,P自增1;
基于得到的信源数P,根据公式估算噪声功率步骤6,P个信源信号中,含有一个期望信号,其余的P‑1个信源信号被视为干扰信号;从v1,…,vM中选取P‑1个特征向量构造干扰子空间矩阵Φ,并基于干扰子空间矩阵Φ构造投影矩阵Σ:H H
1)若信干比大于0,则Φ=[vM‑P+1,…,vM‑1];若信干比小于0,则Φ=[vM‑P+2,…,vM];
H
2)Σ=ΦΦ ;
H
步骤7,将GG矩阵向干扰子空间投影,得到近似的干扰信号协方差矩阵:步骤8,基于干扰信号协方差矩阵 估计P‑1个干扰信号功率,其中第j个干扰信号功率为其中, 表示预估的第j个干扰信号导向矢量;
步骤9,基于估计的噪声功率 和P‑1个干扰信号功率 重构干扰加噪声协方差矩阵其中,IM表示M×M的单位矩阵;
步骤10,将重构的干扰加噪声协方差矩阵 和预估的期望信号导向矢量 带入最小方差无失真响应波束形成器中得到权值向量w:
2.如权利要求1所述的方法,其特征在于,协方差矩阵估计处理具体为:定义 和 分别表示待进行协方差矩阵估计处理的实部采样协方差矩阵和虚部采样协方差矩阵;
计算 和 的收缩系数ρr和ρi:
其中,Tr表示矩阵的迹;
计算收缩矩阵 和
分别计算实部和虚部的估计矩阵 和
将估计矩阵 和 相加,得到协方差矩阵估计处理结果: j为虚数单位;
其中,
1)对于采样自协方差矩阵 其对应的 和 分别为:
2)对于采样互协方差矩阵 其对应的 和 分别为:其中,N表示采样点数, 表示取实部, 表示取虚部。