%radon transform
clear all;
%
 N=800;
 n=1:N;
 fs=200;
 t=n/fs;
x1=exp(j*2*pi*(5*t+0.5*5*t.^2));
x2=exp(j*2*pi*(5*t+0.5*15*t.^2));
x=x1+x2;

%****************************************RWT
tfr=tfrwv(x);
%hough transmate
figure(2);
htl(abs(tfr));
[wh,rho,theta]=htl(abs(tfr));
% colormap([0,0,0]);
% mesh(wh);
% xlabel('極半徑');
% ylabel('角度');
%找出峰值點的坐標
b=max(max(wh)); 
[u,a]=find(wh=b);