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')
