www.pudn.com > Esprit_try.zip > Esprit_try.m, change:2012-07-03,size:670b


clear all 
close all 
derad = pi/180; 
w0 = 2*pi; 
M = 8; 
d = 0.5; 
dd = 0:d:(M-1)*d; 
doa = 3; 
theta = [10 40 60]; 
snr = 10; 
N = 500; 
A = exp(-j*w0*dd.'*sin(theta*derad)); 
mu = w0*d*sin(theta*derad); 
Phim = diag(exp(j*mu));%相位延迟矩阵 
Aphim = A*Phim; 
AA = [A;Aphim]; 
S = randn(doa,N); 
Z = AA*S;%构造数据矩阵 
Z = awgn(Z,snr,'measured'); 
Rzz = Z*Z'/N;%由阵列输出采样值的自相关函数估计Rzz 
[E,Delt] = eig(Rzz); 
D = diag(Delt).'; 
[D,I] = sort(D); 
D = fliplr(D); 
E = fliplr(E(:,I)); 
Es = E(:,1:doa); 
Ex = Es(1:M,:); 
Ey = Es(M+1:2*M,:); 
F = pinv(Ex)*Ey; 
[Ef,Lamda] = eig(F); 
L = diag(Lamda).'; 
Mu = imag(log(L)); 
Doatheta = asind(Mu/(w0*d))