smldpe.m

依据GSM 05.03建议

function sig_err_rate=smldpe(snr_in_dB)
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(E/snr)/2;%噪声功率
%信息序列
input=round(rand(1,10000));
%经信道加高斯燥声
u=2*input-1;
u=sign(u+sgma*randn(size(u)));
u=(u+1)/2;
%误码率
sig_err_rate=(length(find(input~=u)))/length(input);