clear?all;??
close?all;?
clc
M=8;?????%陣元數
P=3;?????%信號源數
K=256;???%抽樣點數
n=4;?????%分成4個子陣
m=5;?????%每個子陣包括5個陣元
Z=[?-30*pi/180??10*pi/180?25*pi/180];???%3個信號的波束來向
w=[pi/5?pi/5?pi/5];???????????%信號頻率?
%SNR=[10?20?5];?????????????????%3個信號的信噪比
SNR=20;
S_=[1?0.9?0.8];?????????????????%相干信號衰減幅度
%%---------------------------------------------------------------
T_1=zeros（mm）;
T=zeros（mm）;
for?l=1:100?????????????????????%蒙特卡羅
R=zeros（mm）;
for?i=1:n
?A=zeros（mP）;???????????????????%空間導向矢量
?for?ii=1:m
????for?jj=1:P
????????A（iijj）=exp（-j*pi*（ii-1）*sin（Z（jj）））;
????end
?end
?D=zeros（PP）;
?for?ii=1:P
????D（iiii）=exp（j*pi*sin（Z（ii）））;
?end
?S=zeros（PK）;???%信號矩陣
?for?ii=1:P
????for?jj=1:K
??????S（iijj）=S_（ii）*sqrt（10^（SNR/10））*exp（j*5*pi*（jj-1））;
????end
?end
?%接收到的數據矩陣
?X=A*D^（i-1）*S+randn（mK）+j*randn（mK）;
?%接收信號的協方差矩陣
??R=R+X*X‘;
end
R=R/n;?????%前向平滑數據協方差矩陣

a=eye（5）;
b=a（:5:-1:1）;??%構造置換矩陣
RXX=（R+b*R.‘*b）/2;?????%前后向平滑接收數據矩陣
end;