%產生線性調頻信號
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');