music1.txt

用C语言实现对经典的MUSIC算法的仿真

clear;
close al
%----产生数据-------------
signal_a=sqrt(10^(snr/10)); 
t=[0:1:(L-1)];
f1 = 0.21;S1=sin(2*pi*f1*t);
f2 = 0.33;S2=sin(2*pi*f2*t);
S=[S1;S2];
A=exp(i*2*pi*ddd*(0:M-1)'*sin(signal_doa*pi/180));
X=A*S*signal_a;
R=X*X'/L;
[U D]=eig(R);
US=U(:,1:M-N);
UN=U(:,M-N+1:M);
ddoa=-90:0.1:90;
steer=exp(i*2*pi*ddd*(0:M-1)'*sin(ddoa*pi/180));
for ij=1:length(ddoa)
aa=steer(:,ij);
p1(ij)=1/(aa'*UN*UN'*aa); 
end
%----画图得到角度-------------
figure;box on;grid on;
plot(ddoa,20*log10(abs(p1)));
legend('MUSIC')