Matemática Aplicada I

Leonardo Flores Gonzalez - Rolando Astete Chuquichaico

Universidad Nacional de Ingenierı́a

Sesión 3
 Reducción de sistemas de orden superior a sistemas de primer orden

Podemos restringirnos a los sistemas de primer orden. Esto es debido a que cualquier
sistema de orden superior es equivalente a uno de primer orden.
Dado el sistema

x(n) = f (t, x, x0 , x00 , . . . , x(n−1) ) (∗)

introducimos nuevas variables dependientes para x y cada una de sus derivadas hasta
de orden (n − 1) :

x1 = x, x2 = x0 , x3 = x00 , · · · , xn = x(n−1)

de esta forma x01 = x0 = x2 , x02 = x00 = x3 , x03 = x000 , · · · , etc. El sistema

resultante es
x01 = x2
x02 = x3
..
.
x0n−1 = xn
x0n = f (t, x, x1 , x2 , . . . , xn )
es un sistema de dimensión n pero de primer orden.

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 2 / 25

 Reducción de sistemas de orden superior a sistemas de primer orden

En notación matricial consideremos x = (x1 , x2 , . . . , xn ) y F = (F1 , F2 , . . . , Fn ) donde

F1 (t, x) = x2 , F2 (t, x) = x3 , . . . , Fn−1 (t, x) = xn y Fn (t, x) = f (t, x). De esta
manera el sistema resultante se puede escribir como

x0 = F(t, x) (∗∗)

Si x = x(t)
 es solución de la ecuación diferencial (∗), entonces
