Concepto en matemáticas
En matemáticas, el algoritmo de Samuelson-Berkowitz calcula de manera eficiente el polinomio característico de una matriz cuyas entradas pueden ser elementos de cualquier anillo conmutativo unitario . A diferencia del algoritmo de Faddeev-LeVerrier , no realiza divisiones, por lo que puede aplicarse a una gama más amplia de estructuras algebraicas.
Descripción del algoritmo
El algoritmo de Samuelson-Berkowitz aplicado a una matriz produce un vector cuyos elementos son el coeficiente del polinomio característico de . Calcula este vector de coeficientes recursivamente como el producto de una matriz de Toeplitz y el vector de coeficientes de una submatriz principal .
Sea una matriz particionada de modo que
La primera submatriz principal de es la matriz . Asociada con la matriz de Toeplitz
definida por
Si es ,
si es , y en general
Es decir, todas las superdiagonales de constan de ceros, la diagonal principal consta de unos, la primera subdiagonal consta de y la enésima subdiagonal consta de .
El algoritmo se aplica entonces de forma recursiva a , produciendo la matriz de Toeplitz multiplicada por el polinomio característico de , etc. Finalmente, el polinomio característico de la matriz es simplemente . El algoritmo de Samuelson-Berkowitz establece entonces que el vector definido por
contiene los coeficientes del polinomio característico de .
Debido a que cada uno de ellos puede calcularse independientemente, el algoritmo es altamente paralelizable .
Referencias
- Berkowitz, Stuart J. (30 de marzo de 1984). "Sobre el cálculo del determinante en tiempos paralelos pequeños utilizando un pequeño número de procesadores". Information Processing Letters . 18 (3): 147–150. doi :10.1016/0020-0190(84)90018-8.
- Soltys, Michael; Cook, Stephen (diciembre de 2004). "La complejidad de la demostración del álgebra lineal" (PDF) . Anales de lógica pura y aplicada . 130 (1–3): 277–323. CiteSeerX 10.1.1.308.6521 . doi :10.1016/j.apal.2003.10.018.
- Kerber, Michael (mayo de 2006). Cálculo sin división de subresultados utilizando matrices de Bezout (PS) (Reporte técnico). Saarbrücken: Max-Planck-Institut für Informatik. Tecnología. Informe MPI-I-2006-1-006.