Articulo de referencia

Polinomio característico

En álgebra lineal , el polinomio característico de una matriz cuadrada es un polinomio invariante bajo la semejanza de matrices y cuyas raíces son los autovalores . Entre sus co...

En álgebra lineal , el polinomio característico de una matriz cuadrada es un polinomio invariante bajo la semejanza de matrices y cuyas raíces son los autovalores . Entre sus coeficientes se encuentran el determinante y la traza de la matriz. El polinomio característico de un endomorfismo de un espacio vectorial de dimensión finita es el polinomio característico de la matriz de dicho endomorfismo sobre cualquier base (es decir, el polinomio característico no depende de la elección de la base ). La ecuación característica , también conocida como ecuación determinante , [ 1 ] [ 2 ] [ 3 ] es la ecuación que se obtiene al igualar el polinomio característico a cero.

En la teoría espectral de grafos , el polinomio característico de un grafo es el polinomio característico de su matriz de adyacencia . [ 4 ]

Motivación

Los valores y vectores propios desempeñan un papel fundamental en el álgebra lineal , ya que, dada una transformación lineal , un vector propio es un vector cuya dirección no cambia con la transformación, y el valor propio correspondiente es la medida del cambio de magnitud resultante del vector.

Más precisamente, supongamos que la transformación está representada por una matriz cuadrada.A.{\displaystyle A.}Entonces un vector propiov{\displaystyle \mathbf {v} }y el valor propio correspondienteλ{\displaystyle \lambda }debe satisfacer la ecuación Av=λv,{\displaystyle A\mathbf {v} =\lambda \mathbf {v},} o, equivalentemente (ya que λv=λIv{\displaystyle \lambda \mathbf {v} =\lambda I\mathbf {v} }), (λIA)v=0{\displaystyle (\lambda IA)\mathbf {v} =\mathbf {0} } dóndeI{\displaystyle I} es la matriz identidad yv0{\displaystyle \mathbf {v} \neq \mathbf {0} } (aunque el vector cero satisface esta ecuación para cadaλ,{\displaystyle \lambda ,}no se considera un vector propio).

De ello se deduce que la matriz(λIA){\displaystyle (\lambda IA)}debe ser singular y su determinante det(λIA)=0{\displaystyle \det(\lambda I-A)=0} debe ser cero.

En otras palabras, los autovalores de A son las raíces de det(incógnitaIA),{\displaystyle \det(xI-A),} que es un polinomio mónico en x de grado n si A es una matriz n × n . Este polinomio es el polinomio característico de A.

Definición formal

Considere unnorte×norte{\displaystyle n\times n}matrizA.{\displaystyle A.}El polinomio característico deA,{\displaystyle A,}denotado porpagA(t),{\displaystyle p_{A}(t),}es el polinomio definido por [ 5 ]pagA(t)=det(tIA){\displaystyle p_{A}(t)=\det(tI-A)} dóndeI{\displaystyle I}denota elnorte×norte{\displaystyle n\times n}matriz identidad .

Algunos autores definen el polinomio característico comodet(AtI).{\displaystyle \det(A-tI).}Ese polinomio difiere del definido aquí por un signo.(1)norte,{\displaystyle (-1)^{n},}por lo que no hay diferencia para propiedades como tener como raíces los valores propios deA{\displaystyle A}; sin embargo, la definición anterior siempre da un polinomio mónico , mientras que la definición alternativa es mónica solo cuandonorte{\displaystyle n}es par.

Ejemplos

Para calcular el polinomio característico de la matriz A=(2110).{\displaystyle A={\begin{pmatrix}2&1\\-1&0\end{pmatrix}}.} Se calcula el determinante de lo siguiente: tIA=(t211t0){\displaystyle tI-A={\begin{pmatrix}t-2&-1\\1&t-0\end{pmatrix}}} y se descubrió que(t2)t1(1)=t22t+1,{\displaystyle (t-2)t-1(-1)=t^{2}-2t+1\,\!,}el polinomio característico deA.{\displaystyle A.}

Otro ejemplo utiliza funciones hiperbólicas de un ángulo hiperbólico φ . Para la matriz, tome A=(aporrear(φ)sinh(φ)sinh(φ)aporrear(φ)).{\displaystyle A={\begin{pmatrix}\cosh(\varphi )&\sinh(\varphi )\\\sinh(\varphi )&\cosh(\varphi )\end{pmatrix}}.} Su polinomio característico es det(tIA)=(taporrear(φ))2sinh2(φ)=t22t aporrear(φ)+1=(tmiφ)(tmiφ).{\displaystyle \det(tI-A)=(t-\cosh(\varphi ))^{2}-\sinh ^{2}(\varphi )=t^{2}-2t\ \cosh(\varphi )+1=(t-e^{\varphi })(t-e^{-\varphi }).}

