mimo_channel.m

计算信道容量的matlab程序

function [corrRT]=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
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)%;