1.m

来自「对于一个使用相关器的二进制通信系统」· M 代码 · 共 20 行

echo on
SNRindB1=0:1:12;
SNRindB2=0:0.1:12;
for i=1:length(SNRindB1),
    % simulated error rate
    smld_err_prb(i)=smldPe54(SNRindB1(i));
    echo off;
end;
echo on;
for i=1:length(SNRindB2),
    SNR=exp(SNRindB2)*log(10)/10);
    % theoretical error rate
    theo_err_prb(i)=Qfunct(sqrt(SNR));
    echo off;
end;
echo on;
%plotting commands
semilogy(SNRindB1,smld_err_prb,'*');
hold
semilogy(SNRindB2,smld_err_prb);