rs_decode_root.m

Reed Solomon码的编译码器MATLAB仿真程序

function root = rs_decode_root(sigma_x)
%本函数实现错误位置的根的计算
% sigma_x为错误位置多项式，root为错误位置的根
table=[1 2 4 8 16 5 10 20 13 26 17 7 14 28 29 31 27 19 3 6 12 24 21 15 30 25 23 11 22 9 18];
l = length(table);
root1 = zeros(1,length(sigma_x)-1);
indexroot1 = 1;
for i=1:l
    Addall = 1;
    for j=1:length(sigma_x)-1
        temp = table(i);
        for k=1:j-1
            temp = rs_symbolmultiple(temp, table(i));
        end
        Addall = rs_symboladd( Addall, rs_symbolmultiple(sigma_x(j+1) ,temp));
    end
    if Addall == 0
        root1(indexroot1) = rs_symbolrev(table(i));
        indexroot1 = indexroot1 + 1;
    end
end
j=0;
for i=1:length(root1)
    if root1(i)~=0
       j=j+1;
   end
end
root=zeros(1,j);
root=root1(1:j);