n=4;k=2;
msg=randint（4k[01]）;
cy=cyclpoly（nk）;
code=encode（msgnk‘cyclic/binary‘）;
nosiycode=rem（code+randerr（4n[0?1;.7?.3]）2）;
newmsg=decode（nosiycodenk‘cyclic‘）;
[numberratio]=biterr（newmsgmsg）;
disp（[‘The?bit?error?rate?is‘num2str（ratio）]）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2010-06-30?09:25??97288394xunhuan\
?????文件?????????264??2007-04-01?13:19??97288394xunhuan\xuanhuan.m
?????文件????????6125??2007-04-01?10:54??97288394xunhuan\xuanhuan2.m