radon_ambiguity.m

％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

%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;

%N=length(x);
% ambifunb(x');
%*****************************************RAT

naf=ambifunb(x);
htl(abs(naf));
% [wh,rho,theta]=htl(abs(naf));
colormap([0,0,0]);

% xlabel('极半径');
% ylabel('角度');


%**************************************%找出峰值点的坐标,计算初始频率和调频斜率(正确)
%找出峰值点的坐标
b=max(max(wh)); 
[u,a]=find(wh>=0.8*b);