论文---基于MATLAB的地震数据的分析.docVIP

地震数据的MATLAB分析 地震带给人类的损失是巨大的，汶川大地震依旧在我们的记忆深处清晰存在。大地震的每次不约而至，都对国家和人民造成了巨大的损失。地震预测是世界性难题，全世界的地震科学家不断在探索，尽最大努力减少其破坏性。地震台站提供的地震观测资料的可靠性和准确性，是地震学家进行地震预测的基础[1]。但是地震波信号变化的不平稳性和复杂性给地震的分析和预测带来了很大困难，并且在地震波的原始记录中往往还掺杂着来自外界的各种干扰，如仪器、环境噪声、爆破、采矿、火车的震动等，这些都给地震波的分析带来了严重影响，甚至导致分析结果的错误。为保证地震分析的准确性，可先对原始记录进行频谱分析，选择性能优良的滤波器对其进行优化处理，把干扰信号尽量滤除，然后再对处理后的地震波数据进行分析处理，则会得到良好的效果。 地震记录的数字化使得利用计算机对地震信号进行分析处理得以实现。 在数字信号的分析处理中，Fourier变换和数字滤波器的应用极为普遍，语音、雷达、地震、图像、机械振动、地质勘探等众多领域都广泛采用数字滤波器。MATLAB是一种集数值分析、矩阵运算、信号处理和图像显示于一体、功能极其强大的高性能软件，其工具箱中包含了各种经典和现代数字信号处理技术，很容易实现Fourier变换和各种数字滤波器的设计，在地震数据的分析处理中起着重要作用。 在地震波的原始记录数据中往往夹杂不同频率范围的噪声干扰信号，为了显示出地震波数据中的优势频率和干扰频率，保证地震分析的准确性，应首先采用频谱分析，再针对干扰波的频率范围，设计合适的滤波器参数。 在对有限长信号序列进行频谱分析时，离散Fourier 变换（DFT）应用非常广泛，它可以很好地反映序列的频谱特性。 设是一个长度为的有限长序列，则的点离散傅里叶变换为 在MATLAB中可调用函数来进行快速傅里叶变换。其中，为时域内的输入信号序列，N为序列长度，为频率域的输出信号，即的频谱特征。 FIR数字滤波器的传递函数为 （2） 其中，是FIR滤波器的单位脉冲响应，其长度为，非零区间为。（1）式表明是的次多项式，它在平面上有个零点。原点是阶重极点。因此，肯定是稳定的。 由FIR滤波器的传递函数，可得到其频率响应表达式 （3） 如果是实序列，并且满足下列中心对称条件[6] 或 （4） 则FIR滤波器设计在逼近平直幅频特性的同时，还能获得严格的线性相位特性。 FIR滤波器虽然具有相位滞后的缺点，但是其相位滞后和群延迟在整个频带上是相等且不变的。一个阶的线性相位FIR滤波器群延迟为，即滤波后的信号简单地延长常数个时间步长，这一特性使通带频率内信号通过滤波器后仍然保持原有波形形状而无相位失真。 FIR数字滤波器的设计方法很多，本文采用目前最优秀的设计方法——切比雪夫逼近法。切比雪夫逼近法是一种等波纹逼近法，它使误差在整个频带均匀分布，对同样的技术指标，这种逼近法需要的滤波器阶数低，而对同样的滤波器阶数，这种逼近法的最大误差最小。设希望设计的滤波器幅度特性为，实际设计的滤波器幅度特性为，误差加权系数为，则切比雪夫最佳一致逼近准则是使加权函数[7] （5） 的最大值达到最小，即满足 （6） 切比雪夫逼近理论可以解决的存在性、唯一性及如何构造等一系列问题，克服了通带和阻带的边缘不易精确确定和Gibbs现象等多种缺点，使设计出的滤波器具有明显优点。因为在一定意义上对做最佳逼近，可获取较好的通带和阻带性能，并能准确地指定通带和阻带的边缘频率，此时又具有了IIR滤波器的优点，所以这是一种很好的设计方法。 MATLAB提供了大量设计FIR数字滤波器的函数，下面针对切比雪夫逼近法，介绍其设计步骤和方法。 （1）根据得出的地震波原始记录的频谱分析波形，确定FIR数字滤波器的技术指标； （2）用函数[8]估算设计的等波纹逼近法的参数：最低滤波器阶数、频率向量、幅度向量和加权向量。其中，为采样频率；可以是模拟频率或归一化频率，但必须以0开始，以（用归一化频率为1）结束，而且省略了0和两个频点；为各逼近段允许的幅度响应偏差（波纹振幅）；a为滤波器在各个频段上的幅度值，一般对通带取值为1，对阻带取值为0。 （3）用函数[8]完成等纹波FIR滤波器的设计，并用函数或完成对输入信号的滤波，得到输出信号。 地震数据的MATLAB分析： 下面以汕头台受到距台站300m处的汽车干扰的波形记录图的数字地震记录资料为例，进行基