S转换

S转换(s-transform)是一种时频分析的工具。
和其他时频分析工具一样，透过S转换，我们可以同时从时域以及频域观察一个信号的能量分布。S转换的特别之处在它既保持与傅立叶变换的直接关系，又可在不同频率有不同的解析度。此外，S转换与小波转换（wavelet transform）有密切的关系，或可视为连续小波转换（continuous wavelet transform）的变形。S转换的清晰度略优于加伯转换（Gabor transform），而不如韦格纳分布（Wigner distribution function）、科恩克莱斯分布、改良式韦格纳分布（Modified Wigner distribution function）。
定义.
一个信号x(t)的S转换为
其中窗函数为高斯窗函数
另种表示-频谱表示式.
藉著折积定理 
S转换能以频域 formula_4 表示，
这里可将formula_6看成formula_7与formula_8的卷积，
将formula_9 以及 formula_10分别取傅立叶变换可得
窗函数(window function).
对一讯号formula_12进行S转换，写为
formula_13
其中 formula_14为S转换之窗函数，常见之窗函数为高斯函数，即
formula_2 
事实上，S转换之窗函数并不局限于高斯函数，考虑以下一般化型态的S转换(generalized S transform):
formula_16
其中 formula_17为 一般化之窗函数(generalized window function), 参数 P 为一系列数值而组成之矩阵，可控制window function，举例来说
formula_18 ，此时 formula_19 为矩阵 P 内之唯一数值
通过选取较小的 formula_19 值，可得到时域上较窄的窗函数，提高时间轴的解析度，然而，频率轴上的解析度将因此变差。
改善时频两域解析度的问题，可选择非对称的窗函数(asymmetric window)，在正方向(forward direction)上有较陡峭之斜率，反方向(backward direction)上有较缓之斜率，借由牺牲终端时间的解析度(一般而言较不重要)换取初始时间上更好的解析度。
逆S转换(inverse S-transform).
S转换可以沿著时间轴方向积分，将可以得到formula_12的频谱formula_4。推导如下，
利用Gaussian window所包含面积等于1的特性，
因此，沿著时间轴t积分，
这表示S频谱是可逆的，同时也提供一个简单的逆转换。
滤波应用(Filtering).
S转换如同其他时频分析转换，皆可以设计波器来达到消除杂讯留下讯号的功用，<br>
利用逆S转换，我们可以设计一个S域的滤波器"U(t,f)"，对x(t)进行讯号处理<br>
离散时间S转换.
S转换相较于加伯转换，虽在清晰度有较好的改善，但也有其缺点，就是运算复杂度变高，积分的范围会随著formula_28的增加而增加。<br>
因此，这里利用上面推导的频谱表示式来推导离散时间S转换<br>
频谱表示式
令formula_30<br>
如果要使用FFT的方式来实作，必须另加条件
首先先对formula_12做傅立叶变换得到formula_4<br>
接著带入频谱表示式中，<br>
当 m=0 时，S转换就定义成
流程.
Step1 : 计算formula_39,这个步骤只需要计算一次。<br>
Step2 : 计算formula_40for formula_41<br>
Step3 : 将formula_42移动至formula_43<br>
Step4 : 将Step2,Step3的结果相乘得到
Step5 : 对B[m,p]取逆离散傅立叶变换(IDFT)可得到，formula_45在formula_41的行向量<br>
Step6 : 重复Step2~5直到formula_45全部定义完成。
S转换特性.
S转换与加伯转换(Gabor Transform)很相似，
唯一的不同就在于S转换的Gaussian Window的宽度会随著formula_50改变。<br>
原因就是formula_28在高频时，Gaussian Window宽度变小，时域解析度好；反之，formula_28在低频时，Gaussian Window宽度变宽，频域解析度好。<br>
但是，当formula_53时，Gaussian Window会无穷无尽的变宽，就丧失时频分析只做局部分析的精神。<br>
一种解决的方式是:使Gaussian Window宽度不再因formula_54改变\，产生频宽剧烈的变化，<br>
S转换一般式
S转换是一种运算量高的时频分析工具，尤其在低频部分，Gaussian Window宽度变宽，频域解析度比加伯转换来的好，所以S转换对于低频讯号分析比较有优势<br>
例如:声音讯号，人耳对高频的部分没有太特别的感觉，但在低频部分却比较敏感，如:中央Do = 262Hz,高八度Do = 512Hz 可以很清楚的听出两个不同的音，
但10000Hz 和 10170Hz对人来说差别不大，再说人耳对3KHz以内的声音最敏感，所以能分析低频讯号就显得重要。<br>
此时，就可以使用S转换，来强调低频讯号，而牺牲高频讯号。
与韦格纳分布的比较.
韦格纳分布是时频分析工具中，具有高清晰度的一个，但最大的缺点是有交叉项（cross-term）的问题。若一个信号是由数个信号成份组合而成，那么使用韦格纳分布来分析时就会受到两两信号成份之间的交叉项干扰，这将会产生一些不必要的噪声。一个信号x的韦格纳分布为
交叉项是在积分中两个x项相乘时产生的。S转换的计算原理与韦格纳分布不同，是直接对formula_60进行转换，不会有交叉项的问题。
与加伯转换的比较.
加伯转换的定义为
formula_61
我们知道加伯转换是短时距傅立叶转换的一种特殊形式，其中只要把短时距傅立叶转换的窗函数用高斯函数来替代就成了加伯转换；S转换则可视为一种窗函数会随f变化的加伯转换；随著频率的升高，高斯函数在时域上的宽度会越来越窄，使得时域上的解析度会增加，反之牺牲频域上的解析度。
加伯转换和S转换原理相同，两者唯一不同的地方就是窗函数的formula_62和强度formula_50，基本上都是由短时距傅立叶转换延伸而来；两者共有的好处是不会像韦格纳分布一样会有交叉项；又S转换在低频时的频率解析度会优于加伯转换。
与小波转换的关系.
连续小波转换可以视为将一个信号对小波做相关（correlation）：
而S转换可以视为连续小波转换乘上一个相位项：
而其用的母小波为：
S转换，加伯转换和短时距傅立叶转换的比较.
假设我们想要对一讯号formula_12做时频分析
formula_68
formula_69
formula_70
如果使用一个宽度为1秒的矩形函数来做短时距傅立叶转换会得到
如果使用一个加伯转换会得到时频分析图
如果用S转换会得到时频分析图
由上三图可知，S转换是一个在低频时频域解析度高，高频时时域解析度高的时频分析；举例来说，对于100Hz与300Hz和1000Hz与1200Hz这两组声音，哪一组对于人耳会有明显的差异呢? 答案是100Hz和300Hz，低频时的频率差异对于人耳较明显，当频率越高时，人耳就越难分辨出频率的差异；同样的道理，S转换便符合我们的需求，低频时讯号变化慢拥有低的时域解析度和拥有高的频域解析度，高频时因为讯号变化很快则拥有高的时域解析度和低的频域解析度。