讨论一下如何自己写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]