En las matemáticas no entiendes las cosas. Te acostumbras a ellas.
– John Von Neumann

Introducción

Estamos listos para comenzar a desarrollar los distintos métodos de resolución de sistemas lineales de primer orden.

En esta entrada desarrollaremos un método relativamente sencillo, pero muy limitado, ya que en general se utiliza cuando sólo tenemos un sistema lineal de dos ecuaciones diferenciales. Este método se conoce como método de eliminación de variables y, como su nombre lo indica, lo que se intenta hacer es eliminar las variables dependientes de $t$ hasta quedarnos con sólo una, esto produce que el resultado sea una sola ecuación diferencial de orden superior (la ecuación correspondiente a la única variable dependiente que nos queda), la cual es posible resolver aplicando alguno de los métodos vistos en la unidad anterior, la solución de dicha ecuación diferencial servirá para obtener el resto de funciones solución del sistema lineal.

Es importante mencionar que para que este método sea práctico y sencillo se requiere que los coeficientes de las ecuaciones que conforman al sistema lineal sean constantes y como el problema se reduce a resolver una ecuación de orden superior es conveniente usar este método sólo cuando tenemos dos ecuaciones diferenciales en el sistema, ya que esto involucrará resolver una ecuación diferencial de segundo orden con coeficientes constantes.

Desarrollemos el método de manera general.

Método de eliminación de variables

Los sistemas de ecuaciones diferenciales que estamos estudiando son de la forma

\begin{align*}
y_{1}^{\prime}(t) &= a_{11}(t)y_{1} + a_{12}(t)y_{2} + \cdots + a_{1n}(t)y_{n} + g_{1}(t) \\
y_{2}^{\prime}(t) &= a_{21}(t)y_{1} + a_{22}(t)y_{2} + \cdots + a_{2n}(t)y_{n} + g_{2}(t) \\
&\vdots \\
y_{n}^{\prime}(t) &= a_{n1}(t)y_{1} + a_{n2}(t)y_{2} + \cdots + a_{nn}(t)y_{n} + g_{n}(t) \label{1} \tag{1}
\end{align*}

Este método lo desarrollaremos para un sistema lineal de dos ecuaciones diferenciales lineales de primer orden tanto homogéneas como no homogéneas. De manera general desarrollemos el caso no homogéneo, el caso homogéneo será un caso particular.

Consideremos el siguiente sistema de ecuaciones diferenciales en su forma normal.

$$\begin{align*}
y_{1}^{\prime}(t) &= a_{11}(t)y_{1} + a_{12}(t)y_{2} + g_{1}(t) \\
y_{2}^{\prime}(t) &= a_{21}(t)y_{1} + a_{22}(t)y_{2} + g_{2}(t)
\end{align*} \label{2} \tag{2}$$

Debido a que se trata de un sistema pequeño regresemos a nuestra notación usual de derivada y sean $x$ y $y$ las variables dependientes de la variable independiente $t$. Así mismo, usemos una distinta notación para los coeficientes $a_{i, j}$, $i, j \in \{1, 2\}$, de tal manera que el sistema lineal (\ref{2}) lo podamos escribir de la siguiente forma.

$$\begin{align*}
\dfrac{dx}{dt} &= ax + by + g_{1}(t) \\
\dfrac{dy}{dt} &= cx + dy + g_{2}(t)
\end{align*}\label{3} \tag{3}$$

Con $a$, $b$, $c$ y $d$ constantes. El método que desarrollaremos es para sistema de la forma (\ref{3}).

De la primer ecuación del sistema despejamos a la variable $y$.

$$y = \dfrac{1}{b} \left( \dfrac{dx}{dt} -ax -g_{1} \right) \label{4} \tag{4}$$

Sustituyamos en la segunda ecuación.

$$\dfrac{d}{dt} \left[ \dfrac{1}{b} \left( \dfrac{dx}{dt} -ax -g_{1} \right) \right] = cx + d \left[ \dfrac{1}{b} \left( \dfrac{dx}{dt} -ax -g_{1} \right) \right] + g_{2}$$

Derivemos en el lado izquierdo y operemos en el lado derecho de la ecuación.

\begin{align*}
\dfrac{1}{b} \left[ \dfrac{d^{2}x}{dt^{2}} -a \dfrac{dx}{dt} -\dfrac{dg_{1}}{dt} \right] &= cx + \dfrac{1}{b} \left( d \dfrac{dx}{dt} -adx -dg_{1} \right) + g_{2} \\
\dfrac{d^{2}x}{dt^{2}} -a \dfrac{dx}{dt} -\dfrac{dg_{1}}{dt} &= bcx + d \dfrac{dx}{dt} -adx -dg_{1} + bg_{2}
\end{align*}

