rs.m

matlab在通信原理中的应用

l=1000;
m=4;
n=2^m-1;
k=n-4;
rand('state',9876);
msg=randint(l,1);
field=gftuple([-1:n-1]',m);
[code,added]=rsenco(msg,field,k);
msg1=[msg;zeros(added,1)];
noi=rand(length(code)/m,1)<.03;
noi=(noi*ones(1,m))';noi=noi(:);
code_noi=rem(code+noi,2);
[dec,err,ccode,err_c]=rsdeco(code_noi,field,k);
err_c=reshape(err_c,n,length(err_c)/n)';
num_err_codesyms=sum(err_c(:,1)>0)
num_err_decbits=sum(abs(dec-msg1))