1.一种音频大地电磁数据工频噪声干扰的滤波方法，其特征在于，对有噪声干扰的音频大地电磁数据作带阻滤波，得到无工频噪声干扰数据，包括如下步骤： (1)判断音频大地电磁数据的原始时间序列文件的信号道功率谱曲线是否存在噪声干扰，得到第一判断结果；具体包括如下步骤：读取采样率为150Hz的原始时间序列文件中的数据，对数据的每个信号道均做功率谱分析，得到每个信号道的功率谱曲线，根据每个信号道的功率谱曲线判断所述信号道是否存在工频噪声干扰，如果信号道的功率谱曲线出现毛刺信号或增益异常，则判断该信号道存在工频噪声干扰；如果信号道的功率谱曲线是一条近似平行横轴且连续光滑的曲线，则判断该信号道不存在工频噪声干扰； (2)所述第一判断结果为是，对存在工频噪声干扰的信号道作带阻滤波前，读取采样率为2400Hz的原始时间序列文件的每个信号道的数据并保存在二维数组存储区中对应的行数中，所述二维数组存储区的分配方法如下：二维数组存储区＝zeros(采样道数，采样率×采集总时间)，其中，zeros指的是二维数组存储区初始化全部全部置0，所述采样道数是指二维数组存储区的行数，所述采样率×采集总时间为二维数组存储区的列数，即每行数据的总长度，所述采样率、采样道数和采集总时间通过读取采样率为2400Hz的原始时间序列文件的标签信息得到； 以FIR数字滤波器作带阻滤波，对存在工频噪声干扰的信号道作带阻滤波，得到滤波后的数据，具体包括如下步骤：设定凯瑟窗参数，调用matlab中的凯瑟窗函数构建凯瑟窗WinFun，WinFun＝(kaiser(N,β))'，β为依赖于滤波器长度N的参数，凯瑟窗函数公式0≤n≤N-1，I0是修正后的零阶贝赛尔函数，Δω＝ω(p)-ω(s)，凯瑟窗参数包括定频率带阻滤波器阻带低频边缘ω(p)，带阻滤波器阻带高频边缘ω(s)，和衰减系数As，ω(p)、ω(s)和As分别设定为40Hz～50Hz，50Hz～60Hz和0～100； 将二维数组存储区中要滤波的信号道数据与WinFun作卷积运算，得到滤波后的数据，完成对存在工频噪声干扰的信号道的滤波； (3)将滤波后的数据保存在一维数组中； (4)对将滤波后的数据进行格式处理，得到格式文件，具体包括如下步骤：将二维数组中被滤波的信号道数据替换成一维数组数据，然后将二维数组中被滤波的信号道数据替换成一维数组数据后的二维数组数据转换成二进制数据，文件保存后缀为TS3，文件名与采样率为150Hz的原始时间序列文件相同。

请联系平台

请联系平台