[发明专利]基于张量子空间投影的高频地波雷达射频干扰抑制方法

权利要求书

1.一种基于张量子空间投影的高频地波雷达射频干扰抑制方法，其特征在于：依赖于射频干扰的距离相关性和方向特性，通过张量分解从多个维度对射频干扰子空间进行联合估计，实现过程包括以下步骤，

步骤1，数据分块处理，包括通过基于雷达数据进行多普勒滑窗处理完成数据分块；

步骤2，射频干扰检测；

步骤3，选择第一个数据块，判断当前选择的数据块是否存在干扰频点，若含有射频干扰，则进入步骤4进行后续处理；否则，不做进一步处理，选取下一个数据块继续进行判断；

步骤4，构造训练张量，包括设高频地波雷达有M个天线通道，利用M个天线通道的无海洋回波的远距离元数据构造三阶张量，作为训练张量张量大小为P×Q×M，P为选取的远距离单元个数，Q等于该数据块的多普勒范围长度；

步骤5，采用高阶奇异值分解对训练张量进行张量分解，如下式：

式中，为核心张量，U_n为对应张量的第n模式展开矩阵X_(n)的左奇异矩阵，×_n表示张量和矩阵的n模式乘积，其中n＝1,2,3；

步骤6，根据步骤5所得结果估计射频干扰子空间；实现方式如下，

取U_n，n＝2,3，中前d_n个大奇异值对应的左奇异向量分别张成X_(n)的射频干扰子空间，即和它们共同决定该数据块中的射频干扰子空间，和对应的子空间投影矩阵分别为

SVD分解后，所得U_n中奇异值是从大到小排序的，取前d_n个奇异值，提取训练张量中的主成分，即射频干扰；

步骤7，距离滑窗和子空间投影处理，包括沿距离轴对数据块进行滑窗处理，得到若干与训练张量大小相等的处理张量，向射频干扰子空间进行投影得到该数据块中的射频干扰分量，从而得到无射频干扰的数据；

步骤8，判断当前是否已经处理最后一个数据块，如果否则返回重复步骤3-7，直至循环处理完所有的数据块，得到无射频干扰的M个天线通道的距离多普勒谱。

2.根据权利要求1所述的一种基于张量子空间投影的高频地波雷达射频干扰抑制方法，其特征在于：步骤1中，进行多普勒滑窗处理，是沿多普勒轴进行滑窗，相邻滑窗之间设置重叠区间，将距离多普勒张量数据分割成若干个子张量数据块，分别进行后续的抗射频干扰处理。

3.根据权利要求2所述的一种基于张量子空间投影的高频地波雷达射频干扰抑制方法，其特征在于：步骤7中，距离滑窗的窗长和步长都为P个距离单元。

4.根据权利要求3所述的一种基于张量子空间投影的高频地波雷达射频干扰抑制方法，其特征在于：步骤7中的子空间投影处理实现方式为，记某处理张量为计算无射频干扰的数据对各处理张量分别处理后，综合各处理张量的射频干扰分量，得到数据块中的射频干扰分量。