clear?all
close?all
ima=sqrt（-1）;
var=sqrt（1/exp（1.0））;
snr=10;
zy=16;???????%陣元數
k=1/2;?????????%陣元間距與波長的關系
theta=-90:0.01:90;??????%范圍????????
theta1=10;???????????%來波方向
theta2=15;???????????%干擾方向
theta3=30;
theta4=40;
theta5=50;
L=200;??????????????%快拍數
?%產生信號是兩個獨立信號源
?for?i=1:L;
??????amp1=randn（1）
??????amp2=randn（1）
??????amp3=randn（1）
??????amp4=randn（1）
??????amp5=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]‘）+...
?????amp3*exp（ima*2*pi*1/2*sin（theta