apsk_xingzuo.m

实现了高斯噪声下APSK的BER仿真

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),':');