Las ciencias matemáticas exhiben particularmente orden, simetría
y límites; y esas son las más grandes formas de belleza.
– Aristóteles

Introducción

Ya sabemos resolver sistemas lineales homogéneos con coeficientes constantes, en esta entrada estudiaremos el caso no homogéneo.

Como hemos visto en las dos unidades anteriores, el método de variación de parámetros ha sido fundamental para resolver ecuaciones diferenciales en el caso no homogéneo. Éste mismo método es el que desarrollaremos en esta entrada para resolver sistemas lineales no homogéneos con coeficientes constantes.

Cabe mencionar que en esta entrada utilizaremos bastante el concepto de exponencial de una matriz y el de matriz fundamental de soluciones.

Sistemas lineales homogéneos

El sistema lineal que hemos estudiado es

$$\begin{array}{}\text{(1)}& {\mathbf{Y}}^{\prime }=\mathbf{A}\mathbf{Y}\end{array}$$

Donde $\mathbf{A}$ es una matriz con componentes constantes

$$\begin{array}{}\text{(2)}& \mathbf{A}=\left(\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& & & ⋮\\ a_{n1}& a_{n2}& \cdots & a_{nn}\end{array}\right)\end{array}$$

Ahora sabemos que la solución general del sistema lineal homogéneo ($\text{1}$) es de la forma

$$\begin{array}{}\text{(3)}& {\mathbf{Y}}_c(t)=c_1{\mathbf{Y}}_1(t)+c_2{\mathbf{Y}}_2(t)+\cdots +c_n{\mathbf{Y}}_n(t)\end{array}$$

En donde las funciones ${\mathbf{Y}}_1(t),{\mathbf{Y}}_2(t),\cdots ,{\mathbf{Y}}_n(t)$ son soluciones linealmente independientes del mismo sistema. Usamos el subíndice $c$ debido a que, como antes, al resolver el caso no homogéneo será necesario resolver primero el sistema homogéneo asociado y la solución general de dicho sistema será la solución complementaria del sistema no homogéneo.

Recordemos que la matriz que tiene por columnas a las funciones ${\mathbf{Y}}_i(t)$, $i=1,2,\cdots ,n$ de ($\text{3}$) corresponde a la matriz fundamental de soluciones.

$$\begin{array}{}\text{(4)}& \hat{\mathbf{Y}}(t)=\left(\begin{array}{cccc}{\mathbf{Y}}_1(t)& {\mathbf{Y}}_2(t)& \cdots & {\mathbf{Y}}_n(t)\end{array}\right)=\left(\begin{array}{cccc}{y}_{11}(t)& y_{12}(t)& \cdots & y_{1n}(t)\\ y_{21}(t)& y_{22}(t)& \cdots & y_{2n}(t)\\ ⋮& & & ⋮\\ y_{n1}(t)& y_{n2}(t)& \cdots & y_{nn}(t)\end{array}\right)\end{array}$$

Si definimos el vector de constantes

$$\begin{array}{}\text{(5)}& \mathbf{C}=\left(\begin{array}{c}{c}_1\\ c_2\\ ⋮\\ c_n\end{array}\right)\end{array}$$

podemos escribir la solución ($\text{3}$) como

$$\begin{array}{}\text{(6)}& {\mathbf{Y}}_c(t)=\hat{\mathbf{Y}}(t)\mathbf{C}\end{array}$$

Recordemos este resultado para más adelante.

Sistemas lineales no homogéneos

El sistema lineal no homogéneo que intentaremos resolver es de la forma

$$\begin{array}{}\text{(7)}& {\mathbf{Y}}^{\prime }=\mathbf{A}\mathbf{Y}+\mathbf{G}\end{array}$$

Donde $\mathbf{G}=\mathbf{G}(t)$ es una matriz de $n\times 1$ con componentes dependientes de $t$.

$$\begin{array}{}\text{(8)}& \mathbf{G}(t)=\left(\begin{array}{c}{g}_1(t)\\ g_2(t)\\ ⋮\\ g_n(t)\end{array}\right)\end{array}$$

Dada la forma de la solución general de un sistema lineal homogéneo ($\text{3}$), parecería natural pensar que el sistema lineal no homogéneo tiene por solución una función de la forma

$$\begin{array}{}\text{(9)}& {\mathbf{Y}}_p(t)=u_1(t){\mathbf{Y}}_1(t)+u_2(t){\mathbf{Y}}_2(t)+\cdots +u_n(t){\mathbf{Y}}_n(t)\end{array}$$

En donde $u_i(t)$, $i=1,2,\cdots ,n$ son funciones escalares de $t$ derivables y las funciones ${\mathbf{Y}}_i(t)$, $i=1,2,\cdots ,n$ forman una matriz fundamental de soluciones $\hat{\mathbf{Y}}(t)$. Si definimos el vector

$$\begin{array}{}\text{(10)}& \mathbf{U}(t)=\left(\begin{array}{c}{u}_1(t)\\ u_2(t)\\ ⋮\\ u_n(t)\end{array}\right)\end{array}$$

Entonces la solución propuesta ($\text{9}$) adquiere la forma

$$\begin{array}{}\text{(11)}& {\mathbf{Y}}_p(t)=\hat{\mathbf{Y}}(t)\mathbf{U}(t)\end{array}$$

El método de variación de parámetros nos permitirá obtener la forma del vector $\mathbf{U}(t)$, una vez obtenida podremos formar la solución general del sistema lineal no homogéneo ($\text{7}$) que, como siempre, será la superposición de la solución complementaria del sistema lineal homogéneo asociado ${\mathbf{Y}}_c(t)$ más la solución particular del sistema lineal no homogéneo ${\mathbf{Y}}_p(t)$, esto es

