decode44.m

正交码的编吗程序

function decodeout=decode44(H,r,M)
   for L=1:M
       s(L)=qammod(L-1,M);
   end
   z1=[];z2=[];z3=[];
   for t=1:4
       z1(t)=conj(H(t,1))*r(t,1)+H(t,2)*conj(r(t,2))+2*(-conj(H(t,3))*r(t,3)-H(t,3)*conj(r(t,3))+conj(H(t,4))*r(t,3)-H(t,4)*conj(r(t,3))+conj(H(t,3))*r(t,4)-H(t,3)*conj(r(t,4))-conj(H(t,4))*r(t,4)-H(t,4)*conj(r(t,4)));
       z2(t)=conj(H(t,2))*r(t,1)-H(t,1)*conj(r(t,2))+2*(conj(H(t,3))*r(t,3)-H(t,3)*conj(r(t,3))+conj(H(t,4))*r(t,3)+H(t,4)*conj(r(t,3))-conj(H(t,3))*r(t,4)-H(t,3)*conj(r(t,4))-conj(H(t,4))*r(t,4)+H(t,4)*conj(r(t,4)));
       z3(t)=2*(conj(H(t,3))*r(t,1)+conj(H(t,4)*r(t,1))+conj(H(t,3))*r(t,2)-conj(H(t,4))*r(t,2)+H(t,1)*conj(r(t,3))+H(t,2)*conj(r(t,3))+conj(H(t,1))*r(t,4)-conj(H(t,2))*r(t,4));
   end
   Z1=sum(z1);Z2=sum(z2);Z3=sum(z3);
    for m=1:M
      d1(m)=(norm(H(1,:)).^2+norm(H(2,:)).^2+norm(H(3,:)).^2+norm(H(4,:)).^2-1)*abs(s(m)).^2+abs(Z1-s(m)).^2;
      d2(m)=(norm(H(1,:)).^2+norm(H(2,:)).^2+norm(H(3,:)).^2+norm(H(4,:)).^2-1)*abs(s(m)).^2+abs(Z2-s(m)).^2;
      d3(m)=(norm(H(1,:)).^2+norm(H(2,:)).^2+norm(H(3,:)).^2+norm(H(4,:)).^2-1)*abs(s(m)).^2+abs(Z3-s(m)).^2;
     
 end
      [y1,i1]=min(d1);
      s1=s(i1);
    
      [y2,i2]=min(d2);
      s2=s(i2);
    
      [y3,i3]=min(d3);
      s3=s(i3);
  
      det_out=[s1;s2;s3];
     decodeout=qamdemod(det_out,M);