MUSIC算法的DOA估计的matlab仿真程序

clear all;
close all;
clc;
%%设置信号为窄带信号，幅度为1
S=[111]
imag1=sqrt（-1）;
element_num=8;%阵元数为8
d_lamda=1/2;%阵元间距d与波长lamda的关系
theta=linspace（-pi/2pi/22000）;
theta0=0/180*pi;%来波方向 （我觉得应该是天线阵的指向）
%theta0=0;%来波方向
theta1=12/180*pi;
theta2=30/180*pi;
L=20;
kk_sum=0;
SNR=8;
SNR1=16;
SNR2=24;
w_l0=exp（-imag1*2*pi*d_lamda*sin（theta0）*[0:element_num-1]‘）;
w_l1=exp（-imag1*2*pi*d_lamda*sin（theta1）*[0:element_num-1]‘）;
w_l2=exp（-imag1*2*pi*d_lamda*sin（theta2）*[0:element_num-1]‘）;
A=[w_l0w_l1w_l2];
S=randn（length（S）L）
x=A*S;
X=awgn（xSNR‘measured‘）;
X1=awgn（xSNR1‘measured‘）;
X2=awgn（xSNR2‘measured‘）;
R=X*X‘/L;
R1=X1*X1‘/L;
R2=X2*X2‘/L;
[VD]=eig（R）;
[V1D1]=eig（R1）;
[V2D2]=eig（R2）;

Sp=diag（D）‘;
Sp1=diag（D1）‘;
Sp2=diag（D2）‘;
[Sp_randI]=sort（Sp‘descend‘）;
[Sp_rand1I1]=sort（Sp1‘descend‘）;
[Sp_rand2I2]=sort（Sp2‘d