$$\begin{array}{}\text{(12)}& \mathbf{Y}(t)={\mathbf{Y}}_c(t)+{\mathbf{Y}}_p(t)\end{array}$$

Variación de parámetros

Comencemos a desarrollar el método de variación de parámetros, como mencionamos antes, el objetivo es encontrar la forma explícita del vector ($\text{10}$) para formar la solución particular del sistema lineal no homogéneo.

Consideremos la función propuesta ($\text{11}$) y derivémosla.

$$\begin{array}{}\text{(13)}& {\mathbf{Y}}_p^{\prime }(t)={\hat{\mathbf{Y}}}^{\prime }(t)\mathbf{U}(t)+\hat{\mathbf{Y}}(t){\mathbf{U}}^{\prime }(t)\end{array}$$

Si sustituimos ($\text{11}$) y ($\text{13}$) en el sistema lineal no homogéneo ($\text{7}$), se tiene

$$\begin{array}{}\text{(14)}& {\hat{\mathbf{Y}}}^{\prime }(t)\mathbf{U}(t)+\hat{\mathbf{Y}}(t){\mathbf{U}}^{\prime }(t)=\mathbf{A}[\hat{\mathbf{Y}}(t)\mathbf{U}(t)]+\mathbf{G}(t)\end{array}$$

Como $\hat{\mathbf{Y}}(t)$ es una matriz fundamental de soluciones sabemos que satisface el sistema homogéneo, es decir,

$$\begin{array}{}\text{(15)}& {\hat{\mathbf{Y}}}^{\prime }(t)=\mathbf{A}\hat{\mathbf{Y}}(t)\end{array}$$

Si sustituimos en ($\text{14}$) la ecuación queda como

$$\mathbf{A}\hat{\mathbf{Y}}(t)\mathbf{U}(t)+\hat{\mathbf{Y}}(t){\mathbf{U}}^{\prime }(t)=\mathbf{A}\hat{\mathbf{Y}}(t)\mathbf{U}(t)+\mathbf{G}(t)$$

O bien,

$$\begin{array}{}\text{(16)}& \hat{\mathbf{Y}}(t){\mathbf{U}}^{\prime }(t)=\mathbf{G}(t)\end{array}$$

La matriz fundamental es no singular, de manera que siempre existe su inversa, esto nos permite establecer que

$$\begin{array}{}\text{(17)}& {\mathbf{U}}^{\prime }(t)={\hat{\mathbf{Y}}}^{-1}(t)\mathbf{G}(t)\end{array}$$

Esta ecuación es matricial y sabemos que es posible integrar sobre matrices, así que integremos la ecuación anterior con el objetivo de hallar la forma de $\mathbf{U}$.

$$\begin{array}{}\text{(18)}& \mathbf{U}(t)=\int {\hat{\mathbf{Y}}}^{-1}(t)\mathbf{G}(t)dt\end{array}$$

Ahora que conocemos la forma de $\mathbf{U}(t)$, sustituimos en la solución propuesta ($\text{11}$), de forma que una solución particular del sistema lineal no homogéneo es

$$\begin{array}{}\text{(19)}& {\mathbf{Y}}_p(t)=\hat{\mathbf{Y}}(t)\int {\hat{\mathbf{Y}}}^{-1}(t)\mathbf{G}(t)dt\end{array}$$

Por lo tanto, de ($\text{6}$) y el resultado ($\text{19}$) concluimos que la solución general del sistema lineal no homogéneo es

$$\begin{array}{}\text{(20)}& \mathbf{Y}(t)=\hat{\mathbf{Y}}(t)\mathbf{C}+\hat{\mathbf{Y}}(t)\int {\hat{\mathbf{Y}}}^{-1}(t)\mathbf{G}(t)dt\end{array}$$

Si $\hat{\mathbf{Y}}(t)$ es la matriz fundamental de soluciones $e^{\mathbf{A}t}$ y considerando que ${\hat{\mathbf{Y}}}^{-1}(t)=e^{-\mathbf{A}t}$, el resultado anterior queda como

$$\begin{array}{}\text{(21)}& \mathbf{Y}(t)=e^{\mathbf{A}t}\mathbf{C}+e^{\mathbf{A}t}\int e^{-\mathbf{A}t}\mathbf{G}(t)dt\end{array}$$

Problema con valores iniciales

Consideremos el problema con valores iniciales

$$\begin{array}{}\text{(22)}& {\mathbf{Y}}^{\prime }=\mathbf{A}\mathbf{Y}+\mathbf{G};\mathbf{Y}(t_0)={\mathbf{Y}}_0\end{array}$$

De nuestro desarrollo anterior consideremos la relación ($\text{17}$).

$${\mathbf{U}}^{\prime }(t)={\hat{\mathbf{Y}}}^{-1}(t)\mathbf{G}(t)$$

En esta ocasión integremos de $t_0$ a $t$ y usemos el teorema fundamental del cálculo.

$$\begin{array}{}\text{(23)}& \mathbf{U}(t)=\mathbf{U}(t_0)+{\int }_{t_0}^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds\end{array}$$

Notemos que si aplicamos la condición inicial $\mathbf{Y}(t_0)={\mathbf{Y}}_0$ sobre la función ($\text{11}$), se obtiene

$$\begin{array}{}\text{(24)}& {\mathbf{Y}}_p(t_0)=\hat{\mathbf{Y}}(t_0)\mathbf{U}(t_0)={\mathbf{Y}}_0\end{array}$$

De donde,

$$\begin{array}{}\text{(25)}& \mathbf{U}(t_0)={\hat{\mathbf{Y}}}^{-1}(t_0){\mathbf{Y}}_0\end{array}$$

Sustituimos este resultado en la ecuación ($\text{23}$).

