中文说明:
%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)
English Description:
%radon transform