pro_ele.m

产生GF（2^13）的元素

%%%%%%%生成有限域（2^13）的元素,低位在前，高位在后%%%%%%%%%%%%
d=zeros(8192,13);
d(1,1)=1;
for i=1:8191    %%%%i表示元素a^(i-1)%%%%
tmp=d(i,13);
d(i+1,13)=d(i,12);
d(i+1,12)=d(i,11);
d(i+1,11)=d(i,10);
d(i+1,10)=d(i,9);
d(i+1,9)=d(i,8);
d(i+1,8)=d(i,7);
d(i+1,7)=d(i,6);
d(i+1,6)=d(i,5);
d(i+1,5)=xor(tmp,d(i,4));
d(i+1,4)=xor(tmp,d(i,3));
d(i+1,3)=d(i,2);
d(i+1,2)=xor(tmp,d(i,1));
d(i+1,1)=tmp;
end

%%%%生成a^(i-1)的多项式表示，以方便乘法器设计。%%%%
%%%%%ps：多项式的最后多了一个+，用“+;”替代“；”。%%%%
fid=fopen('polynomial','w');
for i=161:161
    fprintf(fid,'%c','p');
    fprintf(fid,'%c','1');
    fprintf(fid,'%c','(');
    fprintf(fid,'%d',i);
    fprintf(fid,'%c',')');
    fprintf(fid,'%c','=');
    for j=1:13
       if d(i,j)~=0
         fprintf(fid,'%c','x');
         fprintf(fid,'%c','^');
         fprintf(fid,'%d',j-1);
         fprintf(fid,'%c','+');
       end
       if j==13
         fprintf(fid,'%c\r\n',';');
       end
    end
end