function output=stbc_decode_TX4RX4(H,R)
%------------------------------------------------------------------------
%此函數(shù)為四發(fā)四收條件下的譯碼函數(shù)
%H,R與兩發(fā)兩收中的H，R相對應
%------------------------------------------------------------------------

Nc=size(H,1);
H11=H(:,1);
H12=H(:,2);
H13=H(:,3);
H14=H(:,4);
H21=H(:,5);
H22=H(:,6);
H23=H(:,7);
H24=H(:,8);
H31=H(:,9);
H32=H(:,10);
H33=H(:,11);
H34=H(:,12);
H41=H(:,13);
H42=H(:,14);
H43=H(:,15);
H44=H(:,16);
R11=R(:,1);
R12=R(:,2);
R13=R(:,3);
R14=R(:,4);
R15=R(:,5);
R16=R(:,6);
R17=R(:,7);
R18=R(:,8);
R21=R(:,9);
R22=R(:,10);
R23=R(:,11);
R24=R(:,12);
R25=R(:,13);
R26=R(:,14);
R27=R(:,15);
R28=R(:,16);
R31=R(:,17);
R32=R(:,18);
R33=R(:,19);
R34=R(:,20);
R35=R(:,21);
R36=R(:,22);
R37=R(:,23);
R38=R(:,24);
R41=R(:,25);
R42=R(:,26);
R43=R(:,27);
R44=R(:,28);
R45=R(:,29);
R46=R(:,30);
R47=R(:,31);
R48=R(:,32);
for i=1:1:Nc
    X1(i,1)=R11(i)*conj(H11(i))+R12(i)*conj(H21(i))+R13(i)*conj(H31(i))+R14(i)*conj(H41(i))+conj(R15(i))*H11(i)+conj(R16(i))*H21(i)+conj(R17(i))*H31(i)+conj(R18(i))*H41(i)+...
            R21(i)*conj(H12(i))+R22(i)*conj(H22(i))+R23(i)*conj(H32(i))+R24(i)*conj(H42(i))+conj(R25(i))*H12(i)+conj(R26(i))*H22(i)+conj(R27(i))*H32(i)+conj(R28(i))*H42(i)+...
            R31(i)*conj(H13(i))+R32(i)*conj(H23(i))+R33(i)*conj(H33(i))+R34(i)*conj(H43(i))+conj(R35(i))*H13(i)+conj(R36(i))*H23(i)+conj(R37(i))*H33(i)+conj(R38(i))*H43(i)+...
            R41(i)*conj(H14(i))+R42(i)*conj(H24(i))+R43(i)*conj(H34(i))+R44(i)*conj(H44(i))+conj(R45(i))*H14(i)+conj(R46(i))*H24(i)+conj(R47(i))*H34(i)+conj(R48(i))*H44(i);
    X2(i,1)=R11(i)*conj(H21(i))-R12(i)*conj(H11(i))-R13(i)*conj(H41(i))+R14(i)*conj(H31(i))+conj(R15(i))*H21(i)-conj(R16(i))*H11(i)-conj(R17(i))*H41(i)+conj(R18(i))*H31(i)+...
            R21(i)*conj(H22(i))-R22(i)*conj(H12(i))-R23(i)*conj(H42(i))+R24(i)*conj(H32(i))+conj(R25(i))*H22(i)-conj(R26(i))*H12(i)-conj(R27(i))*H42(i)+conj(R28(i))*H32(i)+...
            R31(i)*conj(H23(i))-R32(i)*conj(H13(i))-R33(i)*conj(H43(i))+R34(i)*conj(H33(i))+conj(R35(i))*H23(i)-conj(R36(i))*H13(i)-conj(R37(i))*H43(i)+conj(R38(i))*H33(i)+...
            R41(i)*conj(H24(i))-R42(i)*conj(H14(i))-R43(i)*conj(H44(i))+R44(i)*conj(H34(i))+conj(R45(i))*H24(i)-conj(R46(i))*H14(i)-conj(R47(i))*H44(i)+conj(R48(i))*H34(i);
    X3(i,1)=R11(i)*conj(H31(i))+R12(i)*conj(H41(i))-R13(i)*conj(H11(i))-R14(i)*conj(H21(i))+conj(R15(i))*H31(i)+conj(R16(i))*H41(i)-conj(R17(i))*H11(i)-conj(R18(i))*H21(i)+...
            R21(i)*conj(H32(i))+R22(i)*conj(H42(i))-R23(i)*conj(H12(i))-R24(i)*conj(H22(i))+conj(R25(i))*H32(i)+conj(R26(i))*H42(i)-conj(R27(i))*H12(i)-conj(R28(i))*H22(i)+...
            R31(i)*conj(H33(i))+R32(i)*conj(H43(i))-R33(i)*conj(H13(i))-R34(i)*conj(H23(i))+conj(R35(i))*H33(i)+conj(R36(i))*H43(i)-conj(R37(i))*H13(i)-conj(R38(i))*H23(i)+...
            R41(i)*conj(H34(i))+R42(i)*conj(H44(i))-R43(i)*conj(H14(i))-R44(i)*conj(H24(i))+conj(R45(i))*H34(i)+conj(R46(i))*H44(i)-conj(R47(i))*H14(i)-conj(R48(i))*H24(i);
    X4(i,1)=R11(i)*conj(H41(i))-R12(i)*conj(H31(i))+R13(i)*conj(H21(i))-R14(i)*conj(H11(i))+conj(R15(i))*H41(i)-conj(R16(i))*H31(i)+conj(R17(i))*H21(i)-conj(R18(i))*H11(i)+...
            R21(i)*conj(H42(i))-R22(i)*conj(H32(i))+R23(i)*conj(H22(i))-R24(i)*conj(H12(i))+conj(R25(i))*H42(i)-conj(R26(i))*H32(i)+conj(R27(i))*H22(i)-conj(R28(i))*H12(i)+...
            R31(i)*conj(H43(i))-R32(i)*conj(H33(i))+R33(i)*conj(H23(i))-R34(i)*conj(H13(i))+conj(R35(i))*H43(i)-conj(R36(i))*H33(i)+conj(R37(i))*H23(i)-conj(R38(i))*H13(i)+...
            R41(i)*conj(H44(i))-R42(i)*conj(H34(i))+R43(i)*conj(H24(i))-R44(i)*conj(H14(i))+conj(R45(i))*H44(i)-conj(R46(i))*H34(i)+conj(R47(i))*H24(i)-conj(R48(i))*H14(i);
            
 end
output=[X1 X2 X3 X4];