相干目标music算法

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