function output=decode_ZF23_BPSK(H,L,R1,R2,R3)

output=[];
G=pinv(H);

for i=1:L,
    a=R1(i);
    b=R2(i);
    c=R3(i);
    r=[a b c].';
  
    ZF = G*r;
    
    QAM_input_I = real(ZF);

    BPSK_Demodu_I = [0 1];      %f(m)=(m+1)/2 + 1, so I=-1 ---> 1, I=1 ---> 2
    idx = find(QAM_input_I>1);
    QAM_input_I(idx) = 1;
    idx = find(QAM_input_I<-1);
    QAM_input_I(idx) = -1;

    output1 = BPSK_Demodu_I(round((QAM_input_I+1)/2) + 1);
    output=[output output1];
    
end