多重信号分类music）算法计算波达角（不同信噪比的比较分析（matlab）

clear all;
close all;
derad=pi/180;
radeg=180/pi;
twpi=2*pi;
kelm=8; %阵元数
dd=0.5;%阵元间距
d=0:dd:（kelm-1）*dd;
iwave=3;%信号源数
theta=[30 60 120]; %波达方向
snr=[-5051015];
n=200;%采样数
A=exp（j*twpi*d.‘*cos（theta*derad））; %方向矢量
S=randn（iwaven）;
X=A*S;%接收信号
l_snr=length（snr）;
SP_y=zeros（180l_snr）;
for q=1:l_snr
X1=awgn（Xsnr（q）‘measured‘）;%添加噪声
Rxx=X1*X1‘/n;%计算协方差
InvS=inv（Rxx）;
[EVD]=eig（Rxx）; %特征分解
EVA=diag（D）‘;
[EVAI]=sort（EVA）; %特征值降序排序
EVA=fliplr（EVA）; %左右翻转 从大到小顺序
EV=fliplr（EV（:I））; % 对应特征向量排序

%构造MUSIC谱函数
for iang=1:180
    %angle（iang）=（iang-181）/2;
    phim=derad*iang;
    a=exp