clear all;
clf;
echo on
SNRindB1=-1:1:4;
for k=1:length(SNRindB1),
    [pban1_1,psan1_1]=bperuncode2pskm2(SNRindB1(k),10000);
    smld_bit_err_prban1_1(k)=pban1_1;
    smld_symbol_err_prban1_1(k)=psan1_1;
    [pban2_1,psan2_1]=bper2pskG2m2(SNRindB1(k),10000);
    smld_bit_err_prban2_1(k)=pban2_1;
    smld_symbol_err_prban2_1(k)=psan2_1;
    [pban3_1,psan3_1]=bper4pskG3m2(SNRindB1(k),10000);
    smld_bit_err_prban3_1(k)=pban3_1;
    smld_symbol_err_prban3_1(k)=psan3_1;
    [pban4_1,psan4_1]=bper4pskG4m2(SNRindB1(k),10000);
    smld_bit_err_prban4_1(k)=pban4_1;
    smld_symbol_err_prban4_1(k)=psan4_1;
end;
SNRindB2=5:1:17;
for k=1:length(SNRindB2),
    [pban1_2,psan1_2]=bperuncode2pskm2(SNRindB2(k),50000);
    smld_bit_err_prban1_2(k)=pban1_2;
    smld_symbol_err_prban1_2(k)=psan1_2;
    [pban2_2,psan2_2]=bper2pskG2m2(SNRindB2(k),50000);
    smld_bit_err_prban2_2(k)=pban2_2;
    smld_symbol_err_prban2_2(k)=psan2_2;
    [pban3_2,psan3_2]=bper4pskG3m2(SNRindB2(k),50000);
    smld_bit_err_prban3_2(k)=pban3_2;
    smld_symbol_err_prban3_2(k)=psan3_2;
    [pban4_2,psan4_2]=bper4pskG4m2(SNRindB2(k),50000);
    smld_bit_err_prban4_2(k)=pban4_2;
    smld_symbol_err_prban4_2(k)=psan4_2;
end;
SNRindB=[SNRindB1 SNRindB2];
smld_bit_err_prban1=[smld_bit_err_prban1_1 smld_bit_err_prban1_2];
smld_bit_err_prban2=[smld_bit_err_prban2_1 smld_bit_err_prban2_2];
smld_bit_err_prban3=[smld_bit_err_prban3_1 smld_bit_err_prban3_2];
smld_bit_err_prban4=[smld_bit_err_prban4_1 smld_bit_err_prban4_2];
smld_symbol_err_prban1=[smld_symbol_err_prban1_1 smld_symbol_err_prban1_2];
smld_symbol_err_prban2=[smld_symbol_err_prban2_1 smld_symbol_err_prban2_2];
smld_symbol_err_prban3=[smld_symbol_err_prban3_1 smld_symbol_err_prban3_2];
smld_symbol_err_prban4=[smld_symbol_err_prban4_1 smld_symbol_err_prban4_2];
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([-1,17,10^(-6),0.1]);