mimo_channel.m

来自「nv_encd.m 卷积编码程序 viterbi.m 卷积译码程序 其它的」· M 代码 · 共 35 行

function f=mimo_channel(Nr,Nt,t)
Np=1;
s=35;
%mm=0
fd=200;
%rand('state',0);
hq=zeros(Nr*Nt,128);
for i=1:Nt*Nr
    for l=1:Np
        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,:)=h1+j*h2;
    end
end

corrR=mimo_corr(30,0,0.5,Nr);
corrT=mimo_corr(5,0,10,Nt);
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;