

B=1;
hiper=0;
Nt=1000;
L=16;
M=64;
flag=1;
recep=0;
w=[0 2 4 6 8 10 12 14 16 18];
Channel=1;
FrameLength=1;
v=0;
[variances,Lc]=CIRpowers(Channel);
hfr=[];
for ih=1:Lc+1
  hfr=[hfr;genh(FrameLength,v)];
end
    hfr=diag(variances.^0.5)*hfr;     
    h=hfr;  