musicfromula.m

一美国佬编的经典MUSIC算法

%MUSIC ALGORITHM
%DOA Estimation by ULA
clear; 
close all; 
N_x=1024; % Length of Signal 
N=8; % Size of Rx Matrix 
A=[4 3 3];
l=1.8;%波长
d=0.4*l;%阵元间距
w=[200 200 200]';%信号频率
xx=pi/4;yy=1.5*pi/4;zz=3*pi/4; %三个信号的入射角
% B=[exp(-j*2*[0:7]*pi*d*sin(xx)/l);exp(-j*2*[0:7]*pi*d*sin(yy)/l);exp(-j*2*[0:7]*pi*d*sin(zz)/l)]';
B=[1 exp(-j*2*pi*d*sin(xx)/l) exp(-j*2*2*pi*d*sin(xx)/l) exp(-j*2*3*pi*d*sin(xx)/l) exp(-j*2*4*pi*d*sin(xx)/l) exp(-j*2*5*pi*d*sin(xx)/l) exp(-j*2*6*pi*d*sin(xx)/l) exp(-j*2*7*pi*d*sin(xx)/l);...
1 exp(-j*2*pi*d*sin(yy)/l) exp(-j*2*2*pi*d*sin(yy)/l) exp(-j*2*3*pi*d*sin(yy)/l) exp(-j*2*4*pi*d*sin(yy)/l) exp(-j*2*5*pi*d*sin(yy)/l) exp(-j*2*6*pi*d*sin(yy)/l) exp(-j*2*7*pi*d*sin(yy)/l);...
1 exp(-j*2*pi*d*sin(zz)/l) exp(-j*2*2*pi*d*sin(zz)/l) exp(-j*2*3*pi*d*sin(zz)/l) exp(-j*2*4*pi*d*sin(zz)/l) exp(-j*2*5*pi*d*sin(zz)/l) exp(-j*2*6*pi*d*sin(zz)/l) exp(-j*2*7*pi*d*sin(zz)/l)]'; %阵列流型
M=3; % Number of Signals 
xxx=10*exp(j*w*[0:N_x-1]) %仿真信号
x=randn(8,N_x)+B*xxx; %加噪声
R=1/1024*x*x';
[V D]=eig(R);
[lambda,index] = sort(abs(diag(D)));
UU=V(:,1:5);
for theta=0:.005:pi
AA=[1 exp(-j*2*pi*d*sin(theta)/l) exp(-j*2*2*pi*d*sin(theta)/l) exp(-j*2*3*pi*d*sin(theta)/l) exp(-j*2*4*pi*d*sin(theta)/l) exp(-j*2*5*pi*d*sin(theta)/l) exp(-j*2*6*pi*d*sin(theta)/l) exp(-j*2*7*pi*d*sin(theta)/l)];%方向矢量
WW=AA*UU*UU'*AA';
P=1/WW;%角谱
hold on
plot(theta,P)
end
grid
xlabel('rad')