x(t) = x(t), x0 (t), x00 (t) . . . , x(n−1) (t) es solución del sistema.(∗∗). Y
recı́procamente si x(t) = (x1 (t), x2 (t), . . . , xn (t)) es solución del sistema (∗∗),
entonces x = (x(t) es solución de la ecuación (∗)
Para el problema de condiciones iniciales
x0 = F(t, x)

donde
x(t0 ) = x0
x(t0 ) = x0 = (x01 , x02 , . . . , x0n )
y
x = f (t, x, x0 , x00 , . . . , x(n−1) )
 (n)

x(t0 ) = x01 , x0 (t0 ) = x02 , . . . , x(n−1) (t0 ) = x0n

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 3 / 25

 Reducción de sistemas de orden superior a sistemas de primer orden

Example
Expresar como un sistema de primer orden:
1 y 000 + 3y 00 − 4y = 0,
solución
y 000 = 4y − 3y 00 sean: x1 = y, x2 = y 0 , x3 = y 00
0 0 0
x
1 0 = y = x 2 , x2 = y 00 = x3 , x03 = y 000
 x1 = x2
x0 = x3
 20
x3 = 4x1 − 3x3
2 y 000 + 3y 00 + 6y 0 + 3y = x
solución
y 000 = x − 3y − 6y 0 − 3y 00 sean: x1 = y, x2 = y 0 , x3 = y 00
0 0 0 00 0 000
x
1 =0 y = x2 , x2 = y = x3 , x3 = y
 x1 = x2
x0 = x3
 20
x3 = x − x1 − 6x2 − 3x3
3 x000 + xx00 = cos t
solución
x000 = cos t − xx00 sean: x1 = x, x2 = x0 , x3 = x00
x1 = x = x2 , x2 = x00 = x3 , x03 = x000
0 0 0

x0 G.

Leonardo Flores = Rolando
x Astete Ch. Matemática Aplicada I Mayo, 2020 4 / 25
 Proposición

Las soluciones del sistema lineal X 0 = AX forman un espacio vectorial

La función X(t) = n
P
k=1 ck Xk (t) se dice que es una combinación lineal de las
funciones Xk (t) : k = 1, 2, · · · , n
0
Sabemos que si XP k (t) : k = 1, 2, · · · , n es un conjunto de soluciones LI de X = AX
n
entonces X(t) = k=1 ck Xk (t) también es solución. (Principio de superposición de
soluciones)
Definition
Se llama Sistema fundamental de soluciones del sistema X 0 = AX en I a cualquier
conjunto de n soluciones L.I. en I.

Theorem
Siempre existe un sistema fundamental de souciones del sistema X 0 = AX. en I.

Demostración.
Sea t0 ∈< a, b > por el teorema de existencia y unicidad de soluciones se sabe que
existen Xk (t), k = 1, 2, · · · , n, funciones que son soluciones del sistema y cada una
verifica X1 (t0 ) = (1, 0, · · · , 0) , X2 (t0 ) = (0, 1, 0, · · · , 0) , · · · , Xn (t0 ) = (0, 0, · · · , 1)
Las funciones Xk (t), k = 1, 2, · · · , n, son LI, ya que si existen n constantes
ck : k Flores
Leonardo = 1, 2,G. ·Rolando
· · , n, Astete
tal queCh. ck Xk (t) = 0, ∀t ∈<
Matemática a, b I>, se tiene que
Aplicada Mayo, 2020 5 / 25
 Sistema de ecuaciones diferenciales de primer orden

Theorem

Si {X1 , X2 , · · · , Xn } es un sistema fundamental de soluciones del sistema X 0 = AX

en I, entones toda solución ϕ(t) del sistema se puede expresar de la forma
ϕ(t) = n
P
k=1 k Xk (t) donde ck : k = 1, 2, · · · , n son constantes
c

Definition

Se denomina matriz fundamental del sistema lineal homogéneo X 0 (t) = AX(t) a una
matriz Φ(t) cuyas columnas son n soluciones linealmente independientes del sistema,
(son un sistema fundamental de soluciones en I)
···
 
x11 x12 x1n
 x21 x22 ··· x2n 
Φ(t) = (X1 (t), X2 , (t), . . . , Xn (t)) = 
 .. .
..
. 
.. 
donde los vectores
.
xn1 xn2 ··· xnn
solución son
     
x11 x12 x1n
 x21   x22   x2n 
X1 =  .  ; X2 =  .  ; · · · Xn =  . 
     
 ..   ..   .. 
xn1
Leonardo Flores xn2
G. Rolando Astete Ch. xnn I
Matemática Aplicada Mayo, 2020 6 / 25
 Sistema de ecuaciones diferenciales de primer orden

Example
e2t te2t e2t + te2t te2t
   
Comprobar que ϕ1 (t) = ϕ2 (t) =
−e2t −te + e2t
2t
−te2t e − te2t
2t

x01
    
3 1 x1
son matrices fundamentales del sistema =
x02 −1 1 x2

Example
2e2t 3et
 
Dada la matriz ϕ(t) = Determinar un sistema de ecuaciones
−e2t 4et
diferenciales para que ϕ(t) sea una matriz fundamental.

Example
4e2t −e−3t
 
Comprobar que ϕ(t) = es una matriz fundamental del sistema
e2t e−3t
 0    
x1 1 4 x1
=
x02 1 −2 x2

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 7 / 25

 Sistema de ecuaciones diferenciales de primer orden

Observación
Se denomina matriz fundamental principal a la matriz Φ tal que para un valor t0 en
el intervalo I verifica que Φ(t0 ) = In×n

Observación
La solución general de un sistema X 0 = AX se puede expresar en términos de la
matriz fundamental del sistema X(t) = Φ(t)C donde C es un vector de constantes
indeterminadas . Una matriz queda caracterizada por verificar:
1 Φ0 (t) = A(t)Φ(t) pues es una matriz solución, y
2 su determinante es distinto de cero

Propiedades de la matriz fundamental

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 8 / 25

 Ejemplo

Example
x01
! ! !
2 6 6 x1
0
Resolver el sistema X = x02 = −3 −7 −6 x2
x03 3 6 5 x3

Solución  
 2−λ 6 6 
|A − λI|
 −3 −7 − λ −6
=

 
 3 6 5−λ 
= −λ3 + 3λ + 2 = − (λ − 2) (λ + 1)2 = 0

λ1 = 2 ;m = 1
los valores propios son , para determinar los vectores
λ2 = −1 ; m = 2
propios resolvemos el sistema (A −  λI) v = 0 para cada valor propio, para:
1
λ1 = 2 ⇒ (A − 2I) v = 0 ⇒ v1 =  −1 
1
   
−2 −2
λ2 = −1 ⇒ (A + I) v = 0 ⇒ v2 =  1  , y v3 =  0 
0 1
−2e−t −2e−t
 2t     
e
las soluciones son X1 =  −e2t  , X2 =  e−t  , y X3 =  0 
2t −t
e 0 e
Leonardo Flores G. Rolando Astete Ch.  Matemática
 Aplicada
 I   Mayo, 2020 9 / 25
 Sistema de ecuaciones diferenciales de primer orden

Propiedades
1 Formula de Jacobi. Si ϕ(t) es una matriz solución entonces:
(det ψ)0 (t) = traza (A(t)) . det ψ(t)
R t
traza(A(s))ds
det(ψ(t)) = det ψ(t0 )e t0

2 Sea Φ(t) es una matriz fundamental si y solo si existe algún t0 ∈< a, b > tal que
det (Φ(t0 )) es distinto de cero.
3 Si Φ(t) es una matriz fundamental entonces cualquier otra matriz fundamental
es de la forma Φ · P donde P es una matriz regular de orden n. ( P.Φ(t) puede
no ser matriz fundamental)
4 Si Φ(t) es una matriz fundamental se verifica que:
Z
X(t) = XG (t) + XP (t) = Φ(t)C + Φ(t) Φ−1 (t)b(t)dt

(Fórmula de variación de constantes)

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 10 / 25

 Sistema de ecuaciones diferenciales no homogéneos de primer orden

Propiedades
1 Si
 Φ 0es una matriz fundamental, y si se considera el problema de valor inicial:
X (t) = A(t)X(t) + b(t)
donde (t0 , X0 ) ∈ I × Rn entonces la solución única
X(t0 ) = X0
es: Z t
X(t) = Φ(t) · Φ−1 (t0 )X0 + Φ(t) Φ−1 (t)b(t)dt
t0

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 11 / 25

 Sistema de ecuaciones diferenciales no homogéneos de primer orden

Example
    
1 1 x 2t
X 0 (t) = + ,
4 1 y 2t + 1
Solución  
 1−λ 1  2
|A − λI| =  4 1−λ 
 λ − 2λ − 3 = (λ − 3) (λ + 1) = 0

λ1 = −1 ; m = 1
los valores propios son , para determinar los vectores
λ2 = 3 ;m = 1
propios resolvemos el sistema (A − λI) v= 0 para cada valor propio, para:
− 21
λ1 = − ⇒ (A + I) v = 0 ⇒ v1 =
1
 1 
λ2 = 3 ⇒ (A − 3I) v = 0 ⇒ v2 = 2 ,
1
 1 −t   1 3t 
−2e 2
e
las soluciones son X1 = −t , X2 = ,
e e3t
 1 −t   1 3t 
−2e 2
e
solución general de la gomogénea X(t) = c1 −t + c2
e e3t
 1 −t 1 3t 
−2e 2
e
matriz fundamental Φ(t) =
e−t e3t
R −1
para la solución particular
Leonardo Flores G. Rolando Astete Ch.
XPMatemática
(t) = Φ(t) Φ (t)F (t)dt
Aplicada I Mayo, 2020 12 / 25
 Ejemplos

Example
e2t te2t e2t + te2t te2t
   
Comprobar que ϕ1 (t) = ϕ2 (t) =
−e2t −te + e2t
2t
−te2t e − te2t
2t

x01
    
3 1 x1
son soluciones fundamentales del sistema =
x02 −1 1 x2
 
3 1
Solución: la matriz A = para compromar que ϕ1 (t) es solución del
−1 1
sistema 
2e2t e2t + 2te2t

ϕ01 (t) = 2t
−2e −2te2t + e2t
  2t
te2t 2e2t 2te2t + e2t
   
3 1 e
2t 2t 2t =
−1 1 −e −te + e −2e2t −2te2t + e2t
Similar para ϕ2 (t).

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 13 / 25

 Ejemplos

Example
2e2t 3et
 
Dada la matriz ϕ(t) = 2t t Determine un sistema de ecuaciones
e 2e
diferenciales para que ϕ(t) sea una matriz fundamental.

2e2t 3et
   
Solución: Sea X1 (t) = 2t , X2 (t) = t son L.I. y son soluciones del
e 2e
 0
x1 = 2x1
sistema buscada por lo tanto ϕ(t) es matriz fundamental del sistema
x0 = x
 0   2 2 
x1 2 0 x1
=
x02 0 1 x2

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 14 / 25

 Ejemplos

Example
Demostrar que z1 = (a + bi) e(α+βi)t , z1 = (a − bi) e(α+βi)t son soluciones complejas
x1 = Real (a + bi) e(α+βi)t
LI de un sistema lineal de ecuaciones diferenciales son
x2 = Im (a + bi) e(α+βi)t
soluciones reales y son LI

Example
4e2t −e−3t
 
Comprobar que ϕ(t) = es una matriz fundamental del sistema
e2t e−3t
 0    
x1 1 4 x1
=
x02 1 −2 x2

Example
e3t (1 − 2t)e−3t
 
Dada la matriz ϕ(t) = encontrar un sistema de ecuaciones
e3t −2te−3t
diferenciales para el que ϕ(t) es una matriz.

Example
2t
 
0
Leonardo Flores G. Rolando Astete Ch. −2eMatemática
0 Aplicada I Mayo, 2020 15 / 25
 Ejemplos

Example
et cos t et sent
   
1 Verificar que X1 = t , X2 = t son
e (cos t − sent) e (cos t + sent)
soluciones del sistema
2 Verificar que X(t) = cX1 (t) + c2 X2 (t) también es solución del sistema
 
2
3 Halle la solución que verifica la condición X(0) =
3

Example
x01 = 4x1 + 2x2

Resolver el sistema
x02 = 3x1 − x2

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 16 / 25

 Autovalores y autovectores de una matriz

Ecuación caracterı́stica, Autovalores y autevectores

sea el siistema X 0 (t) = AX)t) se supone que el sistema tiene una solución de la forma
Para una matriz A , su polinomio caracterı́stico, en el plano complejo viene dado
por:
det(A − λI) = − (λ − λ1 )m1 (λ − λ2 )m2 · · · (λ − λk )mk
donde: λ1 , λ2 , . . . , λk son raices diferentes. mi :se llama multiplicidad
algebraica de λi
Se llama subespacio propio de λi , al espacio Si formado por los vectores v tales
que
(A − λi I) v = 0
Al número ri , la dimensión del espacio propio de λi , se le dice multiplicidad
geométrica de λi .

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 17 / 25

 Autovalores y autovectores de una matriz

