clear;
SNRindB1=0:2:20;
M=16;
k=log2(M);
for i=1:length (SNRindB1),
   smld_err_prb(i)=APSK_modulation(SNRindB1(i));       %仿真誤符號率
end;
semilogy (SNRindB1,smld_err_prb,'*');
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('16APSK系統的蒙特卡羅仿真')