f6_3.m

来自「《时频分析技术》这本书中的MATLAB CODE」· M 代码 · 共 49 行

%产生三分量信号
N=60;
%%信号分量1
[sig1,ifl1]=fmsin(N,0.15,0.35,50,1,0.35,1);
%%信号分量2
[sig2,ifl2]=fmlin(N,0.3,0.1);
%%信号分量3
[sig3,ifl3]=fmconst(N,0.4);
sig=[sig1;zeros(8,1);sig2+sig3];
%时域波形
figure(1);
plot(real(sig),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');
%获取坐标
[tfr,t,f]=tfrwv(sig);
%Wigner-Ville分布
figure(2);
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
[tfr,rtfr,hat]=tfrrspwv(sig);
%重排平滑伪Wigner-Ville分布
figure(3);
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
%计算重排谱图
[tfr,rtfr,hat]=tfrrsp(sig);
figure(4)
%显示重排谱图
contour(t,2*f,rtfr);
xlabel('时间 t');
ylabel('频率 f');
%计算重排M-H分布
[tfr,rtfr,hat]=tfrrpmh(sig);
figure(5)
%显示重排M-H分布
contour(t,2*f,rtfr);
xlabel('时间 t');
ylabel('频率 f');
%计算重排Page分布
[tfr,rtfr,hat]=tfrrppag(sig);
figure(6)
%显示重排Page分布
contour(t,2*f,rtfr);
xlabel('时间 t');
ylabel('频率 f');