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: