Multiply polynomials over Galois field


c = gfconv(a,b)
c = gfconv(a,b,p)
c = gfconv(a,b,field)


The gfconv function multiplies polynomials over a Galois field. (To multiply elements of a Galois field, use gfmul instead.) Algebraically, multiplying polynomials over a Galois field is equivalent to convolving vectors containing the polynomials' coefficients, where the convolution operation uses arithmetic over the same Galois field.

c = gfconv(a,b) multiplies two GF(2) polynomials, a and b, which can be either polynomial strings or numeric vectors. The polynomial degree of the resulting GF(2) polynomial c equals the degree of a plus the degree of b.

c = gfconv(a,b,p) multiplies two GF(p) polynomials, where p is a prime number. a, b, and c are row vectors that give the coefficients of the corresponding polynomials in order of ascending powers. Each coefficient is between 0 and p-1.

c = gfconv(a,b,field) multiplies two GF(pm) polynomials, where p is a prime number and m is a positive integer. a, b, and c are row vectors that list the exponential formats of the coefficients of the corresponding polynomials, in order of ascending powers. The exponential format is relative to some primitive element of GF(pm). field is the matrix listing all elements of GF(pm), arranged relative to the same primitive element. See Representing Elements of Galois Fields for an explanation of these formats.


The command below shows that


over GF(3).

gfc = gfconv([1 1 0 0 1],[0 1 1],3)

The output is

gfc =

     0     1     2     1     0     1     1

The code below illustrates the identity


for the case in which p = 7, r = 5, and = 3. (The identity holds when p is any prime number, and r and s are positive integers.)

p = 7; r = 5; s = 3;
a = gfrepcov([r s]); % x^r + x^s

% Compute a^p over GF(p).
c = 1;
for ii = 1:p
   c = gfconv(c,a,p);

% Check whether c = x^(rp) + x^(sp).
powers = [];
for ii = 1:length(c)
   if c(ii)~=0
      powers = [powers, ii];
if (powers==[r*p+1 s*p+1] | powers==[s*p+1 r*p+1])
   disp('The identity is proved for this case of r, s, and p.')

See Also

| | | |

Introduced before R2006a

Was this topic helpful?