圆阵目标方位估计，music方法matlab

close all;
clear all;
clc;
jay = sqrt（-1）;
A1 = 0.1;                    % 信号幅值
f1 = 1000;                  % 信号频率
fs = 10000;                 % 采样频率
ph1 = 0;                   % 正弦信号的初始相位
T = 1;                   % 采样时间
t = 0:1/fs:T-1/fs;
s1 = A1*exp（jay*（2*pi*f1*t+ph1））;  % 单频信号

M = 8;                     % 阵元数
theta1 = -20;              % 入射角度
Q1=-90;                    % 扫描扇区
Q2=90;
theta_step=0.5;           % 角度扫描步长
Angles=Q1:theta_step:Q2;
SNR = 60;                   % 信噪比
d_lamda = 0.5;             % 阵元间距与信号波长的比值
% % % % % % % % % % % % % % 圆阵
r=1.5;
c=1500;
for kk=0:M-1
    a11（kk+1:）=exp（-jay*2*pi*r*f1*cos（2*pi/M*kk-theta1*pi/180）/c）
end
y=a11*s1;
X = （awgn（y.‘SNR‘measured‘））.‘;
R=1/s