www.pudn.com > xunhuan.rar > xuanhuan.m

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)])