clear all;
clf;
echo on
SNRindB1=5:1:15;
for k=1:length(SNRindB1),
    [pban1_1,psan1_1]=bperuncode8psk(SNRindB1(k),10000);
    smld_bit_err_prban1_1(k)=pban1_1;
    smld_symbol_err_prban1_1(k)=psan1_1;
    [pban2_1,psan2_1]=bper8pskG2(SNRindB1(k),10000);
    smld_bit_err_prban2_1(k)=pban2_1;
    smld_symbol_err_prban2_1(k)=psan2_1;
    [pban3_1,psan3_1]=bper16QAMH3(SNRindB1(k),10002);
    smld_bit_err_prban3_1(k)=pban3_1;
    smld_symbol_err_prban3_1(k)=psan3_1;
    [pban4_1,psan4_1]=bper16QAMH4(SNRindB1(k),10002);
    smld_bit_err_prban4_1(k)=pban4_1;
    smld_symbol_err_prban4_1(k)=psan4_1;
end;
SNRindB2=16:1:20;
for k=1:length(SNRindB2),
    [pban1_2,psan1_2]=bperuncode8psk(SNRindB2(k),50000);
    smld_bit_err_prban1_2(k)=pban1_2;
    smld_symbol_err_prban1_2(k)=psan1_2;
    [pban2_2,psan2_2]=bper8pskG2(SNRindB2(k),50000);
    smld_bit_err_prban2_2(k)=pban2_2;
    smld_symbol_err_prban2_2(k)=psan2_2;
    [pban3_2,psan3_2]=bper16QAMH3(SNRindB2(k),50001);
    smld_bit_err_prban3_2(k)=pban3_2;
    smld_symbol_err_prban3_2(k)=psan3_2;
    [pban4_2,psan4_2]=bper16QAMH4(SNRindB2(k),50001);
    smld_bit_err_prban4_2(k)=pban4_2;
    smld_symbol_err_prban4_2(k)=psan4_2;
end;
SNRindB3=23:3:41;
for k=1:length(SNRindB3),
    [pban1_3,psan1_3]=bperuncode8psk(SNRindB3(k),100000);
    smld_bit_err_prban1_3(k)=pban1_3;
    smld_symbol_err_prban1_3(k)=psan1_3;
    [pban2_3,psan2_3]=bper8pskG2(SNRindB3(k),100000);
    smld_bit_err_prban2_3(k)=pban2_3;
    smld_symbol_err_prban2_3(k)=psan2_3;
    [pban3_3,psan3_3]=bper16QAMH3(SNRindB3(k),100002);
    smld_bit_err_prban3_3(k)=pban3_3;
    smld_symbol_err_prban3_3(k)=psan3_3;
    [pban4_3,psan4_3]=bper16QAMH4(SNRindB3(k),100002);
    smld_bit_err_prban4_3(k)=pban4_3;
    smld_symbol_err_prban4_3(k)=psan4_3;
end;
SNRindB=[SNRindB1 SNRindB2 SNRindB3];
smld_bit_err_prban1=[smld_bit_err_prban1_1 smld_bit_err_prban1_2 smld_bit_err_prban1_3];
smld_bit_err_prban2=[smld_bit_err_prban2_1 smld_bit_err_prban2_2 smld_bit_err_prban2_3];
smld_bit_err_prban3=[smld_bit_err_prban3_1 smld_bit_err_prban3_2 smld_bit_err_prban3_3];
smld_bit_err_prban4=[smld_bit_err_prban4_1 smld_bit_err_prban4_2 smld_bit_err_prban4_3];
smld_symbol_err_prban1=[smld_symbol_err_prban1_1 smld_symbol_err_prban1_2 smld_symbol_err_prban1_3];
smld_symbol_err_prban2=[smld_symbol_err_prban2_1 smld_symbol_err_prban2_2 smld_symbol_err_prban2_3];
smld_symbol_err_prban3=[smld_symbol_err_prban3_1 smld_symbol_err_prban3_2 smld_symbol_err_prban3_3];
smld_symbol_err_prban4=[smld_symbol_err_prban4_1 smld_symbol_err_prban4_2 smld_symbol_err_prban4_3];
echo on;
semilogy(SNRindB,smld_bit_err_prban1,'k');
hold on;
semilogy(SNRindB,smld_bit_err_prban1,'k+');
hold on;
semilogy(SNRindB,smld_bit_err_prban2,'k');
hold on;
semilogy(SNRindB,smld_bit_err_prban2,'k*');
hold on;
semilogy(SNRindB,smld_bit_err_prban3,'k');
hold on;
semilogy(SNRindB,smld_bit_err_prban3,'kx');
hold on;
semilogy(SNRindB,smld_bit_err_prban4,'k');
hold on;
semilogy(SNRindB,smld_bit_err_prban4,'ko');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban1,'k');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban1,'k+');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban2,'k');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban2,'k*');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban3,'k');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban3,'kx');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban4,'k');
% hold on;
% semilogy(SNRindB,smld_symbol_err_prban4,'ko');
grid on;
xlabel('SNR(dB)');
ylabel('Bit Error Probability');
axis([5,40,10^(-6),1]);