SNR=[-5:1:10];
snr = 10.^(SNR/10); 
result_qpsk = fopen('result_qpsk_viterbi.txt','r');
BER_qpsk = fscanf(result_qpsk,'%g');
fclose(result_qpsk);
theo_qpsk=0.5*erfc(sqrt(snr));
theo_16qam=3*qfunc(sqrt(4*snr/5))/4;
semilogy(SNR,theo_qpsk,'k-')
hold on
semilogy(SNR,theo_16qam,'m-')
hold on
semilogy(SNR,BER_qpsk,'r-s')
xlabel('SNR [dB]')
ylabel('BER')
hold on
grid on
result_16qam = fopen('result_16qam_Viterbi.txt','r');
BER_16qam = fscanf(result_16qam,'%g');
fclose(result_16qam);
semilogy(SNR,BER_16qam,'b-s')
hold on
result_qpsk = fopen('result_qpsk_ViterbiSoft.txt','r');
BER_qpsk = fscanf(result_qpsk,'%g');
fclose(result_qpsk);
semilogy(SNR,BER_qpsk,'r-o')
hold on
result_16qam = fopen('result_16qam_ViterbiSoft.txt','r');
BER_16qam = fscanf(result_16qam,'%g');
fclose(result_16qam);
semilogy(SNR,BER_16qam,'b-o')
hold on
title('BER with Convolution coding in AWGN')
legend('Theoretical QPSK','Theoretical 16QAM','QPSK Hard','16QAM Hard','QPSK Soft','16QAM Soft','location','southwest')