MATLAB维特比算法卷积码译码

%%% 《通信原理》P355卷积码产生及译码
clear;close all;
%%  Convolutional code
%   Random sequence generate
Len=50;
x=round（rand（1Len））;
%   （313） PCCC
c1=1;
c2=[1 0 1];
c3=[1 1 1];
o1=mod（conv（xc1）2）;
o2=mod（conv（xc2）2）;
o3=mod（conv（xc3）2）;
o1=o1（1:Len）;
o2=o2（1:Len）;
o3=o3（1:Len）;
O=[o1;o2;o3];
encode=reshape（O1[]）;
%%  Mapping
%   2^（3-1）=4 status: a-00 b-01 c-10 d-11
%   start: 00 send 2 zeros
%   former status     input bit       output       next status
M=cell（81）;
M{1}={‘a‘0[0 0 0]‘a‘};
M{2}={‘a‘1[1 1 1]‘b‘};
M{3}={‘b‘0[0 0 1]‘c‘};
M{4}={‘b‘1[1 1 0]‘d‘};
M{5}={‘c‘0[0 1 1]‘a‘};
M{6}={‘c‘1[1 0 0]‘b‘};
M{7}={‘d‘0[0 1 0]‘c‘};
M{8}={‘d‘1[1 0 1]‘d‘};
%%  Decode
%   simulate errors
errNum=10;
err_pos=randi（3*Len[1errNum]）;
recv=encod