juanji.m

来自「matlab在通信原理中的应用」· M 代码 · 共 10 行

len=1000;
msg=randint(2*len,1);
trel=poly2trellis([5,4],[27 33 0;0 5 13]);
code=convenc(msg,trel);
ncode=rem(code+randerr(3*len,1,[0 1;.96 .04]),2);
decoded=vitdec(ncode,trel,34,'cont','hard');
[number,ratio]=biterr(decoded(68+1:end),msg(1:end-68));
disp(number);
disp(ratio);