%加高斯噪聲后的誤碼率
function[sig_err_rate]=PrimePe(snr_in_dB,signal)
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);
%經信道加高斯燥聲
u=2*signal-1;
u=sign(u+sgma*randn(size(u))+eps);
u=(u+1)/2;
%誤碼率
sig_err_rate=(length(find(signal~=u)))/length(signal);