Propiedades

El polinomio característicopagA(t){\displaystyle p_{A}(t)}de unnorte×norte{\displaystyle n\times n}matrizA{\displaystyle A}es mónico (su coeficiente principal es1{\displaystyle 1}) y su grado esnorte.{\displaystyle n.}El hecho más importante sobre el polinomio característico ya se mencionó en el párrafo motivacional : los valores propios deA{\displaystyle A}son precisamente las raíces depagA(t){\displaystyle p_{A}(t)}(esto también se cumple para el polinomio mínimo deA,{\displaystyle A,}pero su grado puede ser menor quenorte{\displaystyle n}). Todos los coeficientes del polinomio característico son expresiones polinómicas en las entradas de la matriz. En particular, su coeficiente constante det0{\displaystyle t^{0}}esdet(A)=(1)nortedet(A),{\displaystyle \det(-A)=(-1)^{n}\det(A),}el coeficiente detnorte{\displaystyle t^{n}}es 1, y el coeficiente detnorte1{\displaystyle t^{n-1}}es tr(− A ) = −tr( A ) , donde tr( A ) es la traza deA.{\displaystyle A.}(Los signos que se dan aquí corresponden a la definición formal dada en la sección anterior; para la definición alternativa, estos serían en cambiodet(A){\displaystyle \det(A)}y (−1) n – 1 tr( A ) respectivamente. [ 6 ] )

Para un2×2{\displaystyle 2\times 2}matrizA,{\displaystyle A,}El polinomio característico viene dado por: t2tr(A)t+det(A).{\displaystyle t^{2}-\operatorname {tr} (A)t+\det(A).}

Utilizando el lenguaje del álgebra exterior , el polinomio característico de unnorte×norte{\displaystyle n\times n}matrizA{\displaystyle A}puede expresarse como pagA(t)=k=0nortetnortek(1)ktr(kA){\displaystyle p_{A}(t)=\sum _{k=0}^{n}t^{n-k}(-1)^{k}\operatorname {tr} \left(\textstyle \bigwedge ^{k}A\right)} dóndetr(kA){\textstyle \operatorname {tr} \left(\bigwedge ^{k}A\right)}es el rastro de lak{\displaystyle k}el poder exterior deA,{\displaystyle A,}que tiene dimensión(nortek).{\textstyle {\binom {n}{k}}.}Esta traza se puede calcular como la suma de todos los menores principales deA{\displaystyle A}de tamañok.{\displaystyle k.}El algoritmo recursivo de Faddeev-LeVerrier calcula estos coeficientes de forma más eficiente .

Cuando la característica del campo de los coeficientes es0,{\displaystyle 0,}Cada una de esas trazas puede calcularse alternativamente como un único determinante, el de lak×k{\displaystyle k\times k}matriz, tr(kA)=1k¡|trAk100trA2trAk20trAk1trAk21trAktrAk1trA| .{\displaystyle \operatorname {tr} \left(\textstyle \bigwedge ^{k}A\right)={\frac {1}{k!}}{\begin{vmatrix}\operatorname {tr} A&k-1&0&\cdots &0\\\operatorname {tr} A^{2}&\operatorname {tr} A&k-2&\cdots &0\\\vdots &\vdots &&\ddots &\vdots \\\operatorname {tr} A^{k-1}&\operatorname {tr} A^{k-2}&&\cdots &1\\\operatorname {tr} A^{k}&\operatorname {tr} A^{k-1}&&\cdots &\operatorname {tr} A\end{vmatrix}}~.}

El teorema de Cayley-Hamilton establece que al reemplazart{\displaystyle t}porA{\displaystyle A}en el polinomio característico (interpretando las potencias resultantes como potencias matriciales y el término constante)do{\displaystyle c}comodo{\displaystyle c}veces la matriz identidad) produce la matriz cero . En términos informales, cada matriz satisface su propia ecuación característica. Esta afirmación es equivalente a decir que el polinomio mínimo deA{\displaystyle A}divide el polinomio característico deA.{\displaystyle A.}

Dos matrices similares tienen el mismo polinomio característico. Sin embargo, lo contrario no es cierto en general: dos matrices con el mismo polinomio característico no tienen por qué ser similares.

