qam16.m

在高斯信道下

function q_code=qam16(qa_source)
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 a=1:length(qa_source),
   qam_code(a,:)=map(qa_source(a),:);
   q_code(a,1)=qam_code(a,1)+(-1)^(1/2)*qam_code(a,2);
end;