[发明专利]一种基于子空间误差修正的线性扫频干扰抑制方法及系统

权利要求书

1.一种基于子空间误差修正的线性扫频干扰抑制方法，其特征在于：包括以下步骤：

步骤1：对用阵列天线接收到的信号白化处理并进行时频分析得到优化后的接收信号时频分布图；

步骤2：从优化后的时频分布图中遴选自源点；

步骤3：通过自源点估计线性扫频频率参数即初始频率f_i和调频斜率g_i，i为第i个线性扫频干扰信号，i＝2,…,P,P为干扰信号的数量；

步骤4：对所估计的频率参数进行修正；

步骤4.1：求取根据所估计的频率参数构建的LFM干扰信号和接收信号的相关值；

所述相关值是指：对用估计参数构建的LFM干扰信号

与接收信号的某一路x(t)求内积，

式中T为求内积的时间长度，其离散形式为：

L为基矢量序列长度；

步骤4.2：以估计值为中心，选取周边点；

设定初始频率和调频斜率的步进值δ_f和δ_g，在以初始频率f为横轴，调频斜率g为纵轴的平面上，以为中心，选取其左右上下四个点，即

步骤4.3：求取各个点对应的内积

将周围的四个点对应的频率参数值代入(6)式当中，计算相应的内积值

步骤4.4：更新中心点；

比较中心点和周围四个点对应内积的大小，用五个点中对应内积最大的点取代中心点，返回步骤4.2，如果中心点对应的内积已经是最大，则停止运算，输出该内积最大值所对应的点坐标，该坐标为修正后的初始频率和调频斜率；

步骤5：使用修正后的频率参数构建干扰子空间，用子空间投影消除干扰。

2.根据权利要求1所述的一种基于子空间误差修正的线性扫频干扰抑制方法，其特征在于：步骤1得到时频分布图的具体方法为：

步骤1.1：求取接收信号矢量的白化矩阵W；

步骤1.2：用白化矩阵将接收信号矢量白化并用Wigner-Ville变换求取时频分布，得到优化后的接收信号时频分布图。

3.根据权利要求2所述的一种基于子空间误差修正的线性扫频干扰抑制方法，其特征在于：所述步骤2中遴选自源点的方法为：

步骤2.1：对白化后的接收信号矢量z(n)求取空时频分布矩阵序列，设时频分析中所取的数据段长度为L，如下式

n表示时间点，ω表示归一化频率，l表示偏移量，D_zz(n,ω)表示白化后的接收信号矢量z(n)的空时频分布矩阵序列；

步骤2.2：求得时频分布图上所有时频点对应的空时频分布矩阵的迹并归一化，即Trace(E{D_zz(n,ω)})/L，设ζ＝Trace(E{D_zz(n,ω)})/L；

步骤2.3：设置门限值ε，上述空时频分布矩阵的迹均已归一化，其值域为(0,1)，0.1≤ε≤0.3,当空时频分布矩阵的迹ζ＞ε时，则此矩阵对应的时频点(n,ω)属于自源点。