Sensor=8; Snap=100; 
Doa=[10 20 50 60]; 
S_number=length(Doa); 
Lambda=0.1; > 波長 
d=0.5*Lambda; > 陣元間距 
S_x=[0:Sensor-1]*d; A=exp(-j*2*pi/Lambda*S_x'*sin(Doa*pi/180)); 
Snr=10; 
Snr=sqrt(10.^(Snr/10)); > ？？？？ 
Sig=sqrt(0.5)*(randn(S_number,Snap)+j*randn(S_number,Snap)); > 信號 
Noise=sqrt(0.5)*(randn(Sensor,Snap)+j*randn(Sensor,Snap)); > 噪聲 
X=Snr*A*Sig+Noise; 
R=X*X'/Snap; 
[U S V]=svd(R); >特征值分解 
En=U(:,S_number+1:Sensor); 
theta=[-90:0.1:90]*pi/180; 
for num=1:length(theta) 
A_search=exp(-j*2*pi/Lambda*S_x'*sin(theta(num))); 
P_music(num)=abs(1/(A_search'*En*En'*A_search)); 
end 
h=[-90:0.1:90]; 
plot(h,10*log10(P_music/max(P_music))); 
axis([-100 100 -50 +10]) 
xlabel('????????????/??¨????????','fontsize',10) 
ylabel('Pmusic','fontsize',10) 
title({'????????music??????¨'}); 
grid