matlab地震频谱分析,《基于MATLAB的地震数据的分析》.doc

《基于MATLAB的地震数据的分析》.doc

基于MATLAB的地震数据的分析

孙玉柱 冯光 房桂梅

摘要：地震波原始数据中存在的干扰信号，会影响震相分析的准确性。为了滤除干扰信号，对地震波原始信号进行了频谱分析，给出了一种基于MATLAB的FIR数字滤波器的优化设计方案，将其用于地震波数据的分析中，并进行了仿真分析。仿真结果表明，FIR数字滤波器对地震波原始信号进行滤波处理后，提高了震相分析的准确性，得到了理想的效果，达到了预期的目的。

4 基于MATLAB的地震数据的分析

下面以汕头台受到距台站300m处的汽车干扰的波形记录图的数字地震记录资料为例，进行基于MATLAB的地震数据的分析。该地震台采用的是频带范围为0.05~20Hz 的FBS-3A型宽频带数字地震记录仪和EDAS-C24型数据采集器，系统的采样频率为50Hz。

将该地震波原始数据保存在MATLAB的work文件夹中，并对其进行地频谱分析，如图1所示。从图1中可以很清楚的看出地震波原始记录中地震信号的优势频率为0.25Hz，频段范围在0~1Hz之间；干扰的优势频率为12.5Hz，频段范围在10~15Hz之间。

图1 地震波原始数据的时域波形图和频谱分析图

Figure 1 the original time-domain and spectrum analysis of the earthquake data

为了滤除干扰信号，最大限度的保留其中的有用信号，选取的FIR滤波器为带阻滤波器，其参数选定为：通带上截止频率Fp1=7，阻带下截止频率Fs1=7.1，阻带上截止频率Fs1=18.9，通带下截止频率Fp2=