阻带响应加权通带总体响应误差约束空域矩阵滤波器设计方法

摘要

本发明属于阵列信号处理技术领域，涉及到传感器阵列的数据处理，特别涉及到空域矩阵滤波器设计方法。其特征是通过对最优空域矩阵滤波器权值参数迭代的方式，获得恒定的左右阻带响应，且通带响应误差满足设定的约束。加权系数矩阵是由矩阵滤波器阻带响应的包络获得。左右阻带响应可以相等，也可通过调节左右阻带响应比例，设定左右阻带响应值的差值。

主权项

一种阻带响应加权通带总体响应误差约束空域矩阵滤波器设计方法，其特征在于以下步骤，步骤1：k＝0，将探测空域离散化，获得V_P，V_S，Θ_S1，Θ_S2；通过ξ，计算初始最优空域矩阵滤波器设置左阻带、右阻带响应比例γ(θ_s)_，令w₀(θ_s)＝1；步骤2：计算|E_k(θ_s)|＝|H_ka(θ_s)|,s＝1,…,S；求|E_k(θ_s)|的局部极大值点，获取局部极大值的横坐标及其相应的纵坐标步骤3：利用和两点间连线的延长线，计算在横坐标θ₁处的取值z₁，设置(θ₁,max(|E_k(θ₁)|,z₁))为包络加权起始点；步骤4：利用和两点间连线的延长线，计算在横坐标θ_S处的取值z_S，设置(θ_S,max(z_S,|E_k(θ_S)|))为包络加权终点；步骤5：计算(θ₁,max(|E_k(θ₁)|,z₁))、(θ_S,max(z_S,|E_k(θ_S)|))共Q+2个点之间的线段，并取α_k(θ_s)为θ_s在相应线段上的取值；步骤6：计算下列各式${\beta }_k\left({\theta }_s\right)=S{\alpha }_k\left({\theta }_s\right)/\underset{s=1}{\overset{S}{\Sigma }}{\alpha }_k\left({\theta }_s\right)$w_k+1(θ_s)＝β_k(θ_s)γ(θ_s)w_k(θ_s)R_k+1＝diag[w_k+1(θ₁),w_k+1(θ₂),…,w_k+1(θ_S)]$C_P=\frac{1}{\mathrm{NP}}{V}_P{V}_P^H,C_S=\frac{1}{\mathrm{NS}}{V}_S{R}_{k+1}{V}_S^H$$\mathrm{tr}[{\hat{\mu }}_{k+1}{C}_P{({\hat{\mu }}_{k+1}{C}_P+C_S)}^{-1}{C}_S{({\hat{\mu }}_{k+1}{C}_P+C_S)}^{-1}{C}_P+{({\hat{\mu }}_{k+1}{C}_P+C_S)}^{-1}{C}_P]=1-\xi$$H_{k+1}={\hat{\mu }}_{k+1}{C}_P{({\hat{\mu }}_{k+1}{C}_P+C_S)}^{-1}$β_k(θ_s)为第k次迭代的乘积因子，w_k+1(θ_s)为第k+1个加权系数，R_k+1为第k+1次的加权系数矩阵，为第k+1次迭代所对应的Lagrange成子，H_k+1为第k+1个空域矩阵滤波器；判断H_k+1是否满足如下终止条件之一：(a)k+1＝K；此时，此时迭代K次，算法终止；(b)迭代后，空域矩阵滤波器对阻带上所有方位的实际响应值小于常数算法终止；(c)迭代后，空域矩阵滤波器对阻带上所有方位的响应变化率都小于常数值算法终止；步骤7：若迭代终止条件满足，则H_k+1即为最终的空域矩阵滤波器；否则，令k:＝k+1，重复步骤2～6。