马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
Zhlong和破凰曾给出用复morlet小波分析信号一个例子,但是用实morlet小波分析结果相差很远,请问大家这是什么原因?
SampFreq = 30;
t=1/SampFreq:1/SampFreq:4;
sig = sin(12*pi*t);
sig(1:end/2) = sig(1:end/2) + sin(6*pi*t(1:end/2));
sig(end/2+1:end) = sig(end/2+1:end) + sin(18*pi*t(end/2+1:end));
fmax = 0.5; % 最高分析频率(归一化频率)
fmin = 0.005; % 最低分析频率(归一化频率)
fb = 4 ; % 取cmor4-2小波进行实验,带宽参数为4
fc = 2; % 中心频率2Hz
totalscal = 512; % 所取尺度的数目
FreqBins = linspace(fmin,fmax,totalscal);% 将频率轴在分析范围内等间隔划分
Scales = fc./ FreqBins;% 计算相应的尺度参数
RealFreqBins = FreqBins * SampFreq;% 尺度所对应的实际频率
%-------------------------------------------------
MWT=cwt(sig,Scales,'cmor4-2');
%-------------------------------------------------
pcolor(t,RealFreqBins,abs(MWT));
colormap jet;shading interp;colorbar;
得到正确的结果:
如果用
MWT=cwt(sig,1:512,'morl');
pcolor(t, 1:512,MWT); colormap jet;shading interp;colorbar;
如果用
MWT=cwt(sig,Scales,'morl');
为什么用morl 实小波变换得不到正确结果呢?
[ 本帖最后由 swt 于 2008-11-28 19:47 编辑 ] |