pruebaguarda.m

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

function PruebaGuarda(n_mod_tipo,SUI,muestras,BW,figura);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%                                                                        %
%%   Nombre: PruebaGuarda.m                                               %
%%                                                                        %
%%      Funci髇 con la que llamamos a la rutina que simula el sistema     %
%%       con diferentes par醡etros. En este caso, hacemos una variaci髇   %
%%       del tama駉 del prefijo c韈lico para ver c髆o influye en la       %
%%       probabilidad de error de lo transmitido.                         %
%%                                                                        %
%%      Devolvemos una gr醘ica de las diferencias de la simulaci髇        %
%%      con los diferentes tama駉s del prefijo.                           %
%%                                                                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure(figura);

% Los diferentes tama駉s del cyclic prefix que puedo tener
v_G = [1/4 1/8 1/16 1/32];          
v_EbN0_dB=[1:15];
codifico = 1;

% Hago diferentes simulaciones para valores distintos de G
for G = v_G   
    v_ber=[];
    canal = canalSUI(SUI,G,BW);
    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('Guarda',G,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. ( Canal AWGN, BW=',num2str(BW),...
        'MHz y Modulaci髇 de ',modula,' )']);
else
    title(['BER de los s韒bolos recibidos. ( SUI=',num2str(SUI),',BW=',num2str(BW),...
        'MHz y Modulaci髇 de ',modula,' )']);
end

% C醠culo te髍ico de la BER, en funci髇 de la modulaci髇 usada.(De momento s髄o AWGN)
figure(figura);
BERteorica(v_EbN0_dB,n_mod_tipo,SUI);
etiqueta = legend('Guarda=1/4','Guarda=1/8','Guarda=1/16','Guarda=1/32','Te髍ica',...
    'Location','SouthWest');