music的MATLAB算法

%MUSIC ALGORITHM
%DOA Estimation by ULA
clear all;
close all;
clc
N_x=1024; % Length of Signal
N=8; % Size of Rx Matrix

l=1.8;%波长
d=0.5*l;%阵元间距
M=3; % Number of Signals
w=[pi/6 pi/10 pi/8]‘;%信号频率
xx=pi/4;yy=0*pi/4;zz=-pi/4; %三个信号的入射角小于5度就难以辨别了
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*s