一种复合高斯背景下的知识辅助STAP检测方法与流程

技术特征：

1.一种复合高斯背景下的知识辅助STAP检测方法，其特征在于包括下述步骤：

(1)在复合高斯背景下，根据杂波协方差矩阵的先验分布，利用MAP准则推导出协方差矩阵其中，表示协方差矩阵R的均值，v为协方差矩阵所满足复逆Wishart分布中表示R和间距离的参数，q_k和β_k分别表示复合高斯分布中纹理分量所满足逆伽玛分布的形状参数和尺度参数，k＝1,...K，表示第k个距离单元的数据，K表示估计协方差矩阵所选取的训练样本数，N表示距离单元数据的维数，(·)^H表示矩阵共轭转置，(·)^-1表示矩阵求逆运算；

(2)根据阵列雷达的STAP信号模型，建立杂波协方差矩阵结构的先验模型，得到协方差矩阵估计的均值其中，R_c表示根据先验信号模型建立的杂波协方差矩阵结构，α_p、θ_p和f_p分别表示阵列雷达从杂波块接收到杂波信号的复幅度、到达角和多普勒频移，N_c表示一个杂波距离环上划分的杂波块的数目，表示一个杂波块所对应到达角为θ_p和多普勒频移为f_p的空时导向矢量，τ表示服从复合高斯分布的训练样本数据的纹理分量值；

(3)由NSCM方法计算样本协方差矩阵作为杂波协方差矩阵的初始估计，利用IAA方法通过杂波协方差矩阵计算出协方差矩阵均值所需的杂波块功率，将其代入协方差矩阵表达式中计算出协方差矩阵，此过程迭代进行，直至收敛，具体步骤如下：

a1)初始化迭代次数iter＝1，R_NSCM表示归一化样本协方差矩阵，

a2)当iter＝t时，计算其中，

${\mathrm{\τ}}^{\left(t\right)}=\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\frac{{\mathrm{\β}}_k+z_k^H{\hat{R}}_{\left(t\right)}^{-1}{z}_k}{q_k+N+1},{\mathrm{\σ}}_p^{\left(t\right)}=|\frac{e^H\left({\mathrm{\θ}}_p,f_p\right){\hat{R}}_{\left(t\right)}^{-1}z}{e^H\left({\mathrm{\θ}}_p,f_p\right){\hat{R}}_{\left(t\right)}^{-1}e\left({\mathrm{\θ}}_p,f_p\right)}{|}^2,{\stackrel{\‾}{R}}_{\left(t\right)}=\underset{p=1}{\overset{Nc}{\mathrm{\Σ}}}{\mathrm{\σ}}_p^{\left(t\right)}e\left({\mathrm{\θ}}_p,f_p\right)e^H\left({\mathrm{\θ}}_p,f_p\right)/{\mathrm{\τ}}^{\left(t\right)},$σ_p＝|α_p|²表示一个杂波距离环上划分的杂波块的功率值；

a3)当与上次迭代估计的变化量小于设定极小值ε＝0.01时，迭代停止，得到否则令iter＝t+1，返回步骤a2)继续迭代；

(4)将得到的杂波协方差矩阵用于NAMF检测器中，NAMF检测器的检测统计量将待检测单元数据z和协方差矩阵估计值代入到NAMF检测器中得到检测统计量，将其与阈值λ进行比较，如果检测统计量大于阈值，则判定为有目标，否则判定没有目标。