line.m

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

n=5;k=3;
msg=randint(n,k,[0,1]);
genmat=[randint(k,(n-k),[0,1]),eye(k)];
parmat=gen2par(genmat);
code=encode(msg,n,k,'linear',genmat);
nosiycode=rem(code+randerr(5,5,[0 1;.7 .3]),2);
trt=zeros(2^(n-k),n);
[newmsg,err]=decode(nosiycode,n,k,'linear',genmat,trt);
err_words=find(err~=0);
[number,ratio]=biterr(newmsg,msg);
disp(['The bit error rate is',num2str(ratio)])