Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

T3.2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

Taller Global 3.

2 Ecuaciones Diferenciales 2024-2

Resultados de aprendizaje

- Resolver sistemas de ecuaciones mediante el método de valores propios.

Contenidos para desarrollar en el taller

Un sistema de n ecuaciones diferenciales lineales de primer orden es de la forma


x′1 = a11 (t)x1 + a12 (t)x2 + · · · + a1n (t)xn + f1 (t)

x′2 = a21 (t)x1 + a22 (t)x2 + · · · + a2n (t)xn + f2 (t)


. (1)
x′3 = a31 (t)x1 + a32 (t)x2 + · · · + a3n (t)xn + f3 (t)
..
.
x′n = an1 (t)x1 + an2 (t)x2 + · · · + ann (t)xn + fn (t)

Si definimos A(t) = [aij (t)]n y los vectores columnas X = [xi ]n×1 y f = [fi (t)]n×1 entonces (1) tiene la
siguiente forma

X ′ = A(t)X + f (t) (2)

Una solución de la ecuación (2) en el intervalo abierto I es un vector de funciones columna X(t) =
[xi (t)]n×1 tal que las funciones que conforman X satisfacen el sistema (1) en I.
Al igual que en una dimensión, si las funciones aij y fi son continuas en I, entonces existe una única
solución X que satisface las condiciones iniciales x1 (a) = b1 , x2 (a) = b2 ... xn (a) = bn .
Si f = 0 entonces diremos que la ecuación (2) es homogénea.
Para la ecuación homogénea X ′ = A(t)X se sigue cumpliendo el principio de superposición, esto
quiere decir que si X1 , ..., Xn son soluciones de la ecuación lineal homogénea en un intervalo abierto
I entonces la combinación lineal

X(t) = c1 X1 (t) + c2 X2 (t) + · · · + cn Xn (t), ci ∈ R ∀i ∈ {1, ..., n}

es también una solución de la ecuación lineal homogénea en I.


Al igual que en el caso de una ecuación de orden n, existe un Wronskiano que nos dice si n soluciones
(que ahora son vectores) de la ecuación homogénea X ′ = A(t)X son linealmente dependientes o no.
Si X1 , ..., Xn son tales soluciones entonces se define el Wronskiano como
 
x11 (t) x12 (t) · · · x1n (t) x1j
x21 (t) x22 (t) · · · x2n (t)  x2j 
W (t) = .. , donde Xj =  ..  , ∀j ∈ {1, ..., n}
 
.. ..
. . .  . 
xn1 (t) xn2 (t) · · · xnn (t) xnj
Si {X1 , ..., Xn } es linealmente dependiente en I, entonces W = 0 en cada punto de I.
Si {X1 , ..., Xn } es linealmente independiente en I, entonces W ̸= 0 en cada punto de I.
Sean X1 , ..., Xn soluciones L.I. de la ecuación lineal homogénea X ′ = A(t)X en I, donde A es continua.
Si X(t) es cualquier solución de X ′ = A(t)X en I, entonces existen valores c1 , ..., cn tales que

X(t) = c1 X1 (t) + c2 X2 (t) + · · · cn Xn (t)

En el caso de una ecuación no homogénea X ′ = A(t)X + f (t) tenemos que la solución general está
dada por
X(t) = c1 X1 (t) + c2 X2 (t) + · · · cn Xn (t) + Xp (t),
donde Xp es una solución particular de la ecuación no homogénea X ′ = A(t)X + f (t).

Sistema de ecuaciones diferenciales homogéneo de coeficientes constantes


Método de los valores propios

Considere la ecuación (2) con A(t) = A ∈ Mn (R) para todo t ∈ R, es decir, con coeficientes constantes.
Para resolver la ecuación

X ′ (t) = AX, (3)

será fundamental el cálculo de valores y vectores propios. Esto se puede ver en el siguiente Teorema:

