function s=qampsk_modulate(R,Q_modu)%輸入為編碼后的，輸出為復數點Q=4or16or8
if Q_modu==4
   table=exp(j*[pi/4 ,3*pi/4 ,-pi/4,-3*pi/4]);%mapping表
   table=table([1,3,2,4]);
   s=table([ 2 1]*R+1);%s為mapping后的復數點  
   
elseif Q_modu==8
   table=exp(j*[0 ,pi/4 ,2*pi/4 ,3*pi/4 ,pi ,-3*pi/4,-2*pi/4 ,-pi/4]);
   table=table([1,2,4,3,8,7,5,6]);
   s=table([4 2 1]*R+1);%s為mapping后的復數點  


elseif Q_modu==16
     table=1/sqrt(10)*[-3+3*j,-1+3*j,3+3*j,1+3*j,-3+j,-1+j,3+j,1+j,-3-3*j,-1-3*j,3-3*j,1-3*j,-3-j,-1-j,3-j,1-j];
     s=table([8 4 2 1]*R+1);%s為mapping后的復數點 
%  elseif Q_modu==64%利用最大距離
%      table=1/sqrt(42)*[7+7*j,7+5*j,7+j,7+3*j,7-j,7-3*j,7-7*j,7-5*j,...
%                        5+7*j,5+5*j,5+j,5+3*j,5-j,5-3*j,5-7*j,5-5*j,...
%                        1+7*j,1+5*j,1+j,1+3*j,1-j,1-3*j,1-7*j,1-5*j,...
%                        3+7*j,3+5*j,3+j,3+3*j,3-j,3-3*j,3-7*j,3-5*j,...
%                        -1+7*j,-1+5*j,-1+j,-1+3*j,-1-j,-1-3*j,-1-7*j,-1-5*j,...
%                        -3+7*j,-3+5*j,-3+j,-3+3*j,-3-j,-3-3*j,-3-7*j,-3-5*j,...
%                        -7+7*j,-7+5*j,-7+j,-7+3*j,-7-j,-7-3*j,-7-7*j,-7-5*j,...
%                         -5+7*j,-5+5*j,-5+j,-5+3*j,-1-j,-1-3*j,-1-7*j,-1-5*j];
%          s=table([32 16 8 4 2 1]*R+1);%s為mapping后的復數點 
elseif Q_modu==64%利用reim
     table=1/sqrt(42)*[-7-7*j,-7-5*j,-7-j,-7-3*j,-7+7*j,-7+5*j,-7+j,-7+3*j,...
                        -5-7*j,-5-5*j,-5-j,-5-3*j,-5+7*j,-5+5*j,-5+j,-5+3*j,...
                        -1-7*j,-1-5*j,-1-j,-7-3*j,-1+7*j,-1+5*j,-1+j,-1+3*j,...
                        -3-7*j,-3-5*j,-3-j,-3-3*j,-3+7*j,-3+5*j,-3+j,-3+3*j,...
                         7-7*j,7-5*j,7-j,7-3*j,7+7*j,7+5*j,7+j,7+3*j,...
                         5-7*j,5-5*j,5-j,5-3*j,5+7*j,5+5*j,5+j,5+3*j,...
                          1-7*j,1-5*j,1-j,1-3*j,1+7*j,1+5*j,1+j,1+3*j,...
                         3-7*j,3-5*j,3-j,3-3*j,3+7*j,3+5*j,3+j,3+3*j];
         s=table([32 16 8 4 2 1]*R+1);%s為mapping后的復數點 
% 
 end