function mod_out=modulation(mod_in,mod_mode)

%%根據輸入的調制方式，對輸入序列MOD_IN進行調制，采用64QAM，
%%完成對星座圖的映射，輸出為Y.轉化的方法為：先寫出十進制情況下從0 到N-1
%%（N為星座圖的點數）所對應的星座坐標；再將輸入的二進制序列轉化為相應的
%%十進制，以查表的方法查出對應點的復數坐標，即為調制映射后的結果。



%% Input: 
%%-------------------------------------------------------------------------
%%  mod_in:輸入的二進制序列(The sequence to be modulated)
%%-------------------------------------------------------------------------
%% Output:
%%-------------------------------------------------------------------------
%% mod_out:星座圖映射后得到的調制復數結果(The output after modulation)
%%-------------------------------------------------------------------------
%% Global Variable:
%%  g_RT (the vector which contains the modulation mode)
%%-------------------------------------------------------------------------
%% Z :選擇調制方式的參數 (the parameter to choose the modulation mode)
%%
%% R :輸入二進制序列重新排列（按一定要求）后的結果,例如：對64QAM，要把輸入序列調整為
%% 4行，length(g_MOD_IN_16QAM )/4  列的矩陣。(Reshape the input binary sequence to be
%% matrix of n-row,m-column .For example,16QAM,will be reshaped into 4-row,
%%length(g_MOD_IN_16QAM))/4 column )
%%
%% B2D :二進制向十進制轉化后的結果 (convert the binary sequence to  dec )
%%
%% Temp:星座圖陣列 (the  constellation)
%%-------------------------------------------------------------------------
%%********************************************************
%system_parameters  
switch (mod_mode)
     mod_out=zeros(1,length(mod_in)/6);
   R=reshape(mod_in,6,length(mod_in)/6);
   B2D=bi2de(R','left-msb')+1;                                         %將二進制轉為十進制,注意加1，因為matlab沒有a(0)項，而是從a(1)開始      left-msb:表明左邊為最高位    --yzh
   Temp=[-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...
          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 ];
   for i=1:length(mod_in)/6
       mod_out(i)=Temp(B2D(i))/sqrt(42);
   end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
otherwise
    disp('Error! Please input again');
end