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

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

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)

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


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/12338234.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-24
下一篇 2023-05-24

发表评论

登录后才能评论

评论列表(0条)

保存