pam.m

来自「PAM误码率的程序,可以直接运行,有图片」· M 代码 · 共 22 行

echo on
SNRindB1=0:1:12;
SNRindB2=0:0.1:12;
for i=1:length(SNRindB1),
    smld_err_prb(i)=smldPe58(SNRindB1(i));
    echo off;
end;
echo on;
for i=1:length(SNRindB2),
    SNR_per_bit=exp(SNRindB2(i)*log(10)/10);
    theo_err_prb(i)=(3/2)*(0.5*erfc(sqrt(((4/5)*SNR_per_bit)/2)));
    echo off;
end;
echo on;
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,theo_err_prb);
echo off;
grid
xlabel('Eb/No in dB')
title('PAM通信系统的蒙特卡洛仿真')
gtext('（注：“—”理论误码率；“*”误比特率；')