function?[XYZ]=StokesDistribution（exeyM）
ax=abs（ex）;ay=abs（ey）;
Normalization8=sqrt（1+（sqrt（3）+1）^2/2）;%8qam歸一化因子
if?M==8
????ax=abs（ax）/Normalization8;ay=abs（ay）/Normalization8;%because?of?Procedural?problems8qam?special?dispose
elseif?M==32
????ax=abs（ax）/sqrt（20）;ay=abs（ay）/sqrt（20）;%32qam
else?
????ax=abs（ax）/sqrt（（M-1）*2/3）;ay=abs（ay）/sqrt（（M-1）*2/3）;%power?Norm