广义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大时窗变化反而告宏很小裤亮，会比较有用的。

matalb求得结果与手动计算得结果是等并轮价的。

1. 手动计算得结果

f1 = x0*cos(8*t)*exp(-6*t) + sin(8*t)*exp(-6*t)*(v0+6*x0)/8

f1 = expand(f1, 'ArithmeticOnly', true)

2. matalb求得结果

f2 = x0*exp(-6*t)*(cos(8*t) + sin(8*t)*((v0 + 12*x0)/(8*x0) - 3/4))

f2 = expand(f2, '绝谈信ArithmeticOnly'侍灶, true)

3.

simplify(f1-f2)

ans = 0

[st,t,f] = st(timeseries,minfreq,maxfreq,samplingrate,freqsamplingrate)

后面两桥差判掘个参数就是时域和频域敏冲皮的采样率

---