err_pict.m

信道解码viterbi算法的matlab实现

%ce55_cnv.m 二进制双极性（2，1，3）卷积码通信系统的蒙特卡罗仿真绘图
%仿真序列长度10^5，运行时间约２小时
echo on
SNRindB1=0:16;
SNRindB2=0:0.1:16;
smld_err_prb=zeros(1,length(SNRindB1));
smld_err_prb_cnv=zeros(1,length(SNRindB1));
SNR=0;
theo_err_prb=zeros(1,length(SNRindB2));
for i=1:length(SNRindB1)
   smld_err_prb(i)=smldPe55(SNRindB1(i));
   smld_err_prb_cnv(i)=smldPe55_cnv(SNRindB1(i));
end
for  i=1:length(SNRindB2)
   SNR=exp(SNRindB2(i)*log(10)/10);
   theo_err_prb(i)=(1/2)*erfc(sqrt(2*SNR)/sqrt(2));
   %Qfunct y=(1/2)*erfc(x/sqrt(2));
   %theo_err_prb(i)=Qfunct(sqrt(2*SNR));
end
%绘图函数
figure;
semilogy(SNRindB1,smld_err_prb,'b*-');
axis([0,16,1e-5,1]);
xlabel('Eb/N0 in dB');
ylabel('Pe');
title('（２，１，３）卷积编码系统仿真结果与未编码系统的比较');
hold on
semilogy(SNRindB1,smld_err_prb_cnv,'ro-');
semilogy(SNRindB2,theo_err_prb,'b:');