ricefcc 发表于 2012-12-4 22:33

讨论一下如何自己写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的就不一致呢?
请教下大家。

ricefcc 发表于 2012-12-5 09:47

感觉国内论坛就这里搞小波的人最多了

ricefcc 发表于 2012-12-5 21:42

木有人吗?急呢~~~~~~~~~~~~~

ricefcc 发表于 2012-12-8 09:59

每天顶一顶{:3_47:}
页: [1]
查看完整版本: 讨论一下如何自己写CWT