Nt=2;
Nr=3;
t=0.1


%function f=mimo_channel(Nr,Nt,t)
%功能：產(chǎn)生帶有相關(guān)性MIMO信道的信道沖激響應(yīng)
%輸入?yún)?shù)說明：（1）Nr接收天線的陣元的個(gè)數(shù)。（2）Nt發(fā)送天線陣元的個(gè)數(shù)。（3）時(shí)間變量
%輸出參數(shù)說明：mimo_channel MIMO信道的信道沖激響應(yīng)矩陣
%程序：?jiǎn)螐叫诺览肑akes模型仿真。移動(dòng)速度是3km。
% function f=mimo_channel(Nr,Nt,t)
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,3); % correlation at Rx d=0.5lanbuda
corrT=mimo_corr(5,0,5,Nt); % correlation at Tx d=5lanbuda
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;