四种传统算法进行DOA估计

clear all;
close all;
ima=sqrt（-1）;
var=sqrt（1/exp（1.0））;
zy=8;       %阵元数
k=1/2;         %阵元间距与波长的关系
theta=-90:0.01:90;      %范围        
theta1=5;           %来波方向
theta2=15;           %干扰方向
L=100;              %快拍数
u0=[1 0 0 0 0 0 0 0]‘;%MEM

for i=1:L;
   %产生信号是幅度随机的，保证信号之间是不相关的既是两个独立信号源
%     amp1=randn（1）;
%     amp2=randn（1）;
%     x（:i）=amp1*exp（ima*2*pi*1/2*sin（theta1*pi/180）*[0:zy-1]‘）+...
%       amp2*exp（ima*2*pi*1/2*sin（theta2*pi/180）*[0:zy-1]‘）+...
%      var*randn（zy1）;
 %产生信号是两个相干信号源
 x（:i）=exp（ima*2*pi*1/2*sin（theta1*pi/180）*[0:zy-1]‘）+...
      exp（ima*2*pi*1/2*sin（theta2*pi/180）*[0:zy-1]‘）+...
     var*randn（zy1）;
 
end    
% 产生数据协方差矩阵MVM/MEM
Rx=1/L*x* x‘;
R=inv（Rx）;
% 对数据协方差矩阵进行特征值分解MNM
[v d]=eig（Rx）;
Un=v（:[3 4 5 6 7 8]）;
Us=v（:[1 2]）;
h=Us（1:）;
g=h.‘;
Es=Us（[2 3 4 5 6 7 8]:）;
y=Un（1:）;
c=y.‘;
En=Un（[2 3 4 5 6 7 8]:）;
a=En*y‘;
b=c‘*c;
w3=a/b;
w33=[1;w3];
% 产生Rfb rfb 矩阵 FBLP
  x1=x（1:）;x2=x（2:）;x3=x（3: