f4_1.m

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

%产生信号
sig=[fmlin(128,0.05,0.45);fmlin(128,0.35,0.15)];
%时域波形
plot(real(sig),'LineWidth',2);
hold on;
plot(imag(sig),'ro--');
legend('实部','虚部');
xlabel('时间 t');
ylabel('幅值 A');
%计算短时傅立叶变换
[tfr,t,f]=tfrstft(sig);
%幅值图
figure(2)
plot(t,abs(tfr(1:128,:)),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');
%相位图
figure(3)
plot(t,angle(tfr(1:128,:)),'LineWidth',2);
xlabel('时间 t');
ylabel('相位 phi');
%时频表示
figure(4)
contour(t,f(1:256),abs(tfr));
xlabel('时间 t');
ylabel('频率 f');