f6_4.m

来自「工具箱提供了各种常用的信号处理」· M 代码 · 共 17 行

%产生线性调频信号
sig=fmlin(64,0.1,0.4);
%获取坐标
[tfr,t,f]=tfrsp(sig);
%计算重排Morlet小波尺度图
[tfr,rtfr,hat]=tfrrmsc(sig,1:64,64,2.1,1);
%显示Morlet小波尺度图
figure(1);
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
figure(2)
%显示重排Morlet小波尺度图
contour(t,f,rtfr);
xlabel('时间 t');
ylabel('频率 f');