f7_3.m

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

%两线性调频信号叠加
N=64;
sig=sigmerge(fmlin(N,0,0.4),fmlin(N,0.3,0.5),1);
%计算信号的Wibner-Ville分布
[tfr,t,f]=tfrwv(sig);
%绘图
%%绘制等高线图
figure(1);
contour(t,f,abs(tfr).^2);
xlabel('时间 t');
ylabel('频率 f');
%%三维图
figure(2);
mesh(t,f,abs(tfr).^2);
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');

%进行Hough变换
[HT,rho,theta]=htl(tfr,N,N,1);
%%绘制等高线图
figure(3);
contour(rho,theta,HT');
xlabel('rho');
ylabel('theta');
%%三维图
figure(4);
mesh(rho,theta,HT');
xlabel('rho');
ylabel('theta');
zlabel('幅值 A');