Theorem
Para cada valor propio λi se cumple

1 ≤ ri ≤ mi

Theorem
La matriz A es diagonalizable si y solo si ri = mi

Corolario
(Corolario)Si mi = 1 para todo valor propio λi (esto es, si hay exactamente n raı́ces
diferentes del polinomio caracterı́stico), entonces A es diagonalizable .

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 18 / 25

 Autovalores y autovectores de una matriz

Estudiaremos el caso en el cual la matriz no es dicgonalizable, que por teorema

será para ri < mi , esto es, cuando se obtiene menos vectores psropios que la
multiplicidad de λi . En este caso, si deseamos calcular la forma de Jordan,
necesitamos obtener de alguna manera mi − ri vectores. A estos vectores que se
necesitan se denominan vectores propios generalizados .

Dado un vector propio v, el vector propio generalizado v1 se obtiene resolviendo

el sistema
(A − λi I) v1 = v
Si paraa completar el espacio necesitamos otro vector generalizado v2 , es posible
obtenerlo resolviendo el sistema

(A − λi I) v2 = v1

y ası́ sucesivamente.
Los vectores v1 , v2 , . . . , vs ası́ obtenidos forman una cadena, que llamaremos la
cadena de Jordan de largo s asociada al vector propio v
Usaremos en lo siguiente la notación v ,→ v1 ,→ v2 ,→ . . . ,→ vs para denotar a la
cadena de Jordan obtenida a partir de v

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 19 / 25

 Sistema de ecuaciones diferenciales no homogéneos de primer orden

