% AWGN 信道,維特比譯碼,(8,4)漢明碼糾錯,分組交織
Eb_to_Nj_in_dB_Q1=[0,1,2,3,4,5,6,7,8,9,10,11,12];
Pb_Q1=[0.48235,0.48226,0.47974,0.47756,0.47465,0.4707,0.44928,0.41736,0.27968,8.4056e-2,1.0397e-2,8.42e-4,6.65e-5];
Eb_to_Nj_in_dB_Q4=[0,1,2,3,4,5,6,7,8,9,10];
Pb_Q4=[0.47565,0.48111,0.47063,0.45224,0.41212,0.32544,0.13835,0.0268,2.8215e-3,2.145e-4,2.4e-5];
Eb_to_Nj_in_dB_Q8=[0,1,2,3,4,5,6,7,8,9];
Pb_Q8=[0.4675,0.45944,0.44475,0.39577,0.26715,9.013e-2,2.9861e-2,8.0725e-3,8.68e-4,6.45e-5];
Eb_to_Nj_in_dB_Q16=[0,1,2,3,4,5,6,7,8,9];
Pb_Q16=[0.45146,0.43302,0.23829,0.1871,0.1007,4.4843e-2,9.576e-3,2.5685e-3,2.125e-4,1.55e-5];

% Draw Pb vs.Eb/Nj curves
semilogy(Eb_to_Nj_in_dB_Q1,Pb_Q1,':.',Eb_to_Nj_in_dB_Q4,Pb_Q4,':o',Eb_to_Nj_in_dB_Q8,Pb_Q8,':+',Eb_to_Nj_in_dB_Q16,Pb_Q16,':x');
grid on;
title('AWGN信道下(8,4)漢明碼性能曲線');
xlabel('Eb/Nj (dB)');
ylabel('誤碼率 Pb');