前向平滑与前后向平滑算法

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;