DOA算法中的bartlett算法研究

%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