M=16; %16QAM的各种参数
N=2;
m=[0.5 1 2];
K=[1 2 4];
L=[1 2 3 4 6];
a=[4*(1-1/(sqrt(M)))/pi -4*((1-1/(sqrt(M)))^2)/pi];
phi=[1.5/(M-1) 1.5/(M-1)];
theta=[0.5*pi 0.25*pi];
snr_db=0:25;
for i=1:length(snr_db),
snr=exp(snr_db(i)*log(10)/10);
prb1(i)=Pe(N,K(1),L(2),m(1),snr,theta,phi,a); %Pe(N,K,L,m,snr,theta,phi,a)
echo off ;
end;
echo on;
semilogy(snr_db,prb1,'*');
xlabel('SNR per bit expressed in dB');
ylabel('SER');
axis([0 25 10^(-12) 1]);
grid on;
%hold
%semilogy(SNRindB2,theo_err_prb);
评论2