chirp导频信号.m

chirp导频信号

clear all;
close all;
N=512;
for n=1:N
%sig(n)=real(exp(j*pi*1*n+j*pi*0.0012*n.^2));
 sig(n)=exp(j*pi*0.76*n+j*pi*0.001*n.^2);
end
[y, tmin, tmax] = wvd(sig,1,N);
figure(1)
contour(y,5);
a=0:179;%Radon Trannsform
figure(2)
[R,xp]=radon(y,a);
mesh(R);
title('chirp导频信号的rwt')
figure(3)
[H,rho,theta]=htl(y);%Hough Transform
mesh(H);
title('chirp导频信号的wht')