f6_12.m

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

%产生两分量信号
N=64;
sig=sigmerge(fmlin(N,0,0.4),fmlin(N,0.3,0.5),1);

%时域波形
figure(1);
subplot(2,1,1);
plot(1:N,real(sig),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');

%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%显示Wigner-Ville分布等高线
figure(2);
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');

%Hough变换
[WH,rho,theta]=htl(tfr,N,N);
figure(3);
mesh(rho,theta,WH);
xlabel('极半径 rho');
ylabel('角度 theta');