lanyanlin 发表于 2010-11-12 21:33

小波变换的频率成分提取

本帖最后由 lanyanlin 于 2010-11-12 21:34 编辑

想求出信号的实际频率的0.04Hz-0.15Hz范围内的小波系数和能量，请问下面的代码对吗？

s=importdata('interval.txt');
wavename='cmor2-1';
SampFreq = 1;
dataLength=length(s);
fmax = 0.15;      %指定频率范围上限                                    
fmin = 0.04;         %指定频率范围下限                        
fc = centfrq(wavename);   %中心频率                                                         
totalscal =128;                                                                     
FreqBins = linspace(fmin,fmax,totalscal);      
Scales = fc./ FreqBins;                              
RealFreqBins = FreqBins * SampFreq;            
MWT=cwt(s,Scales,wavename);%指定频率范围内的小波系数
h=mean(abs(MWT).^2,1);          %指定频率范围内的能量和

hahaer 发表于 2010-11-12 21:56

个人感觉没有什么问题
只要注意时间-尺度变换中尺度和频率成倒数关系就可以了

c2019 发表于 2010-11-13 15:53

好东西，顶

smallqiao 发表于 2010-11-17 23:11

回复 2 # hahaer 的帖子

不是简简单单的倒数关系吧？
求详解

ice_fire 发表于 2011-4-21 13:51

很简单吗

ice_fire 发表于 2011-4-21 14:49

f = scal2frq(a,wname,delta);
a 尺度wname = 'coif3'; delta = 0.1;采样周期f代表伪频率

yu婧 发表于 2011-5-26 16:12

过来学习。。

杨德昌 发表于 2011-5-26 18:26

学习一下！

hjhnh999 发表于 2011-5-27 16:13

不是很懂啊！
页: [1]
查看完整版本: 小波变换的频率成分提取