La matrizA{\displaystyle A}y su transpuesta tienen el mismo polinomio característico.A{\displaystyle A}es similar a una matriz triangular si y solo si su polinomio característico puede factorizarse completamente en factores lineales sobreK{\displaystyle K}(Lo mismo ocurre con el polinomio mínimo en lugar del polinomio característico). En este casoA{\displaystyle A}es similar a una matriz en forma normal de Jordan .

Polinomio característico de un producto de dos matrices

SiA{\displaystyle A}yB{\displaystyle B}son dos cuadradosnorte×norte{\displaystyle n\times n}matrices y luego polinomios característicos deAB{\displaystyle AB}yBA{\displaystyle BA}coincidir: pagAB(t)=pagBA(t).{\displaystyle p_{AB}(t)=p_{BA}(t).\,}

Prueba

Siλ{\displaystyle \lambda }es un valor propio generalizado distinto de cero deAB{\displaystyle AB}de multiplicidad algebraicak{\displaystyle k}, yv{\displaystyle v}pertenece al núcleo de(BAλ)k{\displaystyle (BA-\lambda )^{k}}, entoncesAv{\displaystyle Av}pertenece al núcleo de(ABλ)k{\displaystyle (AB-\lambda )^{k}}, por lo tanto, los autoespacios generalizados no nulos deAB{\displaystyle AB}yBA{\displaystyle BA}tienen la misma dimensión. Por lo tanto, dado queAB{\displaystyle AB}yBA{\displaystyle BA}son ambosnorte×norte{\displaystyle n\times n}, los autoespacios generalizados restantes, con autovalor 0, tienen la misma dimensión. Por lo tantoAB{\displaystyle AB}yBA{\displaystyle BA}tienen el mismo polinomio característico, porque todos los autovalores generalizados son iguales, con las mismas multiplicidades algebraicas.

En términos más generales, siA{\displaystyle A}es una matriz de ordenmetro×norte{\displaystyle m\times n}yB{\displaystyle B}es una matriz de ordennorte×metro,{\displaystyle n\times m,}entoncesAB{\displaystyle AB}esmetro×metro{\displaystyle m\times m}yBA{\displaystyle BA}esnorte×norte{\displaystyle n\times n}matriz, y uno tiene pagBA(t)=tnortemetropagAB(t).{\displaystyle p_{BA}(t)=t^{n-m}p_{AB}(t).\,}

Para probar esto, uno puede suponernorte>metro,{\displaystyle n>m,}mediante intercambio, si es necesario,A{\displaystyle A}yB.{\displaystyle B.}Luego, al limitarA{\displaystyle A}en la parte inferior pornortemetro{\displaystyle n-m}filas de ceros yB{\displaystyle B}a la derecha, por,nortemetro{\displaystyle n-m}columnas de ceros, uno obtiene dosnorte×norte{\displaystyle n\times n}matricesA{\displaystyle A^{\prime }}yB{\displaystyle B^{\prime }}de tal manera queBA=BA{\displaystyle B^{\prime }A^{\prime }=BA}yAB{\displaystyle A^{\prime }B^{\prime }}es igual aAB{\displaystyle AB}bordeado pornortemetro{\displaystyle n-m}filas y columnas de ceros. El resultado se deduce del caso de matrices cuadradas, comparando los polinomios característicos deAB{\displaystyle A^{\prime }B^{\prime }}yAB.{\displaystyle AB.}

Polinomio característico de A k

Siλ{\displaystyle \lambda }es un valor propio de una matriz cuadradaA{\displaystyle A}con vector propiov,{\displaystyle \mathbf {v} ,}entoncesλk{\displaystyle \lambda ^{k}}es un valor propio deAk{\displaystyle A^{k}}porque Akv=Ak1Av=λAk1v==λkv.{\displaystyle A^{k}{\textbf {v}}=A^{k-1}A{\textbf {v}}=\lambda A^{k-1}{\textbf {v}}=\dots =\lambda ^{k}{\textbf {v}}.}

Se puede demostrar que las multiplicidades también coinciden, y esto se generaliza a cualquier polinomio en lugar deincógnitak{\displaystyle x^{k}}: [ 7 ]

Teorema SeaA{\displaystyle A}ser un cuadradonorte×norte{\displaystyle n\times n}matriz y dejarF(t){\displaystyle f(t)}sea ​​un polinomio. Si el polinomio característico deA{\displaystyle A}tiene una factorización pagA(t)=(tλ1)(tλ2)(tλnorte){\displaystyle p_{A}(t)=(t-\lambda _{1})(t-\lambda _{2})\cdots (t-\lambda _{n})} entonces el polinomio característico de la matrizF(A){\displaystyle f(A)}es dado por pagF(A)(t)=(tF(λ1))(tF(λ2))(tF(λnorte)).{\displaystyle p_{f(A)}(t)=(t-f(\lambda _{1}))(t-f(\lambda _{2}))\cdots (t-f(\lambda _{n})).}

