% function out=gf28(a,b)   %伽羅華域乘法
temp1=4;
temp2=4;
i=8;
for j=0:7  %將temp1中的值化為二進制,存在一個1行的矩陣g中,矩陣每列存一位
g(i-j)=rem(temp1,2);
temp1=fix(temp1/2);
end
for j=0:7 %將temp2中的值化為二進制,存在一個1行的矩陣h中,矩陣每列存一位
h(i-j)=rem(temp2,2);
temp2=fix(temp2/2);
end
s=conv(g,h);  %多項式乘法
[h l]=size(s);
for i=1:l    %奇數=1;偶數=0
    if(rem(s(i),2)==0)s(i)=0;
    else if(rem(s(i),2)~=0)s(i)=1;
        end
    end
end
m=[1 0 0 0 1 1 0 1 1];  %不可月多項式M
[k r]=deconv(s,m);    %多項式除法;l為余數,h為商
[h l]=size(r);  
for i=1:l  %奇數=1;偶數=0
    if(rem(r(i),2)==0)r(i)=0;
    else if(rem(r(i),2)~=0)r(i)=1;
        end
    end
end
c=0;
for i=0:l-1   %結果放入c中
    c=c+r(l-i)*(2^i);
end
out=c; %輸出結果
out
————————————————