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