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

数字信号处理学生所在学院：测试与光电工程学生姓名：黎德任课教师：李志农教师所在学院：测试与光电工程2014年10月2014 级—种典型的时频分析方法一S变换黎德测试与光电工程学院摘 要：目前提出的时频分析方法有很多，并在机械、电子、生物、医学等很多领域得到了广泛 的应用本文主要介绍一种典型时频分析方法，S变换关键词：时频分析方法；s变换s变换是小波变换和短时傅里叶变换结合发展起来的一种时频分析方法，它采用与 频率(即尺度)有关的高斯窗函数，也即采用频率作为自变量的高斯函数，其时频分辨 率随着频率的不同而变化，具有线性、多分辨率、逆变换唯一等特点，而且它获得的二 维时频谱与傅里叶变换保持着直接的联系，已经在盲信号分离、医学图像处理、地震波 分析、故障检测等很多领域得到了广泛应用信号S变换的分辨率与频率有关，且其结果具有无损可逆性［1］,然而S变换中的基本 小波函数是固定不变的，因此其在实际应用中受到了限制•高静怀等提出了根据“实际需 要”恰当地选择或构造基本小波函数的广义S变换方法⑵.Pinnegar给出了既可调节窗 口标准差又非对称的另外一种广义S变换方法^Pinnegar在研究信号频率随时间变化 的问题时，提出了多尺度和非零相位复数窗的广义S变换方法⑷.S变换方法在国内外得到了广泛的应用:在地震勘探中用于时频分析［“、地震相和主 频分析XI、人工地震探测［9-13］、地震波的衰减分析［14］、地震层特征分析［15 ］和油气探测与 预测in ］等.尤其是应用S变换在时变滤波中去除噪声具有明显的效果。

1S变换的定义S变换是Stockwell在研究众多常用的时频分析方法的基础上提出的，它是一种随着 频率发生变化，时频分辨率也会发生变化的时频分析方法［18］由于S变换具有很好的时 频特性，因而在很多领域中都被广泛应用［19-20］s变换在分析高频信号或低频信号时具 有良好的特性，对于高频信号，窗口会比较窄，满足高时间分辨率要求，而对于低频信号， 窗口会比较宽，满足高频率分辨率要求S变换的定义如下：S G ,f）=+『h（t）2以上公式成为h(t)的S变换，t表示的是时间,e _ 2 碑dtT用来控制高斯窗口在时间轴上的位置反变换描述为:h（t）二了 fslr,f 1/Tei 2 呻df2 S变换的原理2.1 一维连续S变换S变换由于1996年提出，是一种可逆的局部时频分析方法/其思想 是对连续小波变换和短时傅立叶变换的发展信号x（t）的S变换S（t,f）定义如下：S(t, f )= J x(tbG — t, f )e(-j2顶)dt式中®G-1 ,f）为高斯窗口-8 (1)（Gaussian Window）； T为控制高斯窗口在t轴位置的参数由式中可以看出，S变换不同于短时傅里叶变换之处在于高斯窗口的高度和宽度随 频率而变化，这样就克服了短时傅里叶变换窗口高度和宽度固定的缺陷。

信号x（t）可以由其S变换S（，f）很好的重构，其S逆变换为：exp(j2时t )dfxC)二了 f S (c, f )dT—8」一8S变换可以看作是对连续小波变换的一种相位修正，并可以从连续小波变换推导而 来信号乙）的连续小波变换可以定义如下：W (T, d )= T x(t b(t -t , d )dt—8式中d,T分别为伸缩参数和时移参数；W（t-T,d）为母小波的伸缩时移变换如果选 取变换核为一高斯窗和以复向量的乘积，即：① C, f )=2 丿e(-j2f注意，此时伸缩参数d为频率f的倒数于是，信号x（t）的S变换式（1）可以表示为式（5）作为变换核进行连续小波变换, 再乘上一个相位校正因子，即：S(T, t)= ej 2沪 W G, d )2.2 —维离散S变换对（1）右边先做传统的Fourier变换，再做Fourier反变换，最后进行变量代换, 将S变换转换成信号x（t）的傅里叶反变换X f）的函数，即：S C, f )=『X (v + f )e-2k2v2丿ej2kvtdv, f 丰 0(7)-8对（7）进行离散化，可以得到S变换的离散表示形式:m=0m + nNT-2k 2 m 2j 2Kmke n2 e nS 册，o]=N、mN N丿其中:NT1 * h(kT)e 卜曽｝Nk=1式中 k,m,n=0,l,2, N—l（n丰 0）。

