资 源 简 介
％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)
