clear all;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 系統參數設置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N=10240;
M=16;
K=4;
v=sqrt(-1);
alpha=2.75;        %16APSK性能最優時內外半徑比,r1為內半徑
P=12;              % 信號功率設定
jiao_du=0:pi/100:2*pi;
x=randint(1,N,M)+1;  %產生一組隨機數
%%%%%%%%%%%%%%%%%%%%%%%%%%  16APSK信源產生模塊 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r1=sqrt(P/(4+12*alpha*alpha));
r2=r1*alpha;
s_apsk=[r1*exp(v*(2*pi*1/4-pi/4));
   r1*exp(v*(2*pi*2/4-pi/4));
   r1*exp(v*(2*pi*3/4-pi/4));
   r1*exp(v*(2*pi*4/4-pi/4));
   r2*exp(v*(2*pi*1/12-pi/12));
   r2*exp(v*(2*pi*2/12-pi/12));
   r2*exp(v*(2*pi*3/12-pi/12));
   r2*exp(v*(2*pi*4/12-pi/12));
   r2*exp(v*(2*pi*5/12-pi/12));
   r2*exp(v*(2*pi*6/12-pi/12));
   r2*exp(v*(2*pi*7/12-pi/12));
   r2*exp(v*(2*pi*8/12-pi/12));
   r2*exp(v*(2*pi*9/12-pi/12));
   r2*exp(v*(2*pi*10/12-pi/12));
   r2*exp(v*(2*pi*11/12-pi/12));
   r2*exp(v*(2*pi*12/12-pi/12))].'; 
for n=1:N
    x_apsk(n)=s_apsk(x(n));
end
figure(1);
plot(real(x_apsk),imag(x_apsk),'.');title('16APK信號星座圖');%grid;
hold on;
plot(r1*cos(jiao_du),r1*sin(jiao_du),':');
hold on;
plot(r2*cos(jiao_du),r2*sin(jiao_du),':');