%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