function[enc_comp]=cm_dpske(E,M,mapping,sequence)
%[enc_comp]=cm_dpske(E,M,mapping,sequence)
%            CM_DPSKE以微分方式編碼一個序列
%            E是平均能量，M是星座點的數目
%            而mapping是一個定義星座點如何分布的向量
%            sequence是未編碼的二進制數據序列
k=log2(M);
N=length(sequence);
%如果N不可為k整除，填充零使其可被k整除
remainder=rem(N,k);
if (remainder~=0),
  for i=N+1:N+k-remainder,   
     sequence(i)=0;
  end;
  N=N+k-remainder;
end;
theta=0;
for i=1:k:N,
   index=0;
   for j=i:i+k-1,
      index=2*index+sequence(j);
   end;
   index=index+1;
   theta=mod(2*pi*mapping(index)/M+theta,2*pi);
   enc_comp((i+k-1)/k,1)=sqrt(E)*cos(theta);
   enc_comp((i+k-1)/k,2)=sqrt(E)*sin(theta);
end;