dftxinhao1.m

来自「雷达信号带限外推。包括burg外推的方法」· M 代码 · 共 19 行

clear;
N=16;
n=0:N-1;
k=0:N-1;
F1=0.320;
F2=0.560;
sga=0.01;
r1=randn(1,N);
r2=randn(1,N);
Z1=sga*sqrt(-2*log(r1)).*(cos(2*pi*r2));
Z2=sga*sqrt(-2*log(r1)).*(sin(2*pi*r2));
Z=Z1+j*Z2;
Xn=3*exp(j*F1*pi*n)+2*exp(j*F2*pi*n)+Z;
nk=n'*k;
WN=exp(-j*2*pi/N);
WNnk=WN.^nk;
Xk=Xn*WNnk;
figure(1),plot(n,Xn);
figure(2),plot(k,abs(Xk));