DOA 估计 music算法matlab

clear all;clc;
%=========参数初始化=========%
m=8;%阵元数
f=10e5;%载波频率
theta1=30;%信号1方向
theta2=35;%信号2方向
th=[theta1;theta2];
nn=1024;%每个阵元采样点数
Ts=[0:1/nn:1-1/nn];
s1=exp（-1i*2*pi*f*Ts）;%信号1
s2=s1;%信号2
S=[s1;s2];
degrad=pi/180;
n=0;
for SN=-10:5:30
%信号1的信噪比
%SN=20;%信号2的信噪比
SN=[500
    0SN];

%--------------------------------------------

%生成阵元接收信号及噪声
%--------------------------------------------
S2=（10.^（SN/10））*S;
wsum=0;
for k=1:50
nr=randn（mnn）;%noise
ni=randn（mnn）;
U=sqrt（2）*（nr+1i*ni）; %每路噪声能量为1.
or=[];
%--------------------------------------------

%协方差矩阵及前后向平滑
%--------------------------------------------
tmp2=[0:m-1]‘; %阵元位置
A=exp（-1i*pi*tmp2*sin（th‘*degrad））; %方向矩阵
X=A*S2+U;
X1=X（1:7:）;X2=X（2:8:）;
Rxx1=X1*X1‘/nn;%第一个子阵
Rxx2=X2*X2‘/nn;%第二个子阵
RXX=（Rxx1+Rxx2）/2;
%--------------------------------------------

%奇异值分解
%--------------------------------------------
[QSSW]=svd（RX