%兩線性調(diào)頻信號(hào)疊加
N=64;
sig=sigmerge(fmlin(N,0,0.4),fmlin(N,0.3,0.5),1);
%計(jì)算信號(hào)的Wibner-Ville分布
[tfr,t,f]=tfrwv(sig);
%繪圖
%%繪制等高線圖
figure(1);
contour(t,f,abs(tfr).^2);
xlabel('時(shí)間 t');
ylabel('頻率 f');
%%三維圖
figure(2);
mesh(t,f,abs(tfr).^2);
xlabel('時(shí)間 t');
ylabel('頻率 f');
zlabel('幅值 A');

%進(jìn)行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');