mul_2_5.asv

修正的euclid算法matlab程序

%d=zeros(31,5);
%d(1,:)=[1,0,0,0,0];
%for i=2:31
%    tmp=d(i-1,5);
%    d(i,5)=d(i-1,4);
%    d(i,4)=d(i-1,3);
%    d(i,3)=xor(tmp,d(i-1,2));
%    d(i,2)=d(i-1,1);
%    d(i,1)=tmp;
%end

%syms a1 a2 a3 a4 a5 b1 b2 b3 b4 b5 x
%y1=a1+a2*x+a3*x^2+a4*x^3+a5*x^4;
%y2=b1+b2*x+b3*x^2+b4*x^3+b5*x^4;
%y=collect(expand(y1*y2));
%y=a5*(x^3+x^2+1)*b5+(a4*b5+a5*b4)*(x^4+x^2)+(a3*b5+a4*b4+a5*b3)*(x^3+x)+(a2*b5+a4*b3+a3*b4+a5*b2)*(x^2+1)+(a5*b1+a1*b5+a4*b2+a2*b4+a3*b3)*x^4+(a4*b1+a2*b3+a3*b2+a1*b4)*x^3+(a1*b3+a3*b1+a2*b2)*x^2+(a2*b1+a1*b2)*x+a1*b1;
%y=a5*(x^3+x^2+1)*b5+(a4*b5+a5*b4)*(x^4+x^2)+(a3*b5+a4*b4+a5*b3)*(x^3+x)+(a2*b5+a4*b3+a3*b4+a5*b2)*(x^2+1)+(a5*b1+a1*b5+a4*b2+a2*b4+a3*b3)*x^4+(a4*b1+a2*b3+a3*b2+a1*b4)*x^3+(a1*b3+a3*b1+a2*b2)*x^2+(a2*b1+a1*b2)*x+a1*b1
%y=collect(expand(y));
%y=mod(y,2);
%y=(a4*b2+a4*b5+a3*b3+a2*b4+a5*b4+a5*b1+a1*b5)*x^4+(a4*b1+a3*b5+a1*b4+a2*b3+a3*b2+a4*b4+a5*b5+a5*b3)*x^3+(a3*b4+a5*b2+a5*b4+a4*b5+a1*b3+a5*b5+a4*b3+a2*b5+a3*b1+a2*b2)*x^2+(a4*b4+a5*b3+a2*b1+a1*b2+a3*b5)*x+a1*b1+a5*b5+a5*b2+a4*b3+a3*b4+a2*b5;
function product=mul_2_5(a,b)%LSB-MSB
product(1)=mod(a(1)*b(1)+a(5)*b(5)+a(5)*b(2)+a(4)*b(3)+a(3)*b(4)+a(2)*b(5),2);
product(2)=mod(a(4)*b(4)+a(5)*b(3)+a(2)*b(1)+a(1)*b(2)+a(3)*b(5),2);
product(3)=mod(a(3)*b(4)+a(5)*b(2)+a(5)*b(4)+a(4)*b(5)+a(1)*b(3)+a(5)*b(5)+a(4)*b(3)+a(2)*b(5)+a(3)*b(1)+a(2)*b(2),2);
product(4)=mod(a(4)*b(1)+a(3)*b(5)+a(1)*b(4)+a(2)*b(3)+a(3)*b(2)+a(4)*b(4)+a(5)*b(5)+a(5)*b(3),2);
product(5)=mod(a(4)*b(2)+a(4)*b(5)+a(3)*b(3)+a(2)*b(4)+a(5)*b(4)+a(5)*b(1)+a(1)*b(5),2);