de_16qam.m

QPSK通过瑞利信道的MATLAB仿真

%2005-3-25
function [symbol_bit,symbol_symbol]=DE_16QAM(xxxxx,HHH)

LLLLL=length(xxxxx);
sixteen_qam=[0 0 0 0
             0 0 0 1
             0 0 1 0
             0 0 1 1
             0 1 0 0
             0 1 0 1
             0 1 1 0
             0 1 1 1
             1 0 0 0
             1 0 0 1
             1 0 1 0
             1 0 1 1
             1 1 0 0
             1 1 0 1
             1 1 1 0
             1 1 1 1];
sixteen_QAM=1/sqrt(10)*[3+3j 3+j 3-3j 3-j 1+3j 1+j 1-3j 1-j -3+3j -3+j -3-3j -3-j -1+3j -1+j -1-3j -1-j];
  
for i=1:LLLLL
    for ii=1:16
        SYMBOL(ii)=abs(xxxxx(i)-sixteen_QAM(ii))^2+(-1+sum(abs(HHH.^2)))*(abs(sixteen_QAM(ii)))^2;
    end
    [value position]=min(SYMBOL);
    symbol_bit(1,(4*i-3):4*i)=sixteen_qam(position,:);
    symbol_symbol(1,i)=sixteen_QAM(1,position);
end