mimo_channel.m

MIMO-OFDM 信道的自适应仿真程序

function f=mimo_channel(Nr, Nt,t)
s=35;      
% mm=O;
fd=5.56;    
rand('state',0);
for i=1 :Nt*Nr
    for l=1:1
       h1=0;
       h2=0;
        for k=l: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+l/(sqrt(2)*sqrt(s-1/2))*cos(2*pi*fd*t+sita(s));
       h(i,1)=h1 +j*h2;
   end
end
h

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