Exposición Ríos Grupo2 Algoritmo de Putzer

UNIVERSIDAD REGIONAL
AMAZÓNICA IKIAM
TRONCO COMÚN
MATEMÁTICAS III
GRUPO 2
Algoritmo de Putzer para el cálculo de una matriz exponencial 𝒆𝑨𝒕
Ríos Henry
21 DE FEBRERO DE 2019
El conocimiento de las
matemáticas añade vigor a la
mente, la libera del prejuicio,
credulidad y superstición.
-John Arbuthnot.
Objetivos
 Reconocer la importancia del Algoritmo de Putzer para
sistemas lineales con coeficientes constantes.
 Proponer un ejercicio resuelto didácticamente de un
Sistema Lineal de Ecuaciones Diferenciales de Primer
Orden, asociado al Problema de Valores Iniciales.
Cálculo directo de la función
exponencial matricial 𝒆 𝑨𝒕
 Existen algunos métodos para el cálculo de 𝑒 A𝑡 directamente que pueden

ser usados en cualquier tipo de matriz A cuadrada de entradas
constantes.
 Estos métodos hacen uso del teorema de Cayley-Hamilton en algún
momento de su aplicación y no requieren pasar por el cálculo de una
forma canónica de Jordan.
Método directo. Proposición 1
 Sea la matriz A: [2x2] con valores propios λ1 ≠ λ2 y defínanse la matrices
𝐴−λ2 𝐼 𝐴−λ1 𝐼
Q1 = λ , Q2 = λ
1 −λ2 2 −λ1
Entonces
m𝐴 = λ1 Q1 + λ2 Q 2 (1)
Q 1 Q1 = Q 1 2 = Q 1 , Q 2 Q 2 = Q 2 2 = Q 2 , Q 1 Q 2 = Q 2 Q 1 = 0 (2)
𝐴𝑘 = λ1 𝑘 Q1 + λ2 𝑘 Q 2 , ∀ 𝑘 ∈ N (3)
∞
(𝐴𝑡)𝑘
෍ = 𝑒 λ1 𝑡 Q1 + 𝑒 λ2 𝑡 Q 2 (4)
𝑘!
𝑘=0
Proposición 2
 Para una matriz A como en la proposición anterior excepto con valores

propios iguales λ1 = λ2 sea ahora la matriz
𝑄 = 𝐴 − λ0 𝐼 , donde λ0 = λ1 = λ2
entonces
I. La matriz Q es nilpotente de grado 2.
II. 𝑒 A𝑡 = (𝐼 + 𝑡 𝐴 − λ0 𝐼 )
Algoritmo de Putzer
 Sea el problema de valores iniciales
𝑑𝑥
𝑥 = 𝐴𝑥 , 0, 𝑥0 , 0 ≤ 𝑡 < ∞; 𝑥 ∈ ℝ𝑛 , 𝑥 = ∈ ℝ𝑛
𝑑𝑡
donde A: [n × n] es una matriz real constante. Se propone obtener la solución
de este problema en la forma
𝑥 𝑡 = 𝑒 A𝑡
sin calcular la forma canónica de Jordan asociada a la matriz A.
Algoritmo de Putzer
 El polinomio característico de A se define como

𝑝 λ = det λ𝐼 𝑛 − 𝐴 = λ𝑛 + 𝑐𝑛−1 λ𝑛−1 + ⋯ + 𝑐0
 Se define el operador diferencial lineal de coeficientes constante

𝐿 = 𝑝 𝐷 = 𝐷𝑛 + 𝑐𝑛−1 D𝑛−1 + ⋯ + 𝑐1 𝐷 + 𝑐0
de manera que la función 𝑧 𝑡 es la solución del problema de valores iniciales

𝐿𝑧 = 0, 𝑧 0 = 𝑧ሶ 0 = ⋯ 𝑧 𝑛−2 = 0, 𝑧 𝑛−1 = 1
Algoritmo de Putzer
 La obtención de la solución de un problema lineal como este no ofrece

ningún problema. Se definen también el vector
𝑐1 𝑐2 . . . 𝑐𝑛−1 1
𝑧(𝑡) 𝑐2 𝑐3 . . . 1
𝑧(𝑡) . .
𝑍 𝑡 = .. 𝑦 𝑙𝑎 𝑚𝑎𝑡𝑟𝑖𝑧 𝐶 = .. .
. .
𝑧 (𝑛−1) 𝑐𝑛−1 1
1
Preposición 3
 La matriz exponencial
