一种典型的时频分析方法—S变换.docx

数字信号处理

学号：

学生所在学院：测试与光电工程

学生姓名：黎德

任课教师：李志农

教师所在学院：测试与光电工程

2014年10月

2014级
一种典型的时频分析方法—S变换
黎德
测试与光电工程学院

摘要：目前提出的时频分析方法有很多，并在机械、电子、生物、医学等很多领域得到了广泛的应用。本文主要介绍一种典型时频分析方法，S变换。
关键词：时频分析方法；S变换

S变换是小波变换和短时傅里叶变换结合发展起来的一种时频分析方法，它采用与频率（即尺度）有关的高斯窗函数，也即采用频率作为自变量的高斯函数，其时频分辨率随着频率的不同而变化，具有线性、多分辨率、逆变换唯一等特点，而且它获得的二维时频谱与傅里叶变换保持着直接的联系，已经在盲信号分离、医学图像处理、地震波分析、故障检测等很多领域得到了广泛应用。
信号S变换的分辨率与频率有关,且其结果具有无损可逆性[1],然而S变换中的基本小波函数是固定不变的,因此其在实际应用中受到了限制.高静怀等提出了根据“实际需要”恰当地选择或构造基本小波函数的广义S变换方法[2].Pinnegar给出了既可调节窗口标准差又非对称的另外一种广义S变换方法[3].Pinnegar在研究信号频率随时间变化的问题时,提出了多尺度和非零相位复数窗的广义S变换方法[4].
S变换方法在国内外得到了广泛的应用:在地震勘探中用于时频分析[5-6]、地震相和主频分析[8-9]、人工地震探测[9-13]、地震波的衰减分析[14]、地震层特征分析[15]和油气探测与预测[16-17]等.尤其是应用S变换在时变滤波中去除噪声具有明显的效果。

1S变换的定义
S变换是Stockwell在研究众多常用的时频分析方法的基础上提出的,它是一种随着频率发生变化,时频分辨率也会发生变化的时频分析方法[18]。由于S变换具有很好的时频特性，因而在很多领域中都被广泛应用[19-20]。s变换在分析高频信号或低频信号时具有良好的特性,对于高频信号,窗口会比较窄,满足高时间分辨率要求,而对于低频信号,窗口会比较宽,满足高频率分辨率要求。
S变换的定义如下：

以上公式成为h(t)的S变换，t表示的是时间，τ用来控制高斯窗口在时间轴上的位置。
反变换描述为：

2S变换的原理
2.1 一维连续S变换
S变换由R.G.Stockwell于1996年提出，是一种可逆的局部时频分析方法/其思想是对连续小波变换和短时傅立叶变换的发展。
信号x(t)的S变换S(t,f)定义如下：
(1)
(2)
式中为高斯窗口（GaussianWindow）；为控制高斯窗口在轴位置的参数。由式中可以看出，S变换不同于短时傅里叶变换之处在于高斯窗口的高度和宽度随频率而变化，这样就克服了短时傅里叶变换窗口高度和宽度固定的缺陷。
信号可以由其S变换很好的重构，其S逆变换为：
(3)
S变换可以看作是对连续小波变换的一种相位修正，并可以从连续小波变换推导而来。信号的连续小波变换可以定义如下：
(4)
式中分别为伸缩参数和时移参数；为母小波的伸缩时移变换。如果选取变换核为一高斯窗和以复向量的乘积，即：
(5)
注意，此时伸缩参数d为频率f的倒数。
于是，信号x(t)的S变换式（1）可以表示为式（5）作为变换核进行连续小波变换，再乘上一个相位校正因子，即：
(6)
2.2一维离散S变换
对（1）右边先做传统的Fourier变换，再做Fourier反变换，最后进行变量代换，将S变换转换成信号x(t)的傅里叶反变换X（f）的函数，即：
(7)
对（7）进行离散化，可以得到S变换的离散表示形式：
(8)
(9)
其中： (10)
式中k,m,n=0,1,2,……N-1(n0)。
信号S变换时，频矩阵中元素的幅值为：
(11)
显然，时间序列x(t)的S变换结果是一个复时间频矩阵，记为S矩阵，其行对应频率，列对应时间，将S矩阵各个元素求模后得到的矩阵记为S模矩阵[21]，其列向量表示信号某一时刻的幅值随频率变化的分布。因此S摸矩阵某位置元素的大小就是相应频率和时间处信号S变换的幅值，S变换的结果二维等高线图等直观的表示。

3S变换算法
由于使用了褶积定理和快速傅里叶变换，S变换的计算速度显得比较快。在实际应用中，可以按以下步骤计算S变换。
计算h(t)的快速傅里叶变换（FFT）的谱后，将谱频移n得到频移谱，n为频率采样点数。
直接计算得到的快速傅里叶变换。
按频率采样点计算。
计算的快速傅里叶反变换（IFFT），得到S变换谱。

4S变换与其他典型信号分析方法对比

4.1S变换与短