f5_12.m

来自「《Matlab时频分析及其应用》的详细代码」· M 代码 · 共 22 行

%生成两分量信号
sig=atoms(128,[32,0.15,20,1;96,0.32,20,1]);
%时域波形
figure(2)
plot(real(sig),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');
%计算Margenau-Hill分布
[tfr,t,f]=tfrmh(sig);
%显示
%%等高线图
figure(3);
contour(t,f,abs(tfr));
xlabel('时间 t');
ylabel('频率 f');
%%三维图
figure(4)
mesh(t,f,abs(tfr));
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');