mimo_channel.m

该程序实现了mimo的带有相关性的Rayleigh信道模型

function f=mimo_channel(Nr,Nt,fd,t)
% 接收天线个数 发射天线个数 最大多普勒频移 时间变量
% 输出 mimo_channel MIMO信道的信道冲击响应。
% 单径信道利用 Jakes 模型仿真。
% 移动速度是3km。

s=35;
%mm=0;
rand('state',0);

for i=1:Nt*Nr
    for l=1:1  % l:多径可分辨独立径数
        h1=0;h2=0;
        for k=1:s-1
            sita(k)=2*pi*rand;
            h1=h1+2/sqrt(s-1/2)*sin(pi*k/(s-1))*cos(2*pi*fd*cos(pi*k/(2*s-1))*t+sita(k));
            h2=h2+2/sqrt(s-1/2)*cos(pi*k/(s-1))*cos(2*pi*fd*cos(pi*k/(2*s-1))*t+sita(k));
        end
        sita(s)=rand;
        h1=h1+1/sqrt(s-1/2)*cos(2*pi*fd*t+sita(s));
        h2=h2+1/sqrt(s-1/2)*cos(2*pi*fd*t+sita(s));
        h(i,l)=h2+j*h1;
    end
end
corrR=mimo_corr(30,0,0.5,Nr);
corrT=mimo_corr(5,0,5,Nt);
corrRT=kron(corrR,corrT);
hr=transpose(chol(corrRT));
h=hr*h;
for p=1:Nr
    for q=1:Nt
        hh(p,q)=h(Nr*(q-1)+p);
    end
end
f=hh;