xxf 发表于 2007-3-22 11:33

求助:关于连续小波变换的问题

有两个问题没搞懂:
1  =wavread('E:\Heart Sound Signal\www.qdheart.com\normal\normal',20000);
   x=decimate(signal,4);
   Fs=fs/4;%采样频率为2000Hz
   coefs=cwt(x,,'morl','3Dplot');
   colormap('jet');
   for m=1:200
      A(m)=sum(coefs(m,:));
   end
   figure(2)
   plot(A)
   红色的部分是我将频率相同的点幅值相加,画一个频率-幅值图,这样对吗?
2 实际频率=采样频率*小波中心频率/尺度
  书上说小波中心频率可以自己定义,可我怎么没找到相关的函数呢
  比如上面的morlet小波,怎么来定义呢?
      =cmorwavf(0,128,129,40,60);
      书上说这样可以得到复morlet小波,可没法对他调用啊。
      括号里的数是我自己写的不知道有没有道理

[ 本帖最后由 zhlong 于 2007-6-4 17:21 编辑 ]

Anny2006 发表于 2007-3-22 11:44

红色部分可以改成:
   A=sum(coefs,2);
不过这样算的只是简单的将同一尺度上的小波系数求和,我认为应该是平方求和
   A=sum(coefs.*coefs,2);
我有发帖,呵呵,有空可以看下,有时间可以讨论一下!
对了,你这样画出来的是尺度--幅值图,不是频率-幅值图!必须用scal2frq将尺度转化为频率:@)

xxf 发表于 2007-3-22 15:22

信号经过小波变换后在某一尺度上会出现负值,所以才要加平方的吗?

xiaopanpan 发表于 2008-11-10 19:21

连续小波变换所得并不代表各尺度对应频率下的幅值,我认为此系数表示信号与该尺度下的小波函数的相似程度。还请大家考证

bkly 发表于 2008-11-12 09:29

同意楼上。小波系数实际代表了,原始信号在不同尺度小波基上的投影。也就是相似程度。带有原信号很多信息,但并非原始信号的直接划分。
如果要求幅值,还是应当重构后采用
页: [1]
查看完整版本: 求助:关于连续小波变换的问题