chirp导频信号2.m

来自「chirp导频信号 共同参考」· M 代码 · 共 32 行

clear all;
close all;
N=256;
for n=1:N
% sig(n)=real(exp(j*pi*1*n+j*pi*0.0012*n.^2));
sig(n)=exp(j*pi*1*n+j*pi*0.0012*n.^2);
end
n=1:N;
% figure(1)
% plot(n,sig);
% title('时域chirp导频信号')
[y, tmin, tmax] = wvd(sig,1,N);
figure(1)
mesh(y);
title('chirp导频信号的wvd')
a=0:179;%Radon Trannsform
figure(2)
[R,xp]=radon(y,a);
mesh(R);
title('chirp导频信号的rwt')
figure(3)
[H,rho,theta]=htl(y);
mesh(H);
title('chirp导频信号的wht')
l=0:N-1;
for i=1:length(l);
   dcft(i,:)=zct(sig,l(i),N);
end;
K=0:length(dcft)-1;
figure(4);
mesh(l,K,dcft/sqrt(N));
title('chirp导频信号的mcft')