$$\begin{array}{}\text{(26)}& \mathbf{U}(t)={\hat{\mathbf{Y}}}^{-1}(t_0){\mathbf{Y}}_0+{\int }_{t_0}^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds\end{array}$$

Aquí debemos tener cuidado, si sustituimos la función ($\text{26}$) en ($\text{11}$), se obtiene

$$\begin{array}{}\text{(27)}& \mathbf{Y}(t)=\hat{\mathbf{Y}}(t){\hat{\mathbf{Y}}}^{-1}(t_0){\mathbf{Y}}_0+\hat{\mathbf{Y}}(t){\int }_{t_0}^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds\end{array}$$

Pero ${\hat{\mathbf{Y}}}^{-1}(t_0){\mathbf{Y}}_0$ es una matriz de constantes, digamos

$${\hat{\mathbf{Y}}}^{-1}(t_0){\mathbf{Y}}_0=\mathbf{C}$$

Entonces el resultado ($\text{27}$) queda como

$$\begin{array}{}\text{(28)}& \mathbf{Y}(t)=\hat{\mathbf{Y}}(t)\mathbf{C}+\hat{\mathbf{Y}}(t){\int }_{t_0}^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds\end{array}$$

Este resultado se parece a la ecuación ($\text{20}$), es decir, a pesar de que sustituimos ($\text{26}$) en ($\text{11}$) esperando obtener la solución particular ${\mathbf{Y}}_p(t)$, en realidad estamos obteniendo la solución general, la solución general del problema de valores iniciales.

Si consideramos nuevamente que $\hat{\mathbf{Y}}(t)=e^{\mathbf{A}t}$, el resultado ($\text{27}$) se reduce significativamente.

$$\begin{array}{}\text{(29)}& \mathbf{Y}(t)=e^{\mathbf{A}t}{e}^{-\mathbf{A}{t}_0}{\mathbf{Y}}_0+e^{\mathbf{A}t}{\int }_{t_0}^t{e}^{-\mathbf{A}s}\mathbf{G}(s)ds\end{array}$$

O bien,

$$\begin{array}{}\text{(30)}& \mathbf{Y}(t)=e^{\mathbf{A}(t-t_0)}{\mathbf{Y}}_0+{\int }_{t_0}^t{e}^{\mathbf{A}(t-s)}\mathbf{G}(s)ds\end{array}$$

Por otro lado, si $t_0=0$, de ($\text{27}$) se obtiene que

$$\mathbf{Y}(t)=\hat{\mathbf{Y}}(t){\hat{\mathbf{Y}}}^{-1}(0){\mathbf{Y}}_0+\hat{\mathbf{Y}}(t){\int }_0^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds$$

Pero recordemos que

$$\begin{array}{}\text{(31)}& e^{\mathbf{A}t}=\hat{\mathbf{Y}}(t){\hat{\mathbf{Y}}}^{-1}(0)\end{array}$$

Entonces la solución anterior queda como

$$\begin{array}{}\text{(32)}& \mathbf{Y}(t)=e^{\mathbf{A}t}{\mathbf{Y}}_0+\hat{\mathbf{Y}}(t){\int }_0^t{\hat{\mathbf{Y}}}^{-1}(s)\mathbf{G}(s)ds\end{array}$$

Si nuestro propósito es determinar la solución general de un sistema lineal no homogéneo ($\text{7}$), primero resolvemos el sistema lineal homogéneo asociado para obtener la solución complementaria en la forma ($\text{3}$). Con las funciones ${\mathbf{Y}}_i(t)$, $i=1,2,\cdots ,n$ obtenidas formamos una matriz fundamental $\hat{\mathbf{Y}}(t)$, se calcula su inversa y se sustituyen las matrices correspondientes en la solución particular ($\text{19}$). Una vez obtenidas ambas soluciones, la solución general del sistema lineal no homogéneo será

$$\mathbf{Y}(t)=c_1{\mathbf{Y}}_1(t)+c_2{\mathbf{Y}}_2(t)+\cdots +c_n{\mathbf{Y}}_n(t)+{\mathbf{Y}}_p(t)={\mathbf{Y}}_c(t)+{\mathbf{Y}}_p(t)$$

Sin embargo, si lo que tenemos es un problema de valores iniciales, debemos nuevamente obtener la solución del sistema lineal homogéneo asociado ya que eso es lo que nos permite formar la matriz fundamental de soluciones $\hat{\mathbf{Y}}(t)$, una vez obtenida esta función calculamos su inversa y se sustituyen las matrices correspondientes en la ecuación ($\text{27}$), esto nos dará la solución completa del problema de valores iniciales, es decir, no es necesario aplicar las condiciones iniciales en la solución complementaria para obtener los valores de las constantes $c_1,c_2,\cdots ,c_n$.

Para concluir con esta entrada realicemos dos ejemplos, en el primero de ellos obtendremos la solución general de un sistema lineal no homogéneo y en el segundo ejemplo resolveremos un problema con valores iniciales. Con estos ejemplos se espera que el método quede claro.

Ejemplo: Obtener la solución general del siguiente sistema lineal no homogéneo.

$${\mathbf{Y}}^{\prime }=\left(\begin{array}{cc}0& 2\\ -1& 3\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}1\\ -1\end{array}\right)e^t$$

Solución: Resolvamos primero el sistema homogéneo asociado.

$${\mathbf{Y}}^{\prime }=\left(\begin{array}{cc}0& 2\\ -1& 3\end{array}\right)\mathbf{Y}$$

En este caso la matriz $\mathbf{A}$ es

$$\mathbf{A}=\left(\begin{array}{cc}0& 2\\ -1& 3\end{array}\right)$$

Determinemos los valores y vectores propios de esta matriz.

$$|\mathbf{A}-\lambda \mathbf{I}|=\left|\begin{array}{cc}-\lambda & 2\\ -1& 3-\lambda \end{array}\right|=0$$

