mimo_ channel.m

mimo通信系统matlab仿真。 qpsk调制

%带有相关性的MIMO信道仿真程序
% 摘自 杨大成《移动传播环境》附录
% 尹霆 2006 Mar 27

function f= mimo_channel(Nr,Nt,t)

%功能： 产生带有相关性的MIMO信道的信道冲激响应
%输入参数说明：
%  Nr :接收天线阵元个数
%  Nt: 发送天线阵元个数
%  t : 时间变量

% 输出参数说明：
% MIMO信道的信道冲激响应
% 单径信道利用Jakes模型仿真
% 移动速度是3km

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


for i=1:Nt*Nr
    for l=1:1
        h1=0;
        h2=0;
        for k=1:s-1
            sita(k)=2*pi*rand;
            h1=h1+sqrt(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+sqrt(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(2)*sqrt(s-1/2))*cos(2*pi*fd*t+sita(s));
        h2=h2+1/(sqrt(2)*sqrt(s-1/2))*cos(2*pi*fd*t+sita(s));
        h(i,l)=h1+j*h2;
    end
end

corrR=mimo_corr(30,0,0.5,Nr);
% correlation at Rx d= 0.5 lanbuda
corrT=mimo_corr(5,0,5,Nt);
% correlation at Tx d=5 lanbuda
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;