Reordenando los términos se tiene lo siguiente.

$$\dfrac{d^{2}x}{dt^{2}} -(a + d) \dfrac{dx}{dt} + (ad -bc) x = \dfrac{dg_{1}}{dt} -dg_{1} + bg_{2} \label{5} \tag{5}$$

Si definimos

$$p = -(a + d), \hspace{1cm} q = (ad -bc) \hspace{1cm} y \hspace{1cm} g(t) = \dfrac{dg_{1}}{dt} -dg_{1} + bg_{2}$$

entonces el resultado (\ref{5}) se puede escribir como

$$\dfrac{d^{2}x}{dt^{2}} + p \dfrac{dx}{dt} + q x = g(t) \label{6} \tag{6}$$

Con $p$ y $q$ constantes. En esta forma es claro que tenemos una ecuación diferencial lineal de segundo orden con coeficientes constantes, basta resolver la ecuación usando los métodos desarrollados en la unidad anterior para obtener la función $x(t)$. Una vez obtenida la solución de (\ref{6}) sustituimos en el despeje inicial que hicimos para $y(t)$ (\ref{4}) y resolvemos, con ello estaremos obteniendo la solución del sistema lineal (\ref{3}).

Caso homogéneo

El caso homogéneo es un caso particular del desarrollo anterior, pues el sistema a resolver es

$$\begin{align*}
\dfrac{dx}{dt} &= ax + by \\
\dfrac{dy}{dt} &= cx + dy
\end{align*}\label{7} \tag{7}$$

El desarrollo es exactamente el mismo considerando que $g_{1}(t) = 0$ y $g_{2}(t) = 0$.

Despejando a $y$ de la primer ecuación, obtenemos

$$y = \dfrac{1}{b} \left( \dfrac{dx}{dt} -ax \right) \label{8} \tag{8}$$

Sustituyendo en la segunda ecuación y siguiendo el mismo procedimiento obtendremos que la ecuación diferencial de segundo orden homogénea para $x$ es

$$\dfrac{d^{2}x}{dt^{2}} -(a + d) \dfrac{dx}{dt} + (ad -bc)x = 0 \label{9} \tag{9}$$

Si nuevamente definimos

$$p = -(a + d), \hspace{1cm} y \hspace{1cm} q = (ad -bc)$$

entonces podemos escribir

$$\dfrac{d^{2}x}{dt^{2}} + p \dfrac{dx}{dt} + qx = 0 \label{10} \tag{10}$$

Resolvamos un par de ejemplos, comencemos con un sistema lineal homogéneo.

Ejemplo: Resolver el siguiente sistema lineal homogéneo.

\begin{align*}
\dfrac{dx}{dt} &= 2x -y \\
\dfrac{dy}{dt} &= 5x -2y
\end{align*}

Solución: Comencemos por despejar a la variable $y$ de la primer ecuación.

$$y = 2x -\dfrac{dx}{dt}$$

Sustituimos en la segunda ecuación.

$$\dfrac{d}{dt} \left( 2x -\dfrac{dx}{dt} \right) = 5x -2 \left( 2x -\dfrac{dx}{dt} \right)$$

Operando, se tiene

\begin{align*}
2 \dfrac{dx}{dt} -\dfrac{d^{2}x}{dt^{2}} &= 5x -4x + 2 \dfrac{dx}{dt} \\
-\dfrac{d^{2}x}{dt^{2}} &= x
\end{align*}

La ecuación de segundo orden a resolver es

$$\dfrac{d^{2}x}{dt^{2}} + x = 0$$

Por supuesto esta ecuación se puede obtener sustituyendo los coeficientes directamente en la ecuación (\ref{9}).

Resolvamos la ecuación. La ecuación auxiliar es

$$k^{2} + 1 = 0$$

cuyas raíces son $k_{1} = i$ y $k_{2} = -i$.

Recordemos que la forma de la solución para raíces complejas $k_{1} = \alpha + i \beta$ y $k_{2} = \alpha -i \beta$ es

$$x(t) =e^{\alpha t}(c_{1} \cos(\beta t) + c_{2} \sin(\beta t)) \label{11} \tag{11}$$

En nuestro caso $\alpha =0$ y $\beta = 1$, entonces la solución es

$$x(t) = c_{1} \cos(t) + c_{2} \sin(t)$$

Vemos que

$$\dfrac{dx}{dt} = -c_{1} \sin(t) + c_{2} \cos(t)$$

Sustituimos en el despeje de $y$.

\begin{align*}
y(x) &= 2(c_{1} \cos(t) + c_{2} \sin(t)) -(-c_{1} \sin(t) + c_{2} \cos(t)) \\
&= 2c_{1} \cos(t) + 2c_{2} \sin(t) + c_{1} \sin(t) -c_{2} \cos(t)
\end{align*}