La ecuación característica es

$${\lambda }^2-3\lambda +2=0$$

Resolviendo para $\lambda$ se obtiene que los valores propios son

$${\lambda }_1=1y{\lambda }_2=2$$

Determinemos los vectores propios correspondientes a cada valor propio.

Caso 1: ${\lambda }_1=1$.

Buscamos un vector $\mathbf{K}\ne \mathbf{0}$, tal que

$$(\mathbf{A}-\mathbf{I})\mathbf{K}=\left(\begin{array}{cc}-1& 2\\ -1& 2\end{array}\right)\left(\begin{array}{c}{k}_1\\ k_2\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right)$$

La ecuación que se obtiene es

$$-k_1+2k_2=0$$

Es decir, $2k_2=k_1$. Elegimos $k_2=1$, entonces $k_1=2$. Por lo tanto el primer vector propio es

$${\mathbf{K}}_1=\left(\begin{array}{c}2\\ 1\end{array}\right)$$

Caso 2: ${\lambda }_2=2$.

Buscamos un vector $\mathbf{K}\ne \mathbf{0}$, tal que

$$(\mathbf{A}-2\mathbf{I})\mathbf{K}=\left(\begin{array}{cc}-2& 2\\ -1& 1\end{array}\right)\left(\begin{array}{c}{k}_1\\ k_2\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right)$$

La ecuación que se obtiene es

$$-k_1+k_2=0$$

Es decir, $k_1=k_2$. Elegimos $k_1=k_2=1$. Por lo tanto, el segundo vector propio es

$${\mathbf{K}}_2=\left(\begin{array}{c}1\\ 1\end{array}\right)$$

Con estos resultados concluimos que la solución general del sistema lineal homogéneo asociado es

$${\mathbf{Y}}_c(t)=c_1{e}^t\left(\begin{array}{c}2\\ 1\end{array}\right)+c_2{e}^{2t}\left(\begin{array}{c}1\\ 1\end{array}\right)$$

Para determinar la solución particular ${\mathbf{Y}}_p(t)$, formemos, con el resultado anterior, la matriz fundamental de soluciones ($\text{4}$).

$$\hat{\mathbf{Y}}(t)=\left(\begin{array}{cc}2e^t& e^{2t}\\ e^t& e^{2t}\end{array}\right)$$

Como también requerimos de la matriz inversa, verifica que

$${\hat{\mathbf{Y}}}^{-1}(t)=\left(\begin{array}{cc}{e}^{-t}& -e^{-t}\\ -e^{-2t}& 2e^{-2t}\end{array}\right)$$

La matriz $\mathbf{G}$ en este caso es

$$\mathbf{G}(t)=\left(\begin{array}{c}{e}^t\\ -e^t\end{array}\right)$$

Sustituyamos estas matrices en la solución particular ($\text{19}$).

$$\begin{array}{rl}{\mathbf{Y}}_p(t)& =\left(\begin{array}{cc}2e^t& e^{2t}\\ e^t& e^{2t}\end{array}\right)\int \left(\begin{array}{cc}{e}^{-t}& -e^{-t}\\ -e^{-2t}& 2e^{-2t}\end{array}\right)\left(\begin{array}{c}{e}^t\\ -e^t\end{array}\right)dt\\ & =\left(\begin{array}{cc}2e^t& e^{2t}\\ e^t& e^{2t}\end{array}\right)\int \left(\begin{array}{c}2\\ -3e^{-t}\end{array}\right)dt\end{array}$$

Resolviendo la integral (sin considerar constantes de integración), se obtiene

$$\int \left(\begin{array}{c}2\\ -3e^{-t}\end{array}\right)dt=\left(\begin{array}{c}2t\\ 3e^{-t}\end{array}\right)$$

Entonces,

$${\mathbf{Y}}_p(t)=\left(\begin{array}{cc}2e^t& e^{2t}\\ e^t& e^{2t}\end{array}\right)\left(\begin{array}{c}2t\\ 3e^{-t}\end{array}\right)=\left(\begin{array}{c}4t{e}^t+3e^t\\ 2t{e}^t+3e^t\end{array}\right)$$

Esto es,

$${\mathbf{Y}}_p(t)=\left(\begin{array}{c}4\\ 2\end{array}\right)t{e}^t+\left(\begin{array}{c}3\\ 3\end{array}\right)e^t$$

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

$$\mathbf{Y}(t)=c_1{e}^t\left(\begin{array}{c}2\\ 1\end{array}\right)+c_2{e}^{2t}\left(\begin{array}{c}1\\ 1\end{array}\right)+\left(\begin{array}{c}3\\ 3\end{array}\right)e^t+\left(\begin{array}{c}4\\ 2\end{array}\right)t{e}^t$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Realicemos ahora un problema con valores iniciales.

Ejemplo: Resolver el siguiente problema con valores iniciales.

$${\mathbf{Y}}^{\prime }=\left(\begin{array}{ccc}1& 0& 0\\ 2& 1& -2\\ 3& 2& 1\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}0\\ 0\\ e^t\cos (2t)\end{array}\right),\mathbf{Y}(0)=\left(\begin{array}{c}0\\ 1\\ 1\end{array}\right)$$

Solución: Primero debemos obtener la solución del sistema lineal homogéneo asociado

$${\mathbf{Y}}^{\prime }=\left(\begin{array}{ccc}1& 0& 0\\ 2& 1& -2\\ 3& 2& 1\end{array}\right)\mathbf{Y}$$

Ello nos permitirá obtener la matriz fundamental de soluciones. En este caso la matriz $\mathbf{A}$ es

$$\mathbf{A}=\left(\begin{array}{ccc}1& 0& 0\\ 2& 1& -2\\ 3& 2& 1\end{array}\right)$$

