deqam16.m

在高斯信道下

function qa_decode=deqam16(code,M)

b=5;
for bj=0:3,
    b=b-2;
    a=-5;
  for bi=1:4,
    a=a+2;
    map(bi+4*bj,1)=a;
    map(bi+bj*4,2)=b;
  end;
end;
for x=1:length(code),
    for y=1:M,
    result(y)=(real(code(x,1))-map(y,1))^2+(imag(code(x,1))-map(y,2))^2;
    end;
    [metic j]=min(result);
    qa_decode(x)=j;
     
 end;