广义S变换MatlAB程序中窗函数的编码问题

n即f的离散采用意义。W是窗函数w，它的傅里叶变换fft(W)，直接和输入信号的傅里叶变换形式H进行了离散卷积袜纯册，然后逆傅里叶变换得到对信号h的广义S变换。

w0(n,:)=exp(-n.^(2*p)*k.^2/(2*N.^2))*(n.^p)/(N*sqrt(2*pi))是窗函数，k是t的离散取样意义。公式的话就是简单的在S变换里把f变成了f^p（f的p次方)

这里p是要赋值的，它是给定的一个数值。如果p=1那就是原S变换，如果p>1则窗函数随f变化呈指数增长，一般这个效果不值得应用。反而是0<p<1时，则窗函数不对称性和窗宽随f呈对数变化，f小时窗变化剧烈，f大时窗变化反而告宏很小裤亮，会比较有用的。

为了进行时频局部化分析，Stockwell等学者于1996年引入了S变换这一数字信号处理技术。从概念上理解，S变运郑换可以看作是短时傅里叶变换与连续小波变换的结合，它可以像短时傅顷亩里叶变换一样得到局部化的时频信息，同时利雀悄森用长宽可变的时间窗，又可以像小波变换一样具有不同的频率分辨率。

---