%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  ML BPSK
%  2*2系統
%  輸入參數：H 
%  輸出參數：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;
clear;

H=[0.5 0.25;
    0.125 0.5;
    0.25 0.125];

x=[-1;1];
n=[0;0;0];

r=H*x+n;

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=dec2bin(k0-1,2)-48