cnv_awgn_pe.m

依据GSM 05.03建议

%经卷积编译码后的改善误码率
function sig_err_rate=cnv_awgn_pe(snr_in_dB)
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(E/snr)/2;%噪声功率
err_num=0;
%信息序列
for i=1:10
input=round(rand(1,267));
G=[1 0 0 1 1;1 1 0 1 1];
k=1;

ch_out=viterbi_encoder(G,k,input);%编码
%经信道加高斯燥声
u=2*ch_out-1;
u=sign(u+sgma*randn(size(u))+eps);
u=(u+1)/2;

[output,survivor_state,cumulated_metric]=viterbi_decoder(G,k,u);%解码
err_num=err_num+length(find(input~=output));
end;
%误码率
sig_err_rate=err_num/2670;