sylm

Parâmetros

a, b

Dois polinómios simples com coeficientes reais ou complexos.

S

matriz quadrada de números reais ou complexos, de tamanho degree(a*b).

Descrição

sylm(a,b) fornece a matriz de Sylvester associada aos polinômios a e b, i.e. a matriz S tal que:

coeff( a*x + b*y )' = S * [coeff(x)';coeff(y)'].

Se a e b são polinômios primos entre si, então

rank(sylm(a,b))=degree(a)+degree(b)) e as instruções

u = sylm(a,b) \ eye(na+nb,1)
x = poly(u(1:nb),'z','coeff')
y = poly(u(nb+1:na+nb),'z','coeff')

u = sylm(a,b) \ eye(na+nb,1)
x = poly(u(1:nb),'z','coeff')
y = poly(u(nb+1:na+nb),'z','coeff')

computam os fatores de Bezout x e y de graus mínimos tais que a*x+b*y = 1

Exemplos

x = poly(0,"x");
y = poly ([1, 2, 3], "x", "coeff")
sylm(x, y)

--> x = poly(0, "x");
--> y = poly([1, 2, 3], "x","coeff")
 y  =
  1 +2x +3x²

--> sylm(x, y)
 ans  =
   0.   0.   1.
   1.   0.   2.
   0.   1.   3.

Resultant of two polynomials: It can be defined as the determinant of the Sylvester matrix of the two polynomials. It is null if and only if both polynomials have at least one common root.

a = poly([1 2 3 4], "x", "roots")
b = poly([-2 -1 5], "x", "roots")
det(sylm(a, b))

// Simpler case:
det(sylm((%s+1)^5, (%s+3)^3))
(-3 -(-1))^(5*3)     // Expected resultant

--> a = poly([1 2 3 4],"x","roots")
 a  =
  24 -50x +35x² -10x³ +x⁴

--> b = poly([-2 -1 5],"x","roots")
 b  =
  -10 -13x -2x² +x³

--> det(sylm(a,b))
 ans  =
   1036800.0

--> // Simpler case:
--> det(sylm((%s+1)^5,(%s+3)^3))
 ans  =
  -32768.000

--> (-3 -(-1))^(5*3)     // Expected resultant
 ans  =
  -32768.

Ver também

• bezout — Maior Comum Divisor de dois polinômios ou dois inteiros, pelo método Bezout
• diophant — Solves the diophantine (Bezout) equation p1*x1 + p2*x2 = b