function y=rsdecode(input)
%rs譯碼程序 
m=5;       
n=2^m-1;    
t=3;         
k=n-2*t;   

w=size(input);
l=w(1)*w(2)/(n*m);%輸入與編碼輸出的序列應(yīng)該是對(duì)應(yīng)的

%把二進(jìn)制轉(zhuǎn)換到gf域上用于譯碼
y1=reshape(input,m,l*n)';
y2=bin2dec(char(y1+48));
y3=reshape(y2,n,l)';
y4=gf(y3,m);

y5=rsdec(y4,n,k);%譯碼

%轉(zhuǎn)換輸出為二進(jìn)制
y6=double(y5.x)';
y7=dec2bin(y6,m)';
b=zeros(1,l*k);
b1=dec2bin(b,m)';
y8=y7-b1;

y=reshape(y8,1,m*l*k);%并串變換