MATLAB如何将系数缩放为连续小波变换函数的尺度图?

时间:2016-04-26 21:02:58

标签: matlab wavelet

我试图比较小波系数随时间和信号之间的差异。为此,MATLAB的{​​{1}}功能。输出cwt给出了小波系数和代表所述系数的尺度系数,就其对信号能量的总贡献而言。

我需要确定的是,我了解这是如何表现的?#34;引擎盖下#34;这样比较是有效的。我使用参考:http://www.mathworks.com/help/wavelet/ref/cwt.html

测试

我关心的是着色选项。在这个例子中,我们使用' scal'。文档没有告诉我它究竟是如何工作的。我试图使用我指定的采样率的其他着色选项,并且看起来缩放不是恒定的。我想知道它是如何计算的。

[coefs sgram] = cwt(foo)

如果我使用load vonkoch vonkoch=vonkoch(1:510); len = length(vonkoch); cw1 = cwt(vonkoch,1:32,'sym2','plot'); title('Continuous Transform, absolute coefficients.') ylabel('Scale') [cw1,sc] = cwt(vonkoch,1:32,'sym2','scal'); title('Scalogram') ylabel('Scale') cw1提取sc[ms is] = max(abs(cw1(:)))的最大值(绝对因为系数有时是负数),我会注意到它们占据相同的索引并且值是

[mc ic] = max(sc(:))

谁的比例mc = 0.0084 ms = 0.1313 。但是,测试其他两个值,ms/mc = 15.6310cw1(4,4)我无法重现此比率。 sc(4,4)如何运作?如果可能的话,将进一步描述其他着色模式的链接以及它们如何将系数缩放成尺度图表示。

1 个答案:

答案 0 :(得分:0)

来自cwt的比例(不再推荐使用此版本的cwt。请使用更新的cwt代替。)在引擎盖下有":

WSCALOGRAM Scalogram for continuous wavelet transform.
%   SC = WSCALOGRAM(TYPEPLOT,COEFS) computes the scalogram 
%   SC (percentage of energy for each coefficient). COEFS is
%   the matrix of continuous wavelet coefficients (see CWT).
%   The scalogram is obtained by computing:
%       S = abs(coefs.*coefs); SC = 100*S./sum(S(:))