psk.m

产生MPSK信号的matlab源代码

clear all
echo off
SNRindB1=0:2:10;
SNRindB2=0:0.1:10;

for i=1:length(SNRindB1),
   [pb,ps]=cm_sm32(SNRindB1(i));
   smld_bit_err_prb(i)=pb;
   smld_symbol_err_prb(i)=ps;
end;

for i=1:length(SNRindB2),
   SNR=exp(SNRindB2(i)*log(10)/10);
   theo_err_prb(i)=2*Qfunct(sqrt(6*SNR)*sin(pi/8))/3;
 end;

semilogy(SNRindB1,smld_bit_err_prb,'*');
hold on
semilogy(SNRindB1,smld_symbol_err_prb,'O');
semilogy(SNRindB2,theo_err_prb);
axis([0,10,10^(-5),10^0]);
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('8PSK通信系统的蒙特卡洛仿真')
%gtext('（注：“—”理论误码率；“*”误比特率；“o”误符号率）')