Articulo de referencia

Jordan forma normal

[ ⌜ λ 1 1 ⌝ λ 1 1 ⌞ λ 1 ⌟ ⌜ λ 2 1 ⌝ ⌞ λ 2 ⌟ [ λ 3 ] ⋱ ⌜ λ norte 1 ⌝ ⌞ λ norte ⌟ ] {\displaystyle \left[{\begin{array}{lcrlrcclr}{\color {red}\ulcorner }\lambda _{1}&1&{\color {r...

[λ11λ11λ1λ21λ2[λ3]λnorte1λnorte]{\displaystyle \left[{\begin{array}{lcrlrcclr}{\color {red}\ulcorner }\lambda _{1}&1&{\color {red}\urcorner }\\&\lambda _{1}&1\,\,\,\,\,\\{\color {red}\llcorner }&&\lambda _{1}{\color {red}\lrcorner }\\&&&{\color {red}\ulcorner }\lambda _{2}&1\,\,{\color {red}\urcorner }\\&&&{\color {red}\llcorner }&\lambda _{2}{\color {red}\lrcorner }\\&&&&&{\color {red}[}\lambda _{3}{\color {red}]}\\&&&&&&\ddots \\&&&&&&&{\color {red}\ulcorner }\lambda _{n}&1\,\,{\color {red}\urcorner }\\&&&&&&&{\color {red}\llcorner }&\lambda _{n}{\color {red}\lrcorner }\\\end{array}}\right]}
Ejemplo de una matriz en forma normal de Jordan. Todas las entradas de la matriz que no se muestran son cero. Los cuadrados resaltados se conocen como "bloques de Jordan". Cada bloque de Jordan contiene un número λi en su diagonal principal y unos directamente encima de ella. Los λi son los autovalores de la matriz; no es necesario que sean distintos .

En álgebra lineal , una forma normal de Jordan , también conocida como forma canónica de Jordan , [ 1 ] [ 2 ] es una matriz triangular superior de una forma particular llamada matriz de Jordan que representa un operador lineal en un espacio vectorial de dimensión finita con respecto a alguna base . Dicha matriz tiene cada entrada no nula fuera de la diagonal igual a 1, inmediatamente encima de la diagonal principal (en la superdiagonal ), y con entradas diagonales idénticas a la izquierda y debajo de ellas. 

Sea V un espacio vectorial sobre un cuerpo K. Entonces, existe una base respecto de la cual la matriz tiene la forma requerida si y solo si todos los autovalores de la matriz pertenecen a K , o equivalentemente, si el polinomio característico del operador se descompone en factores lineales sobre K. Esta condición siempre se satisface si K es algebraicamente cerrado (por ejemplo, si es el cuerpo de los números complejos ). Las entradas diagonales de la forma normal son los autovalores (del operador), y el número de veces que aparece cada autovalor se denomina multiplicidad algebraica del autovalor. [ 3 ] [ 4 ] [ 5 ]

Si el operador viene dado originalmente por una matriz cuadrada M , entonces su forma normal de Jordan también se denomina forma normal de Jordan de M. Cualquier matriz cuadrada tiene una forma normal de Jordan si el campo de coeficientes se extiende a uno que contenga todos los autovalores de la matriz. A pesar de su nombre, la forma normal para una M dada no es del todo única, ya que es una matriz diagonal por bloques formada por bloques de Jordan , cuyo orden no es fijo; es convencional agrupar los bloques correspondientes al mismo autovalor, pero no se impone ningún orden entre los autovalores, ni entre los bloques correspondientes a un autovalor dado, aunque estos últimos podrían, por ejemplo, ordenarse por tamaño decreciente. [ 3 ] [ 4 ] [ 5 ]

La descomposición de Jordan-Chevalley es particularmente simple con respecto a una base para la cual el operador toma su forma normal de Jordan. La forma diagonal para matrices diagonalizables , por ejemplo matrices normales , es un caso especial de la forma normal de Jordan. [ 6 ] [ 7 ] [ 8 ]

La forma normal de Jordan recibe su nombre de Camille Jordan , quien enunció por primera vez el teorema de descomposición de Jordan en 1870. [ 9 ]

Descripción general

Notación

Algunos libros de texto los tienen en la subdiagonal ; es decir, inmediatamente debajo de la diagonal principal en lugar de en la superdiagonal. Los autovalores siguen estando en la diagonal principal. [ 10 ] [ 11 ]

Motivación

Una matriz A de n × n es diagonalizable si y solo si la suma de las dimensiones de sus autoespacios es n . O, equivalentemente, si y solo si A tiene n autovectores linealmente independientes . No todas las matrices son diagonalizables; las matrices que no son diagonalizables se denominan matrices defectuosas . Considere la siguiente matriz:

A=[5421011111301112].{\displaystyle A=\left[{\begin{array}{*{20}{r}}5&4&2&1\\[2pt]0&1&-1&-1\\[2pt]-1&-1&3&0\\[2pt]1&1&-1&2\end{array}}\right].}

Incluyendo la multiplicidad, los autovalores de A son λ = 1, 2, 4, 4. La dimensión del autoespacio correspondiente al autovalor 4 es 1 (y no 2), por lo que A no es diagonalizable. Sin embargo, existe una matriz invertible P tal que J = P −1 AP , donde

J=[1000020000410004].{\displaystyle J={\begin{bmatrix}1&0&0&0\\[2pt]0&2&0&0\\[2pt]0&0&4&1\\[2pt]0&0&0&4\end{bmatrix}}.}

La matrizJ{\displaystyle J}es casi diagonal. Esta es la forma normal de Jordan de A. La sección Ejemplo a continuación detalla el cálculo.

Matrices complejas

En general, una matriz cuadrada compleja A es similar a una matriz diagonal por bloques.

J=[J1Jpag]{\displaystyle J={\begin{bmatrix}J_{1}&\;&\;\\\;&\ddots &\;\\\;&\;&J_{p}\end{bmatrix}}}

donde cada bloque J i es una matriz cuadrada de la forma

Ji=[λi1λi1λi].{\displaystyle J_{i}={\begin{bmatrix}\lambda _{i}&1&\;&\;\\\;&\lambda _{i}&\ddots &\;\\\;&\;&\ddots &1\\\;&\;&\;&\lambda _{i}\end{bmatrix}}.}

Existe una matriz invertible P tal que P −1 AP = J es tal que los únicos elementos no nulos de J están en la diagonal y la superdiagonal. J se denomina forma normal de Jordan de A. Cada J i se denomina bloque de Jordan de A. En un bloque de Jordan dado, cada elemento de la superdiagonal es 1.

Suponiendo este resultado, podemos deducir las siguientes propiedades:

  • Contando las multiplicidades, los valores propios de J , y por lo tanto de A , son las entradas de la diagonal.
  • Dado un valor propio λ i , su multiplicidad geométrica es la dimensión de ker( A λ i I ), donde I es la matriz identidad , y es el número de bloques de Jordan correspondientes a λ i . [ 12 ]
  • La suma de los tamaños de todos los bloques de Jordan correspondientes a un valor propio λ i es su multiplicidad algebraica . [ 12 ]
  • A es diagonalizable si y solo si, para cada autovalor λ de A , sus multiplicidades geométrica y algebraica coinciden. En particular, los bloques de Jordan en este caso sonMatrices de 1 × 1 ; es decir, escalares.
  • El bloque de Jordan correspondiente a λ tiene la forma λI + N , donde N es una matriz nilpotente definida como N ij = δ i , j 1 (donde δ es la delta de Kronecker ). La nilpotencia de N puede aprovecharse al calcular f ( A ), donde f es una función analítica compleja. Por ejemplo, en principio, la forma de Jordan podría dar una expresión analítica para la exponencial exp( A ).
  • El número de bloques de Jordan correspondientes a λ i de tamaño al menos j es dim ker( Aλ i I ) j − dim ker( Aλ i I ) j −1 . Por lo tanto, el número de bloques de Jordan de tamaño j es
    2oscuroker(AλiI)joscuroker(AλiI)j+1oscuroker(AλiI)j1{\displaystyle 2\dim \ker(A-\lambda _{i}I)^{j}-\dim \ker(A-\lambda _{i}I)^{j+1}-\dim \ker(A-\lambda _{i}I)^{j-1}}
  • Dado un valor propio λ i , su multiplicidad en el polinomio mínimo es el tamaño de su bloque de Jordan más grande.

Ejemplo

Consideremos la matrizA{\displaystyle A}a partir del ejemplo de la sección anterior. La forma normal de Jordan se obtiene mediante alguna transformación de similitud :

PAG1APAG=J;{\displaystyle P^{-1}AP=J;}eso es,APAG=PAGJ.{\displaystyle AP=PJ.}

DejarPAG{\displaystyle P}tienen vectores columnapagi{\displaystyle p_{i}},i=1,,4{\displaystyle i=1,\ldots ,4}, entonces

A[pag1pag2pag3pag4]=[pag1pag2pag3pag4][1000020000410004]=[pag12pag24pag3pag3+4pag4].{\displaystyle A{\begin{bmatrix}p_{1}&p_{2}&p_{3}&p_{4}\end{bmatrix}}={\begin{bmatrix}p_{1}&p_{2}&p_{3}&p_{4}\end{bmatrix}}{\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}={\begin{bmatrix}p_{1}&2p_{2}&4p_{3}&p_{3}+4p_{4}\end{bmatrix}}.}

Vemos que

(A1I)pag1=0{\displaystyle (A-1I)p_{1}=0}
(A2I)pag2=0{\displaystyle (A-2I)p_{2}=0}
(A4I)pag3=0{\displaystyle (A-4I)p_{3}=0}
(A4I)pag4=pag3.{\displaystyle (A-4I)p_{4}=p_{3}.}

Parai=1,2,3{\displaystyle i=1,2,3}tenemospagiker(AλiI){\displaystyle p_{i}\in \ker(A-\lambda _{i}I)}, eso es,pagi{\displaystyle p_{i}}es un vector propio deA{\displaystyle A}correspondiente al valor propioλi{\displaystyle \lambda _{i}}. Parai=4{\displaystyle i=4}multiplicando ambos lados por(A4I){\displaystyle (A-4I)}da

(A4I)2pag4=(A4I)pag3.{\displaystyle (A-4I)^{2}p_{4}=(A-4I)p_{3}.}

Pero(A4I)pag3=0{\displaystyle (A-4I)p_{3}=0}, entonces

(A4I)2pag4=0.{\displaystyle (A-4I)^{2}p_{4}=0.}

De este modo,pag4ker(A4I)2.{\displaystyle p_{4}\in \ker(A-4I)^{2}.}

Vectores comopag4{\displaystyle p_{4}}se denominan autovectores generalizados de A.

Ejemplo: Obtención de la forma normal

Este ejemplo muestra cómo calcular la forma normal de Jordan de una matriz dada.

Consideremos la matriz

A=[5421011111301112]{\displaystyle A=\left[{\begin{array}{rrrr}5&4&2&1\\0&1&-1&-1\\-1&-1&3&0\\1&1&-1&2\end{array}}\right]}

lo cual se menciona al principio del artículo.

El polinomio característico de A es

χ(λ)=det(λIA)=λ411λ3+42λ264λ+32=(λ1)(λ2)(λ4)2.{\displaystyle {\begin{aligned}\chi (\lambda )&=\det(\lambda I-A)\\&=\lambda ^{4}-11\lambda ^{3}+42\lambda ^{2}-64\lambda +32\\&=(\lambda -1)(\lambda -2)(\lambda -4)^{2}.\,\end{aligned}}}

Esto demuestra que los autovalores son 1, 2, 4 y 4, según la multiplicidad algebraica. El autoespacio correspondiente al autovalor 1 se puede encontrar resolviendo la ecuación Av = 1 v . Está generado por el vector columna v = (−1, 1, 0, 0) T. De manera similar, el autoespacio correspondiente al autovalor 2 está generado por w = (1, −1, 0, 1) T. Finalmente, el autoespacio correspondiente al autovalor 4 también es unidimensional (aunque se trata de un autovalor doble) y está generado por x = (1, 0, −1, 1) T. Por lo tanto, la multiplicidad geométrica (es decir, la dimensión del autoespacio del autovalor dado) de cada uno de los tres autovalores es uno. Por lo tanto, los dos autovalores iguales a 4 corresponden a un único bloque de Jordan, y la forma normal de Jordan de la matriz A es la suma directa.

J=J1(1)J1(2)J2(4)=[1000020000410004].{\displaystyle J=J_{1}(1)\oplus J_{1}(2)\oplus J_{2}(4)={\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}.}

Hay tres cadenas de Jordan . Dos tienen longitud uno: { v } y { w }, que corresponden a los autovalores 1 y 2, respectivamente. Hay una cadena de longitud dos que corresponde al autovalor 4. Para encontrar esta cadena, calcula

ker(A4I)2=durar{[1000],[1011]}{\displaystyle \ker(A-4I)^{2}=\operatorname {span} \,\left\{{\begin{bmatrix}1\\0\\0\\0\end{bmatrix}},\left[{\begin{array}{r}1\\0\\-1\\1\end{array}}\right]\right\}}

donde yo estoy elMatriz identidad de 4 × 4. Elija un vector en el espacio generado anteriormente que no esté en el núcleo de A − 4 I ; por ejemplo, y = (1,0,0,0) T . Ahora, ( A − 4 I ) y = x y ( A − 4 I ) x = 0 , por lo que { y , x } es una cadena de longitud dos que corresponde al valor propio 4.

La matriz de transición P tal que P −1 AP = J se forma colocando estos vectores uno al lado del otro de la siguiente manera:

PAG=[vwincógnitay]=[1111110000100110].{\displaystyle P=\left[{\begin{array}{c|c|c|c}v&w&x&y\end{array}}\right]=\left[{\begin{array}{rrrr}-1&1&1&1\\1&-1&0&0\\0&0&-1&0\\0&1&1&0\end{array}}\right].}

Un cálculo muestra que la ecuación P −1 AP = J efectivamente se cumple.

PAG1APAG=J=[1000020000410004].{\displaystyle P^{-1}AP=J={\begin{bmatrix}1&0&0&0\\0&2&0&0\\0&0&4&1\\0&0&0&4\end{bmatrix}}.}

Si hubiéramos intercambiado el orden en que aparecen los vectores de la cadena, es decir, cambiando el orden de v , w y { x , y } juntos, los bloques de Jordan se intercambiarían. Sin embargo, las formas de Jordan son equivalentes.

Autovectores generalizados

Dado un valor propio λ , cada bloque de Jordan correspondiente da lugar a una cadena de Jordan de vectores linealmente independientes p i , i = 1, ..., b , donde b es el tamaño del bloque de Jordan. El generador , o vector principal , p b de la cadena es un vector propio generalizado tal que(AλI)bpagb=0{\displaystyle (A-\lambda I)^{b}p_{b}=0} . El vectorpag1=(AλI)b1pagb{\displaystyle p_{1}=(A-\lambda I)^{b-1}p_{b}} es un vector propio ordinario correspondiente a λ . En general, p i es una preimagen de p i −1 bajoAλI{\displaystyle A-\lambda I} . Por lo tanto, el vector principal genera la cadena mediante la multiplicación porAλI{\displaystyle A-\lambda I} . [ 13 ] [ 2 ] Por lo tanto, la afirmación de que toda matriz cuadrada A puede ponerse en forma normal de Jordan es equivalente a la afirmación de que el espacio vectorial subyacente tiene una base compuesta de cadenas de Jordan.

Una prueba

Damos una demostración por inducción de que cualquier matriz cuadrada de valores complejos A puede ponerse en forma normal de Jordan. Dado que se puede demostrar [ 14 ] que el espacio vectorial subyacente es la suma directa de subespacios invariantes asociados con los valores propios, se puede suponer que A tiene un solo valor propio λ . El caso 1 × 1 es trivial. Sea A una matriz n × n . El rango deAλI{\displaystyle A-\lambda I} , denotado porCorrió(AλI){\displaystyle \operatorname {Ran} (A-\lambda I)} , es un subespacio invariante de A . Además, dado que λ es un valor propio de A , la dimensión deCorrió(AλI){\displaystyle \operatorname {Ran} (A-\lambda I)} , r , es estrictamente menor que n , por lo tanto, por la hipótesis inductiva,Corrió(AλI){\displaystyle \operatorname {Ran} (A-\lambda I)}Tiene una base { p 1 , ..., p r } compuesta por cadenas de Jordan.

A continuación, consideremos el núcleo , es decir, el subespacio ker(AλI){\displaystyle \ker(A-\lambda I)} . Si

Corrió(AλI)ker(AλI)={0},{\displaystyle \operatorname {Ran} (A-\lambda I)\cap \ker(A-\lambda I)=\{0\},}

El resultado deseado se deduce inmediatamente del teorema de rango-nulidad . (Este sería el caso, por ejemplo, si A fuera hermitiana ).

De lo contrario, si

Q=Corrió(AλI)ker(AλI){0},{\displaystyle Q=\operatorname {Ran} (A-\lambda I)\cap \ker(A-\lambda I)\neq \{0\},}

Sea la dimensión de Q sr . Cada vector en Q es un vector propio, por lo queCorrió(AλI){\displaystyle \operatorname {Ran} (A-\lambda I)}Debe contener s cadenas de Jordan correspondientes a s autovectores linealmente independientes. Por lo tanto, la base { p 1 , ..., p r } debe contener s vectores, digamos { p 1 , ..., p s }, que son vectores principales de estas cadenas de Jordan. Podemos "extender las cadenas" tomando las preimágenes de estos vectores principales. (Este es el paso clave). Sea q i tal que

(AλI)qi=pagi para i=1,,s.{\displaystyle \;(A-\lambda I)q_{i}=p_{i}{\mbox{ for }}i=1,\ldots ,s.}

Finalmente, podemos elegir cualquier base para

ker(AλI)/Q{\displaystyle \ker(A-\lambda I)/Q}

y luego elevar a vectores { z 1 , ..., z t } en ker(AλI){\displaystyle \ker(A-\lambda I)}Cada z i forma una cadena de Jordan de longitud 1. Solo necesitamos demostrar que la unión de { p 1 , ..., p r }, { z 1 , ..., z t } y { q 1 , ..., q s } forma una base para el espacio vectorial.

Por el teorema de rango-nulidad ,oscuro(ker(AλI)))=norter{\displaystyle \dim(\ker(A-\lambda I)))=n-r} , entoncest=norters{\displaystyle t=n-r-s} , y por lo tanto el número de vectores en la base potencial es igual a n. Para demostrar la independencia lineal, supongamos que alguna combinación lineal de los vectores es 0. AplicandoAλI,{\displaystyle A-\lambda I,}obtenemos alguna combinación lineal de p i , donde los q i se convierten en vectores principales entre los p i . De la independencia lineal de p i , se deduce que los coeficientes de los vectores q i deben ser cero. Además, ninguna combinación lineal no trivial de los z i puede ser igual a una combinación lineal de p i , porque entonces pertenecería aCorrió(AλI){\displaystyle \operatorname {Ran} (A-\lambda I)}y por lo tanto Q , lo cual es imposible por la construcción de z i . Por lo tanto, los coeficientes de z i también serán 0. Esto deja en la combinación lineal original solo los términos p i , que se suponen linealmente independientes, por lo que sus coeficientes también deben ser cero. Hemos encontrado una base compuesta por cadenas de Jordan, y esto muestra que A puede expresarse en forma normal de Jordan.

Unicidad

Se puede demostrar que la forma normal de Jordan de una matriz A dada es única hasta el orden de los bloques de Jordan.

Conocer las multiplicidades algebraicas y geométricas de los autovalores no es suficiente para determinar la forma normal de Jordan de A. Suponiendo que se conoce la multiplicidad algebraica m ( λ ) de un autovalor λ , la estructura de la forma de Jordan se puede determinar analizando los rangos de las potencias ( AλI ) m ( λ ) . Para ver esto, supongamos que una matriz n × n A tiene solo un autovalor λ . Entonces m ( λ ) = n . El entero más pequeño k 1 tal que

(AλI)k1=0{\displaystyle (A-\lambda I)^{k_{1}}=0}

es el tamaño del bloque de Jordan más grande en la forma de Jordan de A. (Este número k 1 también se llama índice de λ . Véase la discusión en una sección posterior). El rango de

(AλI)k11{\displaystyle (A-\lambda I)^{k_{1}-1}}

es el número de bloques de Jordan de tamaño k 1 . De manera similar, el rango de

(AλI)k12{\displaystyle (A-\lambda I)^{k_{1}-2}}

es el doble del número de bloques de Jordan de tamaño k 1 más el número de bloques de Jordan de tamaño k 1 − 1. El caso general es similar.

Esto puede utilizarse para demostrar la unicidad de la forma de Jordan. Sean J₁ y J₂ dos formas normales de Jordan de A. Entonces , J₁ y J₂ son similares y tienen el mismo espectro, incluyendo las multiplicidades algebraicas de los autovalores. El procedimiento descrito en el párrafo anterior puede utilizarse para determinar la estructura de estas matrices. Dado que el rango de una matriz se conserva mediante una transformación de semejanza, existe una biyección entre los bloques de Jordan de J₁ y J₂ . Esto demuestra la unicidad de la afirmación.

Matrices reales

Si A es una matriz real, su forma de Jordan aún puede ser no real. En lugar de representarla con valores propios complejos y unos en la superdiagonal, como se discutió anteriormente, existe una matriz real invertible P tal que P −1 AP = J es una matriz diagonal por bloques real donde cada bloque es un bloque de Jordan real. [ 15 ] Un bloque de Jordan real es idéntico a un bloque de Jordan complejo (si el valor propio correspondienteλi{\displaystyle \lambda _{i}}es real), o es una matriz de bloques en sí misma, que consta de bloques de 2×2 (para valores propios no reales)λi=ai+ibi{\displaystyle \lambda _{i}=a_{i}+ib_{i}}con multiplicidad algebraica dada) de la forma

doi=[aibibiai]{\displaystyle C_{i}=\left[{\begin{array}{rr}a_{i}&-b_{i}\\b_{i}&a_{i}\\\end{array}}\right]}

y describir la multiplicación porλi{\displaystyle \lambda _{i}}en el plano complejo. Los bloques superdiagonales son matrices identidad de 2×2 y, por lo tanto, en esta representación las dimensiones de la matriz son mayores que la forma de Jordan compleja. El bloque de Jordan real completo viene dado por

Ji=[doiIdoiIdoi].{\displaystyle J_{i}={\begin{bmatrix}C_{i}&I&&\\&C_{i}&\ddots &\\&&\ddots &I\\&&&C_{i}\end{bmatrix}}.}

Esta forma de Jordan real es consecuencia de la forma de Jordan compleja. Para una matriz real, los autovectores no reales y los autovectores generalizados siempre pueden elegirse para formar pares conjugados complejos . Al tomar la parte real y la imaginaria (combinación lineal del vector y su conjugado), la matriz adquiere esta forma con respecto a la nueva base.

Matrices con entradas en un campo

La reducción de Jordan se puede extender a cualquier matriz cuadrada M cuyos elementos pertenezcan a un cuerpo K. El resultado establece que cualquier M se puede escribir como una suma D + N , donde D es semisimple , N es nilpotente y DN = ND . Esto se conoce como la descomposición de Jordan-Chevalley . Siempre que K contiene los autovalores de M , en particular cuando K es algebraicamente cerrado , la forma normal se puede expresar explícitamente como la suma directa de bloques de Jordan.

De forma similar al caso en que K son los números complejos, conocer las dimensiones de los núcleos de ( M λI ) k para 1 ≤ km , donde m es la multiplicidad algebraica del valor propio λ , permite determinar la forma de Jordan de M . Podemos ver el espacio vectorial subyacente V como un K [ x ] -módulo considerando la acción de x sobre V como una aplicación de M y extendiéndola por K -linealidad. Entonces, los polinomios ( xλ ) k son los divisores elementales de M , y la forma normal de Jordan se ocupa de representar M en términos de bloques asociados a los divisores elementales.

La demostración de la forma normal de Jordan se suele llevar a cabo como una aplicación al anillo K [ x ] del teorema de estructura para módulos finitamente generados sobre un dominio ideal principal , del cual es un corolario.

Consecuencias

Se puede observar que la forma normal de Jordan es esencialmente un resultado de clasificación para matrices cuadradas, y como tal, varios resultados importantes del álgebra lineal pueden considerarse consecuencias del mismo.

Teorema de mapeo espectral

Utilizando la forma normal de Jordan, el cálculo directo proporciona un teorema de mapeo espectral para el cálculo funcional polinomial : Sea A una matriz n × n con valores propios λ 1 , ..., λ n , entonces para cualquier polinomio p , p ( A ) tiene valores propios p ( λ 1 ), ..., p ( λ n ).

Polinomio característico

El polinomio característico de A espagA(λ)=det(λIA){\displaystyle p_{A}(\lambda )=\det(\lambda I-A)}. Las matrices similares tienen el mismo polinomio característico. Por lo tanto,pagA(λ)=pagJ(λ)=i(λλi)metroi{\textstyle p_{A}(\lambda )=p_{J}(\lambda )=\prod _{i}(\lambda -\lambda _{i})^{m_{i}}}, dóndeλi{\displaystyle \lambda _{i}}es la i -ésima raíz depagJ{\textstyle p_{J}}ymetroi{\displaystyle m_{i}}es su multiplicidad, porque este es claramente el polinomio característico de la forma de Jordan de A.

Teorema de Cayley-Hamilton

El teorema de Cayley-Hamilton afirma que toda matriz A satisface su ecuación característica: si p es el polinomio característico de A , entoncespagA(A)=0{\displaystyle p_{A}(A)=0}Esto se puede demostrar mediante un cálculo directo en la forma de Jordan, ya que siλi{\displaystyle \lambda _{i}}es un valor propio de multiplicidadmetroi{\displaystyle m_{i}}, luego su bloque JordanJi{\displaystyle J_{i}}claramente satisface(JiλiI)metroi=0{\displaystyle (J_{i}-\lambda _{i}I)^{m_{i}}=0}Como los bloques diagonales no se afectan entre sí,i{\displaystyle i}el bloque diagonal de(AλiI)metroi{\displaystyle (A-\lambda _{i}I)^{m_{i}}}es(JiλiI)metroi{\displaystyle (J_{i}-\lambda _{i}I)^{m_{i}}}; por esopagA(A)=i(AλiI)metroi=0{\textstyle p_{A}(A)=\prod _{i}(A-\lambda _{i}I)^{m_{i}}=0}.

Se puede suponer que la forma de Jordan existe sobre un campo que extiende el campo base de la matriz, por ejemplo sobre el campo de descomposición de p ; esta extensión de campo no cambia la matriz p ( A ) de ninguna manera.

polinomio mínimo

El polinomio mínimo P de una matriz cuadrada A es el único polinomio mónico de grado mínimo, m , tal que P ( A ) = 0. Alternativamente, el conjunto de polinomios que anulan una A dada forman un ideal I en C [ x ], el dominio ideal principal de polinomios con coeficientes complejos. El elemento mónico que genera I es precisamente P.

Sean λ 1 , ..., λ q los autovalores distintos de A , y s i el tamaño del bloque de Jordan más grande correspondiente a λ i . Es claro a partir de la forma normal de Jordan que el polinomio mínimo de A tiene grado Σ s i .

Si bien la forma normal de Jordan determina el polinomio mínimo, lo contrario no es cierto. Esto nos lleva al concepto de divisores elementales . Los divisores elementales de una matriz cuadrada A son los polinomios característicos de sus bloques de Jordan. Los factores del polinomio mínimo m son los divisores elementales de mayor grado que corresponden a autovalores distintos.

El grado de un divisor elemental es igual al tamaño del bloque de Jordan correspondiente y, por lo tanto, a la dimensión del subespacio invariante correspondiente. Si todos los divisores elementales son lineales, A es diagonalizable.

Descomposiciones de subespacios invariantes

La forma de Jordan de una matriz A de n × n es diagonal por bloques y, por lo tanto, proporciona una descomposición del espacio euclidiano n- dimensional en subespacios invariantes de A. Cada bloque de Jordan J i corresponde a un subespacio invariante X i . Simbólicamente, ponemos

donorte=i=1kincógnitai{\displaystyle \mathbb {C} ^{n}=\bigoplus _{i=1}^{k}X_{i}}

donde cada X i es el intervalo de la cadena de Jordan correspondiente, y k es el número de cadenas de Jordan.

También se puede obtener una descomposición ligeramente diferente mediante la forma de Jordan. Dado un valor propio λ i , el tamaño de su bloque de Jordan correspondiente más grande s i se llama índice de λ i y se denota por v ( λ i ) . (Por lo tanto, el grado del polinomio mínimo es la suma de todos los índices). Definimos un subespacio Y i por

Yi=ker(λiIA)v(λi).{\displaystyle Y_{i}=\ker(\lambda _{i}I-A)^{v(\lambda _{i})}.}

Esto da como resultado la descomposición

donorte=i=1lYi{\displaystyle \mathbb {C} ^{n}=\bigoplus _{i=1}^{l}Y_{i}}

donde l es el número de autovalores distintos de A. Intuitivamente, agrupamos los subespacios invariantes de bloque de Jordan que corresponden al mismo autovalor. En el caso extremo en que A es un múltiplo de la matriz identidad, tenemos k = n y l = 1.

La proyección sobre Y i y a lo largo de todos los demás Y j ( ji ) se denomina proyección espectral de A en v i y se suele denotar por P ( λ i  ; A ) . Las proyecciones espectrales son mutuamente ortogonales en el sentido de que P ( λ i  ; A ) P (v j  ; A ) = 0 si ij . Además, conmutan con A y su suma es la matriz identidad. Sustituyendo cada v i en la matriz de Jordan J por uno y poniendo a cero todas las demás entradas se obtiene P (v i  ; J ) , además, si UJU −1 es la transformación de similitud tal que A = UJU −1 entonces P ( λ i  ; A ) = UP ( λ i  ; J ) U −1 . No se limitan a dimensiones finitas. Véase más adelante su aplicación a operadores compactos y en cálculo funcional holomorfo para una discusión más general.

Al comparar las dos descomposiciones, observe que, en general, lk . Cuando A es normal, los subespacios X i en la primera descomposición son unidimensionales y mutuamente ortogonales. Este es el teorema espectral para operadores normales. La segunda descomposición se generaliza más fácilmente para operadores compactos generales en espacios de Banach.

Aquí podría ser interesante observar algunas propiedades del índice, ν ( λ ) . De manera más general, para un número complejo λ , su índice se puede definir como el menor entero no negativo ν ( λ ) tal que

ker(AλI)ν(λ)=ker(AλI)metro,metroν(λ).{\displaystyle \ker(A-\lambda I)^{\nu (\lambda )}=\ker(A-\lambda I)^{m},\;\forall m\geq \nu (\lambda ).}

Entonces ν (v) > 0 si y solo si λ es un valor propio de A. En el caso de dimensión finita, ν (v) ≤ la multiplicidad algebraica de v .

Forma normal plana

La forma de Jordan se utiliza para encontrar una forma normal de matrices, salvo conjugación, de modo que las matrices normales conformen una variedad algebraica de grado fijo bajo en el espacio de matrices ambiente.

Los conjuntos de representantes de clases de conjugación de matrices para la forma normal de Jordan o las formas canónicas racionales en general no constituyen subespacios lineales o afines en los espacios de matrices ambiente.

Vladimir Arnold planteó [ 16 ] un problema: encontrar una forma canónica de matrices sobre un cuerpo para el cual el conjunto de representantes de clases de conjugación de matrices sea una unión de subespacios lineales afines (planos). En otras palabras, mapear inyectivamente el conjunto de clases de conjugación de matrices de vuelta al conjunto inicial de matrices de manera que la imagen de esta incrustación —el conjunto de todas las matrices normales— tenga el grado más bajo posible; es una unión de subespacios lineales desplazados.

Fue resuelto para cuerpos algebraicamente cerrados por Peteris Daugulis. [ 17 ] La construcción de una forma normal plana definida de una matriz comienza considerando su forma normal de Jordan.

funciones matriciales

La iteración de la cadena de Jordan motiva diversas extensiones a contextos más abstractos. Para matrices finitas, se obtienen funciones matriciales; esto puede extenderse a operadores compactos y al cálculo funcional holomorfo, como se describe más adelante.

La forma normal de Jordan es la más conveniente para el cálculo de funciones matriciales (aunque puede que no sea la mejor opción para cálculos computacionales). Sea f ( z ) una función analítica de un argumento complejo. Al aplicar la función a un bloque de Jordan J de n × n con valor propio λ, se obtiene una matriz triangular superior:

F(J)=[F(λ)F(λ)F(λ)2F(norte1)(λ)(norte1)¡0F(λ)F(λ)F(norte2)(λ)(norte2)¡000F(λ)F(λ)0000F(λ)],{\displaystyle f(J)={\begin{bmatrix}f(\lambda )&f'(\lambda )&{\tfrac {f''(\lambda )}{2}}&\cdots &{\tfrac {f^{(n-1)}(\lambda )}{(n-1)!}}\\0&f(\lambda )&f'(\lambda )&\cdots &{\tfrac {f^{(n-2)}(\lambda )}{(n-2)!}}\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&0&f(\lambda )&f'(\lambda )\\0&0&0&0&f(\lambda )\end{bmatrix}},}

de modo que los elementos de la k -ésima superdiagonal de la matriz resultante seanF(k)(λ)k¡{\displaystyle {\tfrac {f^{(k)}(\lambda )}{k!}}}Para una matriz de forma normal de Jordan general, la expresión anterior se aplicará a cada bloque de Jordan.

El siguiente ejemplo muestra la aplicación a la función de potencia f ( z )  = z n : 

[λ110000λ110000λ100000λ210000λ2]norte=[λ1norte(norte1)λ1norte1(norte2)λ1norte2000λ1norte(norte1)λ1norte10000λ1norte00000λ2norte(norte1)λ2norte10000λ2norte],{\displaystyle {\begin{bmatrix}\lambda _{1}&1&0&0&0\\0&\lambda _{1}&1&0&0\\0&0&\lambda _{1}&0&0\\0&0&0&\lambda _{2}&1\\0&0&0&0&\lambda _{2}\end{bmatrix}}^{n}={\begin{bmatrix}\lambda _{1}^{n}&{\tbinom {n}{1}}\lambda _{1}^{n-1}&{\tbinom {n}{2}}\lambda _{1}^{n-2}&0&0\\0&\lambda _{1}^{n}&{\tbinom {n}{1}}\lambda _{1}^{n-1}&0&0\\0&0&\lambda _{1}^{n}&0&0\\0&0&0&\lambda _{2}^{n}&{\tbinom {n}{1}}\lambda _{2}^{n-1}\\0&0&0&0&\lambda _{2}^{n}\end{bmatrix}},}

donde los coeficientes binomiales se definen como(nortek)=i=1knorte+1ii{\textstyle {\binom {n}{k}}=\prod _{i=1}^{k}{\frac {n+1-i}{i}}}Para n entero positivo se reduce a la definición estándar de los coeficientes. Para n negativo la identidad(nortek)=(1)k(norte+k1k){\textstyle {\binom {-n}{k}}=(-1)^{k}{\binom {n+k-1}{k}}}puede ser de utilidad.

Operadores compactos

Un resultado análogo a la forma normal de Jordan se cumple para operadores compactos en un espacio de Banach . Nos restringimos a operadores compactos porque cada punto x en el espectro de un operador compacto T es un autovalor; la única excepción es cuando x es el punto límite del espectro. Esto no es cierto para operadores acotados en general. Para dar una idea de esta generalización, primero reformulamos la descomposición de Jordan en el lenguaje del análisis funcional .

Cálculo funcional holomorfo

Sea X un espacio de Banach, L ( X ) los operadores acotados en X , y σ ( T ) el espectro de TL ( X ). El cálculo funcional holomorfo se define de la siguiente manera:

Fijemos un operador acotado T. Consideremos la familia Hol( T ) de funciones complejas que es holomorfa en algún conjunto abierto G que contiene a σ ( T ). Sea Γ = { γi } una colección finita de curvas de Jordan tales que σ ( T ) se encuentra en el interior de Γ, definimos f ( T ) por

F(T)=12πiΓF(z)(zT)1dz.{\displaystyle f(T)={\frac {1}{2\pi i}}\int _{\Gamma }f(z)(z-T)^{-1}\,dz.}

El conjunto abierto G podría variar con f y no necesariamente ser conexo. La integral se define como el límite de las sumas de Riemann, como en el caso escalar. Aunque la integral tiene sentido para f continua , nos restringimos a funciones holomorfas para aplicar la maquinaria de la teoría clásica de funciones (por ejemplo, la fórmula integral de Cauchy). La suposición de que σ ( T ) se encuentra dentro de Γ asegura que f ( T ) esté bien definida; no depende de la elección de Γ. El cálculo funcional es la aplicación Φ de Hol( T ) a L ( X ) dada por

Φ(F)=F(T).{\displaystyle \;\Phi (f)=f(T).}

Necesitaremos las siguientes propiedades de este cálculo funcional:

  1. Φ extiende el cálculo funcional polinomial.
  2. Se cumple el teorema de mapeo espectral : σ ( f ( T )) = f ( σ ( T )).
  3. Φ es un homomorfismo de álgebra.

El caso de dimensión finita

En el caso de dimensión finita, σ ( T ) = { λ i } es un conjunto discreto finito en el plano complejo. Sea e i la función que es 1 en algún entorno abierto de λ i y 0 en cualquier otro lugar. Por la propiedad 3 del cálculo funcional, el operador

mii(T){\displaystyle e_{i}(T)}

es una proyección. Además, sea ν i el índice de λ i y

F(z)=(zλi)νi.{\displaystyle f(z)=(z-\lambda _{i})^{\nu _{i}}.}

El teorema de mapeo espectral nos dice que

F(T)mii(T)=(Tλi)νimii(T){\displaystyle f(T)e_{i}(T)=(T-\lambda _{i})^{\nu _{i}}e_{i}(T)}

tiene espectro {0}. Por la propiedad 1, f ( T ) se puede calcular directamente en la forma de Jordan, y por inspección, vemos que el operador f ( T ) e i ( T ) es la matriz cero.

Por la propiedad 3, f ( T ) e i ( T ) = e i ( T ) f ( T ). Por lo tanto, e i ( T ) es precisamente la proyección sobre el subespacio

Corriómii(T)=ker(Tλi)νi.{\displaystyle \operatorname {Ran} e_{i}(T)=\ker(T-\lambda _{i})^{\nu _{i}}.}

La relación

imii=1{\displaystyle \sum _{i}e_{i}=1}

implica

donorte=iCorriómii(T)=iker(Tλi)νi{\displaystyle \mathbb {C} ^{n}=\bigoplus _{i}\;\operatorname {Ran} e_{i}(T)=\bigoplus _{i}\ker(T-\lambda _{i})^{\nu _{i}}}

donde el índice i recorre los distintos valores propios de T. Esta es la descomposición en subespacio invariante.

donorte=iYi{\displaystyle \mathbb {C} ^{n}=\bigoplus _{i}Y_{i}}

dado en una sección anterior. Cada e i ( T ) es la proyección sobre el subespacio generado por las cadenas de Jordan correspondientes a λ i y a lo largo de los subespacios generados por las cadenas de Jordan correspondientes a v j para ji . En otras palabras, e i ( T ) = P ( λ i ; T ). Esta identificación explícita de los operadores e i ( T ) a su vez proporciona una forma explícita de cálculo funcional holomorfo para matrices:

Para todo f ∈ Hol( T ),
F(T)=λiσ(T)k=0νi1F(k)k¡(Tλi)kmii(T).{\displaystyle f(T)=\sum _{\lambda _{i}\in \sigma (T)}\sum _{k=0}^{\nu _{i}-1}{\frac {f^{(k)}}{k!}}(T-\lambda _{i})^{k}e_{i}(T).}

Nótese que la expresión de f ( T ) es una suma finita porque, en cada entorno de v i , hemos elegido el desarrollo en serie de Taylor de f centrado en v i .

Polos de un operador

Sea T un operador acotado y λ un punto aislado de σ ( T ). (Como se indicó anteriormente, cuando T es compacto, cada punto de su espectro es un punto aislado, excepto posiblemente el punto límite 0).

El punto λ se denomina polo del operador T de orden ν si la función resolvente R T definida por

RT(λ)=(λT)1{\displaystyle R_{T}(\lambda )=(\lambda -T)^{-1}}

tiene un polo de orden ν en λ .

Demostraremos que, en el caso de dimensión finita, el orden de un valor propio coincide con su índice. Este resultado también es válido para operadores compactos.

Consideremos la región anular A centrada en el autovalor λ con un radio ε suficientemente pequeño tal que la intersección del disco abierto B ε ( λ ) y σ ( T ) sea { λ }. La función resolvente R T es holomorfa en A. Extendiendo un resultado de la teoría clásica de funciones, R T tiene una representación en serie de Laurent en A :

RT(z)=ametro(λz)metro{\displaystyle R_{T}(z)=\sum _{-\infty }^{\infty }a_{m}(\lambda -z)^{m}}

dónde

ametro=12πido(λz)metro1(zT)1dz{\displaystyle a_{-m}=-{\frac {1}{2\pi i}}\int _{C}(\lambda -z)^{m-1}(z-T)^{-1}dz}y C es un círculo pequeño centrado en λ . 

Según la discusión anterior sobre el cálculo funcional,

ametro=(λT)metro1miλ(T){\displaystyle a_{-m}=-(\lambda -T)^{m-1}e_{\lambda }(T)}dóndemiλ{\displaystyle e_{\lambda }}es 1 enBε(λ){\displaystyle B_{\varepsilon }(\lambda )}y 0 en cualquier otro lugar.

Pero hemos demostrado que el entero positivo más pequeño m tal que

ametro0{\displaystyle a_{-m}\neq 0}yal=0lmetro{\displaystyle a_{-l}=0\;\;\forall \;l\geq m}

es precisamente el índice de λ , ν ( λ ). En otras palabras, la función R T tiene un polo de orden ν ( λ ) en λ .

Análisis numérico

Si la matriz A tiene múltiples valores propios, o está cerca de una matriz con múltiples valores propios, entonces su forma normal de Jordan es muy sensible a las perturbaciones. Consideremos, por ejemplo, la matriz

A=[11ε1].{\displaystyle A={\begin{bmatrix}1&1\\\varepsilon &1\end{bmatrix}}.}

Si ε = 0, entonces la forma normal de Jordan es simplemente

[1101].{\displaystyle {\begin{bmatrix}1&1\\0&1\end{bmatrix}}.}

Sin embargo, para ε ≠ 0, la forma normal de Jordan es

[1+ε001ε].{\displaystyle {\begin{bmatrix}1+{\sqrt {\varepsilon }}&0\\0&1-{\sqrt {\varepsilon }}\end{bmatrix}}.}

Este mal condicionamiento dificulta enormemente el desarrollo de un algoritmo numérico robusto para la forma normal de Jordan, ya que el resultado depende críticamente de si dos autovalores se consideran iguales. Por esta razón, la forma normal de Jordan suele evitarse en el análisis numérico ; la descomposición estable de Schur [ 18 ] o los pseudoespectros [ 19 ] son ​​mejores alternativas.

Véase también

Notas

  1. Shilov define el término «forma canónica de Jordan» y, en una nota al pie, indica que « forma normal de Jordan» es sinónimo. Estos términos a veces se abrevian como « forma de Jordan» . (Shilov) El término «forma canónica clásica» también se utiliza en ocasiones en el contexto de este artículo. (James y James, 1976)
  2. ^ Holt y Rumynin (2009 , pág.9 ) 
  3. ^ Beauregard y Fraleigh (1973 , págs. 310-316) 
  4. ^ Préstamo Golub y Van (1996 , pág. 355) 
  5. ^ Nering (1970 , págs. 118-127 ) 
  6. ^ Beauregard y Fraleigh (1973 , págs. 270-274) 
  7. Préstamo Golub y Van (1996 , p. 353) 
  8. ^ Nering (1970 , págs. 113-118) 
  9. ^ Brechenmacher, "Histoire du théorème de Jordan de la décomposition matricielle (1870-1930). Formes de représentation et méthodes de décomposition" , Tesis, 2007
  10. Cullen (1966 , pág. 114) 
  11. Franklin (1968 , pág. 122) 
  12. 1 2 Horn & Johnson (1985 , §3.2.1)
  13. Bronson (1970 , págs. 189, 194) 
  14. Roe Goodman y Nolan R. Wallach, Representaciones e invariantes de grupos clásicos , Cambridge UP 1998, Apéndice B.1.
  15. Horn y Johnson (1985 , Teorema 3.4.5)
  16. Arnold, Vladimir I. (2004), "1998-25", en Arnold, Vladimir I. (ed.), Arnold's Problems , Berlín: Springer-Verlag, p.  127, doi : 10.1007/b138219 , ISBN 3-540-20614-0, MR 2078115 Véase también el comentario de la página 613.
  17. Peteris Daugulis (2012), "Una parametrización de conjuntos de órbitas de conjugación de matrices como uniones de planos afines", Álgebra lineal y sus aplicaciones , 436 (3): 709– 721, arXiv : 1110.0907 , doi : 10.1016/j.laa.2011.07.032 , S2CID 119649768 
  18. Véase Golub & Van Loan (2014), §7.6.5; o Golub & Wilkinson (1976) para más detalles.
  19. Véase Golub & Van Loan (2014), §7.9

Referencias

  • Beauregard, Raymond A.; Fraleigh, John B. (1973), A First Course In Linear Algebra: with Optional Introduction to Groups, Rings, and Fields , Boston: Houghton Mifflin Co. , ISBN 0-395-14017-X
  • Bronson, Richard (1970), Métodos matriciales: una introducción , Nueva York: Academic Press , LCCN 70097490 
  • Cullen, Charles G. (1966), Matrices y transformaciones lineales , Reading: Addison-Wesley , LCCN 66021267 
  • Dunford, N.; Schwartz, JT (1958), Operadores lineales, Parte I: Teoría general , Interscience
  • Finkbeiner II, Daniel T. (1978), Introducción a las matrices y transformaciones lineales (3.ª  ed.), WH Freeman and Company
  • Franklin, Joel N. (1968), Teoría de matrices , Englewood Cliffs: Prentice-Hall , LCCN 68016345 
  • Golub, Gene H.; Van Loan, Charles F. (1996), Matrix Computations (3.ª  ed.), Baltimore: Johns Hopkins University Press , ISBN 0-8018-5414-8
  • Golub, Gene H.; Wilkinson, JH (1976), "Sistemas propios mal condicionados y el cálculo de la forma normal de Jordan", SIAM Review , 18 (4): 578–619 , doi : 10.1137/1018113
  • Holt, Derek; Rumynin, Dmitriy (2009), Álgebra I – Álgebra lineal avanzada (MA251) Apuntes de clase (PDF)
  • Horn, Roger A.; Johnson, Charles R. (1985), Análisis matricial , Cambridge University Press , ISBN 978-0-521-38632-6
  • James, Glenn; James, Robert C. (1976), Diccionario de matemáticas (2.ª  ed.), Van Nostrand Reinhold
  • MacLane, Saunders ; Birkhoff, Garrett (1967), Álgebra , Macmillan Publishers
  • Michel, Anthony N.; Herget, Charles J. (1993), Álgebra aplicada y análisis funcional , Dover Publications
  • Nering, Evar D. (1970), Álgebra lineal y teoría de matrices (2.ª  ed.), Nueva York: Wiley , LCCN 76091646 
  • Shafarevich, IR ; Remizov, AO (2012), Álgebra lineal y geometría , Springer , ISBN 978-3-642-30993-9
  • Shilov, Georgi E. (1977), Álgebra lineal , Dover Publications
  • Artículo sobre la forma canónica de Jordan en mathworld.wolfram.com