用S变换做精细时变滤波.doc

-、引言地震信号的频谱是时变的，其中的噪音更是各个时段各不相同。要想将这些各时段频谱不同的噪音准确干净地滤掉而又不影响有效信号，最有效的办法是构造随时间变化的滤波器。目前，地震信号的滤波是由快速傅里叶变换完成的，由于快速傅里叶变换的局限性及其与时间的无关性，使得用其完成时变滤波有一定困难，也很不准确。有鉴于此，笔者尝试应用s变换对地震记录进行时变滤波，取得了较为理想的效果。二、S变换的基本原理1996年，美国地球物理学家Stockwell在前人成果的基础上提出了一种新的时频分析方法一S变换，它组合了短时傅里叶变换和小波变换的优点，同时避免了其不足，因而很快被应用于地震信号的时频分析领域。Stockwell等在文献中对S变换作了详细的阐述。S变换的一维正变换定义。S变换的定义式可用的傅里叶频移谱H(a+ f)表示为：我们在这张地震记录上用S变换的方法试图去除线性干扰。图2-b