psk.m

本程序实现对多进制PSK数字调制系统的仿真。程序运行需要几分钟

            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %   对M=4的PSK系统进行蒙特卡罗仿真   %
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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)=Qfunct(sqrt(2*SNR));    %理论比特误码率
end;
SNRindB3=0:2:10;
SNRindB4=0:0.1:10;
for i=1:length (SNRindB3),
   [pb , ps]=cm_sm321 (SNRindB3 (i));       %仿真比特和符号误码率
   smld_bit_err_prb1 (i)=pb;
   smld_symbol_err_prb1 (i)=ps;
end;
for i=1:length (SNRindB4),
   SNR=exp(SNRindB4(i)*log(10)/10);        %信噪比
   theo_err_prb1(i)=2*Qfunct(sqrt(6*SNR)*sin(pi/8))/3;    %理论比特误码率
end;
 %%%%%%%随后绘图曲线%%%%%%%

 %%%%%%%随后绘图曲线%%%%%%%
semilogy (SNRindB1,smld_bit_err_prb,'*');
hold
semilogy (SNRindB1,smld_symbol_err_prb,'O');
semilogy (SNRindB2,theo_err_prb);
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('4PSK系统的蒙特卡罗仿真')
%gtext('（注：“—”理论误码率；“*”误比特率；“o”误符号率）')
semilogy (SNRindB3,smld_bit_err_prb1,'+');
semilogy (SNRindB3,smld_symbol_err_prb1,'^');
semilogy (SNRindB4,theo_err_prb1,':');
grid on
xlabel('Eb/No in dB')
ylabel('error probability')
title('8PSK系统的蒙特卡罗仿真')
%gtext('（注：“--”理论误码率；“*”误比特率；“o”误符号率）')