Example
Supongamos que se tiene una matriz A tiene un valor propio λ de multiplicidad
algebraica 5 y multiplicidad geométrica 3 (esto es, setiene 3 vectores propios u, v, w
asociados a λ). Luego tenemos que obtener los 2 vectores restantes (vectores propios
generalizados ). Dichos
  vectores pueden
 ser obtenidos con sólo uno de estos esquemas:
 u ,→ u1  u ,→ u1  u
v ,→ v1 v v ,→ v1
w w ,→ w1 w ,→ w1
  
  
 u ,→ u1 ,→ u2  u  u
v v ,→ v1 ,→ v2 v
w w w ,→ w1 ,→ w2
  

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 20 / 25

 Formas canónicas de Jordan

Theorem
El número ks de la cadena de Jordan de largo s es ks = 2ls − ls−1 − ls+1 donde
l0 = 0
ls = dim ker (A − λI)s , para s > 0
Además, si ls+1 = ls , entonces lp = ls , para todo p ≥ s

Cuando se tiene una matriz que es diagonalizable, en su descomposición A = P DP −1 ,

colocamos sus valores propios en la diagonal de la matriz D, y en la matriz P se van
colocando como columnas los vectores propios correspondientes a cada valor propio.

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 21 / 25

 Sistema de ecuaciones diferenciales no homogéneos de primer orden
