Al transformar un sistema de ecuaciones en su forma diagonal, se simplifica considerablemente el proceso de encontrar la solución general. La utilidad de su uso proviene de las siguientes razones:
- Desacopla el sistema: La diagonalización permite desacoplar un sistema de ecuaciones diferenciales acopladas en un conjunto de ecuaciones diferenciales más simples, cada una involucrando una sola variable.
- Simplifica la solución: Las ecuaciones diferenciales desacopladas son mucho más fáciles de resolver que un sistema acoplado.
- Proporciona una interpretación clara: Los autovalores de la matriz diagonalizada proporcionan información sobre la estabilidad y el comportamiento a largo plazo del sistema.
Pasos para aplicar la Diagonalización en ecuaciones diferenciales
Consideremos un sistema de ecuaciones diferenciales de primer orden de la forma:
x'(t) = Ax(t)
donde x(t) es un vector de funciones desconocidas, A es una matriz cuadrada constante y x'(t) es la derivada de x(t) con respecto a t.
-
Diagonalizar la matriz A:
- Encontrar los autovalores y autovectores de A.
- Formar la matriz P cuyas columnas son los autovectores.
- Formar la matriz diagonal D cuyos elementos diagonales son los autovalores.
- A = PDP^(-1)
-
Cambio de variable:
- Introducir un nuevo vector y(t) definido por x(t) = Py(t).
- Sustituir en el sistema original:
Py'(t) = APy(t) - Multiplicando a ambos lados por P^(-1):
y'(t) = Dy(t)
-
Resolver el sistema desacoplado: El sistema resultante es un sistema de ecuaciones diferenciales desacopladas, donde cada ecuación involucra solo una componente de y(t). La solución general de este sistema es fácil de encontrar.
-
Volver a las variables originales: Una vez encontrada la solución y(t), se puede obtener la solución original x(t) mediante x(t) = Py(t).
Ejemplo
Consideremos el sistema:
x'(t) = 2x(t) + y(t)
y'(t) = x(t) + 2y(t)
En forma matricial:
[x'(t)] [2 1] [x(t)]
[y'(t)] = [1 2] [y(t)]
Diagonalizando la matriz A = [[2, 1], [1, 2]], obtenemos:
- Autovalores: λ₁ = 1, λ₂ = 3
- Autovectores: v₁ = [1, -1], v₂ = [1, 1]
- P = [[1, 1], [-1, 1]], D = [[1, 0], [0, 3]]
Realizando el cambio de variable y resolviendo el sistema desacoplado, encontramos la solución general del sistema original.
Comentarios
Publicar un comentario