esprit_doa.m

已知：信号中心波长为2,天线阵元的间距为1米,快拍数为2000

function doa=esprit_doa(Y,n,d)
%
% The ESPRIT method for DOA estimates 
%
% call doa=esprit_doa(Y,n,d)
%
%      Y    <- the ULA data
%      n    <- the number of sources
%      d    <- sensor spacing in wavelengths
%      doa  -> the vector of DOA estimates

%copyright by kily 2008

[m,N]=size(Y);

% compute the sample covariance matrix
R=Y*Y'/N;

% do the eigendecomposition; use svd because it sorts eigenvalues
[U,D,V]=svd(R);

S=U(:,1:n);

phi = S(1:m-1,:)\S(2:m,:);
w=angle(eig(phi));
doa=asin(w/d/pi/2)*180/pi;