doa算法的matlab实现

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）=