[发明专利]海洋大地电磁探测的海浪磁干扰压制方法

权利要求书

1.一种海洋大地电磁探测的海浪磁干扰压制方法，其特征在于，包括：

S1、设置海浪磁干扰压制执行参数，所述参数包括滤波器阶数、最大迭代次数和最小迭代误差；

S2、选择并读取参与海浪磁干扰压制的海洋大地电磁场数据，将电场信号E作为输入信号，磁场信号H作为参考信号；

S3、基于滤波参数C与电磁场的阻抗关系，计算滤波参数初值C₀，所述滤波参数初值C₀的计算方法为：

式中，为阻抗张量；

S4、基于输入信号E与参考信号H中有效信号的相关性，结合信号与噪声之间随时间变化的统计特性,利用前一时刻的滤波参数进行迭代计算，自适应计算现时刻的滤波参数C，所述现时刻的滤波参数C的计算方法为：

C^H(m+1)＝C^H(m)+2μE(m)ε^*(m)

式中，m表示m时刻；电场E(m)为输入信号；磁场H(m)为参考信号；H′(m)为预测信号；C(m)为滤波系数；ε(m)为磁场H(m)与预测信号H′(m)的残差，其计算方法为ε(m)＝H(m)-H′(m)；上标H表示共轭转置，*表示复共轭；μ为滤波器步长，满足λ_max为相关矩阵R＝E[E(m)E^H(m)]最大的特征值；

S5、基于滤波参数C与电场E、磁场H的关系，重构预测信号H′，所述重构预测信号H′的方法为：

H′(m)＝C(m)E(m)；

S6、基于最小均方差算法原理，计算预测信号H′与参考信号H的残差ε；

S7、判断是否满足最小迭代误差或达到最大迭代次数，若满足则转向S8，若不满足则转向S4；

S8、输出压制干扰的预测信号H′。