Determinemos los valores y vectores propios de esta matriz. La ecuación característica se obtendrá de calcular el determinante

$$|\mathbf{A}-\lambda \mathbf{I}|=\left|\begin{array}{ccc}1-\lambda & 0& 0\\ 2& 1-\lambda & -2\\ 3& 2& 1-\lambda \end{array}\right|=0$$

Desarrollando el determinante obtendremos que

$$(1-\lambda )({\lambda }^2-2\lambda +5)=0$$

Resolviendo para $\lambda$ se obtiene que los valores propios de $\mathbf{A}$ son

$${\lambda }_1=1y{\lambda }_2=1+2i,{\lambda }_3=1-2i$$

De acuerdo a los valores propios obtenidos, la manera de resolver el sistema homogéneo será aplicando la teoría vista en la entrada sobre valores propios complejos.

Determinemos los vectores propios correspondientes a cada valor propio.

Caso 1: ${\lambda }_1=1$.

Buscamos un vector $\mathbf{K}\ne \mathbf{0}$, tal que

$$(\mathbf{A}-\mathbf{I})\mathbf{K}=[\left(\begin{array}{ccc}1& 0& 0\\ 2& 1& -2\\ 3& 2& 1\end{array}\right)-\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right)]\left(\begin{array}{c}{k}_1\\ k_2\\ k_3\end{array}\right)=\left(\begin{array}{ccc}0& 0& 0\\ 2& 0& -2\\ 3& 2& 0\end{array}\right)\left(\begin{array}{c}{k}_1\\ k_2\\ k_3\end{array}\right)=\left(\begin{array}{c}0\\ 0\\ 0\end{array}\right)$$

El sistema que se obtiene es

$$\begin{array}{rl}2k_1-2k_3& =0\\ 3k_1+2k_2& =0\end{array}$$

De este sistema se observa que

$$\begin{array}{rl}{k}_1& =k_3\\ k_2& =-{\displaystyle \frac{3k_1}{2}}\end{array}$$

Elegimos $k_1=2=k_3$, de tal manera que $k_2=-3$, así el primer vector propio es

$${\mathbf{K}}_1=\left(\begin{array}{c}2\\ -3\\ 2\end{array}\right)$$

Caso 2: ${\lambda }_2=1+2i$.

Buscamos un vector $\mathbf{K}$, diferente de cero, tal que

$$[\mathbf{A}-(1+2i)\mathbf{I}]\mathbf{K}=\left(\begin{array}{ccc}-2i& 0& 0\\ 2& -2i& -2\\ 3& 2& -2i\end{array}\right)\left(\begin{array}{c}{k}_1\\ k_2\\ k_3\end{array}\right)=\left(\begin{array}{c}0\\ 0\\ 0\end{array}\right)$$

El sistema de ecuaciones que se obtiene es

$$\begin{array}{rl}-2i{k}_1& =0\\ 2k_1-2i{k}_2-2k_3& =0\\ 3k_1+2k_2-2i{k}_3& =0\end{array}$$

De este sistema se observa que $k_1=0$ y $k_3=-i{k}_2$. Elegimos $k_2=1$, de manera que el segundo vector propio sea

$${\mathbf{K}}_2=\left(\begin{array}{c}0\\ 1\\ -i\end{array}\right)$$

Caso 3: ${\lambda }_2=1-2i$.

Sabemos que este caso es el conjugado del caso anterior, por lo que directamente establecemos que el tercer vector propio es

$${\mathbf{K}}_3=\left(\begin{array}{c}0\\ 1\\ i\end{array}\right)$$

La solución general del sistema lineal homogéneo asociado, en su forma compleja, es

$${\mathbf{Y}}_c(t)=c_1{e}^t\left(\begin{array}{c}2\\ -3\\ 2\end{array}\right)+c_2{e}^{(1+2i)t}\left(\begin{array}{c}0\\ 1\\ -i\end{array}\right)+c_3{e}^{(1-2i)t}\left(\begin{array}{c}0\\ 1\\ i\end{array}\right)$$

Sin embargo esta solución no nos sirve de mucho, pues desearíamos construir la matriz fundamental de soluciones con valores reales. Recordando lo visto en la entrada sobre valores propios complejos, podemos encontrar dos funciones ${\mathbf{W}}_1(t)$ y ${\mathbf{W}}_2(t)$, tal que la solución general sea de la forma

$$\begin{array}{}\text{(33)}& {\mathbf{Y}}_c(t)=c_1{e}^{{\lambda }_{1}t}+c_2{\mathbf{W}}_1(t)+c_3{\mathbf{W}}_2(t)\end{array}$$

Recordemos que las funciones ${\mathbf{W}}_1(t)$ y ${\mathbf{W}}_2(t)$ están dadas por

$$\begin{array}{}\text{(34)}& {\mathbf{W}}_1(t)=e^{\alpha t}[\mathbf{U}\cos (\beta t)-\mathbf{V}\sin (\beta t)]\end{array}$$

y

$$\begin{array}{}\text{(35)}& {\mathbf{W}}_2(t)=e^{\alpha t}[\mathbf{U}\sin (\beta t)+\mathbf{V}\cos (\beta t)]\end{array}$$

Consideremos el caso 2 en el que ${\lambda }_2=1+2i$ y

$${\mathbf{K}}_2=\left(\begin{array}{c}0\\ 1\\ -i\end{array}\right)=\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)+i\left(\begin{array}{c}0\\ 0\\ -1\end{array}\right)$$

De estos resultados obtenemos que $\alpha =1$, $\beta =2$ y

$$\mathbf{U}=\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)y\mathbf{V}=\left(\begin{array}{c}0\\ 0\\ -1\end{array}\right)$$

Sustituyamos en la funciones ($\text{34}$) y ($\text{35}$). Por un lado,

