Exp Matriz
Exp Matriz
Fijado un n
umero real x, podemos escribir
x
e =
+ n
X
x
n=0
n!
=1+x+
x2 x3
+
+
2
3!
e := In +
+
X
Ak
k=1
k!
1 0 0
A := 0 2 0 ,
0 0 3
se cumple que
e 1 0
0
eA = 0 e2 0 .
0
0 e3
Cuando la matriz no es diagonal, el calculo de la matriz exponencial se hace mas laborioso. En estas notas se comentan dos metodos que permiten el calculo de la exponencial
de una matriz en el caso general. El primero es haciendo uso de la matriz fundamental de
un sistema de ecuaciones lineal de orden 1.
Definici
on 3. Se llama matriz fundamental del sistema
dY
= AY,
dx
(1)
a cualquier matriz (x) en la cual si denotamos su columna i-esima por la matriz columna Yi (x), entonces {Yi (x)}1in es un sistema fundamental de soluciones de (1). Las
matrices fundamentales suelen representarse haciendo referencia al sistema fundamental
que las forma, e. d.,
(x) = [Y1 (x) Yn (x)] .
dY
= AY, siendo
dx
4 1
0
3 1 .
A= 1
1
2
0
Teorema 5. Fijada una matriz A de orden n, entonces eAx = (x)(0)1 , donde (x)
es una matriz fundamental de (1).
Ejercicio 6. Obtener eAx siendo1
3 4 5
A := 0 5 4 .
0 0 3
Comentamos ahora el segundo procedimiento para hallar la exponencial de una matriz.
En este caso haremos uso de la diagonalizacion o de la descomposicion de Jordan de la
matriz correspondiente.
Teorema 7. Sea A una matriz de orden 3 y 1 , 2 , 3 sus valores propios (no necesariamente distintos). Si los valores propios son no defectuosos, entonces A es diagonalizable
y puede expresarse
A = BDB1 ,
donde
1 0 0
D = 0 2 0 ,
0 0 3
y B es la matriz cuya columna i-esima esta formada por las coordenadas de un vector
propio asociado al valor propio i . Ademas,
eAx = BeDx B1 .
Ejercicio 8. Obtener eAx usando el resultado anterior, siendo
9
4 0
A := 6 1 0 .
6
4 3
Teorema 9. Sea A una matriz de orden 3 y su u
nico valor propio con defecto 2.
Entonces A no es diagonalizable aunque puede expresarse
A = BJB1 ,
donde
1 0
J = 0 1 ,
0 0
y B es la matriz cuya primera columna esta formada por las coordenadas de un vector
propio asociado al valor propio y el resto de columnas por las coordenadas de dos vectores
propios generalizados en el orden conveniente. Ademas,
eAx = BeJx B1 ,
1
donde
eJx
x2
1 x
= ex 0 1 2!
x .
0 0 1
0
1
2
A := 5 3 7 .
1
0
0
Teorema 11. Sea A una matriz de orden 3, 1 un valor propio simple y 2 un valor
propio con defecto 1. Entonces A no es diagonalizable aunque puede expresarse
A = BJB1 ,
donde3
1 0 0
J = 0 2 1 ,
0 0 2
y B es la matriz cuya primera columna esta formada por las coordenadas de un vector
propio asociado al valor propio 1 , la segunda columna esta formada por las coordenadas
de un vector propio asociado al valor propio 2 y la tercera por las coordenadas de un
vector propio generalizado correspondiente. Ademas,
eAx = BeJx B1 ,
donde
eJx
e 1
0
0
= 0 e2 x e2 x x .
0
0
e2 x
3 4
A := 0 5
0 0
anterior, siendo
5
4 .
3
J1 0
0 J2
J=
0 0
2
...
...
0
0
,
Js
1
3
la submatriz 2
se denomina bloque de Jordan de orden 2
0 2
e 1
0 ...
0
0 eJ2 x . . .
0
,
eJx =
Js x
0
0 e
donde cada
eJr x
x2
x
2!
j
=e 0 1
0 0
xp1
(p 1)!
xp2
,
(p 2)!
1 1 1
3 1
(i) A =
; (ii) A = 2 2 2 .
1
1
3 3 3
En el siguiente teorema, se observa como la matriz exponencial aparece al resolver
sistemas lineales de orden 1.
Teorema 14. Fijada una matriz A de orden n e Y0 una matriz columna n1, el problema
d Y = AY;
dx
Y(0) = Y0 ;
tiene por solucion Y = eAx Y0 .
Ejercicio 15. Resolver, haciendo uso de la matriz exponencial, el sistema
4 1
0
3 1 Y.
Y 0 = 1
1
0
2
El uso de la matriz exponencial es muy extendido a la hora de resolver sistemas lineales
completos. No obstante eso sera tratado en otras notas.
Si Jr es un n
umero real, este puede interpretarse como un bloque de Jordan de orden 1