[发明专利]一种音频大地电磁数据工频噪声干扰的滤波方法

权利要求书

1.一种音频大地电磁数据工频噪声干扰的滤波方法，其特征在于，包括如下步骤：对有噪声干扰的音频大地电磁数据作带阻滤波，得到无工频噪声干扰数据。

2.根据权利要求1所述的滤波方法，其特征在于，包括如下步骤：

(1)判断音频大地电磁数据的原始时间序列文件的信号道功率谱曲线是否存在噪声干扰，得到第一判断结果；(2)所述第一判断结果为是，对存在工频噪声干扰的信号道作带阻滤波，得到滤波后的数据；(3)将滤波后的数据保存在一维数组中；(4)对将滤波后的数据进行格式处理，得到格式文件。

3.根据权利要求2所述的滤波方法，其特征在于，所述步骤(1)具体包括如下步骤：读取采样率为150Hz的原始时间序列文件中的数据，对数据的每个信号道均做功率谱分析，得到每个信号道的功率谱曲线，根据每个信号道的功率谱曲线判断所述信号道是否存在工频噪声干扰，如果信号道的功率谱曲线出现毛刺信号或增益异常，则判断该信号道存在工频噪声干扰；如果信号道的功率谱曲线是一条近似平行横轴且连续光滑的曲线，则判断该信号道不存在工频噪声干扰。

4.根据权利要求2所述的滤波方法，其特征在于，所述步骤(2)具体包括如下步骤：设定凯瑟窗参数，调用matlab中的凯瑟窗函数构建凯瑟窗WinFun，WinFun＝(kaiser(N,β))'；将二维数组存储区中要滤波的信号道数据与WinFun作卷积运算，得到滤波后的数据，完成对存在工频噪声干扰的信号道的滤波。

5.根据权利要求4所述的滤波方法，其特征在于，以FIR数字滤波器作带阻滤波。

6.根据权利要求4所述的滤波方法，其特征在于，所述凯瑟窗参数的设定包括：定频率带阻滤波器阻带低频边缘ω(p)、带阻滤波器阻带高频边缘ω(s)和衰减系数As，所述ω(p)、ω(s)和As分别设定为40Hz～50Hz，50Hz～60Hz和0～100。

7.根据权利要求2所述的滤波方法，其特征在于，所述步骤(2)第一判断结果为是，对存在工频噪声干扰的信号道作带阻滤波前，还包括如下步骤：读取采样率为2400Hz的原始时间序列文件的每个信号道的数据并保存在二维数组存储区中对应的行数中。

8.根据权利要求7所述的滤波方法，其特征在于，所述二维数组存储区的分配方法如下：二维数组存储区＝zeros(采样道数，采样率×采集总时间)，其中，zeros指的是二维数组存储区初始化全部全部置0，所述采样道数是指二维数组存储区的行数，所述采样率×采集总时间为二维数组存储区的列数，即每行数据的总长度，所述采样率、采样道数和采集总时间通过读取采样率为2400Hz的原始时间序列文件的标签信息得到。

9.根据权利要求7所述的滤波方法，其特征在于，所述一维数组的数据长度与所述二维数组的列数相同。

10.根据权利要求2所述的滤波方法，其特征在于，步骤(4)具体包括如下步骤：将二维数组中被滤波的信号道数据替换成一维数组数据，然后将替换后的二维数组数据转换成二进制数据，文件保存后缀为TS3，文件名与采样率为150Hz的原始时间序列文件相同。