function x_f=decode_ML22_BPSK(H,L,R1,R2,R3)

x_f=[];
for i=1:L,
    a=R1(i);
    b=R2(i);
    c=R3(i);
    r=[a b c].';
    
    s=[];
    s0=[-1 1];
    for m1=1:2
        for m2=1:2
           s1=[s0(m1);s0(m2)];
           s=[s s1]; 
        end
    end

    estimate=[];
    for m3=1:4
        estimate1=r-H*s(:,m3);
        estimate=[estimate estimate1];
    end

    [gk k0]=min(sum(abs(estimate).^2,1));
    x_f_t=dec2bin(k0-1,2)-48;
    x_f=[x_f x_f_t];

    
    
    
    
    
    
end