function [H1,H2]=ls_time_2x2(y,chantype,sampletime)
% Pure time domain maximum likelihood estimation
% H1,H2: freq channel response from two Tx to this Rx
% y: received time domain signal on this Rx(column vector)

switch chantype,
    case 1,     % SUI-1 Channel
        tau = [0 0.4 0.9];      % tap delay in us
        load X_pinv_sui1.mat;
    case 2,     % SUI-2 Channel
        tau = [0 0.4 1.1];
        load X_pinv_sui2.mat;
    case 3,     % SUI-3 Channel
        tau = [0.0 0.4 0.9];
        load X_pinv_sui3.mat;
    case 4,     % SUI-4 Channel
        tau = [0 1.5 4];
        load X_pinv_sui4.mat;
    case 5,     % SUI-5 Channel
        tau = [0 4 10];
        load X_pinv_sui5.mat;
    case 6,     % SUI-6 Channel
        tau = [0 14 20];
        load X_pinv_sui6.mat;
    otherwise,
        error('Unrecognized channel model');
end
L=max(round(tau/sampletime))+1;

h=X_pinv*y;
h1=h(1:L).'; h2=h(L+1:end).';

H1=fftshift(fft(h1,2048));
H1(:,1877:2048)=[]; H1(:,1:173)=[];
H2=fftshift(fft(h2,2048));
H2(:,1877:2048)=[]; H2(:,1:173)=[];