psk.m

来自「给出 了进行相位调制的仿真程序」· M 代码 · 共 31 行

            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %   对M=2的PSK通信系统进行蒙特卡罗仿真   %
            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            % this script is much alike with the the script for theory
            % simulation
clear all
echo off
SNRindB1=0:2:20;
SNRindB2=0:0.1:12;
SNRindB3=0:0.2:32;
for i=1:length (SNRindB1),
   [pb , ps]=cm_sm32 (SNRindB1 (i));         %仿真比特和符号误码率,此处为+,-pi/2调制验证理论的正确性
   
   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_prb1(i)=Qfunct(sqrt(2*SNR));    %理论比特误码率
end;

 %%%%%%%随后绘图曲线%%%%%%%
semilogy (SNRindB2,theo_err_prb1);
hold on
semilogy (SNRindB1,smld_symbol_err_prb,'*');
grid
xlabel('Eb/No in dB')
ylabel('error probability')
title('BPSK通信系统的蒙特卡洛仿真')
%gtext('（注：“—”理论误码率；“*”误比特率；“o”误符号率）')
% gtext('（注：“r”30度变相理论误码率；“－”180度变相理论误码率）');