f6_2.m

matlab时频分析配套源码第3章的程序

%产生两分量信号
N=128;
%%信号分量1
[sig1,ifl1]=fmsin(N,0.15,0.45,100,1,0.4,-1);
%%信号分量2
[sig2,ifl2]=fmhyp(N,[1,0.5],[32,0.05]);
sig=sig1+sig2;
%时域波形
plot(real(sig),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');
%获取坐标
[tfr,t,f]=tfrsp(sig);
%计算重排谱图
[tfr,rtfr,hat]=tfrrsp(sig);
%显示谱图
figure(2);
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
%%三维图
figure(3)
%显示重排谱图
contour(t,f,rtfr);
xlabel('时间 t');
ylabel('频率 f');