[发明专利]一种运用子空间投影的线性扫频干扰抑制方法及系统

权利要求书

1.一种运用子空间投影的线性扫频干扰抑制方法，其特征在于：包括以下步骤：

步骤1：通过阵列天线获取接收信号，估计接收信号中LFM干扰的初始频率和调频斜率；

步骤1.1：将接收信号经短时傅里叶STFT变换后映射到时频分布图上，获取干扰信号的自源点；

步骤1.2：通过各线性扫频干扰信号的自源点，估计每一路干扰信号的调频斜率和初始频率；

步骤2：根据所估计的初始频率和调频斜率构建干扰子空间；

步骤2.1：估计子空间基矢量长度；

步骤2.1中所述估计子空间基矢量长度的方法为：

步骤2.1.1:求取所构建的子空间基矢量与准确的子空间基矢量匹配程度的相关值,对基矢量长度为L的干扰信号x(t)和以所估计的初始频率和调频斜率构建的干扰信号做内积，写成离散形式的相关值为相关值为L的函数，令归一化的相关值为T_norm(L)＝T(L)/T(L_ini)，其中L_ini表示初始的基矢量长度，其中，为估计出的第i个线性扫频干扰信号的初始频率，为估计出的调频斜率，x_i^*(n)为对接收信号x(t)离散化后的共轭值；

步骤2.1.2：当T_norm(L)＞T_cor时，令L＝L+ΔL，返回步骤2.1.1；T_norm(L)≤T_cor时，则L为估计的子空间基矢量长度，输出L；其中ΔL为L的步进值，T_cor为门限值，0.9＜T_cor＜1；

步骤2.2：根据所估计的子空间基矢量长度构建子空间基矢量；

根据估计的LFM干扰参数构建干扰信号每隔L长度取一段数据作为子空间基矢量，其中第i个干扰信号矢量表示为

式中表示第i个干扰，表示对x_i截取一段L长度的数据；

步骤2.3：构建干扰子空间及其正交子空间；

用时域扩展后的干扰信号矢量构建成一个矩阵，如下

该组矢量张成了干扰的子空间，其正交子空间表示如下

M表示天线阵阵元数。

步骤3：通过子空间投影法消除干扰；

将接收信号矢量投影到干扰子空间的正交子空间，如下式

式中X为接收信号矢量，通过投影过程消除了接收信号中的干扰分量X_u，剩余有用信号和白噪声分量X_s和V。

2.根据权利要求1所述的一种运用子空间投影的线性扫频干扰抑制方法，其特征在于：步骤1.2中所述估计每一路干扰信号的调频斜率和初始频率的方法为：

步骤1.2.1：寻找第i个线性扫频干扰对应时频脊线的邻域，i＝2,…,P,P为干扰信号的数量；在第i个线性扫频干扰信号对应的时频分布脊线上任取一个自源点，过该点作一条直线，使倾角α在0°～180°中变化，设定直线的邻域半径r_line，所述邻域半径r_line不超过干扰在时频分布图的时频脊线在任意时间点上频宽的一半，将时频分布图上各自源点到该直线的距离小于r_line的自源点集合做为该直线的邻域Ω_i，统计不同倾角下直线的邻域Ω_i中自源点的数量，当某一倾角下直线的邻域Ω_i中自源点的数量比该直线处于其它倾角时自源点的数量都大时，则将该直线倾角作为第i个线性扫频干扰信号在时频分布面上的倾角，根据该倾角计算出第i个线性扫频干扰信号的调频斜率；

步骤1.2.2：去除第i个线性扫频干扰对应时频脊线的邻域与其它线性扫频干扰对应时频脊线的邻域的交集，得到去除交集后的区域Ω′_i,用最小二乘法拟合Ω′_i邻域中的自源点，得到第i个干扰信号的拟合直线方程则为估计出的第i个线性扫频干扰信号的初始频率，为估计出的调频斜率。

3.一种运用子空间投影的线性扫频干扰抑制系统，其特征在于：包括处理器，以及与所述处理器连接的存储器，所述存储器存储运用子空间投影的线性扫频干扰抑制程序，所述运用子空间投影的线性扫频干扰抑制程序被所述处理器执行时实现上述权利要求1或2所述方法的步骤。