Método Matricial en Ecuaciones Diferenciales
Método Matricial en Ecuaciones Diferenciales
ANTÚNEZ DE MAYOLO
FACULTAD DE INGENIERÍA DE MINAS,
GEOLOGÍA Y METALURGIA
DOCENTE:
ROSALES APEÑA, Miguel Dionicio
ESTUDIANTES:
- GONZALES SANCHES, Joaquín
2.1 Matriz.......................................................................................................................5
CONCLUSIONES...................................................................................................23
REFERENCIAS BIBLIOGRAFICAS..................................................................24
INTRODUCION
Este trabajo tiene como objetivo explorar los fundamentos teóricos y las
aplicaciones prácticas de los sistemas de ecuaciones diferenciales que involucran
valores propios. Comenzaremos con una revisión de los conceptos básicos de valores y
vectores propios, seguido de una discusión sobre cómo estos conceptos se integran en la
solución de sistemas de ecuaciones diferenciales. Posteriormente, se presentarán
métodos analíticos y numéricos para resolver estos sistemas, acompañados de ejemplos
que ilustran su aplicación
2.1 Matriz
2.2. Sistema de ecuaciones lineales
, ,
,
Siempre y cuando el determinante de la matriz de coeficientes
Uso de las determinantes para resolver sistemas
En un sistema de n ecuaciones lineales con n incógnitas
Donde ,
Ejemplo:
Utilice la regla de Cramer para resolver el sistema
Solución: La solución requiere que se evalúen los cuatro determinantes
Como
Como el sistema tiene solución única, podemos resolverlo ya sea por la regla de Cramer
2.3 Concepto de valores propios (autovalores) y vectores propios (autovectores)
Debemos recordar las definiciones de valores y vectores propios y algunas propiedades
veamos cómo es que estos conceptos son útiles para resolver sistemas lineales de primer
orden homogéneos.
Cálculo de autovectores:
Debido a que un sistema lineal puede ser visto como una ecuación matricial los
resultados de álgebra lineal sobre valores y vectores propios de matrices pueden ser
aplicados aquí. En esta entrada daremos un breve repaso sobre estos conceptos y
veremos cómo es que estos resultados nos pueden ayudar a determinar la solución
general de algunos sistemas de ecuaciones diferenciales.
det [ A−Iλ ] =0
Donde I representa a la matriz de unidad y debe tener la misma dimensión de orden.
Matriz de Unidad: Es una matriz diagonal en la que los elementos de la diagonal
principal son iguales a 1.
Ejemplo: En una matriz 3x3
( )
1 0 0
I= 0 1 0
0 0 1
λ 1=2
λ 2=2
λ 3=0
Por lo tanto:
λ 1=2→ m1=2
λ 3=0 → mi =1
A) Polinomio Característico:
det [ A−Iλ ] ⃗v 1=⃗0
Sea
()
x
⃗v 1= y
z
() ()
x 0
det [ A−Iλ ] y = 0
z 0
( )
2 0 0
1 1 1
−1 1 1
( ) ( )
2 0 0 1 0 0
1 1 1 −λ 0 1 0 =0
−1 1 1 0 0 1
( )
2−λ 0 0
1 1−λ 1 =0
−1 1 1−λ
Resolvemos la matriz3 x 3
( 2− λ ) ( 1−λ )( 1− λ ) +0+ 0−0−2+ λ+ 0=0
( 2− λ ) (−2 λ+ λ 2 )=0
( 2− λ ) (−2+ λ ) λ=0
λ 1=2; λ2=2 y λ3=0
Hallamos los autovectores base para cada valor propio:
o Para λ 1=2
( )( ) ( )
2−2 0 0 x 0
1 1−2 1 y =0
−1 1 1−2 z 0
( )( ) ( )
0 0 0 x 0
1 −1 1 y = 0
−1 1 −1 z 0
{−1x −x +yy−
+ z=0
z=0
Se obtiene que x= y −z, damos el valor de y=β y z=α ∈ R
Entonces la base de los vectores propios λ 1=2 es
() ( ) ( ) ()
x β−α −1 1
⃗v 1= y = α =β 0 +α 1
z β 1 0
{( ) ( )}
−1 1
H ( λ1 ) = 0 , 1
1 0
o Para λ 3=0
( )( ) ( )
2 0 0 x 0
1 1 1 y =0
−1 1 1 z 0
{x+2y+x=0z =0
Se obtiene que x=0 , y=−z , damos el valor de z=α ∈ R
Entonces la base de los vectores propios λ 1=2 es
()( ) ( )
x 0 0
⃗v 2= y = −α =α −1
z α 1
()
0
H ( λ2 ) = −1
1
( )
a11 a 12 a
⋯ 1n
a21 a 22 a 2n
A=
⋮ ⋱ ⋮
an 1 a 2 n ⋯ ann
'
Y = AY
( )( )( )
'
y 1 (t) a11 a12 a y (t)
' ⋯ 1n 1
y 2 (t) = a2 1 a22 a2 n y 2 (t)
⋮ ⋮ ⋱ ⋮ ⋮
'
y n (t) a n 1 a2 n ⋯ ann y n (t)
λ 1t λ2t λnt
Y ( t )=C1 ⃗v 1 e +C2 ⃗v 2 e +…+C n ⃗v n e
Teorema: Sea λ=α +iβ un valor propio complejo de la matriz de coeficientes A del
sistema lineal homogéneo y sean U y V los vectores que definen al vector propio
complejo ⃗v asociado a 𝜆 , entonces
αt
y 1 ( t )=e ¿
αt
y 2 ( t ) =e ¿
son soluciones linealmente independientes de Y ' = AY en (−∞ , ∞ ).
Ejemplo : Dar la solución general de la siguiente ecuación diferencial homogénea
( )( )( )
x' 1 0 0 x
y' = 0 1 −1 y
z' 0 1 1 z
Solución:
Dado la forma de Y ' = AY , trabajamos con la matriz A
( )
1 0 0
A= 0 1 −1
0 1 1
Hallamos los autovalores de la matriz A
det [ A−Iλ ] =0
( )( )
1 0 0 1 0 0
[ A−Iλ ] = 0 1 −1 − 0 1 0 λ=0
0 1 1 0 0 1
( )
1−λ 0 0
0 1−λ −1 =0
0 1 1−λ
( )( ) ( )
0 0 0 a 0
0 0 −1 b = 0
0 1 0 c 0
Las filas son múltiplos
{−c=0
b=0
El valor de “a” al ser libre lo elegimos como a=1 entonces el primer vector propio
es
() ()
a 1
⃗v 1= b = 0
c 0
o Para λ 2=1+i
Buscamos el vector propio ⃗v 1 ≠ 0
det [ A−Iλ ] ⃗v 1=⃗0
( )( ) ( )
1−1−i 0 0 a 0
0 1−1−i −1 b=0
0 1 1−1−i c 0
( )( ) ( )
−i 0 0 a 0
0 −i −1 b = 0
0 1 −i c 0
Las filas son múltiplos
{
−ia=0
−bi−1 c=0 → bi=1 c
−ic=0
El valor de “a” es 0, C al ser libre lo elegimos como c=1 entonces, b=i.Por lo
tanto, el vector propio es
()
0
⃗v 1= i
1
Al ser complejo, también lo podemos expresar como
() () ()
0 0 0
⃗v 1= i = 0 +i 1
1 1 0
()
0
⃗v 1= i =U +Vi
1
En este resultado se obtiene dos soluciones independientes y forman dos soluciones
del sistema de ecuaciones diferencial, una real y la otra imaginaria
( ) () )
0 0
y 2 ( t ) =et 0 cos (t )− 1 sen(t )
1 0
( )( ) ( )
0 0 0
t t
y 2 ( t ) =e 0 − sen (t) =e −sen(t )
cos (t) 0 cos (t)
(( ) ( ) )
0 0
t
y 3 ( t ) =e 0 sen (t)+ 1 cos (t)
1 0
( )
0
t
y 3 ( t ) =e cos (t )
sen(t )
λ1 t
Y ( t )=C1 y 1 ( t ) +C 2 y 1 ( t ) +C 3 ⃗v 3 e
( ) ( ) ()
0 0 1
t t t
Y ( t )=C1 e −sen (t) +C 2 e cos (t) +C 3 e 0
cos(t ) sen(t ) 0
Construimos la matriz
( 2− λ ) ( 1−λ )−12=0
2
2−2 λ−λ+ λ −12=0
2
10−3 λ+ λ =0
λ 1=5 y λ 2=−2
Luego de hallar los valores propios, calculamos los vectores propios base de cada valor
propio
o Para λ 1=5
(−34 3 x
−4 y
=)( ) ( )
0
0
{−34 x−4
x+3 y =0
y =0
Se obtiene que x= y , damos el valor de y=α ∈ R
Entonces la base de los vectores propios λ 1=5 es
(11)
H ( λ1 ) =
o Para λ 1=−2
( )
−3
⃗v 2=H ( λ2 )= 4
1
λ 1t λ2t
Y ( t )=C1 ⃗v 1 e +C2 ⃗v 2 e
( )
−3
()
1 5t
Y ( t )=C1 e +C2 4 e
1
1
−2 t
{ x 1 '=3 x1−13 x 2
x 2 ' =5 x1 + x 2
x 1 ( 0 )=3 ; x 2 ( 0 )=−10
Solución:
( )(
x1 '
x2 ' 5 1 x2
x
= 3 −13 1 )( )
La matriz tiene la forma de
'
X = AX
2
λ −4 λ+68=0
4 ± √ 16−272
λ= =2 ± 8 i
2
Como los valores propios son complejos, encontramos su auto vector considerando
2+8 i
(3−2−8
5
i −13 x
1−2−8 i y
=)( ) ( )
0
0
(1−85 i −13 x
−1−8 i y)( ) ( )
=
0
0
Ambas filas son múltiplos, usamos la segunda fila
5 x−( 1+8 i ) y =0
5 x=( 1+8 i ) y
( )
1
⃗v = ()x
y
=t 5
(1+8 i)
1
=t
(1+8 i)
5 ( )
La matriz (1+85 i ), es el vector propio del autovalor complejo
La solución del sistema de ecuación diferencial correspondiente es
X ( t )=e (2 +8 i) t (1+58i )
RECORDAR: e (2+ 8 i) t=e2 t (cos ( 8 t )+ isen ( 8 t ))
Por lo tanto:
2t
X 1 ( t )=e (cos ( 8t ) +isen ( 8 t ) ) (1+85 i)
Nota: i× i=−1
X 1 ( t )=e
2t
( cos ( 8 t ) +isen ( 8 t ) +8 i cos ( 8 t )−8 sen ( 8 t )
5 cos ( 8 t )+ 5isen ( 8 t ) )
X 1 ( t )=e
2t
(cos ( 85tcos)−8( 8sent ) ( 8 t ))+ e i ( sen ( 85tsen)+ 8( 8cost ) ( 8 t ))
2t
Tiene la forma de
Y 1 ( t ) =U +Vi
X ( t )=C 1 e
2t
( cos (58cost )−sen
(8 t )
( 8t )
) +C e
( sen ( 8 t )+ 8 cos ( 8 t )
2
2t
5 sen ( 8 t ) )
Nos dan la condición inicial de
x 1 ( 0 )=3 ; x 2 ( 0 )=−10
x (0)= (−10
3
)
X ( 0 )=C1 e
0
(cos5( 0cos)−sen
( 0)
( 0)
) +C e
(2
sen ( 0 ) +8 cos ( 0 )
0
5 sen ( 0 ) )
X ( 0 )=C1 (51)+C (80)=(−10
2
3
)
Se tiene como ecuaciones :
{ C1 +8 C 2=3
5 C 1=−10
5
∴ C1=−2 y C2 =
8
X ( t )=−2 e
2t
( 5 cos ( 8 t )
+ e
8 ) (
cos ( 8 t )−sen ( 8 t ) 5 2 t sen ( 8 t ) +8 cos ( 8 t )
5 sen ( 8t ) )
[
x 1 (t)=e2 t 3 cos ( 8t ) +
183
8
sen ( 8 t )
]
[
x 2 ( t )=e2 t −10 cos ( 8t ) +
25
8
sen ( 8 t )
]
EJERCICIO PROPUESTO: (Rainville et. al, 1998, p. 204)
7.
Primero realizamos la ecuación característica de esta matriz:
De modo que:
Por lo tanto, en correspondencia con el valor propio existe un conjunto de
De modo que:
De esta manera, así obtenemos dos conjuntos distintos de soluciones para el sistema:
6.
Primero hallamos la ecuación característica de la matriz del sistema:
Una solución que adquiera esta forma debe recordarnos la situación, donde resolvimos
ecuaciones lineales con coeficientes contantes. Procederemos de manera parecida,
haciendo uso de la fórmula de Euler:
Primera fila:
Parte real:
Parte imaginaria:
Segunda fila:
Parte real:
Parte imaginaria:
Reemplazando en la expresión:
Para b2:
, , reemplazando tenemos:
tendríamos: , el vector
propio será:
4° Escribir las soluciones:
Procedemos a descomponer la matriz del vector propio, de la siguiente manera:
, reemplazando tenemos,
Ejemplo 2)
=0
Resolviendo la matriz, tenemos de la siguiente manera:
valor de .
manera:
4° La solución seria de la siguiente manera:
, entonces el valor de ,
seguidamente aplicamos lo siguiente:
REFERENCIAS BIBLIOGRAFICAS
Ignacio García & Visiers Bañón (2020). Valores propios. Diagonalización de matrices.
https://www.cartagena99.com/recursos/alumnos/apuntes/Algebra%20Capitulo
%2005%20Diagonalizacion.pdf