信号S变换时，频矩阵中元素的幅值为:nNTnS kT, 一 ,(k, n = 0,1,NTN -1)(11)显然，时间序列x（t）的S变换结果是一个复时间频矩阵，记为S矩阵，其行对应频 率，列对应时间，将S矩阵各个元素求模后得到的矩阵记为S模矩阵［21,］其列向量表示 信号某一时刻的幅值随频率变化的分布因此S摸矩阵某位置元素的大小就是相应频率 和时间处信号S变换的幅值，S变换的结果二维等高线图等直观的表示3 S变换算法由于使用了褶积定理和快速傅里叶变换，S变换的计算速度显得比较快在实际应 用中，可以按以下步骤计算S变换后，将谱H倍］频移n得到a•计算h（t）的快速傅里叶变换（FFT）的谱H［］频移谱H寫，n为频率采样点数b. 直接计算得到g（t, f）-占L e弋的快速傅里叶变换G（m, n）= e-以c. 按频率采样点计算h［雳%滋n）d.计算H罟G滋n）的快速傅里叶反变换（ifft），得到S变换谱Sn八NT4 S变换与其他典型信号分析方法对比4.1 S变换与短时傅里叶变换对比S变换相对于短时傅里叶变换的主要优点是:S变换属于多分辨率时频分析方法而 短时傅里叶变换是单分辨率的时频分析方法。

从下面的两种变换的时频分布网格可以更 清楚的了解到两种变换的时频分析窗口的区别Q短时傅里叶变换 （b）S变换4-1时频分布网络短时傅里叶变换由于使用的窗口比较固定，在整个频率范围内，它的频率分辨率并 没有发生变化；而S变换由于可以自动的调节窗口，在高频处频率分辨率相对变低，有 较高的时间分辨率，这正符合非平稳信号的分析要求4.2 S变换与小波对比通过前面介绍，知道S变换和小波变换都是多分辨率分析方法我们更加关注它们之 间存在的区别：小波变换并不是严格意义上的时频分析，它并没有与频率完全对应，而S 变换则不同,它与频率可以对应：小波变换的效果与所选用的小波函数具有直接的关系， 小波存在多解性，在重构时难免会产生信息丢失，而S变换是无损可逆的4.3 S变换与Wig ner — Ville分布对比对于含有多信号分量的信号wigner — ville分布将会产生交叉项，这会严重影响信 号的分析，而S变换不存在交叉项，没有这类影响S变换和wigner 一 ville分布均可以捕捉到信号，但交叉区域存在大量的干扰，严 重的影响了信号的分析，这是wgigner 一 ville分布的不足S变换是一种时频分辨率 自适应变化的分析方法。

总的来说，S变换更适合于信号的时频分析，尤其对于比较复杂的包含较多信号分 量S变换作为一种良好的时频分析方法，已经被广泛的应用在各种信号处理中，如： 地震波低频伴影观测［22］，面波频散分析［23 ］，时频域补偿［24 ］等5 S变换的应用S变换应用于广泛的信号处理中，以下介绍几种S变换的应用1）电能质量扰动信号分析的应用现代电力系统中，电力电子设备的应用越来越广泛，各种非线性、冲击性、波动性 负载也大量增加，使电力系统所遭受的电能质量污染也日趋严重同时信息科技的发展 则对电能质量及供电可靠性提出更高的要求对电能质量进行监测和分析是发现电力质 量问题并进行治理和改善的前提条件，而如何从大量的电能质量扰动信号中自动提取特 征并进行正确分类则是电能质量监测分析系统首先要解决的问题电能质量的扰动信号主要有8种，分别是电压凹陷、电压隆起、短时中断、电压尖 峰、电压缺口、谐波、间谐波和瞬态振荡S变换是由连续小波变换核短时傅里叶变换 结合发展起来的一种时频分析方法，其引入了宽度与频率成反向变换的高斯窗，具有与 频率相关的分辨率由于S变换具有很好的时频特性，因而S变换非常适合电能质量扰动信号的提取S变换利用S变换模矩阵的时频曲线对信号进行处理分析，准确地定位 扰动发生时刻、确定扰动持续时间。

