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

技术特征：

1.一种通带响应误差加权阻带响应约束空域矩阵滤波器设计方法，其特征在于以下步骤：

步骤1：k＝0，将探测的空域离散化，获得V_P，V_S，Θ_P，Θ_S；通过ε，计算初始最优空域矩阵滤波器令w₀(θ_p)＝1；

步骤2：计算|E_k(θ_p)|＝|Ha(θ_p)-a(θ_p)|,p＝1,…,P；求|E_k(θ_p)|的局部极大值点，获取局部极大值的横坐标及其相应的纵坐标

步骤3：利用和两点间连线的延长线，计算在横坐标θ₁处的取值z₁，设置(θ₁,max(|E_k(θ₁)|,z₁))为包络加权起始点；

步骤4：利用和两点间连线的延长线，计算在横坐标θ_S处的取值z_S，设置(θ_P,max(z_P,|E_k(θ_P)|))为包络加权终点；

步骤5：计算(θ₁,max(|E_k(θ₁)|,z₁))、(θ_P,max(z_P,|E_k(θ_P)|))共Q+2个点之间的线段，并取α_k(θ_p)为θ_p在相应线段上的取值；

步骤6：计算下列各式

$\mathrm{\β}\left({\mathrm{\θ}}_p\right)={\mathrm{P\α}}_k\left({\mathrm{\θ}}_p\right)/\underset{p=1}{\overset{P}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left({\mathrm{\θ}}_p\right)$

w_k+1(θ_p)＝β_k(θ_p)w_k(θ_p)

R_k+1＝diag[w_k+1(θ₁),w_k+1(θ₂),…,w_k+1(θ_S)]

$C_P=\frac{1}{\mathrm{NP}}{V}_P{R}_{k+1}{V}_P^H,$$C_S=\frac{1}{\mathrm{NS}}{V}_S{V}_S^H$

$\mathrm{tr}\left[C_P{(C_P+{\widehat{\mathrm{\λ}}}_{k+1}{C}_S)}^{-1}{C}_S{(C_P+{\widehat{\mathrm{\λ}}}_{k+1}{C}_S)}^{-1}{C}_P\right]=\mathrm{\ϵ}$

${\hat{H}}_{k+1}=C_P{(C_P+{\widehat{\mathrm{\λ}}}_{k+1}{C}_S)}^{-1}$

判断H_k+1是否满足如下终止条件之一：

(a)k+1＝K；此时，此时迭代K次，算法终止；

(b)p＝1,…,P；迭代后，空域矩阵滤波器对通带上所有方位的实际响应误差值小于常数算法终止；

(c)p＝1,…,P；迭代后，空域矩阵滤波器对通带上所有方位的响应误差变化率都小于常数值算法终止；

步骤7：若迭代终止条件满足，则H_k+1即为最终的空域矩阵滤波器；否则，令k:＝k+1，重复步骤2～6。