function [pb,ps]=bper2pskG2m2(snr_in_dB,N)
E=1;                                           %energy per symbol
snr=10^(snr_in_dB/10);
numan=2;                                       %the number of transmit antennas
M=2;
L=log(M)/log(2);
sgman=sqrt(E*numan/(2*snr));                   %noise variance
sgmaa=sqrt(0.5);
for k=1:N,
    temp=rand;
    dsource(k)=1+floor(M*temp);                % a integer number between 1 and M, uniform       
end;
dsource1=[dsource(1:N/2)];
dsource2=[dsource((N/2+1):N)];
mappingbits=[0;
             1];
mapping=[1;
         -1];                                   %2PSK
 for k=1:N,                                     %mapping symbol and bits
     bpsk_sig(k)=mapping(dsource(k));
     dsourcebits(k,:)=mappingbits(dsource(k),:);
 end;
 bpsk_sigkrate=[bpsk_sig(1:(N/2));bpsk_sig((N/2+1):N)];
 dsourcebits1=dsourcebits((1:(N/2)),:);
 dsourcebits2=dsourcebits(((N/2+1):N),:);
 for k=1:(N/2),                                 %each frame 2antennas 2time slots
     [a(1) a(2)]=gngauss(sgmaa);                %path gain-- transmit antenna1;  receive antenna1
     [a(3) a(4)]=gngauss(sgmaa);                %path gain-- transmit antenna2;  receive antenna1
     [a(5) a(6)]=gngauss(sgmaa);                %path gain-- transmit antenna1;  receive antenna2
     [a(7) a(8)]=gngauss(sgmaa);                %path gain-- transmit antenna2;  receive antenna2
     pathgain1=[(a(1)+a(2).*i);                 %path gain-- receive antenna1
               (a(3)+a(4).*i)]; 
     pathgain2=[(a(5)+a(6).*i);                 %path gain-- receive antenna2
               (a(7)+a(8).*i)]; 
     [n(1) n(2)]=gngauss(sgman);                %noise timeslot1  receive antenna1
     [n(3) n(4)]=gngauss(sgman);                %noise timeslot2  receive antenna1
     noise1=[(n(1)+n(2).*i);                     %noise receive antenna1
            (n(3)+n(4).*i)];
     [n(5) n(6)]=gngauss(sgman);                %noise timeslot1  receive antenna2
     [n(7) n(8)]=gngauss(sgman);                %noise timeslot2  receive antenna2
     noise2=[(n(5)+n(6).*i);                     %noise receive antenna2
            (n(7)+n(8).*i)]; 
     C=[bpsk_sigkrate(1,k),bpsk_sigkrate(2,k);
        -conj(bpsk_sigkrate(2,k)),conj(bpsk_sigkrate(1,k))];
     r1(:,k)=C*pathgain1+noise1;                   %received signals receive antenna1
     r2(:,k)=C*pathgain2+noise2;                   %received signals receive antenna2
     for j=1:M,
         s1(k,j)=(abs((r1(1,k)*conj(pathgain1(1,1))+conj(r1(2,k))*pathgain1(2,1)+...
             r2(1,k)*conj(pathgain2(1,1))+conj(r2(2,k))*pathgain2(2,1))-mapping(j))).^2+... 
             (-1+(abs(pathgain1(1,1))).^2+(abs(pathgain1(2,1))).^2+...
             (abs(pathgain2(1,1))).^2+(abs(pathgain2(2,1))).^2).*(abs(mapping(j))).^2;
         s2(k,j)=(abs((r1(1,k)*conj(pathgain1(2,1))-conj(r1(2,k))*pathgain1(1,1))+...
             r2(1,k)*conj(pathgain2(2,1))-conj(r2(2,k))*pathgain2(1,1)-mapping(j))).^2+...
             (-1+(abs(pathgain1(1,1))).^2+(abs(pathgain1(2,1))).^2+...
             (abs(pathgain2(1,1))).^2+(abs(pathgain2(2,1))).^2).*(abs(mapping(j))).^2;
     end;
 end;
 numofsymbolerror=0;
 numofbiterror=0;
 for k=1:(N/2),
     [min_metric1 decis1]=min(s1(k,:));
     [min_metric2 decis2]=min(s2(k,:));
     decisbits1(k,:)=mappingbits(decis1,:);
     decisbits2(k,:)=mappingbits(decis2,:);
     for l=1:L,
         if (decisbits1(k,l)~=dsourcebits1(k,l)),
             numofbiterror=numofbiterror+1;
         end;
         if (decisbits2(k,l)~=dsourcebits2(k,l)),
             numofbiterror=numofbiterror+1;
         end;
     end;
     if (decis1~=dsource1(k)),
         numofsymbolerror=numofsymbolerror+1;
     end;
     if (decis2~=dsource2(k)),
         numofsymbolerror=numofsymbolerror+1;
     end;
 end;
 ps=numofsymbolerror/N;
 pb=numofbiterror/(L*N);