root_music-mrse.m

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

%root_MUSIC;  
clear all; close all;clc;
J=sqrt(-1);
source_number=2;
source_doa=[40 70];
sensor_number=13;
snapshot_number=1000;
 

mot=200;
k=0;
ef1=zeros(1,31);
 
A=exp(-J*(0:sensor_number-1)'*pi*sin(source_doa*pi/180));

for snr=-10:1:20
       k=k+1;
       ef1=zeros(1,k);
       ef2=zeros(1,k);
    for j=1:1:mot 
        
        
s=(rand(source_number,snapshot_number)+J*rand(source_number,snapshot_number))/sqrt(2);

x=A*s;
y=awgn(x,snr);

R=y*y'/snapshot_number;
 
[V,D]=eig(R);
Un=V(:,1:sensor_number-source_number);
Un1=Un(1:source_number,:);
Un2=Un(source_number+1:sensor_number,:);

T=[1,0,0,0,0,0,0,0,0,0,0]';
c=Un1*inv(Un2)*T;


c=[1,c(2,1),c(1,1)];
f=roots(c);
esource_doa=asin(angle(f)/pi)*180/pi;
esource_doa=sort(esource_doa);
disp('esource_doa');
disp(esource_doa);
 

f1=(esource_doa(1,1)-source_doa(1,1))^2;
f2=(esource_doa(2,1)-source_doa(1,2))^2;
ef1(1,k)=f1+ef1(1,k);
ef2(1,k)=f2+ef2(1,k);

end

ef1(1,k)=sqrt(ef1(1,k)/mot);
ef2(1,k)=sqrt(ef1(1,k)/mot);
end

snr=-10:1:20;
plot(snr,ef1,'-diamond',snr,ef2,'-square');
title('doa Estimation Performance');
xlabel('snr');
ylabel('RMSE(度)');
legend('source1','source2');