[未答复] 请教连续小波变化cwt函数调用输入参数和频率相关问题

[复制链接]
KevinPeter 发表于 2022-5-12 19:07:40
本人想cwt计算信号的时频图,程序如下,但有报错(见下面),考虑到matglab函数库可能更新,所以使用了cwt 老版本的调用语句,会出现另一个错误,请问如何调整才能顺利计算时频图。

clear all,clc
fs=1024;
x=0:1/fs:1;

y=cos(2*pi*100.*x);

figure
plot(x,y)
title('original signal')

%%%% CWT version running on Matlab 2019b %%%%%
figure
[wt,ff,coi]=cwt(y,'morse',fs);
pcolor(x,ff,abs(wt));shading interp
hold on
plot(x,coi,'r')
title('new CWT version with COI')

以上语句报错信息:
                           错误使用 cwt
                           cwt : you must give at least 4 input parameters


%%%%%%%%% using old CWT form
wavename='cmor4-4';
totalscal=256;%scale length
fc=centfrq(wavename);
cparam=2*fc*totalscal;
a=1:1:totalscal;
scal=cparam./a;
f=scal2frq(scal,wavename,1/fs);
coefs=cwt(y,scal,wavename,1/fs);
figure
imagesc(x,f,abs(coefs));  %% how to plot the COI with the old version??
title('old CWT version without COI')


以上语句报错信息:
                           错误使用 cwt
                           cwt : Fmin and Fmax must be real scalars

请问如何调整才能顺利计算时频图,感恩感谢!



您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表