[发明专利]一种基于CLEAN的稀疏阵列波达方向DOA估计方法

摘要

本发明涉及一种基于CLEAN的稀疏阵列波达方向DOA估计方法，包括以下步骤：利用最小描述长度(MDL)准则来计算目标个数；计算所有快拍数据形成波束的幅度平方累加结果；在计算结果中查找最大值对应的波束号；利用等信号测角法计算“精确波束号”；计算目标角度；判断所有信号是否均已检测；应用CLEAN思想清除检测到的目标响应，转至进行迭代。该方法能够很好的抑制强目标与高旁瓣作用结果对弱目标角度估计的影响，达到很高的测角精度，同时，该方法耗时很短且耗时对阵元数不敏感。

主权项

一种基于CLEAN的稀疏阵列波达方向DOA估计方法，其特征在于步骤如下：步骤1：利用最小描述长度MDL准则计算CLEAN的稀疏阵列信号中的目标个数；步骤2：T>0时，计算所有快拍数据形成波束的幅度平方累加结果否则结束计算过程；其中：y_s为第s个快拍数据x_s利用FFT形成的波束数据，y_s＝FFT(x_s)；L为快拍数；步骤3：在幅度平方累加结果P中查找最大值对应的波束号；步骤4：采用等信号测角法计算“精确波束号”p_t，$p_t=\frac{(p_1+p_2)}{2}+k\·\mathrm{sign}\·\min (\frac{|P\left(p_1\right)-P\left(p_2\right)|}{|P\left(p_1\right)+P\left(p_2\right)|},0.5)$其中，p₁为P中最大值对应的波束号，p₂为次大值的波束号，k为阵列对应的比幅测角斜率，sign为符号位；$p_2=\left\{\begin{array}{cc}{p}_1-1,& P(p_1-1)>P(p_1+1)\\ p_1+1,& P(p_1-1)\≤P(p_1+1)\end{array}\right.$所述$\mathrm{sign}=\left\{\begin{array}{cc}-1,& p_1\<p_2\\ 1,& p_1\≥p_2\end{array}\right.;$步骤5：计算目标角度，其中，λ为波长，d为阵元间距，M为FFT点数；步骤6：应用CLEAN思想清除检测到的目标响应，同时将T值减一，转至步骤2进行迭代。