clear?all
N=100;%采樣數
snr=[10?20]‘;
n=[0:N-1]‘;
s1=randn（N1）;
s2=randn（N1）;
M=10;%陣元數；
p=2;?%?信號數
st=[-1040]*pi/180;
for?t=1:p
????A（:t）=exp（-j*pi*[0:M-1]‘*sin（st（t）））;
end
noise=randn（MN）;
S=[s1?s2]‘;
Ps1=S*S‘/N;
ps=diag（Ps1）;%單個信號的能量
refp=10.^（snr/10）;
tmp=sqrt（refp./ps）;
S2=diag（tmp）*S;
X=A*S2+noise;
%MUSIC
R=X*X‘/N;
[U?D]=eig（R）;
[dindex]=sort（diag（D））;
U=U（:index）;
for?i=-90:90
????a=exp（-j*pi*[0:M-1]‘*sin（i*pi/180））;
????Pmusic（i+91）=abs（1/（a‘*U（:1:end-2）*U（:1:end-2）‘*a））;
end
plot（[-90:90]10*log10（Pmusic/max（Pmusic）））;hold?on
%root-Muisc
z?=sym（‘z‘‘unreal‘）;
for?i=1:10
pz（i）=z^（（i-1））;
pzt（i）=