Teorema: Sea λ un valor propio de la matriz de coeficientes A del sistema (3). Si v es un vector propio
asociado a λ, entonces
X(t) = veλt
es una solución no trivial del sistema.

Método de los valores propios: La descripción general de este método para resolver el sistema (3),
usando el teorema anterior, es de la siguiente forma

1. Primero se calculan los valores propios de A mediante la ecuación característica det(A − λIn ) = 0.
2. Luego, se intenta encontrar n vectores propios linealmente independiente v1 , v2 , ..., vn asociados con
los valores propios obtenidos en el paso anterior.
3. En el paso 2 no siempre es posible, pero cuando lo es, se obtienen n soluciones linealmente indepen-
dientes.
X1 (t) = v1 eλ1 t , X2 (t) = v2 eλ2 t , . . . , Xn (t) = vn eλn t
En este caso, la solución general de X ′ = AX es una combinación lineal

X(t) = c1 X1 (t) + c2 X2 (t) + · · · + cn Xn (t),

de estas n soluciones.

A partir de los pasos anteriores se le pueden presentar lo siguientes casos:

Valores propios distintos: Si lo valores propios λ1 , λ2 , ..., λn son reales y distintos, entonces cada uno se
sustituye en la ecuación (A − λIn )v = 0 para obtener los vectores propios asociados v1 , ..., vn . En tal caso
la solución del sistema homogéneo (3) es

X(t) = c1 v1 eλ1 t + c2 v2 eλ2 t + · · · cn vn eλn t


Valores propios repetidos: Si A es diagonalizable y tiene un valor propio λ con multiplicidad algebrica
m entonces es posible encontrar m vectores propios linealmente independientes v1 , ..., vm correspondientes
al valor propio λ. Así, una parte de la solución general de la ecuación homogénea (3) es
c1 v1 eλt + c2 v2 eλt + · · · + cm vm eλt .

Valores propios complejos: Si λ = p + qi es un valor propio de A, calculamos el vector propio


v = a + bi asociado a λ. Luego, la solución de la ecuación (3) es
X(t) = veλt = ve(p+qi)t = (a + bi)ept (cos(qt) + isen(qt)),
esto es
X(t) = ept (a cos(qt) − bsen(qt)) + iept (b cos qt + asen(qt)).

Finalmente, debido a que las partes real e imaginaria de una solución de valores complejos son también
soluciones, se obtienen las dos soluciones de valores reales
X1 (t) = Re[X(t)] = ept (a cos(qt) − bsen(qt))
X2 (t) = Im[X(t)] = ept (b cos(qt) + asen(qt))

y por lo tanto, parte de la solución del sistema (3) es


c1 X1 + c2 X2

Ejemplo: Resolvamos el sistema  


1 −2 2
X ′ =  −2 1 −2  X
2 −2 1

Para esto, notemos que


1 − λ −2 2
det(A − λI3 ) = −2 1 − λ −2 = 0 ⇔ −(λ + 1)2 (λ − 5) = 0.
2 −2 1 − λ

Así, tenemos que λ1 = −1 (con multiplicidad algebraica 2) y λ2 = 5. Luego, tenemos que

   
1 0
Para λ1 = −1 se obtienen los vectores propios v1 = 1 y v2 = 1
0 1
 
1
Para λ2 = 5 se obtiene el vector propio v3 = −1

1

Por lo tanto, la solución general es


     
1 0 1
X(t) = c1  1  e−t + c2  1  e−t + c3  −1  e5t , c1 , c2 , c3 ∈ R3 .
0 1 1

Observación: Recuerde que para sympy existe el comando .eigenvects para obtener valores y vectores
propios de una matriz.
Ejercicios para desarrollar en el taller

1. Considere el sistema
dx1
= 4x1 − 3x2
dt
dx2
= 3x1 + 4x2
dt
(a) Determine dos soluciones X1 , X2 del sistema.

(b) Verifique, usando el Wronskiano, que las dos soluciones que obtuvo son linealmente independien-
tes.

(c) Concluya la solución general del sistema.

También podría gustarte