function?Get_Color_Constellation（wavein）
????N?=?4;
????Nbeside?=?20;%左右兩側共2*Nbeside個點
????indexdata?=?zeros（length（wavein）N）;
????angledata?=?linspace（0pi/2N）;
????for?m?=?1:N
????????data?=?real（exp（1j*angledata（m）））*real（wavein）+...
????????????imag（exp（1j*angledata（m）））*imag（wavein）;
????????[~sortdata]?=?sort（data）;?
????????indexdata（:m）?=?sortdata;
????end
????dis?=?zeros（length（indexdata）N）;
????for?kk?=?1:N
????????for?k?=?1:length（indexdata）
????????????indexbeg?=?k-Nbeside;
????????????if?indexbeg?????????????????indexbeg?=?1;
????????????end
????????????indexend?=?indexbeg+Nbeside*2+1;
????????????if?indexend?>?length（data）
????????????????indexend?=?length（data）;
????????????????indexbeg?=?length（data）-2*Nbeside;
????????????end??
????????????datatemp?=?wavein（indexdata（indexbeg:indexendkk））-wavein（indexdata（kkk））;
????????????[t1t2]=sort（abs（datatemp））;
????????????dis（indexdata（kkk）kk）?=?t1（2）;
????????e