% 瑞利衰落信道,維特比譯碼,性能曲線,無糾錯編碼
Eb_to_Nj_in_dB_Q1=[0,3,6,9,12,15,18,21,24,27,30,33];
Pb_Q1=[0.4739,0.4499,0.3806,0.2256,7.04e-2,1.44e-2,2.8e-3,6e-4,1.665e-4,8.2e-5,6.25e-5,5.05e-5];
Eb_to_Nj_in_dB_Q2=[0,3,6,9,12,15,18,21,24,27,30,33];
Pb_Q2=[0.4588,0.4163,0.3011,0.1300,3.05e-2,5.7e-3,1.2e-3,3.07e-4,1.19e-4,8.10e-5,5.4e-5,5e-5];
Eb_to_Nj_in_dB_Q4=[0,3,6,9,12,15,18,21,24,27,30,33];
Pb_Q4=[0.4495,0.3726,0.2177,6.63e-2,1.29e-2,2.52e-3,5.67e-4,2.1e-4,1.01e-4,6.9e-5,5.3e-5,4.6e-5];
% Eb_to_Nj_in_dB_Q8=[0,3,6,9,12,15,18,21,24,27];
% Pb_Q8=[0.4124,0.2937,0.1249,3.06e-2,6.16e-3,1.34e-3,3.39e-4,1.6e-4,7.3e-5,5.9e-5];

% Draw Pb vs.Eb/Nj performance curves
semilogy(Eb_to_Nj_in_dB_Q1,Pb_Q1,':.',Eb_to_Nj_in_dB_Q2,Pb_Q2,':o',Eb_to_Nj_in_dB_Q4,Pb_Q4,':x');%Eb_to_Nj_in_dB_Q8,Pb_Q8,':+',
grid on;
title('瑞利衰落信道下無糾錯碼的性能曲線');
xlabel('Eb/Nj (dB)');
ylabel('誤碼率 Pb');
%legend('單音干擾 Q=1','多音干擾 Q=2','多音干擾 Q=4');