function canal = canalSUI(N_SUI,G,BW)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%                                                                        %
%%     Archivo: canalSUI.m                                                %
%%                                                                        %
%%     Descripci髇: Generamos la respuesta impulsiva del canal de un      %
%%     canal variante en el tiempo de acuerdo con el modelo de canal de   %
%%     Jakes.                                                             %
%%      El canal usado depende de los par醡etros que se le indiquen.      %
%%     Podemos simular los canales SUI 1 al 6, con diferentes anchos de   %
%%     banda.                                                             %
%%                                                                        %
%%     Par醡etros:                                                        %
%%      N_SUI : Canal a simular.        G = Tama駉 del prefijo c韈lico    %
%%      v = Velocidad del sistema.      BW = Ancho de banda del canal     %
%%                                                                        %
%%                                                                        %
%%     Autores: Bertrand Muquet, Sebastien Simoens, Shengli Zhou          %
%%      Octubre 2000                                                      %
%%     Modificaci髇 : Carlos Batll閟 - Abril 2007                         %
%%                                                                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Velocidad del receptor 0.001 m/s
v = 0.001;                      

% Vamos a considerar que nosotros estamos mandando un 鷑ico s韒bolo
FrameLength = 1;                 

% Los par醡etros siguientes son con los que calculamos la duraci髇 del s韒bolo en WiMAX
Nfft = 256;
BW = BW*1e6;

% Factor de correci髇 segun indica la norma
if mod(BW,1.75)==0
    n = 8/7;
elseif mod(BW,1.5)==0
    n = 86/75;
elseif mod(BW,1.25)==0
    n = 144/125;
elseif mod(BW,2.75)==0
    n = 316/275;
elseif mod(BW,2)==0
    n = 57/50;
else 
    n = 8/7;
end

if N_SUI~=0
    Fs = floor(n*BW/8000)*8000;         % Frecuencia de muestreo
    deltaF = Fs / Nfft;                 % Separaci髇 entre frecuencias.
    Tb = 1/deltaF;                      % Tiempo del s韒bolo (s髄o datos)
    Ts = Tb * (1+G);                    % Tiempo del s韒bolo (Datos + Prefijo c韈lico)
    T = 1/(Fs*1e-6);                    % Duraci髇 en microsegundos de cada portadora


    [variances,Lc,Dop]=CIRpowers(N_SUI,T);

    hfr=[];
    for ih=1:Lc+1
      hfr=[hfr;genh(FrameLength,v,Dop,Ts)];     
    end
    hfr=diag(variances.^0.5)*hfr;
    %% hfr tiene un tama駉 de (Lc+1)x(FrameLength)
    %% hfr(:,i) contiene la CIR (Respuesta al impulso del canal)
    %% correspondiente a la transmisi髇 del simbolo i

    % Por 鷏timo normalizo los valores del canal
    canal = hfr ./ norm(hfr);
    
elseif N_SUI == 0
    canal = 1;         % Si eligo un canal AWGN, el canal sera la unidad.
end