基于MATLAB的地数据的分析.docVIP

基于MATLAB的地震数据的分析 孙玉柱 冯光 房桂梅 摘要：地震波原始数据中存在的干扰信号，会影响震相分析的准确性。为了滤除干扰信号，对地震波原始信号进行了频谱分析，给出了一种基于MATLAB的FIR数字滤波器的优化设计方案，将其用于地震波数据的分析中，并进行了仿真分析。仿真结果表明，FIR数字滤波器对地震波原始信号进行滤波处理后，提高了震相分析的准确性，得到了理想的效果，达到了预期的目的。 关键词：MATLAB；FIR数字滤波器；优化；滤波 the Analysis of Earthquake Data Based on MATLAB SUN Yuzhu,FENG Guang,FANG Guimei Abstract: The interference that existed in the earthquake data will affect the accuracy of the seismic phase analysis. In order to filter the disturbance signal, this paper carries out spectrum analysis of the earthquake data, proposes an optimum HYPERLINK / \t _blank design method for FIR digital filter based on MATLAB and applies it to the analysis of earthquake data. After the filter of the noise jamming, the true information of the earthquake wave is clearly reflected. The simulation results manifest that it can improve the accuracy of seismic phase analysis and arrive at the purpose desired. Key words: MATLAB；FIR digital filter；optimization；filter 1 引言 地震带给人类的损失是巨大的，汶川大地震依旧在我们的记忆深处清晰存在。大地震的每次不约而至，都对国家和人民造成了巨大的损失。地震预测是世界性难题，全世界的地震科学家不断在探索，尽最大努力减少其破坏性。地震台站提供的地震观测资料的可靠性和准确性，是地震学家进行地震预测的基础[1]。但是地震波信号变化的不平稳性和复杂性给地震的分析和预测带来了很大困难，并且在地震波的原始记录中往往还掺杂着来自外界的各种干扰，如仪器、环境噪声、爆破、采矿、火车的震动等，这些都给地震波的分析带来了严重影响，甚至导致分析结果的错误。为保证地震分析的准确性，可先对原始记录进行频谱分析，选择性能优良的滤波器对其进行优化处理，把干扰信号尽量滤除，然后再对处理后的地震波数据进行分析处理，则会得到良好的效果。 地震记录的数字化使得利用计算机对地震信号进行分析处理得以实现。 在数字信号的分析处理中，Fourier变换和数字滤波器的应用极为普遍，语音、雷达、地震、图像、机械振动、地质勘探等众多领域都广泛采用数字滤波器。MATLAB是一种集数值分析、矩阵运算、信号处理和图像显示于一体、功能极其强大的高性能软件，其工具箱中包含了各种经典和现代数字信号处理技术，很容易实现Fourier变换和各种数字滤波器的设计，在地震数据的分析处理中起着重要作用。 本文首先介绍了数字滤波器，给出了基于MATLAB的FIR数字滤波器的一种优化设计方案，最后将其用于地震波数据的处理中，并进行了仿真分析。 2快速Fourier变换（FFT）及频谱分析 在地震波的原始记录数据中往往夹杂不同频率范围的噪声干扰信号，为了显示出地震波数据中的优势频率和干扰频率，保证地震分析的准确性，应首先采用频谱分析，再针对干扰波的频率范围，设计合适的滤波器参数。 在对有限长信号序列进行频谱分析时，离散Fourier 变换（DFT）应用非常广泛，它可以很好地反映序列的频谱特性[2]。 设是一个长度为的有限长序列，则的点离散傅里叶变换为[3] （1） DFT是信号分析与处理中的一种重要变换，但是当较大时，其计算量太大。快速Fourier 变换（FFT）是减少 DFT 运算次数的一种快速算法，通过在时域将序列逐次分解为一组子序列，然后利用子序列的DFT来实现整个序列的DFT，从而减少离散Fourier变换的运算量，提高了计算效率。 在MATLAB中可调用函数来进行快速