Esta solución la podemos escribir de dos formas.

$$y(x) = c_{1}(2 \cos(t) + \sin(t)) + c_{2}(2 \sin(t) -\cos(t))$$

o bien,

$$y(x) = (2c_{1} -c_{2})\cos(t) + (c_{1} + 2c_{2})\sin(t)$$

Por lo tanto, la solución general del sistema homogéneo es

$$\begin{pmatrix}
x \\ y
\end{pmatrix} = c_{1} \begin{pmatrix}
\cos(t) \\ 2 \cos(t) + \sin(t)
\end{pmatrix} + c_{2} \begin{pmatrix}
\sin(t) \\ 2 \sin(t) -\cos(t)
\end{pmatrix}$$

o bien,

$$\begin{pmatrix}
x \\ y
\end{pmatrix} = \begin{pmatrix}
c_{1} \\ 2c_{1} -c_{2}
\end{pmatrix} \cos(t) + \begin{pmatrix}
c_{2} \\ c_{1} + 2c_{2}
\end{pmatrix} \sin(t)$$

$\square$

Ahora resolvamos un sistema no homogéneo como ejemplo.

Ejemplo: Resolver el siguiente sistema lineal no homogéneo.

\begin{align*}
\dfrac{dx}{dt} &= 4x -y + t + 1 \\
\dfrac{dy}{dt} &= 2x + y + t + 1
\end{align*}

Solución: En este caso no homogéneo se tiene que

$$g_{1}(t) = t + 1 = g_{2}(t)$$

De la primer ecuación despejamos a $y$.

$$y = 4x + t + 1 -\dfrac{dx}{dt}$$

Sustituimos en la segunda ecuación.

$$\dfrac{d}{dt} \left( 4x + t + 1 -\dfrac{dx}{dt} \right) = 2x + \left( 4x + t + 1 -\dfrac{dx}{dt} \right) + t + 1$$

En el lado izquierdo aplicamos la derivada y en el lado izquierdo operamos.

\begin{align*}
4 \dfrac{dx}{dt} + \dfrac{d}{dt}(t + 1) -\dfrac{d^{2}x}{dt^{2}} &= 6x -\dfrac{dx}{dt} + 2t + 2 \\
4 \dfrac{dx}{dt} + 1 -\dfrac{d^{2}x}{dt^{2}} &= 6x -\dfrac{dx}{dt} + 2t + 2
\end{align*}

Reordenando los términos, se tiene

\begin{align*}
5 \dfrac{dx}{dt} -\dfrac{d^{2}x}{dt^{2}} &= 6x + 2t + 1 \\
-\dfrac{d^{2}x}{dt^{2}} + 5 \dfrac{dx}{dt} -6x &= 2t + 1 \\
\end{align*}

La ecuación diferencial de segundo orden no homogénea a resolver es

$$\dfrac{d^{2}x}{dt^{2}} -5 \dfrac{dx}{dt} + 6x = -2t -1$$

Para obtener la función $x(t)$ primero resolveremos el caso homogéneo y posteriormente aplicaremos el método de coeficientes indeterminados para resolver el caso no homogéneo. Recordemos que la solución general será la superposición de ambos resultados.

$$x(t) = x_{c}(t) + x_{p}(t) \label{12} \tag{12}$$

Para el caso homogéneo la ecuación a resolver es

$$\dfrac{d^{2}x}{dt^{2}} -5 \dfrac{dx}{dt} + 6x = 0$$

La ecuación auxiliar es

$$k^{2} -5k + 6 = 0$$

Resolviendo para $k$ se obtiene que $k_{1} = 2$ y $k_{2} = 3$. Como las raíces son reales y distintas, la forma de la solución es

$$x_{c}(t) = c_{1}e^{k_{1}t} + c_{2}e^{k_{2}t} \label{13} \tag{13}$$

Por lo tanto, la solución complementaria es

$$x_{c}(t) = c_{1}e^{2t} + c_{2}e^{3t}$$

Ahora resolvamos la ecuación no homogénea.

$$\dfrac{d^{2}x}{dt^{2}} -5 \dfrac{dx}{dt} + 6x = -2t -1$$

En este caso la función $g$ es

$$g(x) = -2t -1$$

la cual corresponde a un polinomio de grado $1$, entonces proponemos que la solución particular tiene, de igual manera, la forma de un polinomio de grado $1$, esto es

$$x_{p}(t) = At + B$$

Con $A$ y $B$ constantes por determinar. La primera y segunda derivada están dadas como

$$\dfrac{dx_{p}}{dt} = A \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}x}{dt^{2}} = 0$$

