f7_2.m

来自「《Matlab时频分析技术及其应用》随书源代码」· M 代码 · 共 45 行

%产生线性调频信号
N=64;
sig=fmlin(N,0,0.3);
%加噪声
sign=sigmerge(sig,noisecg(N),1);

%绘制时域波形
figure(1);
plot(real(sign),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');

%计算信号的Wibner-Ville分布
[tfr,t,f]=tfrwv(sign);

%绘图
%%绘制等高线图
figure(2);
contour(t,f,abs(tfr).^2);
xlabel('时间 t');
ylabel('频率 f');

%%三维图
figure(3);
mesh(t,f,abs(tfr).^2);
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');

%进行Hough变换
[HT,rho,theta]=htl(tfr,N,N,1);

%%绘制等高线图
figure(4);
contour(rho,theta,HT');
xlabel('变量 rho');
ylabel('变量 theta');

%%三维图
figure(5);
mesh(rho,theta,HT');
xlabel('变量 rho');
ylabel('变量 theta');
zlabel('幅值 A');