f6_1.asv

来自「《Matlab时频分析技术及其应用》（人民邮电2006年出版）一书的MATALB」· ASV 代码 · 共 22 行

%产生线性调频信号
sig=fmlin(128,0.1,0.4);
%设置窗函数
h=window(@Kaiser,17);
%获取坐标
[tfr,t,f]=tfrsp(sig);
%计算重排谱图
[tfr,rtfr,hat]=tfrrsp(sig,1:128,64,h,1);

%显示谱图
figure(1);
contour(t,2*f(1:length(f)/2),tfr);
xlabel('时间 t');
ylabel('频率 f');
%%三维图

figure(2)
%显示重排谱图
contour(t,2*f(1:length(f)/2),rtfr);
xlabel('时间 t');
ylabel('频率 f');