到达角估计MUSIC算法

% MUSIC算法 
clc
clear all
close all
M=10;                                         %阵元数                                         
SNR=5;                                        %信噪比
theta=[5 30 55 80]*pi/180;                    %信号入射角
N=8;                                          %快拍数
f=3000;                                       %信号频率
fs=f*5;                                       %采样频率
c=1500;                                       %波速
lambda=c/f;                                   %波长
d=lambda;                                   %阵元间距
rannumber=1000;
%%%%%%%%%%%%%%%% 生成发射信号和接收信号 %%%%%%%%%%%%%%%%%

S1=sqrt（10^（20/10））*exp（1i*2*pi*f/fs+2*pi*rand（1rannumber））;     %独立的信号源1
S2=sqrt（10^（15/10））*exp（1i*2*pi*2*f/fs+2*pi*rand（1rannumber））;   %独立的信号源2
S3=sqrt（10^（16/10））*exp（1i*2*pi*3*f/fs+2*pi*rand（1rannumber））; 
S4=sqrt（10^（17/10））*exp（1i*2*pi*4*f/fs+2*pi*rand（1rannumber））; 
S=[S1;S2;S3;S4];

A1=[exp（-2*pi*1i*d*[0:M-1]*cos（theta（1））/lambda）].‘;
A2=[exp（-2*pi*1i*d