pruebacodifico.m

802.16-wimax系统物理层的仿真代码

function PruebaCodifico(n_mod_tipo,G,SUI,muestras,BW,figura)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%                                                                        %
%%    Nombre: PruebaCodifico.m                                            %
%%                                                                        %
%%      Funci髇 con la que llamamos a la rutina que simula el sistema     %
%%       con diferentes par醡etros. En este caso, estamos probando si     %
%%       hay cambios en la simulaci髇 si codificamos la informaci髇 a     %
%%       transmitir o si no lo hacemos.                                   %
%%                                                                        %
%%      Devolvemos una gr醘ica de las diferencias de la simulaci髇        %
%%      si hago codificaci髇 o si no la realizo.                          %
%%                                                                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure(figura);

v_EbN0_dB=[1:18];
codifico = 0;           % Indico que voy a empezar por no codificar, 
                        % para estudiar luego la diferencia de hacerlo o no.

canal = canalSUI(SUI,G,BW);

for i=0:1
    v_ber = [];    
    for SNR = v_EbN0_dB
        n_ber = sistema(SNR,n_mod_tipo,G,SUI,codifico,muestras,BW,canal);
        v_ber = [v_ber n_ber];
    end
    dibujar('Codifico',codifico,v_EbN0_dB,v_ber,1);   
    codifico = 1;
end


switch n_mod_tipo
    case 1
        modula = 'BPSK';
    case 2
        modula = 'QPSK';
    case 4
        modula = '16QAM';
    case 6
        modula = '64QAM';
end

if SUI == 0
    title(['BER de los s韒bolos recibidos. ( G=',num2str(G),',BW=',num2str(BW),...
    'MHz,Canal AWGN y Modulaci髇 de ',modula,' )']);
else
    title(['BER de los s韒bolos recibidos. ( G=',num2str(G),',BW=',num2str(BW),...
        'MHz,SUI=',num2str(SUI),' y Modulaci髇 de ',modula,' )']);
end

% C醠culo te髍ico de la BER, en funci髇 de la modulaci髇 usada.(AWGN)
  figure(figura);
  BERteorica(v_EbN0_dB,n_mod_tipo,SUI);
  etiqueta = legend('Sin codificar','Con Codificaci髇','Te髍ica','Location','SouthWest');