xuanhuan.m

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

n=4;k=2;
msg=randint(4,k,[0,1]);
cy=cyclpoly(n,k);
code=encode(msg,n,k,'cyclic/binary');
nosiycode=rem(code+randerr(4,n,[0 1;.7 .3]),2);
newmsg=decode(nosiycode,n,k,'cyclic');
[number,ratio]=biterr(newmsg,msg);
disp(['The bit error rate is',num2str(ratio)])