$$\begin{array}{rl}{\mathbf{W}}_1(t)& =e^t[\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)\cos (2t)-\left(\begin{array}{c}0\\ 0\\ -1\end{array}\right)\sin (2t)]\\ & =e^t[\left(\begin{array}{c}0\\ \cos (2t)\\ 0\end{array}\right)-\left(\begin{array}{c}0\\ 0\\ -\sin (2t)\end{array}\right)]\end{array}$$

Esto es,

$${\mathbf{W}}_1(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)\\ \sin (2t)\end{array}\right)$$

Por otro lado,

$$\begin{array}{rl}{\mathbf{W}}_2(t)& =e^t[\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)\sin (2t)+\left(\begin{array}{c}0\\ 0\\ -1\end{array}\right)\cos (2t)]\\ & =e^t[\left(\begin{array}{c}0\\ \sin (2t)\\ 0\end{array}\right)+\left(\begin{array}{c}0\\ 0\\ -\cos (2t)\end{array}\right)]\end{array}$$

Esto es,

$${\mathbf{W}}_2(t)=e^t\left(\begin{array}{c}0\\ \sin (2t)\\ -\cos (2t)\end{array}\right)$$

Recordemos que estas funciones también se pueden obtener considerando la identidad de Euler. Del caso 2 la solución que se obtiene es

$${\mathbf{Y}}_{2c}(t)=e^{(1+2i)t}\left(\begin{array}{c}0\\ 1\\ -i\end{array}\right)$$

Así,

$${\mathbf{W}}_1(t)=Re\{{\mathbf{Y}}_{2c}(t)\}y{\mathbf{W}}_2(t)=Im\{{\mathbf{Y}}_{2c}(t)\}$$

Usando la identidad de Euler sobre esta solución obtenemos lo siguiente.

$$\begin{array}{rl}{e}^{(1+2i)t}\left(\begin{array}{c}0\\ 1\\ -i\end{array}\right)& =e^t[\cos (2t)+i\sin (2t)][\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)-i\left(\begin{array}{c}0\\ 0\\ 1\end{array}\right)]\\ & =e^t[\cos (2t)\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)+\sin (2t)\left(\begin{array}{c}0\\ 0\\ 1\end{array}\right)]+i{e}^t[\sin (2t)\left(\begin{array}{c}0\\ 1\\ 0\end{array}\right)-\cos (2t)\left(\begin{array}{c}0\\ 0\\ 1\end{array}\right)]\\ & =e^t[\left(\begin{array}{c}0\\ \cos (2t)\\ 0\end{array}\right)+\left(\begin{array}{c}0\\ 0\\ \sin (2t)\end{array}\right)]+i{e}^t[\left(\begin{array}{c}0\\ \sin (2t)\\ 0\end{array}\right)+\left(\begin{array}{c}0\\ 0\\ -\cos (2t)\end{array}\right)]\end{array}$$

De donde,

$${\mathbf{W}}_1(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)\\ \sin (2t)\end{array}\right)y{\mathbf{W}}_2(t)=e^t\left(\begin{array}{c}0\\ \sin (2t)\\ -\cos (2t)\end{array}\right)$$

De esta forma, la solución general del sistema lineal homogéneo asociado es

$$\begin{array}{}\text{(36)}& {\mathbf{Y}}_c(t)=c_1{e}^t\left(\begin{array}{c}2\\ -3\\ 2\end{array}\right)+c_2{e}^t\left(\begin{array}{c}0\\ \cos (2t)\\ \sin (2t)\end{array}\right)+c_3{e}^t\left(\begin{array}{c}0\\ \sin (2t)\\ -\cos (2t)\end{array}\right)\end{array}$$

Esta solución es de la forma ($\text{3}$) por lo que la matriz fundamental de soluciones, formada por estos vectores linealmente independientes, es

$$\hat{\mathbf{Y}}(t)=\left(\begin{array}{ccc}2e^t& 0& 0\\ -3e^t& e^t\cos (2t)& e^t\sin (2t)\\ 2e^t& e^t\sin (2t)& -e^t\cos (2t)\end{array}\right)$$

Para obtener la solución del problema con valores iniciales usaremos el resultado ($\text{29}$) para $t_0=0$.

$$\begin{array}{}\text{(37)}& \mathbf{Y}(t)=e^{\mathbf{A}t}{\mathbf{Y}}_0+e^{\mathbf{A}t}{\int }_0^t{e}^{-\mathbf{A}s}\mathbf{G}(s)ds\end{array}$$

Es decir, consideraremos a la matriz $e^{\mathbf{A}t}$ como la matriz fundamental de soluciones. También es posible usar la relación ($\text{32}$) usando la matriz $\hat{\mathbf{Y}}(t)$ antes establecida. ¿Por qué son equivalentes ambas formas?.

Determinemos la matriz $e^{\mathbf{A}t}$ usando la relación ($\text{31}$). Si evaluamos $t=0$ en la matriz $\hat{\mathbf{Y}}(t)$ se obtiene la matriz

$$\hat{\mathbf{Y}}(0)=\left(\begin{array}{ccc}2& 0& 0\\ -3& 1& 0\\ 2& 0& -1\end{array}\right)$$

Comprueba que la matriz inversa es

$${\hat{\mathbf{Y}}}^{-1}(0)={\left(\begin{array}{ccc}2& 0& 0\\ -3& 1& 0\\ 2& 0& -1\end{array}\right)}^{-1}=\left(\begin{array}{ccc}{\displaystyle \frac{1}{2}}& 0& 0\\ {\displaystyle \frac{3}{2}}& 1& 0\\ 1& 0& -1\end{array}\right)$$

Sustituyamos en ($\text{31}$).

