MUSIC算法估计波达方向

%% MUSIC算法
%%非相干输入信号
%%四信号输入八阵元的线阵
clc;
close all;
clear all;
i=sqrt（-1）;
j=i;
M=8;                                      %阵元数
K=4;                                      %信源数
lma=2;                                    %信号中心波长
d=lma/2;                                  %天线阵元间距
angle1=-40;                               %四个输入信号的入射角
angle2=0; 
angle3=20; 
angle4=60; 
th=[angle1;angle2;angle3;angle4];
degrad=pi/180;
SNR=[10;10;10;10];
%% 构造信号和高斯白噪声
n=1024;
nn=1:n;
v1=0.015;                                                %四个输入信号的频率
v2=0.05;
v3=0.02;
v4=0.035;
A=exp（-2*pi*d*j*[0:M-1]‘*sin（th‘*degrad）/lma）;            %阵列因子
S=[1*cos（v1*nn）;1*sin（v2*nn）;1*sin（v3*nn）;1*sin（v4*nn）];  %输入信号
U=randn（Mn）+j*randn（Mn）;                                %构造噪声源
Ps=diag（S*S‘/n）;