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)}Que una matriz tenga una factorización en factores lineales no siempre es cierto, a menos que la matriz esté sobre un campo algebraicamente cerrado como los números complejos.

Prueba

Esta demostración solo se aplica a matrices y polinomios sobre números complejos (o cualquier cuerpo algebraicamente cerrado). En ese caso, el polinomio característico de cualquier matriz cuadrada siempre se puede factorizar como pagA(t)=(tλ1)(tλ2)(tλnorte){\displaystyle p_{A}(t)=\left(t-\lambda _{1}\right)\left(t-\lambda _{2}\right)\cdots \left(t-\lambda _{n}\right)} dóndeλ1,λ2,,λnorte{\displaystyle \lambda _{1},\lambda _{2},\ldots ,\lambda _{n}}son los valores propios deA,{\displaystyle A,}posiblemente repetido. Además, el teorema de descomposición de Jordan garantiza que cualquier matriz cuadradaA{\displaystyle A}puede descomponerse comoA=S1US,{\displaystyle A=S^{-1}US,}dóndeS{\displaystyle S}es una matriz invertible yU{\displaystyle U}es triangular superior conλ1,,λnorte{\displaystyle \lambda _{1},\ldots ,\lambda _{n}}en la diagonal (con cada autovalor repetido según su multiplicidad algebraica). (La forma normal de Jordan tiene propiedades más fuertes, pero estas son suficientes; alternativamente, se puede usar la descomposición de Schur , que es menos popular pero algo más fácil de demostrar).

DejarF(t)=iαiti.{\textstyle f(t)=\sum _{i}\alpha _{i}t^{i}.} Entonces 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}}} Para una matriz triangular superiorU{\displaystyle U}con diagonalλ1,,λnorte,{\displaystyle \lambda _{1},\dots ,\lambda _{n},}la matrizUi{\displaystyle U^{i}}es triangular superior con diagonalλ1i,,λnortei{\displaystyle \lambda _{1}^{i},\dots ,\lambda _{n}^{i}}enUi,{\displaystyle U^{i},} y por lo tantoF(U){\displaystyle f(U)}es triangular superior con diagonalF(λ1),,F(λnorte).{\displaystyle f\left(\lambda _{1}\right),\dots ,f\left(\lambda _{n}\right).} Por lo tanto, los valores propios deF(U){\displaystyle f(U)}sonF(λ1),,F(λnorte).{\displaystyle f(\lambda _{1}),\dots ,f(\lambda _{n}).} DesdeF(A)=S1F(U)S{\displaystyle f(A)=S^{-1}f(U)S}es similar aF(U),{\displaystyle f(U),}Tiene los mismos autovalores, con las mismas multiplicidades algebraicas.

Función secular y ecuación secular

Función secular

El término función secular se ha utilizado para referirse a lo que ahora se denomina polinomio característico (aunque en algunos textos todavía se emplea el término función secular). Este término proviene del hecho de que el polinomio característico se utilizaba para calcular las perturbaciones seculares (en una escala temporal de un siglo, es decir, lentas en comparación con el movimiento anual) de las órbitas planetarias, según la teoría de oscilaciones de Lagrange .

ecuación secular

La ecuación secular puede tener varios significados.

  • En álgebra lineal , a veces se utiliza en lugar de ecuación característica.
  • En astronomía es la expresión algebraica o numérica de la magnitud de las desigualdades en el movimiento de un planeta que permanecen después de que se hayan tenido en cuenta las desigualdades de un período corto. [ 9 ]

Para álgebras asociativas generales

La definición anterior del polinomio característico de una matrizAMETROnorte(F){\displaystyle A\in M_{n}(F)}con entradas en un campoF{\displaystyle F}generaliza sin ningún cambio al caso cuandoF{\displaystyle F}es simplemente un anillo conmutativo . Garibaldi (2004) define el polinomio característico para elementos de un álgebra arbitraria de dimensión finita ( asociativa , pero no necesariamente conmutativa) sobre un cuerpo.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.