[发明专利]一种机载雷达非平稳杂波抑制方法

权利要求书

1.一种机载雷达非平稳杂波抑制方法，其特征在于：包括如下步骤：

步骤1取出雷达接收到的经脉冲压缩后第l号待检测距离单元的回波数据x_l，并利用常规的迭代自适应谱估计方法估计该距离单元的空时二维功率谱矩阵其中l＝1，...，L，L为需要进行目标检测的距离单元的个数；

步骤2利用估计得到的空时二维功率谱矩阵以及导向矢量矩阵A重构第l号距离单元的空时二维协方差矩阵并对进行修正，得到修正后的协方差矩阵

步骤3在保证期望目标信号增益不变的前提下，根据线性约束最小方差准则，利用修正后的协方差矩阵计算得到第l号距离单元、第k个多普勒通道的空时滤波器系数w_l，k，k＝1，...，K，K为多普勒通道的个数；

步骤4利用第k个多普勒通道的空时滤波器w_l，k抑制数据x_l中的杂波，得到第l号距离单元、第k个多普勒通道的输出z_l，k；

步骤5令k＝k+1，重复步骤(3)～(4)，直至所有K个多普勒通道处理完毕，输出第l号距离单元杂波抑制后的最终多普勒谱：z_l＝[z_l，1，z_l，2，…，z_l，k]^T，其中[·]^T表示转置；

步骤6令l＝l+1，重复步骤(1)～(5)，直至所有L个距离单元处理完毕，输出L个距离单元杂波抑制后的距离多普勒谱Z＝[z₁，z₂，…，z_L]^T。

2.根据权利要求1所述的一种机载雷达非平稳杂波抑制方法，其中步骤2所述的利用估计得到的空时二维功率谱矩阵以及导向矢量矩阵A重构第l号距离单元的空时二维协方差矩阵并对进行修正，得到修正后的协方差矩阵的过程，包括如下步骤：

2a)首先将整个空时平面均分为K＝k_sK_t个网格点，K_s是空间频率轴被等分的数量，K_t是多普勒频率轴被等分的数量，各点对应的归一化空间频率和归一化多普勒频率分别为f_s，n，n＝1，2，...，K_s和f_d，m，m＝1，2，...，K_t，空时导向矢量表示为：

s(fd,m,fs,n)=[1,ej2πfd,m,...,ej2πfd,m(M-1)]T⊗[1,ej2πfs,n,...,ej2πfs,n(N-1)]T]]>

其中M为脉冲数，N为空域接收通道数，[·]^T表示转置，表示Kronecker积；

2b)利用估计得到的空时二维功率谱矩阵和导向矢量矩阵A重构第l号距离单元的空时二维协方差矩阵为：

其中A=[s(fd,1,fs,1),s(fd,2,fs,1),...,s(fd,Kt,fs,Ks)]]]>表示MN×K维的空时导向矢量矩阵；

2c)在空时平面上找到与(f_dk，f_s0)距离最近的4个网格点，这4个网格点的坐标记为(T_di，f_si)(i＝1，…，4)，由功率矩阵确定这4个点对应的功率，记为p_i(i＝1，…，4)，其中f_s0为雷达主波束方向对应的归一化空间频率；

2d)对进行如下修正，得到修正后的协方差矩阵：

R^CN,l=R^l-AaDaAaH]]>

其中D_a＝diag{p₁，p₂，p₃，p₄}，diag{·}表示对角矩阵，A_a＝[s(f_d1，f_s1)，…，s(f_d4，f_s4)]为相邻四点组成的导向矢量矩阵。