8阵元均匀圆阵music算法

%% 均匀线阵仿真图

clc; clear; close all;
m=8;  %阵元数
azimuth=[-55 0 40]; %方位角
p=length（azimuth）; %信号源个数
snap=1000; %最大快拍数
drl=0.5;  %圆阵半径与波长之比
snr=15; %信噪比dB
snrd=10.^（snr/10）;

%% 构造信号

phi=[10 20 50];
f0=2;
fs=1000;
f=[400 600 500];
for t=1:snap;
    s1（t）=3*cos（2*pi*f0/fs*t）*exp（j*f（1）/fs*t+phi（1）*pi/180）;
    s2（t）=5*exp（j*f（2）/fs*t+phi（2）*pi/180）;
    s3（t）=2*exp（j*f（3）/fs*t+phi（3）*pi/180）;
end
S=[s1;s2;s3];

%% 构造噪声

N=（randn（msnap）+j*randn（msnap））/sqrt（2）;


%% 阵列流型矩阵
k=[0:m-1];
for C=1:3
for R=1:8
    A（RC）=exp（-j*