Sustituimos en la ecuación diferencial.

$$0 -5A + 6(At + B) = -2t -1$$

Reordenando, se tiene

$$6At + 6B -5A = -2t -1$$

Para que se cumpla la igualdad es necesario que ocurra lo siguiente.

\begin{align*}
6A &= -2 \\
6B -5A &= -1
\end{align*}

De la primer igualdad se obtiene que

$$A = -\dfrac{1}{3}$$

Sustituyendo este resultado en la segunda igualdad se obtiene que

$$B = \dfrac{1}{9}$$

Por lo tanto, la solución particular es

$$x_{p}(t) = -\dfrac{1}{3}t + \dfrac{1}{9}$$

Entonces concluimos que la solución general de la ecuación diferencial de segundo orden para $x$ es

$$x(t) = c_{1}e^{2t} + c_{2}e^{3t} -\dfrac{1}{3}t + \dfrac{1}{9}$$

Sustituimos este resultado en la ecuación de $y$.

$$y = 4 \left( c_{1}e^{2t} + c_{2}e^{3t} -\dfrac{1}{3}t + \dfrac{1}{9} \right) + t + 1 -\dfrac{d}{dt} \left( c_{1}e^{2t} + c_{2}e^{3t} -\dfrac{1}{3}t + \dfrac{1}{9} \right)$$

Operando, se tiene

$$y = 4c_{1}e^{2t} + 4c_{2}e^{3t} -\dfrac{4}{3}t + \dfrac{4}{9} + t + 1 -2c_{1}e^{2t} -3c_{2}e^{3t} + \dfrac{1}{3}$$

De donde se obtiene finalmente que la solución $y(t)$ es

$$y(x) = 2c_{1}e^{2t} + c_{2}e^{3t} -\dfrac{1}{3}t + \dfrac{16}{9}$$

Por lo tanto, la solución general del sistema lineal no homogéneo es

$$\begin{pmatrix}
x \\ y
\end{pmatrix} = c_{1} \begin{pmatrix}
1 \\ 2
\end{pmatrix} e^{2t} + c_{2} \begin{pmatrix}
1 \\ 1
\end{pmatrix}e^{3t} -\begin{pmatrix}
\dfrac{1}{3} \\ \dfrac{1}{3}
\end{pmatrix}t + \begin{pmatrix}
\dfrac{1}{9} \\ \dfrac{16}{9}
\end{pmatrix}$$

$\square$

Hemos concluido con esta entrada. Este método resulta sencillo y práctico para resolver sistemas lineales de este tipo, sin embargo está limitado a sistemas pequeños y realmente estamos interesados en resolver sistemas mucho más complejos.

En las siguientes entradas desarrollaremos otros métodos más generales.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

1. Resolver los siguientes sistemas lineales homogéneos.

• $\begin{align*}
  \dfrac{dx}{dt} &= x + 2y \\
  \dfrac{dy}{dt} &= 4x + 3y
  \end{align*}$

• $\begin{align*}
  \dfrac{dx}{dt} &= 2x -y \\
  \dfrac{dy}{dt} &= 3x -2y
  \end{align*}$

• $\begin{align*}
  \dfrac{dx}{dt} &= x -4y \\
  \dfrac{dy}{dt} &= -x + 2y
  \end{align*}$

• $\begin{align*}
  \dfrac{dx}{dt} = 2x -3y \\
  \dfrac{dy}{dt} = 3x + 2y
  \end{align*}$

2. Resolver los siguientes sistemas lineales no homogéneos.

• $\begin{align*}
  \dfrac{dx}{dt} &= 2x -y + 3t \\
  \dfrac{dy}{dt} &= 3x -2y + 2t + 4
  \end{align*}$

• $\begin{align*}
  \dfrac{dx}{dt} &= x + 2y + e^{t} \\
  \dfrac{dy}{dt} &= 3x -2y + 3e^{2t} + 2
  \end{align*}$

Más adelante…

En esta entrada presentamos un método sencillo para resolver sistemas lineales compuestos por dos ecuaciones diferenciales lineales de primer orden con coeficientes constantes tanto homogéneas como no homogéneas.

En la siguiente entrada comenzaremos a desarrollar otros métodos de resolución a sistemas lineales, sin embargo estos métodos suelen ser tratados desde una perspectiva del álgebra lineal, así que será importante hacer una pequeño repaso de algunos conceptos y teoremas de álgebra lineal. Unos de los conceptos más importantes que utilizaremos es el de valores y vectores propios.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales I
• Entrada anterior del curso: Exponencial de una matriz y matriz fundamental de soluciones
• Siguiente entrada del curso: Valores y vectores propios para resolver sistemas lineales

• Video relacionado al tema: Método de eliminación de variables

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»