tls_esprit.m

谱估计源程序：是谱估计中的经典算法esprit,MUSIC等算法的源程序

%%TLS-ESPRIT
clear all; close all;clc;
J=sqrt(-1);
source_number=3;
source_doa=[40 70 80];
sensor_number=10;
snapshort_number=100;
snr=10;

A=exp(-J*(0:sensor_number-1)'*pi*sin(source_doa*pi/180));
%theta=[20,40,60];
%fei=[10,30,50];
%for i=1:3
%A(:,i)=exp(-j*pi*cos(theta(i)*pi/180-2*pi*(0:9)/10)*sin(fei(i)*pi/180));
%end

s=(randn(source_number,snapshort_number)+J*randn(source_number,snapshort_number));
x=A*s;
X=awgn(x,snr);
 
M=[exp(j*pi*sin(source_doa(1)*pi/180)),exp(j*pi*sin(source_doa(2)*pi/180)),exp(j*pi*sin(source_doa(3)*pi/180))];
M=diag(M);
y=A*M*s ;
Y=awgn(y,snr);

 
 Z=[X',Y']';
 Rzz=Z*Z'/snapshort_number;
 [Uz,Sz,Vz]=svd(Rzz);
% % [Uz,Sz]=eig(Rzz);
 

 E=Uz(:,1:source_number);
 Ex=E(1:sensor_number,:);
 Ey=E(sensor_number+1:2*sensor_number,:);
 fy=pinv(Ex)*Ey;
 
[U,S]=eig(fy);
for i=1:3
   doa(1,i)=asin(angle(S(i,i))/pi)*180/pi;
end
disp(doa)