mimo_corr.asv

MIMO（多输入多输出）信道模型建matlab环境

clear all
clc

Nr=4;
Nt=2;
t=10;
fig=mimo_channel(Nr, Nt,t)%;
function 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;
function 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;