Channel_coding.tar.gz_Channel Coding

clc clear all close all % Initialisation des paramètres Nb = 100; G = [1 1 0 1 0 0 0;... 0 1 1 0 1 0 0;... 1 1 1 0 0 1 0;... 1 0 1 0 0 0 1]; n = 7; k = 4; %%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% %%% %%% SANS BRUIT %%% %%% %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Génération d'une suite de bits b = randint(1, Nb); % Ecodage des bits b_encode = encode(b, n, k, 'linear/fmt', G); % Decodage des bits b_decode = decode(b_encode, n, k, 'linear/fmt', G); if (b == b_decode) disp('Pas d''erreur de décodage'); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% %%% %%% AVEC BRUIT %%% %%% %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%% Ne = 100; BER = zeros(1, 8); BER_encode = zeros(1, 8); SNR = 0:1:7; Nb_emis_max = 100000; h = waitbar(0,'Please wait...'); for j=1:1:length(SNR) Nb_emis = 0; erreurs = 0; erreurs_encode = 0; while ((erreurs < Ne) || (erreurs_encode < Ne) && (Nb_emis < Nb_emis_max)) waitbar(j/length(SNR)) % Generation d'une suite de bits b = randint(1, Nb); Nb_emis = Nb_emis + Nb; % Encodage des bits b_encode = encode(b, n, k, 'linear/fmt', G); % Ajout du canal AWGN b_encode_bruite = bsc(b_encode, 1/2*erfc(sqrt(n/k*10^(SNR(j)/10)))); b_bruite = bsc(b, 1/2*erfc(sqrt(10^(SNR(j)/10)))); % Decodage des bits bruités b_decode = decode(b_encode_bruite, n, k, 'linear/fmt', G); % Incrémentation du nombre d'erreurs erreurs_encode = erreurs_encode + sum(b_decode ~= b); erreurs = erreurs + sum(b ~= b_bruite); end BER_encode(j) = erreurs_encode / Nb_emis; BER(j) = erreurs / Nb_emis; end close(h); % Tracé des courbes de probabilité d'erreur dmin = gfweight(G, 'gen'); BER_th_encode = bercoding(n/k*SNR, 'block', 'hard', n, k, dmin); figure, semilogy(SNR, BER_th_encode, 'blue') hold on BER_th = 1/2*erfc(sqrt(10.^(SNR./10))); semilogy(SNR, BER_th, 'm') hold on semilogy(SNR, BER, '--black') hold on semilogy(SNR, BER_encode, '--r') xlabel('SNR (dB)') ylabel('BER') title('Courbe de probabilité d''erreur') h = legend('Courbe théorique avec codage', 'Courbe théorique sans codage', 'Sans codage', 'Avec codage', 1); set(h,'Interpreter','none')