讨论一下如何自己写CWT
本帖最后由 ricefcc 于 2012-12-5 09:45 编辑最近困扰在一个非常基础的问题,希望找到一起搞小波的人探讨:
众所周知,连续小波变换CWT算出的矩阵就是信号与小波函数滤波器的卷积,但是我测试了一下db2小波,算出来的却和CWT函数结果不一样,我的方法是:
s=;
W=cwt(s,1:1:10,'db2');
取W(2,:)就是scale a=2下的连续小波变换结果,这个是matlab算出的
下面是我自己算出的:
用wfilters('db2')找出小波高通滤波器为hir=[-0.1294 -0.2241 0.8365 0.4830]
然后S与hir卷积,取中间8个值,得到ws2,这个ws2应该对应于CWT的尺度为2结果
但是ws2和上面的W(2;:)完全不一样,但是对哈尔小波,结果是一致的,俺验证过——为何db的就不一致呢?
请教下大家。 感觉国内论坛就这里搞小波的人最多了 木有人吗?急呢~~~~~~~~~~~~~ 每天顶一顶{:3_47:}
页:
[1]