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

- Radon transformclear 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)***************************************** RATnaf = ambifunb (x) htl (abs (naf)) [wh, rho, theta ] = htl (abs (naf)) colormap ([0,0,0]) xlabel (polar radius) ylabel (angle)************************************** to find the coordinates of the peak point, calculating the initial slope of the frequency and FM (right) find the peak point of the coordinates b = max (max ( wh)) [u, a] = find (wh