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:
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
La matrizes 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.
donde cada bloque J i es una matriz cuadrada de la forma
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
- 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 partir del ejemplo de la sección anterior. La forma normal de Jordan se obtiene mediante alguna transformación de similitud :
- eso es,
Dejartienen vectores columna,, entonces
Vemos que
Paratenemos, eso es,es un vector propio decorrespondiente al valor propio. Paramultiplicando ambos lados porda
Pero, entonces
De este modo,
Vectores comose 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
lo cual se menciona al principio del artículo.
El polinomio característico de A es
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.
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
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:
Un cálculo muestra que la ecuación P −1 AP = J efectivamente se cumple.
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 . El vector es un vector propio ordinario correspondiente a λ . En general, p i es una preimagen de p i −1 bajo . Por lo tanto, el vector principal genera la cadena mediante la multiplicación por . [ 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 de , denotado por , es un subespacio invariante de A . Además, dado que λ es un valor propio de A , la dimensión de , r , es estrictamente menor que n , por lo tanto, por la hipótesis inductiva, Tiene una base { p 1 , ..., p r } compuesta por cadenas de Jordan.
A continuación, consideremos el núcleo , es decir, el subespacio . Si
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
Sea la dimensión de Q s ≤ r . Cada vector en Q es un vector propio, por lo queDebe 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
Finalmente, podemos elegir cualquier base para
y luego elevar a vectores { z 1 , ..., z t } en 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 , , entonces , 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. Aplicando 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 ay 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
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
es el número de bloques de Jordan de tamaño k 1 . De manera similar, el rango de
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 correspondientees real), o es una matriz de bloques en sí misma, que consta de bloques de 2×2 (para valores propios no reales)con multiplicidad algebraica dada) de la forma
y describir la multiplicación poren 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
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 ≤ k ≤ m , 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 es. Las matrices similares tienen el mismo polinomio característico. Por lo tanto,, dóndees la i -ésima raíz deyes 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 , entoncesEsto se puede demostrar mediante un cálculo directo en la forma de Jordan, ya que sies un valor propio de multiplicidad, luego su bloque Jordanclaramente satisfaceComo los bloques diagonales no se afectan entre sí,el bloque diagonal dees; por eso.
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
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
Esto da como resultado la descomposición
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 ( j ≠ i ) 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 i ≠ j . 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, l ≤ k . 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
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:
de modo que los elementos de la k -ésima superdiagonal de la matriz resultante seanPara 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 :
donde los coeficientes binomiales se definen comoPara n entero positivo se reduce a la definición estándar de los coeficientes. Para n negativo la identidadpuede 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 T ∈ L ( 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
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
Necesitaremos las siguientes propiedades de este cálculo funcional:
- Φ extiende el cálculo funcional polinomial.
- Se cumple el teorema de mapeo espectral : σ ( f ( T )) = f ( σ ( T )).
- Φ 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
es una proyección. Además, sea ν i el índice de λ i y
El teorema de mapeo espectral nos dice que
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
La relación
implica
donde el índice i recorre los distintos valores propios de T. Esta es la descomposición en subespacio invariante.
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 j ≠ i . 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 ),
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
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 :
dónde
- y C es un círculo pequeño centrado en λ .
Según la discusión anterior sobre el cálculo funcional,
- dóndees 1 eny 0 en cualquier otro lugar.
Pero hemos demostrado que el entero positivo más pequeño m tal que
- y
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
Si ε = 0, entonces la forma normal de Jordan es simplemente
Sin embargo, para ε ≠ 0, la forma normal de Jordan es
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
- ↑ 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)
- ^ Holt y Rumynin (2009 , pág.9 )
- ^ Beauregard y Fraleigh (1973 , págs. 310-316)
- ^ Préstamo Golub y Van (1996 , pág. 355)
- ^ Nering (1970 , págs. 118-127 )
- ^ Beauregard y Fraleigh (1973 , págs. 270-274)
- ↑ Préstamo Golub y Van (1996 , p. 353)
- ^ Nering (1970 , págs. 113-118)
- ^ 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
- ↑ Cullen (1966 , pág. 114)
- ↑ Franklin (1968 , pág. 122)
- 1 2 Horn & Johnson (1985 , §3.2.1)
- ↑ Bronson (1970 , págs. 189, 194)
- ↑ Roe Goodman y Nolan R. Wallach, Representaciones e invariantes de grupos clásicos , Cambridge UP 1998, Apéndice B.1.
- ↑ Horn y Johnson (1985 , Teorema 3.4.5)
- ↑ 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.
- ↑ 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
- ↑ Véase Golub & Van Loan (2014), §7.6.5; o Golub & Wilkinson (1976) para más detalles.
- ↑ 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
- Álgebra lineal
- teoría matricial
- Formas normales de la matriz
- Descomposiciones matriciales