mimo_corr.m

MIMO系统的多径信道仿真程序

function f=mimo_corr(anglespread,angle,d,M)
% 功能：通过波束到达角、角度扩展、天线之间的间隔和天线的个数，计算出发送端或接收端的相关矩阵
% 输入参数说明：
% （1）anglespread散射体的角度扩展，表示接收端和发送端散射体的分布情况。这里假设角度分布服从拉普拉斯分布
% （2）angle平均到达角，每个入射波和离去波的到达角的均值
% （3）d天线间隔与波长之比，表示接收端和发送端两个天线之间的间隔与波长之比，本程序假设天线是均匀线阵
% （4）M天线阵元数，表示接收端和发送端天线阵元个数
% 输出参数说明：mimo_corr 发送端或接收端任意两个天线之间的相关系数矩阵

function f=mimo_corr(anglespread,angle,d,M)
L=1000;
anglespread1=720;
c=0;
% clear i;
p=zeros(1,L);
fai=zeros(1,L);
fai1=zeros(1,L);
FAI=zeros(1,L);
matrix1=zeros(M,1);
matrix2=zeros(1,M);
correlation1=zeros(M,M);
correlation2=zeros(M,M);
correlation=zeros(M,M);

for m=1:L
    fai1(1,m)=angle-anglespread1+2*anglespread1*m/L;
    fai(1,m)=2*pi*(angle-anglespread1+2*anglespread1*m/L)/360;
    FAI(1,m)=d*sin(fai(1,m));
end
for m=1:L
    p(1,m)=1/(anglespread*sqrt(2))*exp(-sqrt(2)*abs(fai1(1,m)-angle)/anglespread)*2*anglespread1/L;
end
for m=1:L
    c=p(1,m)+c;
end
c;
for m=1:L
    for n=1:M
        matrix1(n,1)=exp(i*FAI(m)*2*pi*(n-1));
    end
    matrix2=matrix1';
    correlation1=matrix1*matrix2*p(1,m);
    correlation2=correlation1+correlation2;
end
for m=1:M
    for n=1:M
        correlation(m,n)=abs(correlation2(m,n))/c;
    end
end
f=correlation;