[发明专利]一种基于自适应滤波器的干扰信号处理方法及装置

权利要求书

1.一种基于自适应滤波器的干扰信号处理方法，其特征在于，包括：

步骤1：获取天文观测信号与干扰参考信号，进行信号同步；

步骤2：对天文观测信号进行延迟处理，信号延迟大小设置为N/2，N表示FIR数字滤波器的阶数；

步骤3：将天文观测信号对应的干扰参考信号分成两路，对其中一路干扰参考信号进行延迟处理，延迟大小等于FIR滤波器阶数的一半；将另一路干扰参考信号输入所述自适应滤波器，以对所述干扰参考信号进行数字滤波；

步骤4：延迟处理的天文观测信号中减去数字滤波后的干扰参考信号，得到对应的输出信号；

步骤5：将所述输出信号以及延迟处理后的干扰参考信号输入滞后相关器，以利用所述滞后相关器进行互相关运算；

步骤6：对所述滞后相关器输出的信号进行积分，并计算积分的平均值；

步骤7：根据所述积分的平均值，更新所述自适应滤波器的滤波系数；

针对观测信号每个采样的输出，重复上述步骤2至步骤7，以获取输出信号的全部采样数据。

2.根据权利要求1所述的方法，其特征在于，所述滞后相关器通过以下公式进行互相关运算：

其中，f(n)是输出信号，r(n)为干扰参考子信号，C(i)为滞后相关器输出的信号，K是采样点数，i和n的取值为0，1，2，……，K。

3.根据权利要求1所述的方法，其特征在于，通过以下公式对所述滞后相关器输出的信号进行积分，并计算积分的平均值：

其中，C(i)为滞后相关器输出的信号，C_av(k)为积分的平均值，N表示自适应滤波器的阶数。

4.根据权利要求1所述的方法，其特征在于，通过以下公式计算所述自适应滤波器的滤波系数：

h(n)＝h(n-1)+ε×C_av(k)

其中，h(n)为滤波系数，C_av(k)为积分的平均值，n＝0，1，……，N-1，ε 是预设的常数，N表示自适应滤波器的阶数。

5.根据权利要求1所述的方法，其特征在于，经过所述数字滤波后的干扰参考信号为：

其中，y(n)为经过所述数字滤波后的干扰参考信号，h(k)为滤波系数，r(n-k)为经过k个延迟单元的干扰参考信号，N表示自适应滤波器的阶数。

6.根据权利要求1所述的方法，其特征在于，还包括：

预先对所述自适应滤波器进行初始化，以将滤波系数设置为零。

7.一种基于自适应滤波器的干扰信号处理装置，其特征在于，包括：

信号同步模块，获取天文观测信号与干扰参考信号，计算两路信号之间的延迟时间，然后对天文与干扰信号进行同步，完成预处理任务；

天文观测信号处理模块，根据预设的自适应滤波器的阶数，对天文观测信号进行延迟处理；

干扰参考信号处理模块，用于将天文观测信号对应的干扰参考信号分成两路，根据所述自适应滤波器的阶数，对其中一路干扰参考信号进行延迟处理；将另一路干扰参考信号输入所述自适应滤波器，以对所述干扰参考信号进行数字滤波；

输出信号确定模块，用于利用延迟处理后的天文观测信号减去数字滤波后的干扰参考信号，得到对应的输出信号；

互相关运算模块，用于将所述输出信号以及延迟处理后的干扰参考信号输入滞后相关器，以利用所述滞后相关器进行互相关运算；

积分均值模块，用于对所述滞后相关器输出的信号进行积分，并计算积分的平均值；

滤波系数更新模块，用于根据所述积分的平均值，更新所述自适应滤波器的滤波系数。

8.一种计算机可读存储介质，其上存储有计算机指令，所述指令被执行时实现权利要求1-6中任意一项所述方法的步骤。