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