Forma canónica de Jordan

De testwiki
Ir a la navegación Ir a la búsqueda
Ejemplo de matriz en forma normal de Jordan. Los bloques en gris se conocen como bloques de Jordan. Nótese que la λi en diferentes bloques puede ser igual.

En álgebra lineal, la forma canónica de Jordan es la forma de la matriz de un endomorfismo de un espacio vectorial en cierta base asociada a la descomposición en suma directa de subespacios invariantes bajo dicho endomorfismo. Dicha forma canónica consistirá en que la matriz estará formada por "bloques de Jordan" en la diagonal y bloques de ceros fuera de ella.

Introducción

Sea f un endomorfismo sobre un K-espacio vectorial V de dimensión n>1 (f:VV). Si el polinomio característico de f se factoriza completamente sobre el cuerpo K (es decir, K es el cuerpo de descomposición del polinomio característico de la matriz), existe una base donde la aplicación lineal viene dada por una "matriz de m bloques" (mn) con la siguiente forma. Plantilla:Ecuación Donde cada submatriz 𝐀k es un bloque de Jordan. Plantilla:Ecuación Donde λ1,...λm son raíces del polinomio característico (valores propios), y k=1mnk=n,
Cuando f es diagonalizable, vale que m=n y k:nk=1, por lo que la forma canónica de Jordan de la matriz es una matriz diagonal.

Ejemplo

Considérese la situación de una matriz diagonalizable. Una matriz cuadrada es diagonalizable si la suma de las dimensiones de los espacios propios (eigenspaces) es el número de filas o columnas de la matriz. Examinemos la matriz siguiente: Plantilla:Ecuación Tenemos valores propios de A que son solo λ = 5, 5, 5, 5. Ahora bien, la dimensión del núcleo de A5Id4 es 1 (donde Id4 representa la matriz identidad de orden 4), por lo tanto A no es diagonalizable. Sin embargo, podemos construir la forma de Jordan de esta matriz. Dado que la dimensión es 1, sabemos que la forma de Jordan está compuesta de solo un bloque de Jordan, es decir, la forma de Jordan de A es: Plantilla:Ecuación Obsérvese que J puede escribirse como 5Id4+N, donde N es una matriz nilpotente. Puesto que ahora tenemos A similar a dicha matriz simple, podremos realizar cálculos que involucren a A usando la forma de Jordan, lo que en muchos casos puede simplificar el cálculo. Por ejemplo, calcular potencias de matrices es significativamente más sencillo usando la forma de Jordan.

Cálculo de la forma de Jordan

1) Supongamos que se quiere diagonalizar la siguiente matriz

A=(121010113)

Primeros calculamos el polinomio característico y vemos si A es diagonalizable.

pA(λ)=det(λIA)=λ3tr(A)λ2+(A11+A22+A33)λdet(A)==λ33λ2+(3+41)λ(4)=λ33λ2+4=(λ+1)(λ2)2

Hacemos pA(λ)=0 y obtenemos los autovalores 1 y 2, este último de multiplicidad algebraica 2. Para ver si A es diagonalizable buscamos los autovectores, estos conforman las bases de los espacios ker(AλI) (es decir los espacios propios Eλ). Para el primer valor propio obtenemos A(1)I:

