pruebabw.m

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

function PruebaBW (G,SUI,n_mod_tipo,muestras,figura)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%                                                                        %
%%    Nombre: PruebaBW.m                                                  %
%%                                                                        %
%%      Funci髇 con la que llamamos a la rutina que simula el sistema     %
%%       con diferentes par醡etros. En este caso, hacemos una variaci髇   %
%%       del valor del ancho de banda disponible en el sistema y          %
%%       estudiamos su efecto.                                            %
%%                                                                        %
%%      Devolvemos una gr醘ica de las diferencias de la simulaci髇        %   
%%      al usar diferentes mapeos de los bits a transmitir.               %
%%                                                                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure(figura);

% Algunos de los BW disponibles en la norma.
v_BW = [28 20 15 10 2.50 1.25];
v_EbN0_dB=[1:20];
codifico = 1;

% Realizo la simulaci髇 para los diferentes anchos de banda disponibles.
for BW = v_BW         
    canal = canalSUI(SUI,G,BW);
    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('BW',BW,v_EbN0_dB,v_ber,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),',Canal AWGN',...
        ' y Modulaci髇 de ',modula,' )']);
else
    title(['BER de los s韒bolos recibidos. ( G=',num2str(G),',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('BW.-28 MHz','BW.-20 MHz','BW.-15 MHz','BW.-10 MHz',...
   'BW.-2.5 MHz','BW.-1.25 MHz','Te髍ica','Location','SouthWest');