多重信号分类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=10;

A=exp（j*twpi*d.‘*cos（theta*derad））; %方向矢量

n=[10 50 100 200 400];;%采样数
l_n=length（n）;

SP_y=zeros（180l_n）;

for q=1:l_n
    S=randn（iwaven（q））;
X=A*S;%接收信号
X1=awgn（Xsnr‘measured‘）;%添加噪声
Rxx=X1*X1‘/n（q）;%计算协方差
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;