S变换能够对电能质量扰动信号进行简单直观地分 类2） 初至拾取方法的应用初至拾取就是要确定地震道上纯噪音信号和噪音与地震叠加信号之间的分界时刻 初至的拾取在地震资料处理中有着非常重要的作用，它是折射波资料分析、VSP处理和 井间地震资料处理的基础工作，并为其他处理工作提供可靠参数它的准确与否直接影 响到后续工作的好坏初至拾取的方法主要分为两大类，一是自动拾取，另外是手动拾 取两种方法各有优缺点对于信噪比较低的地区（特别是复杂山区低信噪比地区）， 自动拾取的初至一般会产生较大的误差，而手动又不可能一道一道拾取，只能选取几个 初至较清楚的道称为参考道，其它道通过内插得到（称为内插道），在这种情况下，参 考道的初至如果不准确，势必影响内插道的初至因此参考道的初至的准确显得尤为重 要基于此，S变换不但显示了多分辨率的特点，还能将时间-频率的关系直观的在时频 图上表现出来，让处理人员更好的对信号进行处理操作当人工检测无法准确的进行地 震信号的初至识别时，基于S变换的地震信号初至识别方法能够处理比较准确的检测出 地震信号的初至，给人工识别提供一定的技术支持，是地震信号初至识别的一种新的尝 试和方法。

3） 结构损伤信号处理的应用在很多行业中，采用各种数字信号处理方法对实际测量的振动信号进行分析和处理， 提取各种特征，用于参数检测、质量评价、状态监视和故障诊断，因此振动信号的数字 处理方法一直是近10年的主要研究方向之一沏近几年来，盲信号分离和循环统计量 也开始应用于振动信号分析中［26 27］在实际的结构健康监测（SHM ）数据中，除了时域信息 外，还有多种频率成分，利用时频分析技术，可以获得很多的分析参数，能够更好地对 资料进行处理和解释S变换结合了短时傅里叶变换与小波变换的优点，并在一定程度上克服了它们的缺 点，适用于非平稳信号的时频分析与傅里叶变换不同的是，傅里叶变换只显示整个时 域每个频率的总强度，S变换实现每个频点在时域每个时间点的频率信息复合材料超声回波信号降噪的应用复合材料由于具有低热胀系数、高耐热性，在动力结构方面的高比强度、高比模量、 耐磨损和抗老化等特点，被广泛的应用到固体火箭发动机等生产中［28］复合材料在生产 和使用过程中，不可避免地会产生各种缺陷，严重影响结构的强度，甚至导致事故的发 生超声是复合材料常用的无损检测方法［29］但由于超声检测复合材料过程中会产生晶 粒噪声和电噪声［30］，因此必须对超声回波信号进行消噪处理。

超声信号降噪处理的方法有很多，常用的时频分析消噪方法有短时傅里叶变换、小 波变换等Mo短时傅里叶变换虽然一定程度上克服了标准傅里叶变换不具有的局部分析 能力的缺陷，但由于窗函数选定后只能改变窗口在相平面上的位置，不能改变窗口的形 状，因此不适用于超声等非平稳信号小波分析能够对超声信号进行降噪，但小波阈值 的选取和小波阈值函数的确定比较复杂，且小波软阈值降噪后信号会产生较大的失真， 硬阈值降噪后信号不连续，重构信号会产生振荡［32］作为短时傅里叶变换和小波变换的一种扩展,S变换逐渐引起人们的兴趣,S变换 是 由Stock well等学者于首次提出的，。