Para una matriz que no es necesariamente diagonalizable, podemos encontrar su
descomposición en la forma canónica de una manera similar.
Construcción de dicha matriz
1 En la matriz J, se colocan los n valores propios en la diagonal. Ası́, si el

polinomio caracterı́stico es de la matriz es

det(A − λI) = − (λ − λ1 )m1 (λ − λ2 )m2 · · · (λ − λk )mk
entonces se debe colocar m1 veces el valor λ1 , m2 veces el valor λ2 , y ası́
sucesivamente. Los mi valores propios λi , , forman un supra bloque de Jordan.
2 Si para cada valor λi , existen l cadenas, entonces el supra bloque de Jordan de
valor λi estará compuesto por l bloques, y las dimensiones de cada bloque
corresponderá al largo de cada cadena que le corresponde.
3 Finalmente, la matriz P , se van colocando los vectores propios, y si hay cadenas,
se colocan los elementos de la cadena de la manera como fueron obtenidos.
Example
 
8 0 0 0 0 0

 0 8 0 0 0 0 

 0 0 0 0 0 0 
Consideremos la matriz A =  

 0 0 0 0 1 0 

 0 0 0 0 0 0 
0 0 0 0 0 1
Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 22 / 25
 Formas canónicas de Jordan

Example
 
−1 0 −2
Calcule la forma canónica de Jordan de A =  0 −1 −3 ,
2 −1 −3
 1 
2
|A − λI| = −(λ + 2)2 (λ + 1) = 0 para λ = −1 v1 =  1  para
0
 
2
λ2 = −2, m2 = 2 v2 =  3 ,
1

Example
 
3 −4 4
Calcule la forma canónica de Jordan de A =  2 −3 4 
1 −2 3

Example
 
tA 1 4
Calcúlese la exponencial de e donde A =
−1 −3
Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 23 / 25
 Sistema de ecuaciones diferenciales de primer orden

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 24 / 25

 Bibliografia

Ivorra Castillo. () Funciones de variable Compleja

Nagle,R. Kent, Saff, Edward B. y Snider, Arthur D.

a
Ecuaciones Diferenciales y Problemas con Valores en la Frontera.4 edición
México. Pearson-Addison-Wesley, 2005.
Murray R. Spiegel. Ecuaciones diferenciales aplicadas. Prentice Hall
Hispanoamericana1997.
Dennis Zill. Ecuaciones diferenciales con aplicaciones. Grupo Editorial
Iberoamericano 1992
Bertha Zelaya Análisis Numérico Universidad del Paı́s Vasco

Ruel V. Churchill Variable Compleja Mc Graw Hill

Leonardo Flores G. Rolando Astete Ch. Matemática Aplicada I Mayo, 2020 25 / 25