f6_1.m

希望对大家学习有所帮助

%产生线性调频信号
sig=fmlin(128,0.1,0.4);
%设置窗函数
h=window(17,'Kaiser');
%获取坐标
[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');