bch.m

matlab在通信原理中的应用

m=4;n=2^m-1;
params=bchpoly(n);
k=params(3,2);
t=params(3,3);
msg=randint(6,k);
code=bchenco(msg,n,k);
nosiycode=rem(code+randerr(6,n,1:t),2);
[newmsg,err,ccode]=bchdeco(nosiycode,k,t);
if ccode==code
    disp('All errors were corrected.')
end
if newmsg==msg
    disp('The massage was recovered perfectly.')
end
[number,ratio]=biterr(newmsg,msg);
disp(['The bit error rate is',num2str(ratio)])