$$\begin{array}{rl}{e}^{\mathbf{A}t}& =\left(\begin{array}{ccc}2e^t& 0& 0\\ -3e^t& e^t\cos (2t)& e^t\sin (2t)\\ 2e^t& e^t\sin (2t)& -e^t\cos (2t)\end{array}\right)\left(\begin{array}{ccc}{\displaystyle \frac{1}{2}}& 0& 0\\ {\displaystyle \frac{3}{2}}& 1& 0\\ 1& 0& -1\end{array}\right)\\ & =e^t\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2t)+\sin (2t)& \cos (2t)& -\sin (2t)\\ 1+{\displaystyle \frac{3}{2}}\sin (2t)-\cos (2t)& \sin (2t)& \cos (2t)\end{array}\right)\end{array}$$

Por lo tanto, la matriz que consideraremos como matriz fundamental de soluciones es

$$e^{\mathbf{A}t}=e^t\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2t)+\sin (2t)& \cos (2t)& -\sin (2t)\\ 1+{\displaystyle \frac{3}{2}}\sin (2t)-\cos (2t)& \sin (2t)& \cos (2t)\end{array}\right)$$

Como también requerimos de la inversa de esta matriz, verifica que

$$e^{-\mathbf{A}t}=e^{-t}\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2t)-\sin (2t)& \cos (2t)& \sin (2t)\\ 1-{\displaystyle \frac{3}{2}}\sin (2t)-\cos (2t)& -\sin (2t)& \cos (2t)\end{array}\right)$$

En este caso la matriz $\mathbf{G}(t)$ es

$$\mathbf{G}(t)=\left(\begin{array}{c}0\\ 0\\ e^t\cos (2t)\end{array}\right)$$

Sustituyamos todos estos resultados en la solución ($\text{37}$).

$$\mathbf{Y}(t)=e^{\mathbf{A}t}\left(\begin{array}{c}0\\ 1\\ 1\end{array}\right)+e^{\mathbf{A}t}{\int }_0^t{e}^{-s}\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2s)-\sin (2s)& \cos (2s)& \sin (2s)\\ 1-{\displaystyle \frac{3}{2}}\sin (2s)-\cos (2s)& -\sin (2s)& \cos (2s)\end{array}\right)\left(\begin{array}{c}0\\ 0\\ e^s\cos (2s)\end{array}\right)ds$$

Por un lado,

$$e^{\mathbf{A}t}\left(\begin{array}{c}0\\ 1\\ 1\end{array}\right)=e^t\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2t)+\sin (2t)& \cos (2t)& -\sin (2t)\\ 1+{\displaystyle \frac{3}{2}}\sin (2t)-\cos (2t)& \sin (2t)& \cos (2t)\end{array}\right)\left(\begin{array}{c}0\\ 1\\ 1\end{array}\right)=e^t\left(\begin{array}{c}0\\ \cos (2t)-\sin (2t)\\ \cos (2t)+\sin (2t)\end{array}\right)$$

De tarea moral, determina las constantes $c_1$, $c_2$ y $c_3$ aplicando los valores iniciales sobre la solución complementaria ($\text{36}$). ¿Qué relación tiene tu resultado con la operación anterior?.

Por otro lado,

$$e^{-s}\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2s)-\sin (2s)& \cos (2s)& \sin (2s)\\ 1-{\displaystyle \frac{3}{2}}\sin (2s)-\cos (2s)& -\sin (2s)& \cos (2s)\end{array}\right)e^s\left(\begin{array}{c}0\\ 0\\ \cos (2s)\end{array}\right)=\left(\begin{array}{c}0\\ \sin (2s)\cos (2s)\\ {\cos }^2(2s)\end{array}\right)$$

Sustituimos estas matrices en $\mathbf{Y}(t)$.

$$\mathbf{Y}(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)-\sin (2t)\\ \cos (2t)+\sin (2t)\end{array}\right)+e^{\mathbf{A}t}{\int }_0^t\left(\begin{array}{c}0\\ \sin (2s)\cos (2s)\\ {\cos }^2(2s)\end{array}\right)ds$$

Resolvamos la integral.

$$\begin{array}{rl}{\int }_0^t\left(\begin{array}{c}0\\ \sin (2s)\cos (2s)\\ {\cos }^2(2s)\end{array}\right)ds& ={\left(\begin{array}{c}0\\ -{\displaystyle \frac{1}{8}}\cos (4s)\\ {\displaystyle \frac{s}{2}}+{\displaystyle \frac{\sin (4s)}{8}}\end{array}\right)|}_t–{\left(\begin{array}{c}0\\ -{\displaystyle \frac{1}{8}}\cos (4s)\\ {\displaystyle \frac{s}{2}}+{\displaystyle \frac{\sin (4s)}{8}}\end{array}\right)|}_0\\ & =\left(\begin{array}{c}0\\ -{\displaystyle \frac{1}{8}}\cos (4t)\\ {\displaystyle \frac{t}{2}}+{\displaystyle \frac{\sin (4t)}{8}}\end{array}\right)-\left(\begin{array}{c}0\\ -{\displaystyle \frac{1}{8}}\\ 0\end{array}\right)\\ & =\left(\begin{array}{c}0\\ {\displaystyle \frac{1-\cos (4t)}{8}}\\ {\displaystyle \frac{4t+\sin (4t)}{8}}\end{array}\right)\end{array}$$

Entonces,

$$\mathbf{Y}(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)-\sin (2t)\\ \cos (2t)+\sin (2t)\end{array}\right)+e^{\mathbf{A}t}\left(\begin{array}{c}0\\ {\displaystyle \frac{1-\cos (4t)}{8}}\\ {\displaystyle \frac{4t+\sin (4t)}{8}}\end{array}\right)$$

Ahora realicemos el producto del segundo sumando.