Es decir, la multiplicidad algebraica deλ{\displaystyle \lambda }enF(A){\displaystyle f(A)}es igual a la suma de multiplicidades algebraicas deλ{\displaystyle \lambda '}enA{\displaystyle A}encimaλ{\displaystyle \lambda '}de tal manera queF(λ)=λ.{\displaystyle f(\lambda ')=\lambda .} En particular,tr(F(A))=i=1norteF(λi){\displaystyle \operatorname {tr} (f(A))=\textstyle \sum _{i=1}^{n}f(\lambda _{i})}ydet(F(A))=i=1norteF(λi).{\displaystyle \operatorname {det} (f(A))=\textstyle \prod _{i=1}^{n}f(\lambda _{i}).} Aquí un polinomioF(t)=t3+1,{\displaystyle f(t)=t^{3}+1,}por ejemplo, se evalúa en una matrizA{\displaystyle A}simplemente comoF(A)=A3+I.{\displaystyle f(A)=A^{3}+I.}

El teorema se aplica a matrices y polinomios sobre cualquier cuerpo o anillo conmutativo . [ 8 ] Sin embargo, la suposición de quepagA(t){\displaystyle p_{A}(t)} has a factorization into linear factors is not always true, unless the matrix is over an algebraically closed field such as the complex numbers.

Proof

This proof only applies to matrices and polynomials over complex numbers (or any algebraically closed field). In that case, the characteristic polynomial of any square matrix can be always factorized as pA(t)=(tλ1)(tλ2)(tλn){\displaystyle p_{A}(t)=\left(t-\lambda _{1}\right)\left(t-\lambda _{2}\right)\cdots \left(t-\lambda _{n}\right)} where λ1,λ2,,λn{\displaystyle \lambda _{1},\lambda _{2},\ldots ,\lambda _{n}} are the eigenvalues of A,{\displaystyle A,} possibly repeated. Moreover, the Jordan decomposition theorem guarantees that any square matrix A{\displaystyle A} can be decomposed as A=S1US,{\displaystyle A=S^{-1}US,} where S{\displaystyle S} is an invertible matrix and U{\displaystyle U} is upper triangular with λ1,,λn{\displaystyle \lambda _{1},\ldots ,\lambda _{n}} on the diagonal (with each eigenvalue repeated according to its algebraic multiplicity). (The Jordan normal form has stronger properties, but these are sufficient; alternatively the Schur decomposition can be used, which is less popular but somewhat easier to prove).

Let f(t)=iαiti.{\textstyle f(t)=\sum _{i}\alpha _{i}t^{i}.} Then f(A)=αi(S1US)i=αiS1USS1USS1US=αiS1UiS=S1(αiUi)S=S1f(U)S.{\displaystyle {\begin{aligned}f(A)&=\textstyle \sum \alpha _{i}(S^{-1}US)^{i}\\&=\textstyle \sum \alpha _{i}S^{-1}USS^{-1}US\cdots S^{-1}US\\&=\textstyle \sum \alpha _{i}S^{-1}U^{i}S\\&=S^{-1}(\textstyle \sum \alpha _{i}U^{i})S\\&=S^{-1}f(U)S.\end{aligned}}} For an upper triangular matrix U{\displaystyle U} with diagonal λ1,,λn,{\displaystyle \lambda _{1},\dots ,\lambda _{n},} the matrix Ui{\displaystyle U^{i}} is upper triangular with diagonal λ1i,,λni{\displaystyle \lambda _{1}^{i},\dots ,\lambda _{n}^{i}} in Ui,{\displaystyle U^{i},} and hence f(U){\displaystyle f(U)} is upper triangular with diagonal f(λ1),,f(λn).{\displaystyle f\left(\lambda _{1}\right),\dots ,f\left(\lambda _{n}\right).} Therefore, the eigenvalues of f(U){\displaystyle f(U)} are f(λ1),,f(λn).{\displaystyle f(\lambda _{1}),\dots ,f(\lambda _{n}).} Since f(A)=S1f(U)S{\displaystyle f(A)=S^{-1}f(U)S} is similar to f(U),{\displaystyle f(U),} it has the same eigenvalues, with the same algebraic multiplicities.

Secular function and secular equation

Secular function

The term secular function has been used for what is now called characteristic polynomial (in some literature the term secular function is still used). The term comes from the fact that the characteristic polynomial was used to calculate secular perturbations (on a time scale of a century, that is, slow compared to annual motion) of planetary orbits, according to Lagrange's theory of oscillations.

Secular equation

Secular equation may have several meanings.

  • In linear algebra it is sometimes used in place of characteristic equation.
  • In astronomy it is the algebraic or numerical expression of the magnitude of the inequalities in a planet's motion that remain after the inequalities of a short period have been allowed for.[9]
  • In molecular orbital calculations relating to the energy of the electron and its wave function it is also used instead of the characteristic equation.

For general associative algebras

The above definition of the characteristic polynomial of a matrix AMn(F){\displaystyle A\in M_{n}(F)} with entries in a field F{\displaystyle F} generalizes without any changes to the case when F{\displaystyle F} is just a commutative ring. Garibaldi (2004) defines the characteristic polynomial for elements of an arbitrary finite-dimensional (associative, but not necessarily commutative) algebra over a field F{\displaystyle F}y demuestra las propiedades estándar del polinomio característico en esta generalidad.

Complejidad teórica: cálculo mediante multiplicación rápida de matrices

Es posible calcular el polinomio característico de forma rápida mediante el uso de algoritmos rápidos de multiplicación de matrices en el tiempoO(norteω){\displaystyle O({n^{\omega }})}paraω{\displaystyle \omega }ligeramente por encima de 2,37. Los algoritmos respectivos se dan en [ 10 ] .

Véase también

Referencias

  1. ^ Guillemin, Ernst (1953). Introducción a la teoría de circuitos . Wiley. págs.366  , 541. ISBN 0471330663.{{cite book}}: Incompatibilidad de ISBN/Fecha ( ayuda )
  2. Forsythe, George E.; Motzkin, Theodore (enero de 1952). "Una extensión de la transformación de Gauss para mejorar la condición de los sistemas de ecuaciones lineales" (PDF) . Matemáticas de la computación . 6 (37): 18–34 . doi : 10.1090/S0025-5718-1952-0048162-0 . Consultado el 3 de octubre de 2020 .
  3. Frank, Evelyn (1946). "Sobre los ceros de polinomios con coeficientes complejos" (PDF) . Boletín de la Sociedad Matemática Americana . 52 (2): 144– 157. doi : 10.1090/S0002-9904-1946-08526-2 .
  4. "Polinomio característico de un grafo – Wolfram MathWorld" . Consultado el 26 de agosto de 2011 .
  5. Steven Roman (1992). Álgebra lineal avanzada (2.ª ed.). Springer. pág . 137. ISBN   3540978372.
  6. Teorema 4 en estos apuntes de clase
  7. Horn, Roger A.; Johnson, Charles R. (2013). Análisis matricial (2.ª ed.). Cambridge University Press . págs. 108–109, Sección 2.4.2. ISBN  978-0-521-54823-6.
  8. Lang, Serge (1993). Álgebra . Nueva York: Springer. pág. 567, Teorema 3.10. ISBN 978-1-4613-0041-0OCLC 852792828 
  9. "ecuación secular" . Consultado el 21 de enero de 2010 .
  10. Neiger, Vincent; Pernet, Clément (diciembre de 2021). "Cálculo determinista del polinomio característico en el tiempo de la multiplicación de matrices" . Journal of Complexity . 67. arXiv : 2010.04662 . doi : 10.1016/j.jco.2021.101572 .
  • TS Blyth y EF Robertson (1998) Álgebra lineal básica , pág. 149, Springer ISBN 3-540-76122-5.
  • John B. Fraleigh y Raymond A. Beauregard (1990) Álgebra lineal, segunda edición, p. 246, Addison-Wesley ISBN 0-201-11949-8.
  • Garibaldi, Skip (2004), "El polinomio característico y el determinante no son construcciones ad hoc", American Mathematical Monthly , 111 (9): 761–778 , arXiv : math/0203276 , doi : 10.2307/4145188 , JSTOR 4145188 , MR 2104048  
  • Werner Greub (1974) Álgebra lineal , 4.ª edición, págs. 120-125 , Springer, ISBN 0-387-90110-8.
  • Paul C. Shields (1980) Álgebra lineal elemental , 3.ª edición, pág. 274, Worth Publishers ISBN 0-87901-121-1.
  • Gilbert Strang (1988) Álgebra lineal y sus aplicaciones, 3.ª edición, pág. 246, Brooks/Cole ISBN 0-15-551005-3.