％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)