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