demapping16qam.m

这是一个基于16QAM调制的OFDM系统仿真。

function Code=demapping16QAM(bit)
bitlength=length(bit);
A=1/sqrt(10);
for ii=1:bitlength
    if real(bit(ii))>=0
        if imag(bit(ii))>=0%------------1#Quadrand
            if real(bit(ii))>=2*A
                if imag(bit(ii))>=2*A
                    temp(ii*4)=0;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=0;
                else
                    temp(ii*4)=0;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=0;
                end
            else
                if imag(bit(ii))>=2*A
                    temp(ii*4)=0;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=1;
                else
                    temp(ii*4)=0;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=1;
                end
            end
        else        %----------4#Quadrant
             if real(bit(ii))>=2*A
                  if imag(bit(ii))>=-2*A
                    temp(ii*4)=1;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=0;
                  else
                      temp(ii*4)=1;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=0;
                  end
             else
                 if imag(bit(ii))>=-2*A
                     temp(ii*4)=1;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=1;
                 else
                     temp(ii*4)=1;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=0;
                    temp(ii*4-3)=1;
                 end
             end
        end
    else
       if imag(bit(ii))>=0   %-----------2#Quadrant
           if real(bit(ii))>=-2*A
               if imag(bit(ii))>=2*A
                   temp(ii*4)=0;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=1;
               else
                   temp(ii*4)=0;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=1;
               end
           else
               if imag(bit(ii))>=2*A
                   temp(ii*4)=0;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=0;
               else
                   temp(ii*4)=0;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=0;
               end
           end 
       else   %-----------3#Quadrant
           if real(bit(ii))>=-2*A
               if imag(bit(ii))>=-2*A
                   temp(ii*4)=1;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=1;
               else
                   temp(ii*4)=1;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=1;
               end
           else
               if imag(bit(ii))>=-2*A
                   temp(ii*4)=1;
                    temp(ii*4-1)=1;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=0;
               else
                   temp(ii*4)=1;
                    temp(ii*4-1)=0;
                    temp(ii*4-2)=1;
                    temp(ii*4-3)=0;
               end
           end
       end
    end
end
Code=temp;