rs_enc.m

对做射频方面的朋友很有用的

function data_out = rs_enc(data_in,m,n,k,genpoly)
reg_out =gf(zeros(1,n-k+1),m);
% mult_out =gf(zeros(1,n-k),m);
% regout =zeros(1,n-k);
% adderout =zeros(1,n-k);
% multout =zeros(1,n-k);
% datain = [7 10 5 11 12 6 3 8 4 2 1];
%% Genarator polynomial
% genpoly = rsgenpoly(15,11,19,12);
%10          13          12          10
%% input ring
for i=1:k
%     adder = gf_add(datain(i),regout(1),m);
    adder = data_in(i) + reg_out(1);
    for j=1:n-k
        mult_out = adder .* genpoly(j+1);
%         mult_out = gf_mult(adder,genpoly(j+1),m);
        adder_out = mult_out + reg_out(j+1);
%         adder_out = gf_add(mult_out, reg_out(j+1),m);
        reg_out(j) = adder_out;
    end
end
data_out = [data_in reg_out(1:n-k)];