基于虚拟阵列的 MUSIC算法 MATLAB

clc;
clear all;
f=1*10^8;
c=3*10^8;
lamda=c/f;
d=lamda;
n=9;%实际
M=18; %虚拟阵元
thita1=20;
thita2=35;
thita=-90:1:90;
N=2000;
S1=3*exp（j*pi/6*[1:N]）;%来向信号
S2=3*exp（j*pi/4*[1:N]）;
S=[S1;S2];
Rs=S*S‘/N;
 A1=exp（-j*2*pi*d*[0:n-1]*sin（thita1*pi/180）/lamda）‘;
 A2=exp（-j*2*pi*d*[0:n-1]*sin（thita2*pi/180）/lamda）‘;
 A=[A1A2];%实际阵元矩阵
 RS=A*Rs*A‘;
 R1=awgn（RS10）;%实际数据协方差
 jiao=[0:1:45];
 for b=1:46;
     for nl=1:n
      As（nlb）=exp（-j*2*pi*d*（nl-1）*sin（jiao（b）*pi/180）/lamda）‘;%实际的区域阵元
     end
     for ml=1:M
     Ax（mlb）=exp（-j*2*pi*d/2*（ml-1）*sin（jiao（b）*pi/180）/lamda）‘;%虚拟的区域阵元
     end
 end    
 B