$$\begin{array}{rl}{e}^{\mathbf{A}t}\left(\begin{array}{c}0\\ {\displaystyle \frac{1-\cos (4t)}{8}}\\ {\displaystyle \frac{4t+\sin (4t)}{8}}\end{array}\right)& =e^t\left(\begin{array}{ccc}1& 0& 0\\ -{\displaystyle \frac{3}{2}}+{\displaystyle \frac{3}{2}}\cos (2t)+\sin (2t)& \cos (2t)& -\sin (2t)\\ 1+{\displaystyle \frac{3}{2}}\sin (2t)-\cos (2t)& \sin (2t)& \cos (2t)\end{array}\right)\left(\begin{array}{c}0\\ {\displaystyle \frac{1-\cos (4t)}{8}}\\ {\displaystyle \frac{4t+\sin (4t)}{8}}\end{array}\right)\\ & =e^t\left(\begin{array}{c}0\\ \cos (2t)\left({\displaystyle \frac{1-\cos (4t)}{8}}\right)-\sin (2t)\left({\displaystyle \frac{4t+\sin (4t)}{8}}\right)\\ \sin (2t)\left({\displaystyle \frac{1-\cos (4t)}{8}}\right)+\cos (2t)\left({\displaystyle \frac{4t+\sin (4t)}{8}}\right)\end{array}\right)\\ & =e^t\left(\begin{array}{c}0\\ -{\displaystyle \frac{t\sin (2t)}{2}}+{\displaystyle \frac{\cos (2t)-\cos (4t)\cos (2t)-\sin (4t)\sin (2t)}{8}}\\ {\displaystyle \frac{t\cos (2t)}{2}}+{\displaystyle \frac{\sin (2t)+\sin (4t)\cos (2t)-\cos (4t)\sin (2t)}{8}}\end{array}\right)\end{array}$$

Así,

$$\mathbf{Y}(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)-\sin (2t)\\ \cos (2t)+\sin (2t)\end{array}\right)+e^t\left(\begin{array}{c}0\\ -{\displaystyle \frac{t\sin (2t)}{2}}+{\displaystyle \frac{\cos (2t)-\cos (4t)\cos (2t)-\sin (4t)\sin (2t)}{8}}\\ {\displaystyle \frac{t\cos (2t)}{2}}+{\displaystyle \frac{\sin (2t)+\sin (4t)\cos (2t)-\cos (4t)\sin (2t)}{8}}\end{array}\right)$$

Haciendo las operaciones correspondientes se obtiene finalmente que la solución al problema con valores iniciales es

$$\mathbf{Y}(t)=e^t\left(\begin{array}{c}0\\ \cos (2t)-(1+{\displaystyle \frac{1}{2}}t)\sin (2t)\\ (1+{\displaystyle \frac{1}{2}}t)\cos (2t)+{\displaystyle \frac{5}{4}}\sin (2t)\end{array}\right)$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Vemos que este método puede ser bastante largo y complicado, de hecho se puede volver una tarea imposible de hacer a mano si se tienen sistemas con matriz $\mathbf{A}$ de $3\times 3$ o más. Se recomienda, en la medida de lo posible, usar algún programa computacional para llevar a cabo algunas de las operaciones, lo importante es entender como llevar a cabo el método.

Con esto concluimos lo que corresponde al estudio de los distintos métodos para resolver sistemas lineales. Prácticamente hemos concluido con la unidad 3 del curso.

En las siguientes dos entradas de esta unidad trataremos de justificar los teoremas de existencia y unicidad en el caso de los sistemas lineales, esto con el propósito de justificar toda la teoría desarrollada a lo largo de la unidad.

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. Obtener la solución general de los siguientes sistemas lineales no homogéneos.

• ${\mathbf{Y}}^{\prime }=\left(\begin{array}{cc}1& -1\\ 1& 1\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}\cos (t)\\ \sin (t)\end{array}\right)e^t$

• ${\mathbf{Y}}^{\prime }=\left(\begin{array}{cc}1& 2\\ -{\displaystyle \frac{1}{2}}& 1\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}\csc (t)\\ \sec (t)\end{array}\right)e^t$

• ${\mathbf{Y}}^{\prime }=\left(\begin{array}{ccc}1& 1& 0\\ 1& 1& 0\\ 0& 0& 3\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}{e}^t\\ e^{2t}\\ t{e}^{3t}\end{array}\right)$

2. Resolver los siguientes problemas con valores iniciales.

• ${\mathbf{Y}}^{\prime }=\left(\begin{array}{ccc}2& 0& 1\\ 0& 2& 0\\ 0& 1& 3\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}1\\ 0\\ 1\end{array}\right)e^{2t},\mathbf{Y}(0)=\left(\begin{array}{c}1\\ 1\\ 1\end{array}\right)$

• ${\mathbf{Y}}^{\prime }=\left(\begin{array}{ccc}-1& -1& -2\\ 1& 1& 1\\ 2& 1& 3\end{array}\right)\mathbf{Y}+\left(\begin{array}{c}1\\ 0\\ 0\end{array}\right)e^t,\mathbf{Y}(0)=\left(\begin{array}{c}0\\ 0\\ 0\end{array}\right)$

Más adelante…

En la siguiente entrada demostraremos los teoremas de existencia y unicidad para el caso de los sistemas lineales de primer orden con coeficientes constantes homogéneos y no homogéneos y posteriormente, en la última entrada de esta unidad, justificaremos el teorema de existencia y unicidad en el caso general, es decir, para sistemas lineales y no lineales.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales I
• Entrada anterior del curso: Sistemas lineales homogéneos con coeficientes constantes – Valores propios repetidos
• Siguiente entrada del curso: Teorema de existencia y unicidad para sistemas lineales

• Video relacionado al tema: Sistemas de ecuaciones lineales no homogéneas: solución por variación de parámetros

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»