A+I=(221000114)(110001000){x+y=0z=0{x=yz=0

de donde se resuelve que ker(AλI)={(x,x,0):x}. En particular, un vector de la base (el más simple sin contar el nulo) es (1,1,0).

Ahora para el otro autovalor

A2I=(121030111)(101010000){x+z=0y=0{x=zy=0

Es decir que obtuvimos ker(AλI)={(z,0,z):z} pero este espacio es unidimensional luego no alcanzan los autovectores para construir una base de 3 (suponiendo que estamos trabajando en este espacio) y por lo tanto A NO ES diagonalizable. ¿Cómo reducir esta matriz a una forma simple entonces si no la podemos hacer diagonal? Precisamente de esto se trata la forma de Jordan, buscaremos un vector más, linealmente independiente respecto de los anteriores de modo que podamos construir la matriz de pasaje tal que A quede triangular, en vez de diagonal, expresada como bloques de Jordan.

Sea B esta base, debe estar conformada por tres vectores y solo tenemos dos. Hay varias maneras de encontarlo, una es proponer B={(1,1,0),(1,0,1),(a,b,c)} y buscar las coordenadas (a,b,c) tal que se cumpla J=P1APPJ=AP donde J es la matriz en forma canónica de Jordan y P es la matriz cambio de base de B a la base canónica de 3. Por la manera en la que definimos la base B la matriz J tiene que ser

J=(100021002)

es decir que basta efectuar los productos mencionados e igualarlos, queda

PJ=AP(122a+1102b022c+1)=(12a2b+c10b02a+b+3c)

Queda formado entonces el siguiente sistema

{a+2bc=1b=0abc=1{a=c1b=0,(a,b,c)=(c1,0,c)

que son infinitos vectores de la forma (a,b,c)=(1,0,1)t+(1,0,0) (notemos que se trata de una múltiplo del autovector asociado al autovalor 2 pero con una coordenada sumada). En particular para t=1 se obtiene una de las infinitas soluciones P=(110100011).

En general, cualquier matriz de la forma P=(11t110001t) cumple que

P1AP=J=(100021002).

Nota: podemos resolver el sistema PJ=AP triangulando la matriz ampliada (A2I|101). Este algoritmo será analizado en el siguiente ejemplo con más detalle.

2) Tomemos ahora una matriz similar a la anterior

A=(201110113)

de la cual buscamos la forma de Jordan. El polinomio característico es

pA(λ)=det(λIA)=λ36λ2+12λ8=(λ2)3.

Veamos que para pA(λ)=0 obtenemos un único autovalor λ=2, esto significa que A no es diagonalizable ya que la única manera de obtener que rg(A3I)=0 (y por lo tanto dim(E2)=3) es que (A3I) sea la matriz nula, o lo que es lo mismo, que A sea una matriz diagonal cuya única entrada no nula sea 3.

Busquemos entonces los autovectores...

A2I=(001110111)(110001000){xy=0z=0{x=yz=0

Es decir (x,y,z)E2(x,y,z)=y(1,1,0) esto equivale a afirmar que el vector (1,1,0) genera el subespacio E2=ker(A2I).

Para encontrar otro vector linealmente independiente, podemos triangular la matriz A-2I ampliada con las coordenadas del autovector (asociado a este autovalor 2) como columna.

(A2I|110)=(001111011110)(110000110000){xy=0z=1

Si llamamos y=u tenemos vectores de la forma (x,y,z)=(u,u,1)=u(1,1,0)+(0,0,1) ¡que son combinaciones lineales del autovector! Sin embargo todavía nos falta un vector más para construir una base de 3, este se obtiene sustituyendo la solución 𝐯=u(1,1,0)+(0,0,1) en el producto (A2I)[𝐯].

(A2I)(uu1)=(101) y por lo tanto la base buscada es {(1,1,0),(1,0,1),(u,u,1)}. Es importante ponerlos en este orden, de otro modo la matriz en esta base no estará constituida por bloques de Jordan.

Formamos

P=(11u10u011)P1AP=J=(210021002)u

3) Hallar la forma canónica de Jordan de la matriz

A=(1230001200001200002100001)

Hallamos el polinomio característico:

PA(λ)=(1λ)4(2λ)

Sus raíces son λ1=1 y λ2=2 con multiplicidades 4 y 1 respectivamente.
Buscamos los autovectores, comencemos con

λ1=1_

Triangulamos

AI=(0230000200000200001100000)(0100000100000100000100000)x2=x3=x4=x5=0

Por lo tanto el espacio propio asociado a este autovalor es E1=ker(AI)={(t,0,0,0,0):t}. Resulta evidente que la matriz no es diagonalizable, ya que la multiplicidad geométrica es menor que 4, o lo que es equivalente, dim(E1)=1<4.

Por ahora, busquemos el otro autovector.


λ2=2_

A2I=(1230001200001200000100001)(10014001040001200000100000){x1=14x4x2=4x4x3=2x4x5=0

Es decir que, si designamos t=x4 obtenemos E2=ker(A2I)={t(14,4,2,1,0):t}.


Buscamos la base en la cual A tiene la forma de Jordan. Para λ1 tenemos que hallar 3 vectores más que sean linealmente independientes con (1,0,0,0,0), pues la multiplicidad de λ1es 4 y nosotros tenemos un único vector. Una forma de encontrar estos vectores es la siguiente.

Hallar las potencias (Aλ1I)2,(Aλ1I)3,... hasta que la dimensión del último sea la multiplicidad de la raíz (4 en este caso).

C:=(A1I)2=(0046000040000220001100000)

Obtenemos que el rango es 3, luego su nulidad es 2. Resolviendo el sistema CX=0 se obtiene que todas las coordenadas de los vectores de ker(A1I)2 han de valer cero, excepto las dos primeras. Como ker(AλI)ker(AλI)2, sabemos que podemos expandir la base de ker(A1I) para obtener una base de ker(A1I)2. Elegimos entonces el vector (0,1,0,0,0) . Así: ker(A1I)2=(1,0,0,0,0),(0,1,0,0,0).

D:=(A1I)3=(00014600044000220001100000)

El rango de esta matriz es rg(D)=2. Su nulidad es por tanto 3. Resolvemos el sistema DX=0 y observamos que las dos últimas coordenadas han de valer 0. Expandimos la base de ker(A1I)2 para obtener la de ker(A1I)3, por ejemplo con el vector (0,0,1,0,0):

ker(A1I)3=(1,0,0,0,0),(0,1,0,0,0),(0,0,1,0,0)
E:=(A1I)4=(000141400044000220001100000)

En este caso, la nulidad de E es n(E)=4, y como la dimensión de ker(A1I)4 (es decir, la nulidad de E) no puede ser superior a la multiplicidad algebraica del autovalor 1, que es 4, ya hemos llegado a la dimensión máxima. Resolvemos el sistema EX=0 y concluímos que la suma de las últimas dos coordenadas ha de ser nula. Ahora tomamos un vector v4ker(A1Id)4 pero que no pertenezca a ninguno de los anteriores. Por ejemplo, v4=(0,0,0,1,1). Obtenemos así la base de ker(A1Id)4:

ker(A1Id)4=(1,0,0,0,0),(0,1,0,0,0),(0,0,1,0,0),(0,0,0,1,1).

Ahora hay que hallar v3,v2yv1.

v3=(A1Id)v4=(0,0,2,0,0)
v2=(A1Id)2v4=(A1Id)v3=(6,4,0,0,0)
v1=(A1Id)3v4=(A1Id)v2=(8,0,0,0,0).

Como ker(A2Id)=(14,4,2,1,0), ya tenemos los 5 vectores de la nueva base.

La matriz de cambio de base es

P=(86001404004002020001100010)

Para hallar la matriz de Jordan solo hay que hacer las imágenes por A, de los vectores de la base de Jordan, y expresarlos en dicha base:

Av1=v1(1,0,0,0,0)
Av2=v1+v2(1,1,0,0,0)
Av3=v2+v3(0,1,1,0,0)
Av4=v3+v4(0,0,1,1,0)
Av5=2v5(0,0,0,0,2)
J=(1100001100001100001000002)

Se cumple J=P1AP

Véase también

Enlaces externos


Plantilla:Control de autoridades