𝑛−1
𝑒 A𝑡 = ෍ 𝑞𝑗 (𝑡) 𝐴𝑗
𝑗=0
resuelve el problema de valores iniciales (20) y donde q0, q1, . . ., qn−1 son los
elementos del vector
𝑞 = 𝐶𝑍
Ejemplo
 El primer ejemplo es uno ya resuelto usando el algoritmo asociado a la

forma canónica de Jordan. Determinar la matriz 𝑒 A𝑡 donde
2 −7 0
𝐴= 3 −8 0
0 0 1
cuyo polinomio característico tiene tres valores propios distintos
p λ = λ3 + 5λ2 − λ − 5; y λ1 = 1, λ2 = −1, λ3 = −5
Ejemplo
 Este algoritmo de Putzer requiere resolver la ecuación diferencial
𝑧 ′′′ + 5𝑧 ′′ − 𝑧 ′ − 5𝑧 = 0, condiciones iniciales 𝑧 0 = 𝑧 ′ 0 = 0, 𝑧 ′′ 0 = 1
 La solución de este problema es inmediata
1 −5𝑡 1 −𝑡 1
𝑧 𝑡 = 𝑒 − 𝑒 + 𝑒𝑡
24 8 12
Ejemplo
 Con esta función y los coeficientes del polinomio característico se construyen

las matrices
1 −5𝑡 1 −𝑡 1
𝑒 − 𝑒 + 𝑒𝑡
𝑧(𝑡) 24 8 12
5 1 1
𝑍 = 𝑧′(𝑡) = − 𝑒 −5𝑡 + 𝑒 −𝑡 + 𝑒 𝑡
𝑧′′(𝑡) 24 8 12
25 −5𝑡 1 −𝑡 1
24 8 12
y
−1 5 1
𝐶= 5 1 0
1 0 0
Ejemplo
 Las funciones q se obtienen de

1 −5𝑡 5 −𝑡 5
− 𝑒 + 𝑒 + 𝑒𝑡
𝑞0 24 8 12
1 1
𝑞 = 𝐶𝑍 = 𝑞1 = − 𝑒 −𝑡 + 𝑒 𝑡
𝑞2 2 2
1 −5𝑡 1 −𝑡 1
24 8 12
 La función exponencial matricial buscada es

2
𝑒 A𝑡 = ෍ 𝑞𝑗 𝐴𝑗
𝑗=0
= 𝑞0 𝐼3 + 𝑞1 𝐴 + 𝑞2 𝐴2
−5𝑡 −𝑡 −5𝑡 −𝑡
1 −3𝑒 + 7𝑒 7(𝑒 − 𝑒 ) 0
= −3(𝑒 −5𝑡 − 𝑒 −𝑡 ) 7𝑒 −5𝑡 − 3𝑒 −𝑡 0
4
1 0 4𝑒 𝑡
y que es la matriz obtenida arriba.
Referencias
 Morones, R. (2012). Cálculo de la función exponencial de una matriz. Dept.

de Matemáticas, ITAM.
 Valverde Sandoval, O. G. (2005). El método de EJ Putzer en la solución de
sistemas lineales con corrientes constantes.

Exposición Ríos Grupo2 Algoritmo de Putzer

Cargado por

Copyright:

Exposición Ríos Grupo2 Algoritmo de Putzer

Cargado por

Información del documento

Descripción original:

Derechos de autor

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Exposición Ríos Grupo2 Algoritmo de Putzer

Cargado por

Copyright:

UNIVERSIDAD REGIONAL

 Existen algunos métodos para el cálculo de 𝑒 A𝑡 directamente que pueden

 Para una matriz A como en la proposición anterior excepto con valores

 Sea el problema de valores iniciales

 El polinomio característico de A se define como

 Se define el operador diferencial lineal de coeficientes constante

de manera que la función 𝑧 𝑡 es la solución del problema de valores iniciales

 La obtención de la solución de un problema lineal como este no ofrece

 El primer ejemplo es uno ya resuelto usando el algoritmo asociado a la

cuyo polinomio característico tiene tres valores propios distintos

𝑧 ′′′ + 5𝑧 ′′ − 𝑧 ′ − 5𝑧 = 0, condiciones iniciales 𝑧 0 = 𝑧 ′ 0 = 0, 𝑧 ′′ 0 = 1

 La solución de este problema es inmediata

 Con esta función y los coeficientes del polinomio característico se construyen

 Las funciones q se obtienen de

 La función exponencial matricial buscada es

 Morones, R. (2012). Cálculo de la función exponencial de una matriz. Dept.

También podría gustarte