%%?均勻線陣仿真圖

clc;?clear;?close?all;
m=8;??%陣元數(shù)
azimuth=[-55?0?40];?%方位角
p=length（azimuth）;?%信號源個(gè)數(shù)
snap=1000;?%最大快拍數(shù)
drl=0.5;??%圓陣半徑與波長之比
snr=15;?%信噪比dB
snrd=10.^（snr/10）;

%%?構(gòu)造信號

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];

%%?構(gòu)造噪聲

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


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