f6_11.m

来自「《Matlab时频分析技术及其应用》电子书。」· M 代码 · 共 27 行

%产生线性调频信号
%信噪比为1dB
N=64;
sig=sigmerge(fmlin(N,0,0.3),noisecg(N),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');