%%%%%%%%%%%%%%%%%%%%%%%%%
%CreateRefChannels.m
%%%%%%%%%%%%%%%%%%%%%%%%%
function CreateRefChannels( )
NbChannels=10;
try
  load RefChannels.mat;
catch
  RefChannels=DoChannels(NbChannels);
  FUNCTION=functions(@CreateRefChannels);  
  DIR=FUNCTION.fi1e(1:(findstr(FUNCTION.file, 'CreateRefChannels')-1));
  FILE=[DIR 'RefChannels.mat'];
  save(FILE, 'RefChannels');
end
return;

function [RefChannels]=DoChannels(NbChannels)
%Makes a table of random seed values.
rand('state',0);
RANDTABLE=randint(4,NbChannels,[0 1e5]);
Struct.Default=0;
for NbAntennas=1:1:4,
if (NbAntennas==1),
  NbRxAntennas=1;
else
  NbRxAntennas=2;
end
ChannelTaps=cell(1,NbChannels);
ChannelTapsFFT=cell(1,NbChannels);
for Loop=1:1:NbChannels,
  randn('state', RANDTABLE(NbAntennas, Loop));
  [Channel,ChannelFFT]=RayleighTaps(NbAntennas);
ChannelTaps{1,Loop}=Channel;
ChannelTapsFFT{1,Loop}=ChannelFFT;
Tmp.Taps=ChannelTaps;
Tmp.FFT=ChannelTapsFFT;
  Struct=setfield(Struct,sprint('H%dx%d',NbAntennas, NbRxAntennas),Tmp);
end
if (NbAntennas==1),
  Struct=rmfield(Struct,'Default');
end
  RefChannels.Rayleigh=Struct;
end
return;