ch8example20func.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 12 行

function Pe=ch8example20func(P0,SNR)
s1=5; s0=0; % 发送电平（单极性）
A=s1;
if P0>1 | P0<0
    P0=rand; % 输入超限的处理
end
sigma2=A^2./SNR; % 由信噪比求出方差
P1=1-P0;         % 发‘1’概率
C_opt=(s0+s1)./2+sigma2./(s1-s0).*log(P0./P1); % 计算最佳判决门限
Pe0=0.5-0.5*erf((C_opt-s0)./(sqrt(2*sigma2))); % 发0出错率
Pe1=0.5+0.5*erf((C_opt-s1)./(sqrt(2*sigma2))); % 发1出错率
Pe=P0*Pe0+P1*Pe1;                              % 平均错误率