juzhen.asv

雷达信号带限外推。包括burg外推的方法

N=16;
I=64;
K1=8;
K2=;
P=K2-K1+1;
A=[];
for bb=1:1:P
    for aa=1:1:N
        A(aa,bb)=exp(j*(aa-1)*2*pi*(K1+bb-1)/(N+I));
    end
end 
n=0:N-1;
k=0:N-1;
F1=0.325;
F2=0.420;
sga=0.001;
r1=rand(1,N);
r2=rand(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;
D=(A'*A)\(A'*Xn');
m=K1:K2;
plot(m,D);
D