%bartlett?algorithm?
clear?all;
j=sqrt（-1）;
P=4;%P?number?of?signal?source
N=8;%N?array?number
L=100;%L?length?of?snap
f=9*10^8;%載頻
c=3*10^8;%
d=0.5;%陣元間距
doa1=13/180*pi;%DOA?of?source?1
doa2=30/180*pi;%DOA?of?source?2
doa3=50/180*pi;%DOA?of?source?3
doa4=85/180*pi;%DOA?of?source?4
doa=[doa1doa2doa3doa4];
snr_in_db=6;
snr=10.^（snr_in_db./10）;
AM=[snrsnrsnrsnr];%信源的幅值向量
DOA=sin（doa）;
for?k=1:N
???for?m=1:P
??????A（km）=exp（j*2*pi*f*d*（k-1）*DOA（m）/c）;%A?denotes?arry?response?matrix
???end
end
randn（‘state‘sum（100*clock））;
pn=randn（1P*L）;
count=1;
for?k=1:P
???for?m=1:L
??????s（km）=AM（1k）*p