Archivo de la etiqueta: Ecuaciones de segundo orden

Ecuaciones Diferenciales I: Ecuaciones lineales no homogéneas de segundo orden – Método de coeficientes indeterminados

Introducción

En la entrada anterior estudiamos las ecuaciones diferenciales de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0 \label{1} \tag{1}$$

Con $a$, $b$ y $c$ constantes reales. Vimos que, una vez que hemos determinado la ecuación auxiliar , las soluciones están dadas de acuerdo al valor del discriminante $\Delta = b^{2} -4ac$.

En esta entrada estudiaremos el mismo tipo de ecuación pero en el caso no homogéneo, es decir, de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = g(x) \label{2} \tag{2}$$

Con $a$, $b$, $c$ constantes y $g(x) \neq 0$.

Ya sabemos que, en general, para resolver una ecuación diferencial no homogénea de la forma

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n -1} \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{1} \dfrac{dy}{dx} + a_{0}y = g(x) \label{3} \tag{3}$$

debemos hacer dos cosas: encontrar la función complementaria $y_{c}(x)$ que corresponde a la solución general de la ecuación homogénea asociada de (\ref{3}) y posteriormente encontrar la solución particular $y_{p}(x)$ de dicha ecuación. Una vez obtenidas, la solución general de (\ref{3}) será la superposición $y(x) = y_{c}(x) + y_{p}(x)$. En la entrada anterior vimos como resolver el caso homogéneo de este tipo de ecuaciones cuando los coeficientes son constantes por lo que en esta entrada nos ocuparemos de desarrollar un método que nos permita obtener soluciones particulares.

El método para resolver este tipo de ecuaciones se conoce como método de coeficientes indeterminados.

Dentro de este método existen al menos dos caminos desde una perspectiva diferente, uno de ellos es conocido como método de superposición y se desarrolla desde el punto de vista del principio de superposición para ecuaciones no homogéneas, mientras que un segundo método utiliza el concepto de operadores diferenciales anuladores y se conoce como método del anulador. Éste segundo método no lo estudiaremos en este curso pero puedes investigar acerca de él si lo deseas.

Coeficientes indeterminados: método de superposición

La idea fundamental de este método consiste en hacer una conjetura acerca de la forma que debe tener la solución particular $y_{p}(x)$, esta intuición esta motivada por el tipo de funciones que constituyen a la función $g(x)$ de (\ref{3}). Debido a que la solución de cierta manera depende de la forma de $g(x)$, es que este método se ve limitado a cierto tipo de funciones que componen a $g(x)$.

Como lo hemos venido haciendo, desarrollaremos esta teoría para el caso de ecuaciones de segundo orden ya que al aumentar el orden aumenta la complejidad en los cálculos, pero la idea se mantiene.

La ecuación a resolver es (\ref{2}), esto es

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = g(x)$$

Este método se puede aplicar sólo en el caso el que se cumple lo siguiente:

  • Los coeficientes $a$, $b$ y $c$ son constantes y,
  • $g(x)$ es una constante $c$, una función polinomial, una función exponencial $e^{\alpha x}$, una función seno o coseno $(\sin(\beta x)$ o $\cos(\beta x))$ o sumas infinitas y productos de estas funciones.

De acuerdo al segundo punto, este método no es aplicable a funciones como

$$g(x) = \ln(x); \hspace{1cm} g(x) = \dfrac{1}{x}; \hspace{1cm} g(x) = \tan(x); \hspace{1cm} g(x) = \arcsin(x)$$

etcétera. En cambio, es aplicable a funciones como

$$g(x) = P_{n}(x) = a_{n}x^{n} + a_{n -1}x^{n -1} + \cdots + a_{1}x + a_{0}$$

$$g(x) = P_{n}(x) e^{\alpha x}; \hspace{1cm} g(x) = P_{n}(x) e^{\alpha x} \sin(\beta x); \hspace{1cm} y \hspace{1cm} g(x) = P_{n}(x) e^{\alpha x} \cos(\beta x)$$

donde $n$ es un número entero no negativo que indica el grado del polinomio y $\alpha$ y $\beta$ son números reales.

Lo importante a considerar en este método es que el conjunto de funciones que consiste en constantes, polinomios, exponenciales, senos y cosenos tiene la notable propiedad de que las derivadas de sus sumas y productos vuelven a ser funciones de este mismo estilo y debido a que la combinación lineal $a \dfrac{d^{2}y_{p}}{dx^{2}} + b \dfrac{dy_{p}}{dx} + cy_{p}$ debe ser idéntica a la función $g(x)$ es razonable suponer que la solución particular $y_{p}(x)$ tiene la misma forma que $g(x)$.

Ahora que hemos motivado la noción de este método desarrollemos la teoría que nos permite determinar las soluciones a ecuaciones diferenciales lineales de segundo orden no homogéneas con coeficientes constantes (\ref{2}).

Comencemos con el caso en el que la función $g(x)$ es un polinomio.

$g(x)$ es un polinomio

Supongamos que la función $g(x)$ es un polinomio de grado $n$ de la forma

$$g(x) = a_{n}x^{n} + a_{n -1}x^{n -1} + \cdots + a_{1}x + a_{0} = \sum_{k = 0}^{n}a_{k}x^{k} \label{4} \tag{4}$$

con $a_{k}$, $k = 0, 1, 2, \cdots, n$ constantes.

Buscamos la solución particular $y_{p}(x)$ tal que

$$a \dfrac{d^{2}y_{p}}{dx^{2}} + b \dfrac{dy_{p}}{dx} + cy_{p} = \sum_{k=0}^{n}a_{k}x^{k} \label{5} \tag{5}$$

Como mencionamos anteriormente, vamos a suponer que la solución particular tiene la misma forma que $g(x)$ así, en este caso, $y_{p}(x)$ será también un polinomio de grado $n$.

$$y_{p}(x) = A_{n}x^{n} + A_{n -1}x^{n -1} + \cdots + A_{1}x + A_{0} = \sum_{k = 0}^{n}A_{k}x^{k} \label{6} \tag{6}$$

con $A_{k}$, $k = 0, 1, 2, \cdots, n$ constantes. Calculemos la primera y segunda derivada y sustituyamos en la ecuación (\ref{5}).

$$\dfrac{dy_{p}}{dx} = \sum_{k = 1}^{n} kA_{k}x^{k -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{p}}{dx^{2}} = \sum_{k = 2}^{n}k(k -1)A_{k}x^{k -2}$$

Si desglosas las sumas notarás la razón y la importancia de que el valor de $k$ comience en $1$ y en $2$ en la suma de la primera y la segunda derivada respectivamente. Más adelante trabajaremos con series y mencionaremos con mayor detalle este tipo de situaciones.

Sustituyendo en la ecuación diferencial (\ref{5}) se tiene lo siguiente:

$$a \sum_{k = 2}^{n} k(k -1)A_{k}x^{k -2} + b \sum_{k = 1}^{n}k A_{k}x^{k -1} + c \sum_{k = 0}^{n}A_{k}x^{k} = \sum_{k = 0}^{n}a_{k}x^{k} \label{7} \tag{7}$$

Para que se cumpla la igualdad los coeficientes de cada $x^{k}$, $k = 0, 1, \cdots, n$, tienen que ser iguales. Para que resulte más intuitivo el cómo igualar los coeficientes vamos a reescribir la ecuación (\ref{7}) desglosando las sumas.

\begin{align*}
a \left[ 2A_{2} + 6A_{3}x + 12A_{4}x^{2} + \cdots + n(n -1)A_{n}x^{n -2} \right] \\
+ b \left[ A_{1} + 2A_{2}x + 3A_{3}x^{2} + \cdots + (n -1)A_{n -1}x^{n -2} + n A_{n}x^{n -1} \right] \\
+ c \left[ A_{0} + A_{1}x + A_{2}x^{2} + \cdots + A_{n -2}x^{n -2} + A_{n -1}x^{n -1} + A_{n}x^{n} \right] \\
= a_{0} + a_{1}x + a_{2}x^{2} + \cdots + a_{n -2}x^{n -2} + a_{n -1}x^{n -1} + a_{n}x^{n}
\end{align*}

Ahora si resulta directo igualar los coeficientes de cada $x^{k}$, $k = 0, 1, \cdots, n$. Para que se cumpla la ecuación (\ref{7}) se deben cumplir las siguientes relaciones:

\begin{align*}
2aA_{2} + bA_{1} + cA_{0} &= a_{0} \\
6aA_{3} + 2bA_{2} + cA_{1} &= a_{1} \\
12aA_{4} + 3bA_{3} + cA_{2} &= a_{2} \\
\vdots \\
n(n -1)aA_{n} + (n -1)bA_{n -1} + cA_{n -2} &= a_{n -2} \\
nbA_{n} +c A_{n -1} &= a_{n -1} \\
c A_{n} &= a_{n}
\end{align*}

Lo que tenemos es un sistema de $n +1$ ecuaciones. Si $c \neq 0$, de la última relación se obtiene que $A_{n} = \dfrac{a_{n}}{c}$. De la penultima se obtiene que $A_{n -1} = \dfrac{a_{n-1} -nbA_{n}}{c}$ en donde ya se conoce el valor de $A_{n}$, y así sucesivamente de manera que seremos capaces de determinar todos los valores $A_{k}$, $k = 0, 1, \cdots, n$ y por tanto, obtendremos la solución particular $y_{p}(x)$ dada en (\ref{6}).

El paso fundamental en este método es suponer que la solución particular $y_{p}(x)$ tiene la misma forma que la función $g(x)$, pero debemos ser cuidadosos con ello, por ejemplo, si ocurre que $c = 0$ y suponemos que $y_{p}(x)$ es un polinomio de grado $n$, entonces la expresión $a\dfrac{d^{2}y_{p}}{dx^{2}}+b\dfrac{dy_{p}}{dx}$ será una expresión de grado $n -1$, pero estará igualada a $g(x)$ que es una expresión de grado $n$ y esto no debe ocurrir. Para solucionar esto será necesario considerar a $y_{p}(x)$ como un polinomio de grado $n + 1$, esto es

$$y_{p}(x) = x (A_{n}x^{n} + A_{n -1}x^{n -1} + \cdots + A_{1}x + A_{0}) \label{8} \tag{8}$$

de esta manera la expresión $a\dfrac{d^{2}y_{p}}{dx^{2}}+b\dfrac{dy_{p}}{dx}$ sí será de grado $n$ y por tanto podremos igualar coeficientes con $g(x)$ y así obtener los valores $A_{k}$, $k = 1, 2 \cdots, n$, y con ello la solución particular $y_{p}(x)$.

Finalmente, si $b = c = 0$ tendremos la ecuación $a\dfrac{d^{2}y_{p}}{dx^{2}} = g(x)$ la cual se puede resolver integrando un par de veces con respecto a $x$.

Realicemos un ejemplo.

Ejemplo: Determinar la solución particular $y_{p}(x)$ de la ecuación diferencial $9\dfrac{d^{2}y}{dx^{2}} -6\dfrac{dy}{dx} + y = 9 -x^{3}$

Solución: En este caso la función $g(x) = 9 -x^{3}$ corresponde a un polinomio de grado $n = 3$ y como $c = 1 \neq 0$ vamos a suponer que la solución particular es también un polinomio de grado $3$. Importante, aunque en la función $g$ no estén los términos con $x$ y $x^{2}$ en la solución sí debemos considerarlos. De acuerdo a (\ref{6}), la solución particular tiene la siguiente forma:

$$y_{p}(x) = A_{3}x^{3} + A_{2}x^{2} + A_{1}x + A_{0}$$

Vamos a calcular la primera y segunda derivada para sustituir en la ecuación diferencial.

$$\dfrac{dy_{p}}{dx} = 3A_{3}x^{2} + 2A_{2}x + A_{1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{p}}{dx^{2}} = 6A_{3}x + 2A_{2}$$

Sustituyendo:

$$(54A_{3}x + 18A_{2}) -(18A_{3}x^{2} + 12A_{2}x + 6A_{1}) + (A_{3}x^{3} + A_{2}x^{2} + A_{1}x + A_{0}) = 9 -x^{3}$$

Agrupando términos semejantes:

$$A_{3}x^{3} + (-18A_{3} + A_{2})x^{2} + (54A_{3} -12A_{2} + A_{1})x + (18A_{2} -6A_{1} + A_{0}) = 9 -x^{3}$$

La función $g$ la podemos escribir como $g(x) = -1x^{3} + 0x^{2} + 0x + 9$ para que podamos igualar los coeficientes de cada $x^{k}$, $k = 0, 1, 2, 3$.

\begin{align*}
A_{3} &= -1\\
-18A_{3} + A_{2} &= 0 \\
54A_{3} -12A_{2} + A_{1} &= 0 \\
18A_{2} -6A_{1} + A_{0} &= 9
\end{align*}

Resolviendo el sistema de ecuaciones obtenemos que $A_{3} = -1$, $A_{2} = -18$, $A_{1} = -162$ y $A_{0} = -639$. Por lo tanto, la solución particular de la ecuación diferencial es

$$y_{p}(x) = -x^{3} -18x^{2} -162x -639$$

$\square$

Puedes comprobar que efectivamente es solución de la ecuación diferencial. Por otro lado, nota que en el ejemplo hemos obtenido sólo la solución particular, aplicando lo visto en la entrada anterior, de tarea moral determina la solución complementaria $y_{c}(x)$ que corresponde a la solución de la ecuación homogénea asociada para así finalmente obtener la solución general $y(x) = y_{c}(x) + y_{p}(x)$.

Ahora veamos que ocurre cuando la función $g(x)$ es una función exponencial o el producto de un polinomio con una función exponencial.

$g(x)$ es producto de un polinomio y una función exponencial

Antes de desarrollar la teoría para este caso veamos una situación que nos muestra que no siempre la solución particular debe tener la misma forma que la función $g(x)$.

Supongamos que queremos obtener la solución particular de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -5\dfrac{dy}{dx} + 4y = 8e^{x}$$

Se puede suponer razonablemente una solución particular de la forma $y_{p}(x) = Ae^{x}$. Con $A$ una constante por determinar. Si derivamos obtenemos que $\dfrac{dy_{p}}{dx} = Ae^{x} = \dfrac{d^{2}y_{p}}{dx^{2}}$. Al sustituir en la ecuación diferencial da como resultado la expresión contradictoria $0 = 8e^{x}$, esto nos muestra que se hizo una conjetura equivocada para $y_{p}$. Lo que sucede es que la solución complementaria $y_{c}(x) = c_{1}e^{x} + c_{2}e^{4x}$ ya contiene nuestra suposición $Ae^{x}$. Esto significa que $e^{x}$ es una solución de la ecuación diferencial homogénea asociada y un múltiplo constante $Ae^{x}$ cuando se sustituye en la ecuación diferencial necesariamente da cero.

Como lo vemos en este ejemplo, la solución particular que propongamos no debe ser una función que ya este contenida en la solución complementaria. A continuación desarrollaremos un método que nos permite obtener la forma de la solución particular cuando la función $g(x)$ esta constituida por un polinomio y una función exponencial.

Supongamos que la función $g(x)$ es el producto de un polinomio con una función exponencial:

$$g(x) = \left( a_{n}x^{n} + a_{n -1}x^{n -1} + \cdots + a_{1}x + a_{0} \right) e^{\alpha x} = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) e^{\alpha x} \label{9} \tag{9}$$

Con $\alpha$ una constante distinta de cero ya que en caso de serlo regresamos al caso anterior. En esta ocasión la ecuación diferencial a resolver es

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) e^{\alpha x} \label{10} \tag{10}$$

Manteniendo la esencia del método, vamos a suponer que la solución particular tiene la forma de la función $g(x)$. Supongamos que la solución particular es de la forma $y_{p}(x) = e^{\alpha x} h(x)$, con $h(x)$ una función continua y derivable en el intervalo de solución. Lo que haremos será determinar la forma de la función $h(x)$, para ello vamos a calcular las derivadas correspondientes y sustituir en la ecuación diferencial (\ref{10}).

$$\dfrac{dy_{p}}{dx} = e^{\alpha x} \dfrac{dh}{dx} + \alpha e^{\alpha x}h \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{p}}{dx^{2}} = e^{\alpha x} \dfrac{d^{2}h}{dx^{2}} + 2 \alpha e^{\alpha x} \dfrac{dh}{dx} + \alpha^{2}e^{\alpha x}h$$

Sustituyendo en la ecuación diferencial.

$$a \left[ e^{\alpha x} \dfrac{d^{2}h}{dx^{2}} + 2 \alpha e^{\alpha x} \dfrac{dh}{dx} + \alpha^{2}e^{\alpha x}h \right] + b \left[ e^{\alpha x} \dfrac{dh}{dx} + \alpha e^{\alpha x}h \right] + c\left( e^{\alpha x}h \right) = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) e^{\alpha x}$$

Como $e^{\alpha x} \neq 0$ $\forall x \in \mathbb{R}$, podemos cancelar la exponencial de todos los términos y reordenar para obtener la siguiente expresión:

$$a \dfrac{d^{2}h}{dx^{2}} + \left( 2 \alpha a + b \right) \dfrac{dh}{dx} + \left( a \alpha^{2} + b \alpha + c \right) h = \sum_{k = 0}^{n}a_{k}x^{k} \label{11} \tag{11}$$

Lo que obtuvimos es una ecuación diferencial de segundo orden con coeficientes constantes en donde el lado derecho de la ecuación corresponde a una función polinomial, es decir, hemos reducido el problema al primer caso. Podemos aplicar lo visto anteriormente para obtener la forma de la solución $h(x)$. Consideremos los siguientes casos:

  • $a \alpha^{2} + b \alpha + c \neq 0$

Al considerar este caso lo que podemos proponer es que la función $h(x)$ sea un polinomio de grado $n$, esto es, $h(x) = A_{n}x^{n} + A_{n -1}x^{n -1} + \cdots + A_{1}x + A_{0} = \sum_{k = 0}^{n}A_{k}x^{k}$, de esta manera ambos lados de la ecuación (\ref{11}) serán expresiones de grado $n$ y podremos igualar coeficiente a coeficiente tal como lo hicimos en el caso anterior. Por lo tanto, la solución particular a la ecuación diferencial (\ref{10}) es

$$y_{p}(x) = \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) e^{\alpha x} \label{12} \tag{12}$$

  • $a \alpha^{2} + b \alpha + c = 0$ y $2 \alpha a + b \neq 0$

En este caso, por el mismo argumento que en el caso en el que $g(x)$ es sólo un polinomio de grado $n$, debemos asegurarnos de que la combinación lineal $a \dfrac{d^{2}h}{dx^{2}} + \left( 2 \alpha a + b \right) \dfrac{dh}{dx}$ sea una expresión de grado $n$ por lo que será necesario que $h(x)$ sea un polinomio de grado $n + 1$. La propuesta para $h(x)$ es que sea de la forma $h(x) = x \left( A_{n}x^{n} + A_{n -1}x^{n -1} + \cdots + A_{1}x + A_{0} \right) = x \left( \sum_{k = 0}^{n}A_{k}x^{k} \right)$. Por lo tanto, para este caso la solución particular es la siguiente:

$$y_{p}(x) = x \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) e^{\alpha x} \label{13} \tag{13}$$

  • $a \alpha^{2} + b \alpha + c = 0$ y $2\alpha a + b = 0$

Siguiendo el mismo razonamiento, en este caso la solución particular debe ser:

$$y_{p}(x) = x^{2} \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) e^{\alpha x} \label{14} \tag{14}$$

Ahora que conocemos las posibles formas de las soluciones particulares, vamos a resolver la ecuación que nos genero conflicto.

Ejemplo: Hallar la solución particular $y_{p}(x)$ de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -5\dfrac{dy}{dx} + 4y = 8e^{x}$

Solución: Vamos a suponer que la solución tiene la forma $y_{p}(x) = h(x)e^{x}$, con $\alpha = 1$. Las derivadas correspondientes son

$$\dfrac{dy_{p}}{dx} = \dfrac{dh}{dx}e^{x} + he^{x} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{p}}{dx^{2}} = \dfrac{d^{2}h}{dx^{2}} e^{x} + 2\dfrac{dh}{dx} e^{x} + he^{x}$$

Sustituyendo obtenemos

$$\left( \dfrac{d^{2}h}{dx^{2}} e^{x} + 2\dfrac{dh}{dx} e^{x} + he^{x} \right) -5 \left( \dfrac{dh}{dx}e^{x} + he^{x} \right) + 4he^{x} = 8e^{x}$$

Cancelando la exponencial y reordenando los términos obtenemos la siguiente ecuación para $h(x)$

$$\dfrac{d^{2}h}{dx^{2}} -3\dfrac{dh}{dx} = 8$$

Como no aparece explícitamente la función $h(x)$ pero sí su primer derivada esto nos obliga a proponer como solución un polinomio de grado $1$ ya que en el lado derecho de la ecuación tenemos un polinomio de grado cero (una constante). Sea $h(x) = Ax$ la función propuesta. Derivando obtenemos $\dfrac{dh}{dx} = A$ y $\dfrac{d^{2}h}{dx^{2}} = 0$. Sustituyendo $-3A = 8$ de donde $A = -\dfrac{8}{3}$. Por lo tanto, la solución particular de la ecuación dada es

$$y_{p}(x) = -\dfrac{8}{3}xe^{x}$$

$\square$

Realicemos un ejemplo más.

Ejemplo: Hallar la solución particular $y_{p}(x)$ de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} -6y = -5e^{2x}$.

Solución: En este caso $g(x) = -5e^{2x}$ y podemos identificar que $a = 1, b = 1, c = -6$ y $\alpha = 2$. Podemos directamente sustituir en la ecuación (\ref{11}) pero es conveniente no memorizarla sino realizar el procedimiento que nos permite llegar a ella. Proponemos la solución $y_{p}(x) = h(x) e^{2x}$, la primera y segunda derivada están dadas como

$$\dfrac{dy_{p}}{dx} = e^{2x}\dfrac{dh}{dx} + 2e^{2x}h \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = e^{2x}\dfrac{d^{2}h}{dx^{2}} + 4e^{2x}\dfrac{dh}{dx} +4e^{2x}h$$

Sustituyendo en la ecuación diferencial:

$$\left( e^{2x}\dfrac{d^{2}h}{dx^{2}} + 4e^{2x}\dfrac{dh}{dx} +4e^{2x}h \right) + \left( e^{2x}\dfrac{dh}{dx} + 2e^{2x}h \right) -6e^{2x}h = -5e^{2x}$$

Cancelando la exponencial de todos los términos y reordenando obtenemos la expresión correspondiente a (\ref{11})

$$\dfrac{d^{2}h}{dx^{2}} + 5\dfrac{dh}{dx} = -5$$

Notamos que el lado derecho de la ecuación es una constante que, estrictamente hablando, corresponde a un polinomio de grado cero, por otro lado, en la ecuación diferencial no aparece explícitamente la función $h$ pero sí sus derivadas, en particular la primer derivada, esto nos indica que $h$ debe ser un polinomio de grado $1$ para que efectivamente la expresión $\dfrac{d^{2}h}{dx^{2}} + 5\dfrac{dh}{dx}$ sea de grado cero. Tomando en cuenta esto podemos establecer que la forma de $h$ debe ser $h(x) = Ax$, así la solución particular será $y_{p}(x) = Axe^{2x}$. Si regresamos a la teoría desarrollada te darás cuenta que la solución particular tiene justamente la forma de la ecuación (\ref{13}) para $n = 0$ ya que estamos en el caso en el que $a \alpha^{2} + b \alpha + c = 0$ y $2 \alpha a + b \neq 0$.

Determinemos el valor del coeficiente $A$ considerando la función $h(x) = Ax$, calculando las derivadas se tiene que

$$\dfrac{dh}{dx} = A \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}h}{dx^{2}} = 0$$

Sustituyendo en la ecuación diferencial de $h$ tenemos que $0 + 5A = -5$, de donde $A = -1$, así $h(x) = -x$ y por lo tanto, la solución particular de la ecuación diferencial es

$$y_{p}(x) = -xe^{2x}$$

$\square$

De tarea moral determina la solución complementaria del ejemplo para así obtener la solución general $y(x) = y_{c}(x) + y_{p}(x)$.

Finalmente veamos que ocurre si la función $g(x)$ es una función seno, coseno o el producto de un polinomio con alguna de estas funciones trigonométricas.

Caso en el que g(x) es producto de un polinomio y una función seno o coseno

Consideremos nuevamente el polinomio

$$P(x) = a_{n}x^{n} + a_{n -1}x^{n -1} + \cdots + a_{1}x + a_{0} = \sum_{k = 0}^{n}a_{k}x^{k}$$

El tipo de ecuación que queremos resolver es de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b\dfrac{dy}{dx} + cy = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \cos(\beta x) \label{15} \tag{15}$$

o bien

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \sin(\beta x) \label{16} \tag{16}$$

Con $\beta$ una constante. Para resolver este tipo de ecuaciones es conveniente considerar la siguiente propiedad

$$e^{i\beta x} = \cos (\beta x) + i \sin(\beta x) \label{17} \tag{17}$$

Con $i^{2} = -1$. Y consideremos la ecuación diferencial

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) e^{i \beta x} \label{18} \tag{18}$$

Que podemos escribir como

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \cos (\beta x) + i \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \sin(\beta x) \label{19} \tag{19}$$

Supongamos que la función $z(x) = u(x) + i v(x)$ es solución de la ecuación diferencial (\ref{19}). Consideremos el operador polinomial $\mathcal{L} = a \dfrac{d^{2}}{dx^{2}} + b\dfrac{d}{dx} + c$, de manera que podemos reescribir a la ecuación como

$$\mathcal{L}[z] = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \cos(\beta x) + i \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \sin(\beta x) \label{20} \tag{20}$$

pero notemos que

$$\mathcal{L}[z] = \mathcal{L}[u + iv] = \mathcal{L}[u] + i \mathcal{L}[v] \label{21} \tag{21}$$

Comparando las ecuaciones (\ref{20}) y (\ref{21}) deducimos que

$$\mathcal{L}[u] = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \cos(\beta x) \label{22} \tag{22}$$

y

$$\mathcal{L}[v] = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \sin(\beta x) \label{23} \tag{23}$$

Es decir, la parte real $Re(z) = u(x)$ es solución de la ecuación

$$a \dfrac{du^{2}}{dx^{2}} + b \dfrac{du}{dx} + cu = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \cos(\beta x) \label{24} \tag{24}$$

y la parte imaginaria $Im(z) = v(x)$ es solución de la ecuación

$$a \dfrac{dv^{2}}{dx^{2}} + b \dfrac{dv}{dx} + cv = \left( \sum_{k = 0}^{n}a_{k}x^{k} \right) \sin(\beta x) \label{25} \tag{25}$$

Realicemos un par de ejemplos en el que veamos cómo obtener las funciones $u(x)$ o $v(x)$ de la solución $z(x) = u(x) + iv(x)$.

Ejemplo: Determinar la solución particular $y_{p}(x)$ de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -\dfrac{3}{2}\dfrac{dy}{dx} -y = 3 \cos(x)$.

Solución: La ecuación que queremos resolver es de la forma (\ref{15}), pero es conveniente sustituir la función coseno por la función compleja $e^{ix}$ para resolver la ecuación en la forma (\ref{18}) con $\beta = 1$. Haciendo la sustitución, la ecuación que intentaremos resolver será

$$\dfrac{d^{2}y}{dx^{2}} -\dfrac{3}{2}\dfrac{dy}{dx} -y = 3 e^{ix}$$

Ahora tenemos una ecuación en la que el lado derecho corresponde al producto de un polinomio de grado cero con una función exponencial, esto nos indica que podemos suponer una solución de la forma $z(x) = h(x)e^{ix}$. Derivando dos veces tenemos que

$$\dfrac{dz}{dx} = e^{ix}\dfrac{dh}{dx} + ie^{ix}h \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}z}{dx^{2}} = e^{ix}\dfrac{d^{2}h}{dx^{2}} + 2ie^{ix}\dfrac{dh}{dx} -e^{ix}h$$

Sustituyendo en la ecuación diferencial:

$$\left( e^{ix}\dfrac{d^{2}h}{dx^{2}} + 2ie^{ix}\dfrac{dh}{dx} -e^{ix}h \right) -\dfrac{3}{2} \left( e^{ix}\dfrac{dh}{dx} + ie^{ix}h \right) -\left( he^{ix} \right) = 3 e^{ix}$$

Cancelando las exponenciales de todos los términos y reordenando obtenemos:

$$\dfrac{d^{2}h}{dx^{2}} + \left( 2i -\dfrac{3}{2} \right) \dfrac{dh}{dx} + \left( -2 -\dfrac{3}{2}i \right)h = 3$$

Es claro que la forma de la función $h$ debe ser la de un polinomio de grado cero, ya que es el tipo de función que tenemos en el lado derecho de la ecuación. Digamos que $h(x) = A$, ambas derivadas son cero por lo que al sustituir en la ecuación diferencial para $h$ obtenemos que $\left( -2 -\dfrac{3}{2}i \right)A = 3$. Para despejar a la constante $A$ vamos a multiplicar ambos lados de la ecuación por el número complejo conjugado $\left( -2 +\dfrac{3}{2}i \right)$:

\begin{align*}
\left( -2 +\dfrac{3}{2}i \right) \left( -2 -\dfrac{3}{2}i \right)A &= \left( -2 +\dfrac{3}{2}i \right)3 \\
\dfrac{25}{4}A &= -6 +\dfrac{9}{2}i \\
A &= -\dfrac{24}{25} + \dfrac{18}{25}i
\end{align*}

Como $h(x) = A$, entonces $h(x) = -\dfrac{24}{25} + \dfrac{18}{25}i$, sustituyendo en la solución se tiene que

\begin{align*}
z(x) &= \left( -\dfrac{24}{25} + \dfrac{18}{25}i \right) e^{ix} \\
&= \left( -\dfrac{24}{25} + \dfrac{18}{25}i \right) \left( \cos(x) + i \sin(x) \right) \\
&= \left( -\dfrac{24}{25} \cos(x) -\dfrac{18}{25} \sin(x) \right) + i \left( -\dfrac{24}{25} \sin(x) + \dfrac{18}{25} \cos(x) \right)
\end{align*}

Si $z(x) = u(x) +iv(x)$, notamos que

$$Re(z) = u(x) = -\dfrac{24}{25} \cos(x) -\dfrac{18}{25} \sin(x)$$

y

$$Im(z) = v(x) = -\dfrac{24}{25} \sin(x) + \dfrac{18}{25} \cos(x)$$

Ya que en la ecuación original esta presente la función coseno, de acuerdo a $(\ref{24})$, la solución de la ecuación será la función $u(x)$, es decir, la parte real de la función $z(x)$. Por lo tanto, la solución particular es

$$y_{p}(x) = -\dfrac{24}{25} \cos(x) -\dfrac{18}{25} \sin(x)$$

$\square$

Realicemos un ejemplo más.

Ejemplo: Hallar la solución particular $y_{p}(x)$ de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 4y = 12 \sin(2x)$

Solución: La ecuación a resolver es de la forma (\ref{16}) pero es conveniente sustituir la función seno por la función compleja $e^{2ix}$ para resolver la ecuación en la forma (\ref{18}) con $\beta = 2$. Haciendo la sustitución, la ecuación que intentaremos resolver será

$$\dfrac{d^{2}y}{dx^{2}} + 4y = 12 e^{2ix}$$

Como tenemos una ecuación en la que el lado derecho corresponde al producto de un polinomio de grado cero con una función exponencial, entonces podemos suponer una solución de la forma $z(x) = h(x)e^{2ix}$. Derivando dos veces tenemos que

$$\dfrac{dz}{dx} = e^{2ix}\dfrac{dh}{dx} + 2ie^{2ix}h \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}z}{dx^{2}} = e^{2ix}\dfrac{d^{2}h}{dx^{2}} + 4ie^{2ix}\dfrac{dh}{dx} -4e^{2ix}h$$

Sustituyendo en la ecuación diferencial:

$$\left( e^{2ix}\dfrac{d^{2}h}{dx^{2}} + 4ie^{2ix}\dfrac{dh}{dx} -4e^{2ix}h \right) + 4\left( he^{2ix} \right) = 12 e^{2ix}$$

Cancelando las exponenciales de todos los términos y reordenando obtenemos:

$$\dfrac{d^{2}h}{dx^{2}} + 4i \dfrac{dh}{dx} = 12$$

Es claro que la forma de la función $h$ debe ser la de un polinomio de grado $1$, pues dicha función no aparece explícitamente en el lado izquierdo de la ecuación pero sí su primer derivada que es la que determina el grado más alto de la ecuación y queremos que sea de grado cero ya que es el grado de la función que tenemos en el lado derecho. Supongamos que $h(x) = Ax$, su primer derivada es $\dfrac{dh}{dx} = A$ y su segunda derivada es cero, al sustituir en la ecuación diferencial para $h$ obtenemos que $\left( 4i \right)A = 12$, de donde $A = -3i$. Entonces $h(x) = -3ix$ y por tanto $z(x) = -3ixe^{2ix}$. Vemos que

\begin{align*}
z(x) &= \left( -3ix \right) \left( \cos(2x) + i \sin(2x) \right) \\
&= \left( 3x \sin(2x) \right) + i \left(-3x \cos(2x) \right)
\end{align*}

Si $z(x) = u(x) + i v(x)$, notamos que $u(x) = 3x \sin(2x)$ y $v(x) = -3x \cos(2x)$. Debido a que en la ecuación original esta presente la función seno, de acuerdo a $(\ref{25})$, la solución de la ecuación será la función $v(x)$, es decir, la parte imaginaria de la función $z(x)$. Por lo tanto, la solución particular es

$$y_{p}(x) = -3x \cos(2x)$$

$\square$

Como pudiste notar, éste método de resolución requiere, en buena medida, de nuestra intuición para proponer la forma correcta de la solución. A lo largo de esta entrada hemos desarrollado la teoría que nos permite proponer una solución adecuada, sin embargo existen distintas variantes de opciones para la función $g(x)$ y por tanto para la solución particular $y_{p}(x)$. A continuación mostraremos los resultados generales que nos permiten proponer una solución particular dada una función $g(x)$.

Resultados generales

Como mencionamos al inicio de esta entrada, en la ecuación (\ref{2}) la función $g(x)$ puede ser un polinomio, una exponencial, una función seno o coseno o el producto entre estas funciones. De manera general diremos que la función $g(x)$ puede tener la siguiente forma:

$$g(x) = e^{\alpha x} \left[ P_{n}(x) \cos(\beta x) + Q_{m}(x) \sin(\beta x) \right] \label{26} \tag{26}$$

Donde $\lambda = \alpha \pm i\beta$ es raíz de la ecuación auxiliar y $P_{n}(x)$ y $Q_{m}(x)$ son los polinomios en $x$ de grado $n$ y $m$ respectivamente.

$$P_{n}(x) = \sum_{k = 0}^{n}a_{k}x^{k} \hspace{1cm} y \hspace{1cm} Q_{m}(x) = \sum_{k = 0}^{m}b_{k}x^{k}$$

Lo que buscamos es una solución particular de la forma

$$y_{p}(x) = x^{s} e^{\alpha x} \left[ p_{r}(x) \cos(\beta x) + q_{r}(x) \sin(\beta x) \right] \label{27} \tag{27}$$

Donde $r = \max(n, m)$, $p_{r}(x)$ y $q_{r}(x)$ son polinomios en $x$ de grado $r$, cuyos coeficientes están indeterminados, dados como

$$p_{r}(x) = \sum_{k = 0}^{r}A_{k}x^{k} \hspace{1cm} y \hspace{1cm} q_{r}(x) = \sum_{k = 0}^{r}B_{k}x^{k}$$

y $s$ es la multiplicidad de la raíz $\lambda = \alpha \pm i\beta$ de la ecuación auxiliar. La forma de $y_{p}(x)$ se puede resumir de la siguiente manera:

  1. $g(x) = P_{n}(x)$:
  • Si las raíces de la ecuación auxiliar son distintas de cero, $\lambda_{i} \neq 0, i = 1, 2, \cdots, s$, entonces la forma de la solución particular es

$$y_{p}(x) = p_{n}(x) = \sum_{k = 0}^{n}A_{k}x^{k} \label{28} \tag{28}$$

  • Si alguna raíz es igual a cero, $\lambda _{i} = 0$, entonces la forma de la solución particular es

$$y_{p}(x) = x^{s}p_{n}(x) = x^{s} \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) \label{29} \tag{29}$$

  1. $g(x) = P_{n}(x) e^{\alpha x}$:
  • Si $\alpha$ no es raíz, entonces la forma de la solución particular es

$$y_{p}(x) = p_{n}(x) e^{\alpha x} = \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) e^{\alpha x} \label{30} \tag{30}$$

  • Si $\alpha $ es raíz repetida $s$ veces (de orden $s$) entonces la forma de la solución particular es

$$y_{p}(x) = x^{s} p_{n}(x)e^{\alpha x} = x^{s} \left( \sum_{k = 0}^{n}A_{k}x^{k} \right) e^{\alpha x} \label{31} \tag{31}$$

  1. $g(x) = P_{n}(x) \cos(\beta x) + Q_{m}(x) \sin(\beta x)$:
  • Si $\pm i \beta$ no son raíces y $r= \max(n, m)$, entonces la forma de la solución particular es

$$y_{p}(x) = p_{r}(x) \cos(\beta x) + q_{r}(x) \sin(\beta x) \label{32} \tag{32}$$

  • Si $\pm i \beta$ son raíces de orden $s$ y $r = \max(n, m)$, entonces la forma de la solución particular es

$$y_{p}(x) = x^{s} \left[ p_{r}(x) \cos(\beta x) + q_{r}(x) \sin(\beta x) \right] \label{33} \tag{33}$$

  1. $g(x) = e^{\alpha x} \left[ P_{n}(x) \cos(\beta x) + Q_{m}(x) \sin(\beta x) \right]$:
  • Si $\alpha \pm i\beta$ no son raíces, entonces la forma de la solución particular es

$$y_{p}(x) = e^{\alpha x} \left[p_{r}(x) \cos(\beta x) + q_{r}(x) \sin(\beta x) \right] \label{34} \tag{34}$$

  • Si $\alpha \pm i\beta$ son raíces de orden $s$, entonces la forma de la solución particular es

$$y_{p}(x) = x^{s}e^{\alpha x} \left[ p_{r}(x) \cos(\beta x) + q_{r}(x) \sin(\beta x) \right] \label{35} \tag{35}$$

Este método es conocido como método de superposición y la razón es porque los resultados anteriores se aplican a cada sumando de la función $g(x)$. Veamos un ejemplo.

Ejemplo: Determinar la forma de la solución particular de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -9 \dfrac{dy}{dx} +14y = 3x^{2} -5 \sin(2x) + 7xe^{6x}$

Solución: La función $g$ es $g(x) = 3x^{2} -5 \sin(2x) + 7xe^{6x}$. Al primer sumando $g_{1}(x) = 3x^{2}$ le corresponde una solución de la forma

$$y_{p1}(x) = A_{2}x^{2} + A_{1}x + A_{0}$$

Al segundo sumando $g_{2}(x) = -5 \sin(2x)$ le corresponde una solución de la forma

$$y_{p2}(x) = B_{1} \cos(2x) + B_{0} \sin(2x)$$

Y al tercer sumando $g_{3}(x) = 7xe^{6x}$ le corresponde una solución de la forma

$$y_{p3}(x) = (C_{1}x + C_{0})e^{6x}$$

La solución particular estará dada por la superposición de las tres soluciones correspondientes a los tres sumandos de la función $g(x)$:

\begin{align*}
y_{p}(x) &= y_{p_{1}}(x) + y_{p_{2}}(x) + y_{p_{3}}(x) \\
&= A_{2}x^{2} + A_{1}x + A_{0} + B_{1} \cos(2x) + B_{0} \sin(2x) + (C_{1}x + C_{0})e^{6x}
\end{align*}

En esta suposición ningún término duplica un término de la función complementaria $y_{c}(x) = c_{1}e^{2x} + c_{2}e^{7x}$.

De tarea moral determina el valor de los coeficientes para formar la solución general de la ecuación diferencial.

$\square$

Concluyamos la entrada con algunos ejemplos más.

Ejemplo: Determinar la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 6 \dfrac{dy}{dx} -7y = 3e^{2x} -e^{-x}$.

Solución: La ecuación auxiliar es $k^{2} + 6k -7 = 0$, de donde $k_{1} = 1$ y $k_{2} = -7$. En este caso la multiplicidad es de $s = 1$, es decir no hay raíces repetidas. La solución complementaria es

$$y_{c}(x) = c_{1}e^{x} + c_{2}e^{-7x}$$

La solución particular tiene la forma $y_{p}(x) = A_{0}e^{2x} + B_{0}e^{-x}$, derivando obtenemos

$$\dfrac{dy_{p}}{dx} = 2A_{0}e^{2x} -B_{0}e^{-x} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{p}}{dx^{2}} = 4A_{0}e^{2x} + B_{0}e^{-x}$$

Sustituyendo en la ecuación:

$$(4A_{0}e^{2x} + B_{0}e^{-x}) + 6(2A_{0}e^{2x} -B_{0}e^{-x}) -7(A_{0}e^{2x} + B_{0}e^{-x}) = 3e^{2x} -e^{-x}$$

Reordenando obtenemos

$$9A_{0}e^{2x} -12B_{0}e^{-x} = 3e^{2x} -e^{-x}$$

De donde $9A_{0} = 3$ y $-12B_{0} = -1$. Resolviendo obtenemos que $A_{0} = \dfrac{1}{3}$ y $B_{0} = \dfrac{1}{12}$. Por lo tanto la solución particular es

$$y_{p}(x) = \dfrac{1}{3}e^{2x} + \dfrac{1}{12}e^{-x}$$

Así, la solución general a la ecuación diferencial es la superposición de ambas soluciones

$$y(x) = c_{1}e^{x} + c_{2}e^{-7x} + \dfrac{1}{3}e^{2x} + \dfrac{1}{12}e^{-x}$$

$\square$

Ejemplo: Determinar la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -2\dfrac{dy}{dx} = 2e^{2x} + 4 \cos(2x)$

Solución: La ecuación auxiliar es $k^{2} -2k = 0$, de donde $k_{1} = 0$ y $k_{2} = 2$. La solución complementaria es

$$y_{c}(x) = c_{1} + c_{2}e^{2x}$$

Notamos que $\alpha = 2$ es raíz de la ecuación auxiliar, por lo que la solución correspondiente al término $2e^{2x}$ esta dada por (\ref{31}) con multiplicidad $s = 1$, entonces la forma de la solución particular es

$$y_{p}(x) = A_{0}xe^{2x} + B_{1} \cos(2x) + B_{0} \sin (2x)$$

Derivando obtenemos

$$\dfrac{dy_{p}}{dx} = A_{0}e^{2x} + 2A_{0}xe^{2x} -2B_{1}\sin(2x) + 2B_{0} \cos(2x)$$

y

$$\dfrac{d^{2}y_{p}}{dx^{2}} = 2A_{0}e^{2x} + 2A_{0}e^{2x} + 4A_{0}xe^{2x} -4B_{1} \cos(2x) -4B_{0} \sin(2x)$$

Sustituyendo en la ecuación diferencial y reordenando los términos se obtiene la relación

$$2A_{0}e^{2x} + (-4B_{1} -4B_{0}) \cos(2x) + (-4B_{0} + 4B_{1}) \sin(2x) = 2e^{2x} + 4 \cos(2x)$$

De donde

\begin{align*}
2A_{0} &= 2 \\
-4B_{1} -4B_{0} &= 4 \\
-4B_{0} + 4B_{1} &= 0
\end{align*}

Resolviendo obtenemos que $A_{0} = 1$, $B_{1} = -\dfrac{1}{2}$ y $B_{0} = -\dfrac{1}{2}$.

Por lo tanto, la solución particular es

$$y_{p}(x) = xe^{2x} -\dfrac{1}{2} \cos(2x) -\dfrac{1}{2} \sin(2x)$$

Así, la solución general es la superposición de ambas soluciones

$$y(x) = c_{1} + c_{2}e^{2x} + xe^{2x} -\dfrac{1}{2} \cos(2x) -\dfrac{1}{2} \sin(2x)$$

$\square$

Ejemplo: Determinar la solución particular de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -4y = -12e^{-2x} + 15 \cos(x) + 8x$.

Solución: La ecuación auxiliar es $k^{2} -4 = 0$ con $k^{2} = 2$, es decir, $k_{1} = 2$ y $k_{2} = -2$. La solución complementaria es

$$y_{c}(x) = c_{1}e^{2x} + c_{2}e^{-2x}$$

Notamos que $\alpha = -2$ es raíz de la ecuación auxiliar, así que la solución particular es de la forma

$$y_{p}(x) = A_{0}xe^{-2x} + B_{1} \cos(x) + B_{0} \sin(x) + C_{1}x + C_{0}$$

Derivando obtenemos

$$\dfrac{dy_{p}}{dx} = A_{0}e^{-2x} -2A_{0}xe^{-2x} -B_{1} \sin(x) + B_{0} \cos(x) + C_{1}$$

y

$$\dfrac{d^{2}y_{p}}{dx^{2}} = -4A_{0}e^{-2x} + 4A_{0}xe^{-2x} -B_{1} \cos(x) -B_{0} \sin(x)$$

Sustituyendo en la ecuación diferencial y reordenando los términos obtenemos:

$$-4A_{0}e^{-2x} -5B_{1} \cos(x) -5B_{0} \sin(x) -4C_{1}x -4C_{0} = -12e^{-2x} + 15 \cos(x) + 8x$$

De donde

\begin{align*}
-4A_{0} &= -12 \\
-5B_{1} &= 15 \\
-5B_{0} &= 0 \\
-4C_{1} &= 8 \\
-4C_{0} &= 0
\end{align*}

Resolviendo se obtiene que $A_{0} = 3$, $B_{1} = -3$, $B_{0} = 0$, $C_{1} = -2$ y $C_{0} = 0$. Por lo tanto, la solución particular es

$$y_{p}(x) = 3xe^{-2x} -3 \cos(x) -2x$$

Y la solución general es

$$y_{p}(x) = c_{1}e^{2x} + c_{2}e^{-2x} + 3xe^{-2x} -3 \cos(x) -2x$$

$\square$

Tarea moral

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

  1. Hallar la solución general de las siguientes ecuaciones diferenciales.
  • $\dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} -2y = 3e^{2x} -x^{2}$
  • $\dfrac{d^{2}y}{dx^{2}} -2 \dfrac{dy}{dx} + y = 4 \cos(3x) -2 \sin(2x)$
  • $\dfrac{d^{2}y}{dx^{2}} + \dfrac{3}{5} \dfrac{dy}{dx} -\dfrac{1}{5} y = \cos(x) + 5x^{2}$
  • $\dfrac{d^{2}y}{dx^{2}} -3 \dfrac{dy}{dx} -9y = 4 \cos(2x) -5e^{-x}$
  • $\dfrac{d^{2}y}{dx^{2}} -3 \dfrac{dy}{dx} -10y = 50 \cos(5x) -7e^{-2x} + 12e^{x} + 20x$
  1. Resolver el siguiente problema con valor inicial.
  • $\dfrac{d^{2}x}{dt^{2}} + \omega^{2} x = F_{0} \sin(\omega t); \hspace{1cm} x(0) = 0; \hspace{0.4cm} x^{\prime}(0) = 0$.
  1. Considera la ecuación diferencial $a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = e^{\alpha x}$, con $a$, $b$, $c$ y $\alpha$ constantes. La ecuación auxiliar de la ecuación homogénea asociada es $ak^{2} + bk + c = 0$.
  • Si $\alpha$ no es una raíz de la ecuación auxiliar, demuestra que se puede encontrar una solución particular de la forma $y_{p}(x) = Ae^{\alpha x}$, donde $A = \dfrac{1}{a \alpha^{2} + b \alpha + c}$.
  • Si $\alpha$ es una raíz de la ecuación auxiliar de multiplicidad uno, muestra que se puede encontrar una solución particular de la forma $y_{p}(x) = Axe^{\alpha x}$, donde $A = \dfrac{1}{2a \alpha + b}$- Explica cómo se sabe que $\alpha \neq -\dfrac{b}{2a}$.
  • Si $\alpha$ es una raíz de la ecuación auxiliar de multiplicidad dos, demuestra que podemos encontrar una solución particular de la forma $y_{p}(x) = Ax^{2}e^{\alpha x}$, donde $A = \dfrac{1}{2a}$.

Más adelante…

El método desarrollado en esta entrada es bastante práctico para resolver ecuaciones diferenciales de segundo orden no homogéneas con coeficientes constantes, sin embargo es un método limitado sólo a cierto tipo de funciones para $g(x)$.

En la siguiente entrada desarrollaremos un método general que nos permite resolver éste mismo tipo de ecuaciones para cualquier función $g(x)$. Dicho método se conoce como método de variación de parámetros y corresponde a una extension a orden superior del método desarrollado en la unidad anterior.

Entradas relacionadas

Ecuaciones Diferenciales I: Ecuaciones diferenciales homogéneas con coeficientes constantes

Introducción

Continuando con nuestro desarrollo de métodos de resolución de ecuaciones diferenciales de orden superior, en particular de segundo orden, en esta entrada estudiaremos un método de resolución aplicado sólo a ecuaciones diferenciales lineales homogéneas de segundo orden con coeficientes constantes, es decir, de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0$$

Con $a, b$ y $c$ contantes.

Antes de comenzar vamos a motivar el análisis que desarrollaremos a lo largo de la unidad considerando primero las ecuaciones diferenciales de primer orden.

Ecuaciones diferenciales lineales homogéneas de primer orden con coeficientes constantes

En la primer unidad estudiamos las ecuaciones lineales homogéneas de la forma

$$a_{1}(x) \dfrac{dy}{dx} + a_{0}(x) y = 0 \label{1} \tag{1}$$

Si $a_{1}(x) \neq 0$ podemos definir $P(x) = \dfrac{a_{0}(x)}{a_{1}(x)}$ y reescribir la ecuación anterior en su forma canónica como

$$\dfrac{dy}{dx} + P(x) y = 0 \label{2} \tag{2}$$

Cuando estudiamos su método de resolución concluimos que la solución general a esta ecuación es

$$y(x) = k e^{-\int{P(x) dx}} \label{3} \tag{3}$$

Consideremos ahora el caso en el que $a_{1}(x) = a \neq 0$ y $a_{0}(x) = b$ son números constantes, con esto la ecuación (\ref{1}) se puede escribir como

$$a \dfrac{dy}{dx} + b y = 0 \label{4} \tag{4}$$

Si definimos la constante $c = \dfrac{b}{a}$ podremos escribir la ecuación anterior en su forma canónica

$$\dfrac{dy}{dx} + cy = 0 \label{5} \tag{5}$$

En la unidad anterior vimos que esta ecuación se puede resolver ya sea por variables separables o con ayuda de un factor integrante, sin embargo hay un método de resolución que sólo implica hacer un poco de álgebra.

Definiendo $k = -c$ la ecuación (\ref{5}) la podemos escribir de la siguiente manera:

$$\dfrac{dy}{dx} = ky \label{6} \tag{6}$$

La solución se puede intuir rápidamente, buscamos una función $y$ tal que su derivada sea igual a ella misma multiplicada por una constante y la función que satisface esto es

$$y(x) = e^{kx} \label{7} \tag{7}$$

Puedes observar que al derivarla recuperamos la ecuación (\ref{6})

$$\dfrac{dy}{dx} = k e^{kx} = ky$$

Ahora que conocemos la solución, sustituyamos $y = e^{kx}$ y su derivada $\dfrac{dy}{dx} = k e^{kx}$ en la ecuación (\ref{4})

$$a k e^{kx} + b e^{kx} = 0$$

Factorizando la función exponencial tenemos

$$e^{kx} (ak + b) = 0$$

Como $e^{kx} \neq 0$, $\forall x \in \mathbb{R}$ entonces necesariamente

$$ak + b = 0 \label{8} \tag{8}$$

De donde $k = -\dfrac{b}{a}$, sustituyendo en la solución (\ref{7}) tenemos que

$$y(x) = e^{-\frac{b}{a} x}$$

Por lo tanto, la solución general a la ecuación diferencial lineal homogénea de primer orden con coeficientes constantes (\ref{4}) es

$$y(x) = c_{1}e^{-\frac{b}{a} x} \label{9} \tag{9}$$

Ejemplo: Determinar la solución general de la ecuación diferencial $8 \dfrac{dy}{dx} + 16y = 0$.

Solución: Al ser una ecuación con coeficientes constantes, sabemos que la forma de la solución es $y = e^{kx}$. Si sustituimos $y$ y su derivada $\dfrac{dy}{dx}$ en la ecuación diferencial obtenemos que $e^{kx}(8k + 16) = 0$ de donde se debe satisfacer que $8k + 16 = 0$, si despejamos a $k$ se obtiene $k = -\dfrac{16}{8} = -2$, de manera que una solución a la ecuación diferencial es $y = e^{-2x}$, y por tanto su solución general en el intervalo $\delta = (-\infty, \infty)$ es $y(x) = c_{1} e^{-2x}$.

$\square$

Lo interesante es que esta idea de soluciones exponenciales se puede extender a ecuaciones diferenciales homogéneas de orden superior

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n -1} \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{1} \dfrac{dy}{dx} + a_{0} y = 0$$

donde los coeficientes $a_{i}$, $i = 0, 1, 2, …, n$, son constantes reales y $a_{n} \neq 0$.

Este análisis nos sirvió para motivar el siguiente método de resolución de ecuaciones diferenciales de segundo orden.

Ecuaciones diferenciales lineales homogéneas de segundo orden con coeficientes constantes

La ecuación diferencial que queremos resolver es de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0 \label{10} \tag{10}$$

Con $a, b$ y $c$ constantes. Igual que antes, lo que intentamos es encontrar una solución de la forma $y(x) = e^{kx}$, su primera y segunda derivada están dadas por

$$\dfrac{dy}{dx} = k e^{kx} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k^{2} e^{kx}$$

Si sustituimos en la ecuación (\ref{10}), tenemos

\begin{align*}
ak^{2} e^{kx} + bke^{kx} + ce^{kx} &= 0 \\
e^{kx} (ak^{2} + bk + c) &= 0
\end{align*}

Nuevamente $e^{kx} \neq 0, \forall x \in \mathbb{R}$, así que necesariamente

$$ak^{2} + bk + c= 0 \label{11} \tag{11}$$

El problema se ha reducido a obtener las raíces de la ecuación anterior, dicha ecuación tiene un nombre especial.

Definición: La ecuación $ak^{2} + bk + c= 0$ se conoce como la ecuación auxiliar de la ecuación diferencial (\ref{10}).

La dos raíces de la ecuación auxiliar son:

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a}$$

De tus cursos de álgebra recordarás que el discriminate $\Delta = b^{2} -4ac$ puede ser positivo, cero o negativo, en cada caso ocurre lo siguiente:

  • Si $\Delta > 0$ entonces $k_{1}$ y $k_{2}$ son reales y distintos.
  • Si $\Delta = 0$ entonces $k_{1}$ y $k_{2}$ son reales e iguales, y
  • Si $\Delta < 0$ entonces $k_{1}$ y $k_{2}$ son números conjugados complejos.

Estudiemos cada caso y veamos el tipo de solución que se obtiene en cada uno.

Caso 1: Discriminante positivo

Ya vimos que la solución es de la forma $y = e^{kx}$. Si el discriminante es positivo entonces al resolver la ecuación auxiliar obtendremos dos raíces reales y distintas $k_{1}$ y $k_{2}$ de manera que se tendrán dos soluciones a la ecuación diferencial (\ref{10}), dichas soluciones son:

$$y_{1}(x) = e^{k_{1} x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{k_{2} x}$$

Notemos lo siguiente.

\begin{align*}
W(y_{1}, y_{2}) &= \begin{vmatrix}
e^{k_{1} x} & e^{k_{2} x} \\
k_{1} e^{k_{1} x} & k_{2} e^{k_{2} x} \end{vmatrix} \\
&= k_{2} e^{k_{2} x} e^{k_{1} x} -k_{1} e^{k_{1} x} e^{k_{2} x} \\
&= (k_{2} -k_{1}) e^{(k_{1} + k_{2}) x} \neq 0
\end{align*}

Como el Wronskiano es distinto de cero $\forall x \in \mathbb{R}$ esto nos indica que ambas funciones son linealmente independientes en $\mathbb{R}$ y por tanto forman un conjunto fundamental de soluciones de manera que la solución general a la ecuación (\ref{10}), en el caso en el que el discriminante es positivo, es

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} \label{12} \tag{12}$$

Con

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a}$$

Caso 2: Discriminante igual a cero

En este caso como $\Delta = b^{2} -4ac = 0$, entonces $k_{1} = k_{2} = -\dfrac{b}{2a}$, de esta manera sólo obtendremos una solución exponencial

$$y_{1}(x) = e^{k_{1} x} = e^{k_{2} x} = e^{-\frac{b}{2a} x} \label{13} \tag{13}$$

Para obtener la segunda solución vamos a aplicar el método de reducción de orden visto en la entrada anterior, en donde obtuvimos que una segunda solución linealmente independiente es de la forma

$$y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x) dx}}}{y_{1}^{2}(x)} dx} \label{14} \tag{14}$$

En este caso de coeficientes constantes si la ecuación (\ref{10}) la dividimos por la constante $a \neq 0$ obtenemos

$$\dfrac{d^{2}y}{dx^{2}} + \dfrac{b}{a} \dfrac{dy}{dx} + \dfrac{c}{a}y = 0$$

Si definimos las constantes $P = \dfrac{b}{a}$ y $Q = \dfrac{c}{a}$ la ecuación anterior se puede reescribir como

$$\dfrac{d^{2}y}{dx^{2}} + P \dfrac{dy}{dx} + Qy = 0 \label{15} \tag{15}$$

Vamos a sustituir $P(x) = P = \dfrac{b}{a}$ y $y_{1}(x) = e^{-\frac{b}{2a} x}$ en la segunda solución (\ref{14}):

\begin{align*}
y_{2}(x) &= e^{-\frac{b}{2a} x} \int{\dfrac{e^{-\int{\frac{b}{a}dx}}}{\left( e^{-\frac{b}{2a} x} \right)^{2}} dx} \\
&= e^{-\frac{b}{2a} x} \int{\dfrac{e^{-\frac{b}{a}x}}{e^{-\frac{b}{a} x}} dx} \\
&= e^{-\frac{b}{2a} x} \int{dx} \\
&= x e^{-\frac{b}{2a} x}
\end{align*}

Con esto, la segunda solución a la ecuación (\ref{10}) para este segundo caso es:

$$y_{2}(x) = xe^{k_{1}x} = xe^{k_{2}x} = x e^{-\frac{b}{2a} x} \label{16} \tag{16}$$

Usando (\ref{13}) y (\ref{16}) nuevamente notamos que

\begin{align*}
W(y_{1}, y_{2}) &= \begin{vmatrix}
e^{k_{1} x} & xe^{k_{1} x} \\
k_{1} e^{k_{1} x} & e^{k_{1} x} + x k_{1} e^{k_{1} x} \end{vmatrix} \\
&= e^{k_{1} x}(e^{k_{1} x} + x k_{1} e^{k_{1} x}) -x e^{k_{1} x} (k_{1} e^{k_{1} x}) \\
&= e^{2k_{1} x} + x k_{1} e^{2k_{1} x} -x k_{1} e^{2k_{1} x} \\
&= e^{2k_{1} x} \neq 0
\end{align*}

Como el Wronskiano es distinto de cero $\forall x \in \mathbb{R}$, entonces ambas funciones son linealmente independientes en $\mathbb{R}$ y por tanto forman un conjunto fundamental de soluciones de manera que la solución general a la ecuación (\ref{10}), en el caso en el que el discriminante es cero, es

$$y(x) = c_{1}e^{k_{1}x} + c_{2} xe^{k_{1}x} \label{17} \tag{17}$$

Con $k_{1} = -\dfrac{b}{2a}$.

Revisemos el último caso.

Caso 3: Discriminante negativo

Sabemos que

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a}$$

Si el discriminante es menor a cero definimos $-w = b^{2} -4ac$ con $w \in \mathbb{R}^{+}$, tal que

$$k_{1} = \dfrac{-b + \sqrt{-w}}{2a} = \dfrac{-b}{2a} + i\dfrac{\sqrt{w}}{2a}$$

y

$$k_{2} = \dfrac{-b -\sqrt{-w}}{2a} = \dfrac{-b}{2a} -i\dfrac{\sqrt{w}}{2a}$$

Donde $i^{2} = -1$, definimos $\alpha = \dfrac{-b}{2a}$ y $\beta = \dfrac{\sqrt{w}}{2a}$, de esta forma las raíces a la ecuación auxiliar son los números complejos conjugados $k_{1} = \alpha + i\beta$ y $k_{2} = \alpha -i\beta$, donde $\alpha$ y $\beta$ son números reales positivos. Dado que ambas raíces son distintas, similar al caso 1, tendremos dos soluciones dadas por

$$y_{1}(x) = e^{(\alpha + i\beta) x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{(\alpha -i\beta) x}$$

De tarea moral demuestra que $W(y_{1}, y_{2}) = -2i \beta e^{2 \alpha x} \neq 0$, y por tanto $\{ y_{1}, y_{2} \}$ forma un conjunto fundamental de soluciones, es así que la solución general está dada por la combinación lineal

$$y(x) = c_{1}e^{(a + i\beta) x} + c_{2}e^{(a -i\beta) x} \label{18} \tag{18}$$

A pesar de que la función anterior es la solución a la ecuación diferencial (\ref{10}), es común trabajar con una solución que contenga funciones reales en lugar de exponenciales complejas. Con el fin de hallar una expresión distinta para la solución $y(x)$, vamos a considerar la formula de Euler

$$e^{i\theta} = \cos \theta + i\sin \theta \label{19} \tag{19}$$

con $\theta \in \mathbb{R}$. De esta formula y considerando las propiedades de que $\cos (-\beta x) = \cos (\beta x)$ y $\sin (-\beta x) = -\sin (\beta x)$ es que podemos escribir las siguientes expresiones

$$e^{i\beta x} = \cos (\beta x) + i \sin (\beta x) \hspace{1cm} y \hspace{1cm} e^{-i\beta x} = \cos (\beta x) -i \sin (\beta x)$$

De las ecuaciones anteriores puedes observar que se cumplen las siguientes igualdades

$$e^{i\beta x} + e^{-i\beta x} = 2 \cos (\beta x) \hspace{1cm} y \hspace{1cm} e^{i\beta x} -e^{-i\beta x} = 2i \sin (\beta x)$$

Estas ecuaciones nos servirán mas adelante ya que, dada la solución (\ref{18}), si ocurre que $c_{1} = c_{2} = 1$, entonces se obtiene la solución

\begin{align*}
y_{1}(x) &= e^{(a + i\beta )x} + e^{(a -i\beta )x} \\
&= e^{ax} (e^{i\beta x} + e^{-i\beta x}) \\
&= 2 e^{ax} \cos (\beta x)
\end{align*}

Y si ocurre que $c_{1} = 1$ y $c_{2} = -1$, entonces se obtiene la solución

\begin{align*}
y_{2}(x) &= e^{(a + i\beta )x} -e^{(a -i\beta )x} \\
&= e^{ax}(e^{i\beta x} -e^{-i\beta x}) \\
&= 2ie^{ax}\sin (\beta x)
\end{align*}

Con estos resultados vemos que las funciones $g(x) = e^{\alpha x} \cos (\beta x)$ y $h(x) = e^{\alpha x} \sin (\beta x)$ son ahora funciones reales y además de ello son soluciones a la ecuación (\ref{10}).

En la primer entrada de esta segunda unidad demostramos que un múltiplo constante $y(x) = c y_{1}(x)$ de una solución $y_{1}$ de una ecuación diferencial lineal homogénea es también una solución, usando este resultado es que podemos asegurar que las funciones $C_{1}g(x)$ y $C_{2}h(x)$ son también solución. De tarea moral muestra que

$$W(C_{1} e^{\alpha x} \cos (\beta x), C_{2} e^{\alpha x} \sin (\beta x)) = C_{1}C_{2} \beta e^{2 \alpha x} \neq 0$$

Es decir, el Wronskiano de las soluciones es distinto de cero. De esta manera podemos afirmar que ambas funciones forman un conjunto fundamental de soluciones en cierto intervalo $\delta$ y, por lo tanto, podemos concluir que la solución general real a la ecuación (\ref{10}) en el caso en el que el discriminante es negativo es

$$y(x) = C_{1} e^{ax} \cos (\beta x) + C_{2} e^{ax} \sin (\beta x) = e^{ax} (C_{1} \cos (\beta x) + C_{2}\sin (\beta x)) \label{20} \tag{20}$$

Ejemplos

Realicemos una serie de ejemplos en los que tengamos que identificar a que caso pertenecen las ecuaciones diferenciales y así poder obtener su solución.

Ejemplo: Resolver el siguiente PVI: $\dfrac{d^{2}y}{dx^{2}} -4\dfrac{dy}{dx} -5y = 0$, con $y(1) = 0$ y $y^{\prime}(1) = 2$.

Solución: Consideremos la solución $y = e^{kx}$ y sus primeras dos derivadas $\dfrac{dy}{dx}= ke^{kx}$ y $\dfrac{d^{2}y}{dx^{2}}= k^{2}e^{kx}$, sustituyendo en la ecuación diferencial obtenemos lo siguiente:

$$k^{2} e^{kx} -4ke^{kx} -5e^{kx} = e^{kx}(k^{2} -4k -5) = 0$$

Como $e^{kx} \neq 0, \forall x \in \mathbb{R}$ entonces la ecuación auxiliar es $k^{2} -4k -5 = 0$, resolviendo para $k$ tenemos

$$k= \dfrac{4\pm \sqrt{16+20}}{2}= \dfrac{4\pm 6}{2}$$

De donde $k_{1} = 5$ y $k_{2} = -1$. Como ambas raíces son reales y distintas (y $\Delta = 36 > 0$), entonces la ecuación pertenece al caso 1 por lo que podemos afirmar que la solución es de la forma (\ref{12})

$$y(x) = c_{1} e^{5x} + c_{2} e^{-x}$$

La derivada está dada como

$$\dfrac{dy}{dx} = 5 c_{1} e^{5x} -c_{2}e^{-x}$$

Aplicando las condiciones iniciales tenemos

$$y(1) = c_{1}e^{5} + c_{2}e^{-1} = 0 \hspace{1cm} y \hspace{1cm} y^{\prime}(1) = 5 c_{1}e^{5} -c_{2}e^{-1} = 2$$

Si resolvemos el sistema de ecuaciones obtendremos que $c_{1} = \dfrac{e^{-5}}{3}$ y $c_{2} = -\dfrac{e^{1}}{3}$. Por lo tanto, la solución particular del PVI es

$$y(x) = \dfrac{1}{3}e^{5(x -1)} -\dfrac{1}{3}e^{1 -x}$$

$\square$

Ejemplo: Resolver el siguiente PVF: $\dfrac{d^{2}y}{dx^{2}} -10\dfrac{dy}{dx} + 25y = 0$ con $y(0) = 1$ y $y(1) = 0$.

Solución: Consideramos nuevamente la solución $y = e^{kx}$ y sus derivadas $\dfrac{dy}{dx} = ke^{kx}$ y $\dfrac{d^{2}y}{dx^{2}} = k^{2}e^{kx}$, sustituyendo en la ecuación diferencial tenemos

$$k^{2}e^{kx} -10ke^{kx} + 25e^{kx} = e^{kx}(k^{2} -10k + 25) = 0$$

Como $e^{kx} \neq 0, \forall x \in \mathbb{R}$ entonces la ecuación auxiliar es $k^{2} -10k +25 = 0$. Resolviendo para $k$ tenemos

$$k = \dfrac{10 \pm \sqrt{100 -100}}{2} = \dfrac{10}{2} = 5$$

Notamos que $k_{1} = k_{2} = 5$, es decir, son raíces reales e iguales ($\Delta = 0$) por lo que estamos situados en el segundo caso y la solución está dada por la ecuación (\ref{17})

$$y(x) = c_{1} e^{5x} + c_{2}x e^{5x}$$

Aplicamos las condiciones en la frontera

$$y(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} y(1) = c_{1}e^{5} + c_{2}e^{5} = 0$$

Obteniendo que $c_{1} = 1$ y $c_{2} = -1$. Por lo tanto, la solución particular al PVF es

$$y(x) = e^{5x} -x e^{5x}$$

Y además es única al no depender de ningún parámetro libre.

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 2 \dfrac{dy}{dx} + 3y = 0$.

Solución: Considerando la solución $y= e^{kx}$, calculando su primera y segunda derivada y sustituyendo en la ecuación diferencial obtendremos que la ecuación auxiliar es $k^{2} + 2k + 3 = 0$. Resolviendo para $k$, tenemos

$$k = \dfrac{-2 \pm \sqrt{4 -12}}{2} = \dfrac{-2 \pm \sqrt{-8}}{2} = \dfrac{-2 \pm 2i\sqrt{2}}{2} = -1 \pm i \sqrt{2}$$

Las raíces son $k_{1} = -1 + i \sqrt{2}$ y $k_{2} = -1 -i \sqrt{2}$, así mismo, identificamos que $\alpha = -1$ y $\beta = \sqrt{2}$. Al tratarse de raíces complejas notamos que corresponde al caso 3 y su solución compleja esta dada, de acuerdo a la ecuación (\ref{18}), como

$$y(x) = c_{1} e^{(-1 + i \sqrt{2})x} + c_{2} e^{(-1 -i \sqrt{2})x}$$

Sin embargo, una solución real es mucho más práctica de forma que si consideramos la ecuación (\ref{20}) concluimos que la solución general real de la ecuación diferencial es

$$y(x) = e^{-x}(C_{1} \cos (\sqrt{2} x) + C_{2} \sin (\sqrt{2} x))$$

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial $2\dfrac{d^{2}y}{dx^{2}} + 2\dfrac{dy}{dx} + y = 0$.

Solución: Considerando la solución $y= e^{kx}$, calculando su primera y segunda derivada y sustituyendo en la ecuación es como obtendremos que la ecuación auxiliar es $2k^{2} + 2k + 1 = 0$. Resolviendo para $k$, tenemos

$$k = \dfrac{-2 \pm \sqrt{4 -4}}{4} = -\dfrac{2}{4}$$

En este caso $k_{1} = k_{2} = -\dfrac{1}{2}$, es decir, las raíces son reales e iguales de manera que estamos en el caso dos y por lo tanto la solución general es

$$y(x) = c_{1}e^{-\frac{x}{2}} + c_{2}x e^{-\frac{x}{2}}$$

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} -6y = 0$.

Solución: Repitiendo el mismo paso de siempre obtenemos que la ecuación auxiliar es $k^{2} -k -6 = 0$, resolviendo para $k$ tenemos

$$k = \dfrac{1 \pm \sqrt{1 + 24}}{2}$$

Por lo que las raíces son $k_{1} = 3$ y $k_{2} = -2$. Al ser reales y distintas notamos que estamos en circunstancias del caso uno y por tanto la solución general es

$$y(x) = c_{1}e^{3x} + c_{2}e^{-2x}$$

$\square$

Uno más

Ejemplo: Obtener la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 8\dfrac{dy}{dx} + 16y = 0$.

Solución: La ecuación auxiliar en este caso es $k^{2} + 8k + 16= 0$, resolviendo para $k$ se tiene

$$k = \dfrac{8 \pm \sqrt{64 -64}}{2} = \dfrac{8}{2} = 4$$

Las raíces son $k_{1} = k_{2} = 4$, al ser reales e iguales concluimos que la solución general es

$$y(x) = c_{1}e^{4x} + c_{2}x e^{4x}$$

$\square$

Revisemos dos ecuaciones con una forma particular.

Dos ecuaciones particulares

Las ecuaciones diferenciales

$$\dfrac{d^{2}y}{dx^{2}} + m^{2}y = 0 \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} -m^{2}y = 0$$

Con $m \in \mathbb{R}$, son importantes en matemáticas aplicadas. Vamos a obtener la forma de la solución de cada una de ellas.

Para el caso de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + m^{2}y = 0 \label{21} \tag{21}$$

Considerando nuevamente la solución $y(x) = e^{kx}$, la ecuación auxiliar en este caso es $k^{2} + m^{2} = 0$ cuyas raíces son complejas y están dadas por $k_{1} = im$ y $k_{2} = -im$. Si recordamos el caso 3, las raíces son de la forma $k_{1} = \alpha +i\beta$ y $k_{2} = \alpha -i\beta$ así, en este caso tenemos que $\alpha = 0$ y $\beta = m$. Por lo tanto, de la ecuación (\ref{20}) concluimos que la solución a la ecuación (\ref{21}) es

$$y(x) = c_{1} \cos (mx) + c_{2}\sin (mx) \label{22} \tag{22}$$

Ejemplo: Obtener la solución general de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 9y = 0$.

Solución: La ecuación a resolver es de la forma (\ref{21}) con $m = 3$. La ecuación auxiliar en este caso es $k^{2} + 9 = 0$, de donde $k_{1} = i3$ y $k_{2} = -i3$, es decir, $\alpha = 0, \beta = m = 3$. De acuerdo a la ecuación (\ref{22}) la solución general a la ecuación diferencial es

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

$\square$

Por otra parte, para el caso de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -m^{2}y = 0 \label{23} \tag{23}$$

La ecuación auxiliar es $k^{2} -m^{2} = 0$ cuyas raíces son $k_{1} = m$ y $k_{2} = -m$, vemos que en este caso las raíces son reales, entonces la ecuación pertenece al caso 1 y por tanto la solución general está dada de acuerdo a la ecuación (\ref{12}), esto es

$$y(x) = c_{1} e^{mx} + c_{2} e^{-mx} \label{24} \tag{24}$$

Lo interesante de la ecuación (\ref{23}) es que si en su solución (\ref{24}) se elige $c_{1} = c_{2} = \dfrac{1}{2}$ tenemos que

$$y_{1}(x) = \dfrac{1}{2} \left( e^{mx} + e^{-mx} \right) = \cosh (mx)$$

Y si se elige $c_{1} = \dfrac{1}{2}$ y $c_{2} = -\dfrac{1}{2}$ entonces

$$y_{2}(x) = \dfrac{1}{2} \left( e^{mx} -e^{-mx} \right) = \sinh (mx)$$

Puedes comprobar que $W(y_{1}, y_{2}) \neq 0$ lo que prueba que ambas soluciones son linealmente independientes en algún intervalo $\delta$. Este análisis nos permite establecer una forma alternativa de la solución a la ecuación (\ref{23}) dada por

$$y(x) = c_{1} \cosh (mx) + c_{2} \sinh (mx) \label{25} \tag{25}$$

Ejemplo: Resolver el PVI $\dfrac{d^{2}y}{dx^{2}} -3y = 0$ con $y(0) = 1$ y $y^{\prime}(0) = 5$ usando primero la forma de la solución dada en (\ref{24}) y posteriormente resolverlo de nuevo usando la forma de la solución dada en (\ref{25}).

Solución: La ecuación a resolver es de la forma (\ref{23}) con $m = \sqrt{3}$. Considerando la solución $y = e^{kx}$ obtenemos que la ecuación auxiliar es $k^{2} -3 = 0$, de donde $k_{1} = m = \sqrt{3}$ y $k_{2} = -m = -\sqrt{3}$. Usando la ecuación (\ref{24}) concluimos que la solución a la ecuación es

$$y(x) = c_{1} e^{\sqrt{3} x} + c_{2} e^{-\sqrt{3} x}$$

La derivada de la solución es

$$\dfrac{dy}{dx} = \sqrt{3} c_{1} e^{\sqrt{3} x} -\sqrt{3} c_{2} e^{-\sqrt{3} x}$$

Aplicando las condiciones iniciales obtenemos lo siguiente:

$$y(0) = c_{1} + c_{2} = 1 \hspace{1cm} \Rightarrow \hspace{1cm} c_{2} = 1 -c_{1}$$

Por otro lado

$$y^{\prime}(0) = \sqrt{3} c_{1} -\sqrt{3} c_{2} = 5$$

Si sustituimos el valor de $c_{2}$ tenemos

\begin{align*}
\sqrt{3} c_{1} -\sqrt{3} (1 -c_{1}) &= 5 \\
{\sqrt{3}}(c_{1} -1 +c_{1}) &= 5 \\
2 c_{1} -1 &= \dfrac{5}{\sqrt{3}} \\
2c_{1} &= \dfrac{5}{\sqrt{3}} +1 \\
c_{1} &= \dfrac{5}{2\sqrt{3}} + \dfrac{1}{2}
\end{align*}

Sustituyendo en $c_{2}$ se tiene

$$c_{2} = 1 -\left( \dfrac{5}{2\sqrt{3}} -\dfrac{1}{2} \right) = \dfrac{1}{2} -\dfrac{5}{2\sqrt{3}}$$

Por lo tanto, la solución particular del PVI es

$$y(x) = \dfrac{1}{2} \left( 1 + \dfrac{5}{\sqrt{3}} \right) e^{\sqrt{3} x} + \dfrac{1}{2} \left(1 -\dfrac{5}{\sqrt{3}} \right) e^{-\sqrt{3} x}$$

Si por otro lado consideramos la forma de la solución (\ref{25}) obtenemos que

$$y(x) = c_{1} \cosh(\sqrt{3} x) + c_{2} \sinh(\sqrt{3} x)$$

y su derivada

$$\dfrac{dy}{dx} = \sqrt{3} c_{1} \sinh(\sqrt{3} x) + \sqrt{3} c_{2} \cosh(\sqrt{3} x)$$

Aplicando las condiciones iniciales

$$y(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} y^{\prime}(0) = {\sqrt{3}c_{2}} = 5$$

De donde $c_{1} = 1$ y $c_{2} = \dfrac{5}{\sqrt{3}}$. Por lo tanto, la solución alternativa al PVI es

$$y(x) = \cosh (\sqrt{3} x) + \dfrac{5}{\sqrt{3}} \sinh(\sqrt{3} x)$$

$\square$

Para concluir la entrada revisemos brevemente estos mismos resultados para el caso de las ecuaciones de orden superior a dos.

Ecuaciones de orden superior

Es posible aplicar éste método de resolución para el caso de las ecuaciones diferenciales lineales homogéneas de orden superior a dos con coeficientes constantes, sin embargo desarrollar esta teoría desde cero puede ser muy complejo. Ahora que hemos sido muy detallados en el desarrollo de este método para el caso de las ecuaciones de segundo orden, vamos sólo a enunciar los posibles resultados para las ecuaciones de orden superior.

Es importante recordar que estamos estudiando las ecuaciones homogéneas y lineales con coeficientes constantes por lo que una ecuación de orden $n > 2$ con estas características tiene la siguiente forma:

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n -1} \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{2} \dfrac{d^{2}y}{dx^{2}} + a_{1} \dfrac{dy}{dx} + a_{0}y = 0 \label{26} \tag{26}$$

Donde $a_{i}, i = 0, 1, \cdots, n$ son constantes.

Igual que antes, se considera la solución $y(x) = e^{kx}$ de manera que al sustituir esta función y las correspondientes derivadas en la ecuación (\ref{26}) obtenemos la siguiente ecuación auxiliar

$$a_{n} k^{n} + a_{n -1} k^{n -1} + \cdots + a_{2} k^{2} + a_{1} k + a_{0} =0 \label{27} \tag{27}$$

que tendrá $n$ raíces, estas raíces como sabemos pueden ser reales o complejas, iguales o distintas. La solución general para caso está dada de la siguiente manera.

  • Si las raíces son reales y distintas, la solución estará dada por

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} + \cdots + c_{n} e^{k_{n}x} \label{28} \tag{28}$$

  • Si las raíces son reales e iguales, la solución estará dada por

$$y(x) = e^{kx}(c_{1} + c_{2}x + c_{3}x^{2} + \cdots + c_{n}x^{n -1}) \label{29} \tag{29}$$

En el caso de orden superior es posible tener raíces reales y de ellas que unas sean iguales y otras diferentes, en este caso se usan los dos puntos anteriores respectivamente.

Por ejemplo, supongamos que de una ecuación de sexto orden se obtienen seis raíces tales que

$$k_{1} \neq k_{2} = k_{3} \neq k_{4} \hspace{1cm} y \hspace{1cm} k_{1} \neq k_{4} = k_{5} = k_{6}$$

Entonces la solución estaría dada por

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} + c_{3}x e^{k_{3} x} + c_{4} e^{k_{4} x} + c_{5}x e^{k_{5} x} + c_{6}x^{2} e^{k_{6} x}$$

Es decir, el par de raíces $k_{1} \neq k_{2}$ al ser diferentes genera la solución $c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x}$, el par de raíces $k_{2} = k_{3}$ al ser iguales genera la solución $c_{2} e^{k_{2} x} + c_{3}x e^{k_{3} x}$, el par de raíces $k_{1} \neq k_{4}$ genera la solución $c_{1} e^{k_{1} x} + c_{4} e^{k_{4} x}$ y finalmente las raíces $k_{4} = k_{5} = k_{6}$ genera la solución $c_{4} e^{k_{4} x} + c_{5}x e^{k_{5} x} + c_{6}x^{2} e^{k_{6} x}$, esto de acuerdo a los dos puntos anteriores (\ref{28}) y (\ref{29}).

Finalmente

  • Si las raíces son complejas, para cada par conjugado la solución es:

$$y(x) = e^{\alpha x}(C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$

Si hay otro par igual, la función

$$y(x) = e^{\alpha x}x (C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$

será solución y si hay otro par igual, la función

$$y(x) = e^{\alpha x}x^{2} (C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$

será solución y así sucesivamente.

Ejemplo: Resolver el siguiente PVI: $\dfrac{d^{3}y}{dx^{3}} -7 \dfrac{d^{2}y}{dx^{2}} + 4\dfrac{dy}{dx} + 12y = 0$, con $y(0) = 1; y^{\prime}(0) = 0; y^{\prime \prime}(0) = 36$.

Solución: Proponiendo la solución $y = e^{kx}$ la ecuación auxiliar que se obtiene es $k^{3} -7k^{2} +4k +12 = 0$. Factorizando esta ecuación obtenemos que $(k + 1)(k -2)(k -6) = 0$, de donde $k_{1} = -1$, $k_{2} = 2$ y $k_{3} = 6$. Como las raíces son reales y diferentes, entonces la solución a la ecuación diferencial es de la forma (\ref{28}):

$$y(x) = c_{1} e^{-x} + c_{2} e^{2x} + c_{3} e^{6x}$$

Para aplicar las condiciones iniciales calculamos la primera y segunda derivada de la solución.

$$\dfrac{dy}{dx} = -c_{1} e^{-x} + 2c_{2} e^{2x} + 6c_{3} e^{6x} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = c_{1} e^{-x} + 4c_{2} e^{2x} + 36c_{3} e^{6x}$$

Si aplicamos las condiciones iniciales vamos a obtener el siguiente sistema de ecuaciones, respectivamente:

\begin{align*}
c_{1} + c_{2} + c_{3} &= 1 \\
-c_{1} + 2c_{2} + 6c_{3} &= 0 \\
c_{1} + 4c_{2} + 36c_{3} &= 36 \\
\end{align*}

Resolviendo el sistema de ecuaciones obtendremos que $c_{1} = \dfrac{16}{7}$, $c_{2} = -\dfrac{5}{2}$ y $c_{3} = \dfrac{17}{14}$. Por lo tanto, la solución particular de la ecuación diferencial de orden 3 es

$$y(x) = \dfrac{16}{7} e^{-x} -\dfrac{5}{2} e^{2x} + \dfrac{17}{14} e^{6x}$$

$\square$

La dificultad de resolver ecuaciones de orden mayor a 2 realmente radica en que se vuelve más complicado encontrar las raíces de la ecuación auxiliar y resolver el sistema de ecuaciones que se genere en problemas con valores iniciales o con valores en la frontera. Un método para factorizar la ecuación auxiliar y obtener las raíces puede ser la división sintética que seguramente conoces de tus cursos de álgebra, así como el método de Gauss-Jordan para obtener las soluciones de un sistema de ecuaciones lineales.

Tarea moral

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

  1. Obtén la solución general de las siguientes ecuaciones diferenciales con coeficientes constantes.
  • $4 \dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} = 0$
  • $3 \dfrac{d^{2}y}{dx^{2}} + 2 \dfrac{dy}{dx} + y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 13y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -2 \sqrt{3} \dfrac{dy}{dx} + 3y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} + 4 \dfrac{dy}{dx} -y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 5y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -6 \dfrac{dy}{dx} + 13y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -\dfrac{4}{3} \dfrac{dy}{dx} + \dfrac{4}{9} y = 0$
  1. Resuelve los siguientes problemas con valores iniciales:
  • $\dfrac{d^{2}y}{dx^{2}} + y = 0; \hspace{1cm} y \left( \dfrac{\pi}{3} \right) = 0, \hspace{0.5cm} y^{\prime} \left( \dfrac{\pi}{3} \right) = 2$
  • $\dfrac{d^{3}y}{dx^{3}} + 2 \dfrac{d^{2}y}{dx^{2}} -5 \dfrac{dy}{dx} -6y = 0; \hspace{1cm} y(0) = y^{\prime}(0) = 0, \hspace{0.5cm} y^{\prime \prime} (0) = 1$
  1. Resuelve los siguientes problemas con valores en la frontera:
  • $\dfrac{d^{2}y}{dx^{2}} + 4y = 0; \hspace{1cm} y(0) = 0, \hspace{0.5cm} y(\pi) = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -2 \dfrac{dy}{dx} + 2y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y(\pi) = 1$
  1. Resolver el siguiente problema con valores iniciales usando primero la forma de la solución dada en (\ref{24}) y posteriormente resolverlo de nuevo usando la forma de la solución dada en (\ref{25}).
  • $\dfrac{d^{2}y}{dx^{2}} -y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y^{\prime}(1) = 0$

Más adelante…

Continuando con el desarrollo de métodos de resolución de ecuaciones diferenciales de orden superior, en la siguiente entrada desarrollaremos un método que permite resolver ecuaciones que no son homogéneas.

Entradas relacionadas

Ecuaciones Diferenciales I: Método de reducción de orden

Introducción

Hemos comenzado esta segunda unidad estudiando algunas de las propiedades de las soluciones a ecuaciones diferenciales lineales homogéneas y no homogéneas de orden superior. Como mencionamos en la entrada anterior, es momento de comenzar a desarrollar distintos métodos de resolución de ecuaciones diferenciales de orden superior, sin embargo, debido a la complejidad que surge de aumentar el orden, en esta entrada sólo consideraremos las ecuaciones diferenciales de segundo orden.

En esta entrada desarrollaremos el método de reducción de orden, como su nombre lo indica, lo que haremos básicamente es hacer un cambio de variable o una sustitución adecuada que permita que la ecuación de segundo orden pase a ser una ecuación de primer orden y de esta manera aplicar alguno de los métodos vistos en la unidad anterior para resolver la ecuación.

Hay dos distintas formas de reducir una ecuación de segundo orden, la primera de ellas consiste en hacer un cambio de variable $z = \dfrac{dy}{dx}$, esta forma se aplica en ecuaciones tanto lineales como no lineales pero deben satisfacer algunas condiciones y por otro lado, la segunda forma se aplica sólo a ecuaciones lineales homogéneas en las que tenemos conocimiento previo de una solución no trivial. En este segundo caso, considerando que conocemos una solución $y_{1}(x)$, haremos la sustitución $y_{2}(x) = u(x) y_{1}(x)$ para reducir de orden a la ecuación y al resolverla obtendremos la función $u(x)$ y por tanto la segunda solución $y_{2}(x)$ tal que $\{ y_{1}, y_{2} \}$ forme un conjunto fundamental de soluciones de la ecuación diferencial y de esta manera podamos establecer la solución general.

Comencemos por desarrollar la primer forma bajo un cambio de variable.

Ecuaciones reducibles a ecuaciones de primer orden

Hay cierto tipo de ecuaciones de segundo orden que pueden reducirse a una ecuación de primer orden y ser resueltas por los métodos que ya conocemos, vistos en la unidad anterior. Un primer tipo de ecuación son las ecuaciones lineales en las que la variable dependiente $y$ no aparece explícitamente.

Sabemos que una ecuación diferencial lineal no homogénea de segundo orden tiene la siguiente forma:

$$a_{2}(x) \dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x) y = g(x) \label{1} \tag{1}$$

Si la variable dependiente $y$ no se encuentra explícitamente en la ecuación obtenemos la siguiente forma:

$$a_{2}(x) \dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} = g(x) \label{2} \tag{2}$$

Es quizá natural pensar que una forma de resolver la ecuación (\ref{2}) es integrarla dos veces, es esto lo que haremos considerando el siguiente cambio de variable:

$$z = \dfrac{dy}{dx}; \hspace{1cm} \dfrac{dz}{dx} = \dfrac{d^{2}y}{dx^{2}} \label{3} \tag{3}$$

Sean $a_{2}(x) \neq 0$, $P(x) = \dfrac{a_{1}(x)}{a_{2}(x)}$ y $Q(x) = \dfrac{g(x)}{a_{2}(x)}$. Si sustituimos estas funciones y el cambio de variable (\ref{3}) en la ecuación (\ref{2}) lograremos reducirla a una ecuación lineal de primer orden con $z$ la variable dependiente.

$$\dfrac{dz}{dx} + P(x) z = Q(x) \label{4} \tag{4}$$

En la unidad anterior desarrollamos distintos métodos para resolver este tipo de ecuaciones. Una vez que resolvamos la ecuación (\ref{4}) y regresemos a la variable original veremos que dicho resultado nuevamente corresponde a una ecuación de primer orden que podrá ser resuelta una vez más con los métodos vistos anteriormente. Realicemos un ejemplo:

Ejemplo: Reducir de orden a la ecuación diferencial de segundo orden $x \dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} = x$ para $x > 0$ y obtener su solución.

Solución: Primero vamos a dividir toda la ecuación por $x \neq 0$.

$$\dfrac{d^{2}y}{dx^{2}} -\dfrac{1}{x} \dfrac{dy}{dx} = 1$$

Ahora podemos hacer el cambio de variable (\ref{3}) para obtener la forma (\ref{4}).

$$\dfrac{dz}{dx} -\dfrac{1}{x}z = 1 \label{5} \tag{5}$$

Ya no deberíamos tener problema con resolver esta ecuación. El método más eficaz para resolverla es aplicar el método por factor integrante. De la ecuación reducida (\ref{5}) notamos que $P(x) = -\dfrac{1}{x}$ y $Q(x) = 1$. El factor integrante, es este caso, es

$$\mu(x) = e^{\int {P(x)} dx} = e^{-\int \frac{1}{x} dx} = e^{-\ln(x)} = \dfrac{1}{x} \hspace{1cm} \Rightarrow \hspace{1cm} \mu(x) = \dfrac{1}{x}$$

Multiplicamos la ecuación (\ref{5}) por el factor integrante,

$$\dfrac{1}{x} \dfrac{dz}{dx} -\dfrac{z}{x^{2}} = \dfrac{1}{x}$$

e identificamos que

$$\dfrac{d}{dx} \left( \dfrac{z}{x} \right) = \dfrac{1}{x} \dfrac{dz}{dx} -\dfrac{z}{x^{2}}$$

De ambas ecuaciones tenemos que

$$\dfrac{d}{dx} \left( \dfrac{z}{x} \right) = \dfrac{1}{x}$$

Ahora podemos integrar ambos lados de la ecuación con respecto a $x > 0$.

\begin{align*}
\int \dfrac{d}{dx} \left( \dfrac{z}{x} \right) dx &= \int \dfrac{1}{x} dx \\
\dfrac{z}{x} &= \ln (x) + c_{1} \\
z(x) &= x \ln (x) + xc_{1}
\end{align*}

Hemos resuelto la ecuación para la variable $z$, regresemos a la variable original para resolver la nueva ecuación de primer orden.

$$\dfrac{dy}{dx} = x \ln(x) + xc_{1} \label{6} \tag{6}$$

Esta ecuación ahora puede ser resuelta por separación de variables en su versión simple de integración directa, integremos ambos lados de la ecuación con respecto a $x$,

\begin{align*}
\int \dfrac{dy}{dx} dx &= \int x \ln(x) dx + \int xc_{1} dx \\
y(x) &= \int x \ln(x) dx + c_{1} \dfrac{x^{2}}{2}
\end{align*}

Hay que resolver la integral $\int{x \ln(x) dx}$, para ello tomemos $u(x) = \ln(x)$ y $\dfrac{dv}{dx} = x$, así mismo, $\dfrac{du}{dx} = \dfrac{1}{x}$ y $v(x) = \dfrac{x^{2}}{2}$. Entonces

\begin{align*}
\int{x \ln(x) dx} &= \dfrac{x^{2}}{2} \ln(x) -\int{\dfrac{x}{2} dx} \\
&= \dfrac{x^{2}}{2} \ln(x) -\dfrac{x^{2}}{4} + c_{2}
\end{align*}

Sustituyendo en la función $y(x)$.

$$y(x) = \dfrac{x^{2}}{2} \ln(x) -\dfrac{x^{2}}{4} + c_{1} \dfrac{x^{2}}{2} + c_{2}$$

Por lo tanto, la solución general a la ecuación diferencial $x \dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} = x$ es

$$y(x) = \dfrac{x^{2}}{2} \left( \ln(x) -\dfrac{1}{2} \right) + c_{1} \dfrac{x^{2}}{2} + c_{2} \label{7} \tag{7}$$

Verifica por tu cuenta que es la solución general ya que el conjunto $\left\{ y_{1}(x) = \dfrac{x^{2}}{2}, y_{2}(x) = 1 \right\}$ es un conjunto fundamental de soluciones de la ecuación homogénea asociada y $y_{p}(x) = \dfrac{x^{2}}{2} \left( \ln(x) -\dfrac{1}{2} \right)$ es una solución particular de la ecuación no homogénea.

$\square$

Reducción de orden en ecuaciones no lineales

Es posible aplicar un método similar en ecuaciones de segundo orden que pueden ser tanto lineales como NO son lineales, en este caso, a diferencia del caso anterior, la variable dependiente $y$ puede aparecer en la ecuación, sin embargo es necesario que la variable independiente $x$ sea la que no aparezca explícitamente. Este tipo de ecuaciones también pueden reducirse a una ecuación de primer orden pero tomando el siguiente cambio de variable:

$$\dfrac{dy}{dx} = z; \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = z \dfrac{dz}{dy} \label{8} \tag{8}$$

Donde la segunda expresión se deduce de aplicar la regla de la cadena

$$\dfrac{d^{2}y}{dx^{2}} = \dfrac{dz}{dx} = \dfrac{dz}{dy} \dfrac{dy}{dx} = z \dfrac{dz}{dy}$$

Realicemos un ejemplo con una ecuación no lineal.

Ejemplo: Reducir de orden a la ecuación diferencial de segundo orden $\dfrac{d^{2}y}{dx^{2}} -2y \left( \dfrac{dy}{dx}\right)^{3} = 0$ y obtener su solución.

Solución: Es importante notar que es no lineal debido a que la primer derivada es de tercer grado y además esta multiplicada por la función $y$ lo cual no debe ocurrir en el caso lineal.

La ecuación a resolver es $\dfrac{d^{2}y}{dx^{2}} -2y \left( \dfrac{dy}{dx}\right)^{3} = 0$, hacemos el cambio de variable (\ref{8}) y separamos variables

\begin{align*}
z \dfrac{dz}{dy} -2yz^{3} &= 0 \\
\dfrac{dz}{dy} &= 2yz^{2} \\
\dfrac{1}{z^{2}} \dfrac{dz}{dy} &= 2y
\end{align*}

Integramos ambos lados de la ecuación con respecto a $y$

\begin{align*}
\int{\dfrac{1}{z^{2}} \dfrac{dz}{dy} dy} &= \int{2y dy} \\
\int{\dfrac{dz}{z^{2}}} &= 2 \int{y dy} \\
-\dfrac{1}{z} &= y^{2} + c_{1} \\
z &= -\dfrac{1}{y^{2} + c_{1}}
\end{align*}

Regresando a la variable original y separando de nuevo las variables, tenemos

\begin{align*}
\dfrac{dy}{dx} &= -\dfrac{1}{y^{2} + c_{1}} \\
(y^{2} + c_{1}) \dfrac{dy}{dx} &= -1
\end{align*}

Integrando ambos lados de la ecuación con respecto a $x$

\begin{align*}
\int{(y^{2} + c_{1}) \dfrac{dy}{dx} dx} &= -\int{dx} \\
\int{y^{2} dy} + \int{c_{1} dy} &= -\int{dx} \\
\dfrac{y^{3}}{3} + c_{1}y &= -x + c_{2}
\end{align*}

Por lo tanto, la solución implícita a la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -2y \left( \dfrac{dy}{dx}\right)^{3} = 0$ es

$$\dfrac{y^{3}}{3} + c_{1}y = c_{2} -x$$

$\square$

Realicemos un ejemplo más con una ecuación lineal.

Ejemplo: Encontrar la solución general de la ecuación diferencial $4 \dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} = 0$.

Solución: Como la ecuación no contiene explícitamente a la función $y$ ni a la variable independiente $x$, entonces podemos aplicar cualquier cambio de variable, ya sea (\ref{3}) u (\ref{8}). Vamos a resolverla aplicando ambos casos.

Primero consideremos el cambio de variable (\ref{8}):

\begin{align*}
4z \dfrac{dz}{dy} + z &= 0 \\
4 \dfrac{dz}{dy} &= -1 \\
\dfrac{dz}{dy} &= -\dfrac{1}{4}
\end{align*}

Integramos ambos lados de la ecuación con respecto a $y$

\begin{align*}
\int{\dfrac{dz}{dy} dy} &= -\int{\dfrac{1}{4} dy} \\
\int{dz} &= -\dfrac{1}{4} \int{dy} \\
z &= -\dfrac{1}{4} y + c_{1}
\end{align*}

Regresando a la variable original

\begin{align*}
\dfrac{dy}{dx} &= -\dfrac{1}{4} y + c_{1} \\
\dfrac{dy}{dx} + \dfrac{y}{4} &= c_{1}
\end{align*}

Resolvemos esta ecuación por factor integrante.

$$\mu(x) = e^{\int {P(x)} dx} = e^{\int \frac{1}{4} dx} = e^{\frac{x}{4}} \hspace{1cm} \Rightarrow \hspace{1cm} \mu(x) = e^{\frac{x}{4}}$$

Multiplicamos ambos lados de la ecuación por el factor integrante

\begin{align*}
e^{\frac{x}{4}} \dfrac{dy}{dx} + e^{\frac{x}{4}} \dfrac{y}{4} &= e^{\frac{x}{4}} c_{1} \\
\dfrac{d}{dx}\left( y e^{\frac{x}{4}} \right) &= c_{1} e^{\frac{x}{4}}
\end{align*}

Ahora integramos ambos lados con respecto a $x$

\begin{align*}
\int{\dfrac{d}{dx}\left( y e^{\frac{x}{4}} \right) dx} &= \int{c_{1} e^{\frac{x}{4}} dx} \\
y e^{\frac{x}{4}} &= c_{1} \int{e^{\frac{x}{4}} dx} \\
y e^{\frac{x}{4}} &= c_{1} 4 e^{\frac{x}{4}} + c_{2} \\
y(x) &= c_{2} e^{-\frac{x}{4}} + 4c_{1}
\end{align*}

Renombrando a las constantes concluimos que la solución general a la ecuación diferencial $4 \dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} = 0$ es

$$y(x) = k_{1} e^{-\frac{x}{4}} + k_{2}$$

Resolvamos de nuevo la ecuación pero ahora aplicando el cambio de variable (\ref{3}),

\begin{align*}
4 \dfrac{dz}{dx} + z &= 0 \\
\dfrac{1}{z} \dfrac{dz}{dx} &= -\dfrac{1}{4}
\end{align*}

Integramos ambos lados con respecto a $x$

\begin{align*}
\int{\dfrac{1}{z} \dfrac{dz}{dx} dx} &= -\int{\dfrac{1}{4} dx} \\
\int{\dfrac{dz}{z}} &= -\dfrac{1}{4}\int{dx} \\
\ln|z| &= -\dfrac{x}{4} + c_{1} \\
z &= c_{2}e^{-\frac{x}{4}}
\end{align*}

Con $c_{2} = e^{c_{1}}$. Regresando a la variable original

$$\dfrac{dy}{dx} = c_{2}e^{-\frac{x}{4}}$$

Integramos ambos lados con respecto a $x$

\begin{align*}
\int{\dfrac{dy}{dx} dx} = \int{c_{2}e^{-\frac{x}{4}} dx} \\
\int{dy} = c_{2} \int{e^{-\frac{x}{4}} dx} \\
y = -c_{2}4 e^{-\frac{x}{4}} + c_{3}
\end{align*}

Si renombramos las constantes obtenemos nuevamente que

$$y(x) = k_{1} e^{-\frac{x}{4}} + k_{2}$$

$\square$

Ahora pasemos a un método de reducción de orden en el que previamente debemos conocer una solución. Este método es el que usualmente encontrarás como método de reducción de orden.

Reducción de orden conocida una solución

Es posible reducir una ecuación diferencial lineal homogénea de segundo orden

$$a_{2}(x) \dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x) y = 0 \label{9} \tag{9}$$

a una ecuación diferencial de primer orden siempre que se conozca previamente una solución no trivial $y_{1}(x)$. Recordamos de la entrada anterior que una ecuación de la forma (\ref{9}) tiene como solución general la combinación lineal $y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x)$, con $y_{1}$ y $y_{2}$ funciones que forman un conjunto fundamental de soluciones en cierto intervalo $\delta$. Si conocemos $y_{1}$ podremos reducir la ecuación a una de primer orden y resolverla para obtener la solución $y_{2}$ y por tanto obtener la solución general.

Este método también es conocido como método de reducción de orden pues tiene el mismo objetivo que los casos anteriores, reducir de orden a una ecuación diferencial. La idea general del método es la siguiente:

Comenzaremos con el conocimiento previo de una solución no trivial $y_{1}(x)$ de la ecuación homogénea (\ref{9}) definida en un intervalo $\delta$. Lo que buscamos es una segunda solución $y_{2}(x)$ tal que $y_{1}$ y $y_{2}$ formen un conjunto fundamental de soluciones en $\delta$, es decir, que sean soluciones linealmente independientes entre sí. Recordemos que si ambas soluciones son linealmente independientes entonces el cociente $\dfrac{y_{2}}{y_{1}}$ no es contante en $\delta$, es decir, $\dfrac{y_{2}(x)}{y_{1}(x)} = u(x)$ que lo podemos reescribir como $y_{2}(x) = u(x) y_{1}(x)$. Como queremos encontrar $y_{2}$ y previamente conocemos $y_{1}$ entonces debemos determinar la función $u(x)$, dicha función se determina al sustituir $y_{2} = uy_{1}$ en la ecuación diferencial dada, esto reducirá a dicha ecuación a una de primer orden donde la variable dependiente será $u$.

Desarrollemos el método de manera general para encontrar la expresión de $u$ y por tanto de $y_{2}$ y finalmente realicemos un ejemplo.

Método de reducción de orden

Este método se aplica a las ecuaciones diferenciales de la forma

$$a_{2}(x) \dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x) y = 0$$

Si dividimos esta ecuación por $a_{2}(x) \neq 0$ podemos obtener la forma estándar

$$\dfrac{d^{2}y}{dx^{2}} + P(x) \dfrac{dy}{dx} + Q(x)y = 0 \label{10} \tag{10}$$

Con $P(x) = \dfrac{a_{1}(x)}{a_{2}(x)}$ y $Q(x) = \dfrac{a_{0}(x)}{a_{2}(x)}$ ambas continuas en algún intervalo $\delta$. Supongamos además que $y_{1}(x)$ es una solución conocida de (\ref{10}) en $\delta$ y que $y_{1}(x) \neq 0 $ para toda $x \in \delta$. Si se define $y(x) = u(x) y_{1}(x)$, derivando se tiene que

$$\dfrac{dy}{dx} = u \dfrac{dy_{1}}{dx} + y_{1} \dfrac{du}{dx} \label{11} \tag{11}$$

Derivando una segunda ocasión

$$\dfrac{d^{2}y}{dx^{2}} = u \dfrac{d^{2}y_{1}}{dx^{2}} + 2\dfrac{dy_{1}}{dx} \dfrac{du}{dx} + y_{1} \dfrac{d^{2}u}{dx^{2}} \label{12} \tag{12}$$

Sustituyendo (\ref{11}) y (\ref{12}) en la forma estándar (\ref{10}) obtenemos lo siguiente:

\begin{align*}
\dfrac{d^{2}y}{dx^{2}} + P\dfrac{dy}{dx} + Qy &= \left[ u \dfrac{d^{2}y_{1}}{dx^{2}} + 2\dfrac{dy_{1}}{dx} \dfrac{du}{dx} + y_{1} \dfrac{d^{2}u}{dx^{2}} \right] + P \left[ u \dfrac{dy_{1}}{dx} + y_{1} \dfrac{du}{dx} \right] + Q \left[ u y_{1}\right] \\
&= u \left[ \dfrac{d^{2}y_{1}}{dx^{2}} + P\dfrac{dy_{1}}{dx} + Qy_{1} \right] + y_{1} \dfrac{d^{2}u}{dx^{2}} + \left( 2 \dfrac{dy_{1}}{dx} + Py_{1} \right) \dfrac{du}{dx} \\
&= 0
\end{align*}

Como $y_{1}(x)$ es solución sabemos que

$$\dfrac{d^{2}y_{1}}{dx^{2}} + P\dfrac{dy_{1}}{dx} + Qy_{1} = 0$$

Entonces el resultado anterior se reduce a lo siguiente:

$$y_{1} \dfrac{d^{2}u}{dx^{2}} + \left( 2 \dfrac{dy_{1}}{dx} + Py_{1} \right) \dfrac{du}{dx} = 0 \label{13} \tag{13}$$

Consideremos el cambio de variable $w = \dfrac{du}{dx}$ y $\dfrac{dw}{dx} = \dfrac{d^{2}y}{dx^{2}}$, entonces la ecuación (\ref{13}) se puede escribir como

$$y_{1} \dfrac{dw}{dx} + \left( 2 \dfrac{dy_{1}}{dx} + Py_{1} \right) w = 0 \label{14} \tag{14}$$

Esta ecuación es tanto lineal como separable. Separando las variables e integrando, se obtiene

\begin{align*}
\dfrac{1}{w}\dfrac{dw}{dx} + 2\dfrac{1}{y_{1}} \dfrac{dy_{1}}{dx} &= -P \\
\int{\dfrac{dw}{w}} + 2\int{\dfrac{dy_{1}}{y_{1}}} &= -\int{P dx} \\
\ln |w| + 2 \ln|y_{1}| + k &= -\int{P dx} \\
\ln |w y^{2}_{1}| + k &= -\int{P dx} \\
wy^{2}_{1} &= k_{1}e^{-\int{P dx}}
\end{align*}

Despejando a $w$ de la última ecuación, usando $w = \dfrac{du}{dx}$ e integrando nuevamente

\begin{align*}
\dfrac{du}{dx} &= \dfrac{k_{1}e^{-\int{P dx}}}{y^{2}_{1}} \\
\int{du} &= \int{\dfrac{k_{1}e^{-\int{P dx}}}{y^{2}_{1}} dx} \\
u &= k_{1} \int{\dfrac{e^{-\int{P} dx}}{y^{2}_{1}} dx} + k_{2}
\end{align*}

Eligiendo $k_{1} = 1$ y $k_{2} = 0$ obtenemos la expresión para la función $u(x)$,

$$u(x) = \int{\dfrac{e^{-\int{P} dx}}{y^{2}_{1}} dx} \label{15} \tag{15}$$

Si sustituimos en $y(x) = y_{2}(x) = u(x)y_{1}(x)$ obtenemos que la segunda solución a la ecuación diferencial (\ref{10}) es

$$y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x)} dx}}{y^{2}_{1}(x)} dx} \label{16} \tag{16}$$

De tarea moral puedes probar que la función $y_{2}$ satisface la ED y que $y_{1}$ y $y_{2}$ son linealmente independientes en algún intervalo en el que $y_{1}$ no es cero.

Realicemos un ejemplo en el que apliquemos este método.

Ejemplo: Encontrar la solución general a la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 16y = 0$ dada la solución no trivial $y_{1}(x) = \cos(4x)$.

Solución: Apliquemos directamente la expresión (\ref{16}) para obtener la solución $y_{2}(x)$.

La ecuación diferencial a resolver es $\dfrac{d^{2}y}{dx^{2}} + 16y = 0$, si la comparamos con la forma estándar (\ref{10}) notamos que $P(x) = 0$ y $Q(x) = 16$. Sustituyendo en $y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x)} dx}}{y^{2}_{1}(x)} dx}$, tenemos

\begin{align*}
y_{2}(x) &= \cos(4x) \int{\dfrac{e^{0}}{\cos^{2}(4x)} dx} \\
&= \cos(4x) \int{\dfrac{1}{\cos^{2}(4x)} dx}
\end{align*}

Para resolver la integral consideramos el cambio de variable $s = 4x$, $ds = 4 dx$

$$\int{\dfrac{1}{\cos^{2}(4x)} dx} = \dfrac{1}{4} \int{\sec^{2}(s) ds}$$

Y sabemos que $\int{\sec^{2}(s) ds} = \tan(s)$, así

$$y_{2}(x) = \cos(4x) \left( \dfrac{1}{4} \tan(4x) + k_{1} \right)$$

Hacemos $k_{1} = 0$

$$y_{2}(x) = \dfrac{\cos(4x)}{4} \left( \dfrac{\sin(4x)}{\cos(4x)} \right) = \dfrac{\sin(4x)}{4}$$

Como la solución general corresponde a la combinación lineal $y(x) = c_{1} y_{1}(x) + c_{2} y_{2}(x)$, en las constantes $c_{1}$ y $c_{2}$ se pueden englobar todas las constantes que pudieran aparecer, por ello es que podemos tomar $k_{1} = 0$ y además podemos evitar la constante $\dfrac{1}{4}$ de $y_{2}$ y considerar que $y_{2}(x) = \sin(4x)$. Veamos que efectivamente satisface la ecuación diferencial.

$$\dfrac{dy_{2}}{dx} = 4 \cos(4x) \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{d^{2}y_{2}}{dx^{2}} = -16 \sin(4x)$$

Sustituyendo en la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + 16y = -16 \sin(4x) + 16 \sin(4x) = 0$$

Cumple con la ecuación diferencial, lo mismo podemos verificar con la solución dada $y_{1}(x) = \cos(4x)$.

$$\dfrac{dy_{1}}{dx} = -4 \sin(4x) \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{d^{2}y_{1}}{dx^{2}} = -16 \cos(4x)$$

Sustituyendo en la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + 16y = -16 \cos(4x) + 16 \cos(4x) = 0$$

Como ambas soluciones son linealmente independientes, entonces forman un conjunto fundamental de soluciones. Otra forma de verificarlo es mostrando que el Wronskiano es distinto de cero y lo es ya que $W(y_{1}, y_{2}) = 4 \neq 0$.

Por lo tanto, la solución general a la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} + 16y = 0$ corresponde a la combinación lineal

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

$\square$

Con esto concluimos esta entrada sobre un primer método para resolver algunas ecuaciones de segundo orden. En la siguiente entrada continuaremos estudiando un nuevo método.

Tarea moral

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

  1. Obtén la solución general de las siguientes ecuaciones diferenciales lineales.
  • $x \dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} = 0$
  • $(x-1) \dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} = 0$
  1. Resuelve las siguientes ecuaciones diferenciales no lineales.
  • $(y -1)\dfrac{d^{2}y}{dx^{2}} = \left( \dfrac{dy}{dx} \right)^{2} $
  • $\left( \dfrac{dy}{dx} \right)^{2} -2 \dfrac{d^{2}y}{dx^{2}} = 0$
  1. Dada una solución no trivial de las siguientes ecuaciones diferenciales, hallar la segunda solución tal que ambas formen un conjunto fundamental de soluciones y determina la solución general.
  • $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 0; \hspace{1cm} y_{1}(x) = e^{2x}$
  • $\dfrac{d^{2}y}{dx^{2}} -25y = 0; \hspace{1cm} y_{1}(x) = e^{5x}$
  1. Demuestra que la función

$$y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x)} dx}}{y^{2}_{1}(x)} dx}$$

Satisface la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + P(x) \dfrac{dy}{dx} + Q(x)y = 0$$

Siempre que $y_{1}(x)$ sea solución a la misma ecuación.

  1. Usando el inciso anterior, demostrar que

$$\left \{ y_{1}(x), y_{1}(x) \int{\dfrac{e^{-\int{P(x)} dx}}{y^{2}_{1}(x)} dx} \right \}$$

es un conjunto fundamental de soluciones de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + P(x) \dfrac{dy}{dx} + Q(x)y = 0$$

Más adelante…

En esta entrada conocimos un método de reducción de orden basado en un cambio de variable para ecuaciones lineales y no lineales de segundo orden que satisfacen algunas condiciones y desarrollamos el método de reducción de orden para ecuaciones diferenciales lineales homogéneas en el caso en el que previamente tenemos conocimiento de una solución no trivial.

En la siguiente entrada estudiaremos otro método para resolver un tipo particular de ecuaciones diferenciales, éstas son las ecuaciones diferenciales lineales homogéneas con coeficientes constantes, de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0$$

Con $a, b$ y $c$ constantes.

Entradas relacionadas

Ecuaciones Diferenciales I – Videos: Método de la transformada de Laplace

Introducción

Bienvenidos a la última entrada de la segunda unidad del curso, donde revisaremos el método de la transformada de Laplace para resolver problemas de condición inicial de la forma $$ay^{\prime \prime}+by’+cy=f(t) \ ; \ y(0)=y_{0} \ , \ y'(0)=y_{1}$$ con $a$, $b$ y $c$ constantes.

Este método nos permite transformar el problema de resolver la ecuación diferencial por los métodos estudiados en esta misma unidad, por un problema algebraico donde encontraremos la expresión de la transformada de Laplace $\mathscr{L} \{y(t)\}$ de la función solución, y debemos hallar quién es la función $y(t)$ cuya transformada de Laplace es $\mathscr{L} \{y(t)\}$.

Comenzaremos definiendo la transformada de Laplace de una función cuyo dominio es el intervalo $[0, \infty)$, y demostraremos algunas de las propiedades más importantes que cumple esta transformada y que utilizaremos para nuestros propósitos.

Posteriormente resolveremos el problema de condición inicial de manera general, mencionaremos el problema de hallar la transformada inversa de Laplace de una función con ayuda de una tabla de transformadas y transformadas inversas, y revisaremos dos ejemplos particulares donde mostraremos cómo se utiliza el método en la práctica.

Para finalizar consideraremos nuevamente el problema de condición inicial $$ay^{\prime \prime}+by’+cy=f(t) \ ; \ y(0)=y_{0} \ , \ y'(0)=y_{1}$$ donde ahora la función $f(t)$ es una función continua por pedazos. Este tipo de problemas suele aparecer con frecuencia en la física, y con la ayuda de la transformada de Laplace vamos a resolver un ejemplo particular, con ayuda de una función auxiliar y un teorema que enunciaremos y probaremos previamente.

Como te podrás dar cuenta, hicimos un cambio en la notación de la derivada de una función. Durante el curso hemos utilizado la notación de Leibniz $\frac{dy}{dt}$, $\frac{d^{2}y}{dt^{2}}$,…, para denotar a las derivadas de la función $y(t)$. Sin embargo, en esta entrada utilizaremos la notación $y'(t)$, $y^{\prime \prime}(t)$,…, para simplificar la escritura.

Transformada de Laplace y sus propiedades

En el primer video de esta entrada definimos la transformada de Laplace $\mathscr{L} \{y(t)\}$ de una función cuyo dominio es el intervalo $[0, \infty)$, y probamos algunas propiedades que cumple esta transformada y que nos servirán para resolver problemas de condición inicial.

Solución a problemas de condición inicial por método de la transformada de Laplace

En el primer video de esta sección resolvemos el problema de condición inicial $$ay^{\prime \prime}+by’+cy=f(t) \ ; \ y(0)=y_{0} \ , \ y'(0)=y_{1}$$ por el método de la transformada de Laplace.

En el segundo video resolvemos un par de problemas de condición inicial particulares.

Te presentamos una tabla de transformadas y transformadas inversas de Laplace, que aparece en el libro Elementary Differential Equations and Boundary Value Problems, William E. Boyce y Richard C. DiPrima, para que puedas realizar los cálculos presentados en los videos. Esta tabla no es única, por lo que puedes buscar en textos o en internet tablas diferentes según lo requieras.

Tabla de transformadas de Laplace y transformadas inversas
Tabla de transformadas y transformadas inversas de Laplace. Boyce y DiPrima (2012).

Solución a problemas de condición inicial con funciones discontinuas por método de transformada de Laplace

En el último video de esta entrada resolvemos un problema de condición inicial de la forma $$ay^{\prime \prime}+by’+cy=f(t) \ ; \ y(0)=y_{0} \ , \ y'(0)=y_{1}$$ donde $f(t)$ es una función continua por pedazos. Previamente definimos la función auxiliar $$H_{c}(t)= \begin{cases} 1 & 0 \leq t < c \\ 1 & t \geq c \ \end{cases} $$ y probamos un teorema que nos ayudan a resolver el problema de condición inicial.

Tarea moral

  • Encuentra la transformada de Laplace de la función $f(t)=t$.
  • Encuentra la transformada de Laplace de la función $f(t)=\cos{\beta t}$, $\beta$ constante.
  • Prueba que $$\mathscr{L} \{f ^{\prime \prime} (t)\}=s^{2}\mathscr{L}\{f(t)\}-sf(0)-f^{\prime}(0)$$ bajo las hipótesis del último teorema del primer video.
  • Resuelve el problema de condición inicial $$y^{\prime \prime}-2y’+2y=\cos{t} \ \ ; \ \ y(0)=1, \ \ y'(0)=0$$ por el método de la transformada de Laplace.
  • Prueba que bajo las condiciones del primer teorema enunciado en el primer video, se cumplen las siguiente propiedad: $$F^{(n)}(s)=\mathscr{L}\{(-t)^{n}f(t)\}$$ donde $F^{(n)}$ denota a la $n$-ésima derivada de $F$.
  • Resuelve el problema de condiciones iniciales $$y^{\prime \prime}+2y’+2y=f(t) \ \ ; \ \ y(0)=0, \ \ y'(0)=1$$ donde $$f(t)= \begin{cases} 1 & \pi \leq t < 2\pi \\ 0 & 0 \leq t < \pi \ , t \geq 2\pi. \end{cases} $$

Más adelante

Con esta entrada concluimos el estudio a las ecuaciones diferenciales de segundo orden. Como mencionamos en esta entrada, toda la teoría desarrollada en la segunda unidad se puede extender a ecuaciones de orden $n>2$. Sin embargo, a partir de la tercera unidad utilizaremos un método distinto para resolver ecuaciones de orden $n\geq2$.

Lo primero que haremos en la siguiente entrada será transformar una ecuación diferencial que orden $n\geq2$ en un sistema de $n$ ecuaciones diferenciales de primer orden, hablaremos de las ventajas de hacer esta transformación y daremos una introducción a los sistemas de ecuaciones diferenciales de primer orden, que será el objeto de estudio de la tercera unidad del curso.

Entradas relacionadas

Ecuaciones Diferenciales I: Soluciones a ecuaciones diferenciales de orden superior

Introducción

En la entrada anterior comenzamos a estudiar los problemas con valores iniciales (PVI) y problemas con valores en la frontera (PVF), ambos para el caso de las ecuaciones diferenciales lineales de orden superior. Vimos también que si $y_{1}, y_{2}, \cdots, y_{k}$ son $k$ soluciones de una ecuación homogénea de $n$-ésimo orden en un intervalo $\delta$, entonces la combinación lineal $y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{k}y_{k}(x)$, donde las $c_{i}$, $i = 1, 2, \cdots, k$ son constantes, también es solución en el intervalo $\delta$, este resultado es conocido como principio de superposición y nuestro propósito en esta entrada es estudiar las propiedades de todas estas soluciones donde la independencia lineal de las funciones jugará un papel muy importante en la construcción del conjunto fundamental de soluciones y de la solución general.

Es importante tener presente el concepto de conjunto fundamental de soluciones presentado en la entrada anterior. A continuación se muestra de nuevo la definición.

Definición: Cualquier conjunto $\{y_{1}, y_{2}, \cdots, y_{n}\}$ de $n$ soluciones linealmente independientes de una ecuación diferencial lineal homogénea de $n$-ésimo orden en un intervalo $\delta$ es un conjunto fundamental de soluciones en dicho intervalo.

Soluciones a ecuaciones diferenciales lineales de orden superior

Estamos interesados en soluciones linealmente independientes de una ecuación diferencial lineal homogénea de orden superior de la forma

\begin{align}
a_{n}(x) \dfrac{d^{n}y}{dx^{n}} + a_{n-1}(x) \dfrac{d^{n-1}y}{dx^{n-1}} + \cdots + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = 0 \label{1} \tag{1}
\end{align}

Al intentar responder la cuestión de si el conjunto de $n$ soluciones $\{y_{1}, y_{2}, \cdots, y_{n}\}$ de (\ref{1}) es linealmente independiente podemos apelar directamente a la definición de independencia lineal, sin embargo esta pregunta se puede responder de una forma mecánica usando un determinante llamado el Wronskiano.

Definición: Supongamos que cada una de las funciones $f_{1}(x), f_{2}(x), \cdots, f_{n}(x)$ tiene al menos $n -1$ derivadas. El determinante

\begin{align}
W(f_{1}, f_{2}, \cdots, f_{n}) = \begin{vmatrix}
f_{1} & f_{2} & \cdots & f_{n} \\
f_{1}^{\prime} & f_{2}^{\prime} &\cdots &f_{n}^{\prime} \\
\vdots &\vdots & & \vdots \\
f_{1}^{(n -1)} & f_{2}^{(n -1)} & \cdots & f_{n}^{(n -1)}
\end{vmatrix} \label{2} \tag{2}
\end{align}

donde las primas denotan derivadas, se llama el Wronskiano de las funciones.

El Wronskiano es una herramienta que podemos utilizar para determinar si el conjunto de soluciones a la ecuación (\ref{1}) es un conjunto linealmente independiente y la forma de hacerlo es a través del siguiente teorema conocido como criterio para soluciones linealmente independientes.

Teorema: Sean $y_{1}, y_{2}, \cdots, y_{n}$, $n$ soluciones de la ecuación diferencial lineal homogénea de $n$-ésimo orden (\ref{1}) en el intervalo $\delta$. El conjunto de soluciones es linealmente independiente en $\delta$ si, y solo si $W(y_{1}, y_{2}, \cdots, y_{n}) \neq 0$ para toda $x$ en el intervalo $\delta$.

Este teorema nos dice que sólo basta mostrar que el Wronskiano es distinto de cero para garantizar que el conjunto de soluciones $\{y_{1}, y_{2}, \cdots, y_{n}\}$ es linealmente independiente y por tanto formará un conjunto fundamental de soluciones.

Al conjunto de soluciones linealmente independiente $\{y_{1}, y_{2}, \cdots, y_{n}\}$ de la ecuación (\ref{1}) se le denomina fundamental porque, así como cualquier vector en $\mathbb{R}^{3}$ se puede expresar como una combinación lineal de los vectores linealmente independientes $\hat{i}, \hat{j}$ y $\hat{k}$, cualquier solución a una ecuación de la forma (\ref{1}) se puede expresar como una combinación lineal de las $n$ soluciones del conjunto fundamental, podemos decir que las soluciones $\{y_{1}, y_{2}, \cdots, y_{n}\}$ son los bloques básicos para la solución general de la ecuación.

En el siguiente teorema se enuncia la forma general de la solución a la ecuación (\ref{1}).

Teorema: Sea $\{y_{1},y_{2}, \cdots, y_{n}\}$ un conjunto fundamental de soluciones de la ecuación diferencial lineal homogénea de $n$-ésimo orden (\ref{1}) en el intervalo $\delta$. Entonces la solución general de la ecuación en el intervalo es

$$y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x) \label{3} \tag{3}$$

donde $c_{i} , i = 1,2, \cdots, n$ son constantes arbitrarias.

Aterricemos estas ideas generales al caso de las ecuaciones diferenciales lineales homogéneas de segundo orden.

Ecuaciones lineales homogéneas de segundo orden

Una ecuación diferencial lineal homogénea de segundo orden es de la forma

$$a_{2}(x) \dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = 0 \label{4} \tag{4}$$

Sobre esta ecuación desarrollaremos la siguiente teoría. Primero definamos el Wronskiano para el caso $n = 2$.

Definición: Sean $y_{1}(x)$ y $y_{2}(x)$ soluciones de la ecuación (\ref{4}), el Wronskiano de las funciones es

\begin{align}
W(y_{1}, y_{2}) = \begin{vmatrix}
y_{1} & y_{2} \\
\dfrac{dy_{1}}{dx} & \dfrac{dy_{2}}{dx} \\
\end{vmatrix} = y_{1}\dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx} \label{5} \tag{5}
\end{align}

Ahora que conocemos la forma del Wronskiano para $n = 2$ demostremos el teorema de la solución general para el caso $n = 2$.

Teorema: Sean $y_{1}(x)$ y $y_{2}(x)$ soluciones de la ecuación (\ref{4}) en el intervalo $\delta$ y supongamos que $W(y_{1}, y_{2}) \neq 0$ para toda $x \in \delta$, entonces $y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x)$ es la solución general de la ecuación diferencial (\ref{4}).

Demostración: Sea $y(x)$ una solución a la ecuación diferencial (\ref{4}) en el intervalo $\delta$ y sea $x_{0} \in \delta$ tal que $y(x_{0}) = \alpha$ y $\dfrac{dy}{dx}(x_{0}) = \beta$ con $\alpha$ y $\beta$ constantes. Supongamos que existen $c_{1}$ y $c_{2}$ constantes tales que

$$\alpha = c_{1}y_{1}(x_{0}) + c_{2}y_{2}(x_{0}) \label{6} \tag{6}$$

y

$$\beta = c_{1} \dfrac{dy_{1}}{dx}(x_{0}) + c_{2} \dfrac{dy_{2}}{dx}(x_{0}) \label{7} \tag{7}$$

esto debido a que por hipótesis $y_{1}(x)$ y $y_{2}(x)$ son soluciones a la ecuación diferencial y por tanto la combinación lineal también lo será. Aplicando el teorema de existencia y unicidad obtenemos que la solución $y(x)$ tiene que ser de la forma $y(x) = c_{1}y_{1} + c_{2}y_{2}$ por lo que nuestro problema se reduce a demostrar que las constantes $c_{1}$ y $c_{2}$ existen.

Si multiplicamos a la ecuación (\ref{6}) por $\dfrac{dy_{2}}{dx}(x_{0})$ y a la ecuación (\ref{7}) por $y_{2}(x_{0})$ obtenemos lo siguiente, respectivamente:

$$\alpha \dfrac{dy_{2}}{dx}(x_{0}) = c_{1}y_{1}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) + c_{2}y_{2}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) \label{8} \tag{8}$$

y

$$\beta y_{2}(x_{0}) = c_{1} y_{2}(x_{0}) \dfrac{dy_{1}}{dx}(x_{0}) + c_{2} y_{2}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) \label{9} \tag{9}$$

Restémosle a la ecuación (\ref{8}) la ecuación (\ref{9}).

\begin{align*}
\alpha \dfrac{dy_{2}}{dx}(x_{0}) -\beta y_{2}(x_{0}) &= c_{1} y_{1}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) -c_{1} y_{2}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) \\
&= c_{1} \left( y_{1}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) -y_{2}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) \right) \label{10} \tag{10}
\end{align*}

Sabemos que el Wronskiano, en $x = x_{0}$, está definido como

$$W(y_{1}(x_{0}), y_{2}(x_{0})) = y_{1}(x_{0}) \dfrac{dy_{2}}{dx}(x_{0}) -y_{2}(x_{0}) \dfrac{dy_{1}}{dx}(x_{0})$$

Por comodidad denotaremos a $W(y_{1}(x_{0}), y_{2}(x_{0}))$ como $W(x_{0})$. Entonces la ecuación (\ref{10}) se puede escribir de la siguiente manera:

$$\alpha \dfrac{dy_{2}}{dx}(x_{0}) -\beta y_{2}(x_{0}) = c_{1} W(x_{0})$$

Debido a que por hipótesis $W(y_{1}, y_{2}) \neq 0$ para toda $x \in \delta$, en particular lo es en $x = x_{0}$, por tanto podemos despejar a la constante $c_{1}$ y así obtener un valor para dicha constante lo que muestra su existencia.

Para obtener la expresión de $c_{2}$ seguimos el mismo camino, multiplicamos a la ecuación (\ref{6}) por $\dfrac{dy_{1}}{dx}(x_{0})$ y a la ecuación (\ref{7}) por $y_{1}(x_{0})$ y repetimos el mismo procedimiento demostrando que existe un valor para la constante $c_{2}$.

Como hemos encontrado valores para $c_{1}$ y $c_{2}$ entonces existen y por lo tanto la solución general a la ecuación (\ref{4}) es $y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x)$.

$\square$

Ya hemos definido lo que es el conjunto fundamental de soluciones de una ecuación diferencial de orden $n$, para el caso $n = 2$ lo podemos definir de la siguiente manera:

Definición: Decimos que $\{ y_{1}(x), y_{2}(x) \}$ es un conjunto fundamental de soluciones para la ecuación (\ref{4}) si cualquier solución a la ED puede escribirse como combinación lineal de $y_{1}$ y $y_{2}$, o lo que es equivalente, que $y_{1}$ y $y_{2}$ sean linealmente independientes en $\delta$.

Así, si encontramos un conjunto fundamental de soluciones $\{ y_{1}(x), y_{2}(x) \}$, entonces $W(y_{1}, y_{2}) \neq 0$ para toda $x \in \delta$ y por tanto $y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x)$ será la solución general de la ecuación diferencial (\ref{4}).

Del criterio para soluciones linealmente independientes se puede hacer notar que cuando $y_{1}, y_{2}, \cdots, y_{n}$ son $n$ soluciones de la ecuación (\ref{1}) en un intervalo $\delta$, el Wronskiano $W(y_{1}, y_{2}, \cdots, y_{n})$ es siempre igual a cero o nunca es cero en todo $\delta$. Vamos a demostrar este hecho para el caso $n = 2$.

Teorema: Sean $a_{2}(x) \neq 0$, $a_{1}(x)$ y $a_{0}(x)$ funciones continuas en $\delta$ de la ecuación (\ref{4}) y sean $y_{1}(x)$ y $y_{2}(x)$ soluciones a la misma ecuación en $\delta$. Entonces $W(y_{1}, y_{2}) = 0$ o $W(y_{1}, y_{2}) \neq 0$ para todo $x \in \delta$.

Demostración: Como $y_{1}(x)$ y $y_{2}(x)$ son soluciones de la ecuación (\ref{4}), entonces

$$a_{2}(x) \dfrac{d^{2}y_{1}}{dx^{2}} + a_{1}(x) \dfrac{dy_{1}}{dx} + a_{0}(x)y_{1} = 0 \label{11} \tag{11}$$

y

$$a_{2}(x) \dfrac{d^{2}y_{2}}{dx^{2}} + a_{1}(x) \dfrac{dy_{2}}{dx} + a_{0}(x)y_{2} = 0 \label{12} \tag{12}$$

Si multiplicamos a la ecuación (\ref{11}) por $y_{2}$ y a la ecuación (\ref{12}) por $y_{1}$ obtenemos lo siguiente, respectivamente.

$$y_{2}a_{2}(x) \dfrac{d^{2}y_{1}}{dx^{2}} + y_{2} a_{1}(x) \dfrac{dy_{1}}{dx} + y_{2}a_{0}(x)y_{1} = 0 \label{13} \tag{13}$$

y

$$y_{1}a_{2}(x) \dfrac{d^{2}y_{2}}{dx^{2}} + y_{1}a_{1}(x) \dfrac{dy_{2}}{dx} + y_{1}a_{0}(x)y_{2} = 0 \label{14} \tag{14}$$

A la ecuación (\ref{14}) vamos a restarle la ecuación (\ref{13}):

$$a_{2}(x) \left( y_{1} \dfrac{d^{2}y_{2}}{dx^{2}} -y_{2}\dfrac{d^{2}y_{1}}{dx^{2}} \right) + a_{1}(x) \left( y_{1} \dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx}\right) = 0 \label{15} \tag{15}$$

Sabemos que $W(y_{1}, y_{2}) = y_{1} \dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx}$ y notemos lo siguiente:

\begin{align*}
\dfrac{dW}{dx} &= \dfrac{d}{dx} \left( y_{1} \dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx}\right) \\
&= \dfrac{dy_{1}}{dx} \dfrac{dy_{2}}{dx} + y_{1} \dfrac{d^{2}y_{2}}{dx^{2}} -\dfrac{dy_{2}}{dx} \dfrac{dy_{1}}{dx} -y_{2} \dfrac{d^{2}y_{1}}{dx^{2}} \\
&= y_{1} \dfrac{d^{2}y_{2}}{dx^{2}} -y_{2}\dfrac{d^{2}y_{1}}{dx^{2}}
\end{align*}

Es decir,

$$\dfrac{dW}{dx} = y_{1} \dfrac{d^{2}y_{2}}{dx^{2}} -y_{2}\dfrac{d^{2}y_{1}}{dx^{2}} \label{16} \tag{16}$$

En términos del Wronskiano la ecuación (\ref{15}) se puede escribir como

$$a_{2}(x) \dfrac{dW}{dx} + a_{1}(x) W = 0 \label{17} \tag{17}$$

Como $a_{2}(x) \neq 0$ para toda $x \in \delta$, entonces podemos definir la función $P(x) = \dfrac{a_{1}(x)}{a_{2}(x)}$ tal que la ecuación (\ref{17}) se pueda escribir como

$$\dfrac{dW}{dx} + P(x) W = 0 \label{18} \tag{18}$$

La ecuación (\ref{18}) corresponde a una ecuación diferencial lineal homogénea de primer orden y ya sabemos que la solución es de la forma $W(x) = ke^{-\int{P(x)} dx}$, de manera que hay dos posibilidades:

  • Si $k = 0 \hspace{0.5cm} \Rightarrow \hspace{0.5cm} W(x) = 0, \hspace{0.5cm} \forall x \in \delta$
  • Si $k \neq 0 \hspace{0.5cm} \Rightarrow \hspace{0.5cm} W(x) \neq 0, \hspace{0.5cm} \forall x \in \delta$

$\square$

El criterio para soluciones linealmente independientes nos garantiza que si el Wronskiano es distinto de cero entonces el conjunto de soluciones es linealmente independiente en $\delta$, lo opuesto es cierto bajo ciertas condiciones, si el Wronskiano es igual a cero entonces el conjunto de soluciones es linealmente dependiente. Demostremos este hecho.

Teorema: Sean $a_{2}(x) \neq 0$, $a_{1}(x)$ y $a_{0}(x)$ funciones continuas en $\delta$ de la ecuación (\ref{4}) y sean $y_{1}$ y $y_{2} \neq 0$ soluciones a la misma ecuación en $\delta$. Supongamos que $W(y_{1}(x), y_{2}(x)) = 0$, $\forall x \in \delta$. Entonces $y_{1}(x)$ y $y_{2}(x)$ son linealmente dependientes.

Demostración: Por hipótesis $W(y_{1}(x), y_{2}(x)) = 0$, $\forall x \in \delta$, es decir

$$y_{1} \dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx} = 0 \label{19} \tag{19}$$

Consideremos el siguiente resultado:

$$\dfrac{d}{dx} \left( -\dfrac{y_{1}}{y_{2}} \right) = \dfrac{1}{y^{2}_{2}} \left( y_{1} \dfrac{dy_{2}}{dx} -y_{2}\dfrac{dy_{1}}{dx} \right) \label{20} \tag{20}$$

Donde hemos considerado la hipótesis $y_{2} \neq 0$. Si usamos la hipótesis (\ref{19}) obtenemos que $\dfrac{d}{dx} \left( -\dfrac{y_{1}}{y_{2}} \right) = 0$, $\forall x \in \delta$, integrando esta ecuación obtenemos que $-\dfrac{y_{1}}{y_{2}} = -k$, o bien, $y_{1}(x) = k y_{2}(x)$, con $k$ una constante. Esto demuestra que $y_{1}$ y $y_{2}$ son linealmente dependientes.

$\square$

Hay que tener muy presentes las hipótesis de este teorema pues es posible que el Wronskiano sea cero aún cuando las funciones consideradas en un cierto intervalo sean linealmente independientes en él.

Como consecuencia del teorema anterior podemos establecer el criterio para soluciones linealmente independientes en el caso $n = 2$.

Corolario: Dos soluciones $y_{1}(x)$ y $y_{2}(x)$ de la ecuación diferencial (\ref{4}) son linealmente independientes en $\delta$ si, y sólo si $W(y_{1}(x), y_{2}(x)) \neq 0$, $\forall x \in \delta$.

Realicemos un ejemplo.

Ejemplo: En la entrada anterior de tarea moral tenías que verificar que las funciones $y_{1}(x) = e^{-3x}$ y $y_{2}(x) = e^{4x}$ forman un conjunto fundamental de soluciones de la ecuación diferencial $\dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} -12y = 0$ en $\delta = (-\infty, \infty)$. Vamos a demostrar esto mismo usando los teoremas vistos anteriormente.

Solución: Consideremos las soluciones $y_{1}(x) = e^{-3x}$ y $y_{2}(x) = e^{4x}$ y sus respectivas derivadas $\dfrac{dy_{1}}{dx} = -3e^{-3x}$ y $\dfrac{dy_{2}}{dx} = 4e^{4x}$. Calculemos el Wronskiano.

\begin{align*}
W(y_{1}, y_{2}) = \begin{vmatrix}
e^{-3x} & e^{4x} \\
-3e^{-3x} & 4e^{4x} \\
\end{vmatrix} = e^{-3x}(4e^{4x}) -e^{4x}(-3e^{-3x}) = 7e^{-x} \neq 0
\end{align*}

Como $W(y_{1}, y_{2}) = 7 e^{-x} \neq 0$, entonces $\{ y_{1}(x) = e^{-3x}, y_{2}(x) = e^{4x}\}$ forma un conjunto fundamental de soluciones y la solución general está dada por $y(x) = c_{1}e^{-3x} + c_{2}e^{4x}$.

$\square$

Con esto concluimos el estudio de algunas propiedades importantes de las soluciones a la ecuación diferencial lineal homogénea de orden superior, terminemos esta entrada con el estudio del caso no homogéneo.

Ecuaciones NO homogéneas

La ecuación diferencial lineal no homogénea de $n$-ésimo orden es

$$a_{n}(x) \dfrac{d^{n}y}{dx^{n}} + a_{n-1}(x) \dfrac{d^{n-1}y}{dx^{n-1}} + \cdots + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = g(x) \label{21} \tag{21}$$

Nuestro objetivo es obtener la forma general de la solución a la ecuación no homogénea (\ref{21}) y estudiar algunas propiedades de las soluciones.

Definición: Una función $y_{p}(x)$ que es libre de parámetros arbitrarios y que satisface la ecuación diferencial lineal no homogénea de $n$-ésimo orden (\ref{21}) se dice que es una solución particular o integral particular de la ecuación.

Si recordamos al operador polinomial

$$\mathcal{L} = a_{n}(x)D^{n} + a_{n -1}(x)D^{n -1} + \cdots + a_{1}(x)D + a_{0}(x)$$

la definición anterior implica que $\mathcal{L}\{y_{p}\} = g(x)$. Veamos el siguiente resultado.

Teorema: Sean $y_{1}, y_{2}, \cdots, y_{n}$ soluciones a la ecuación diferencial homogénea (\ref{1}) en el intervalo $\delta$ y sea $y_{p}$ la solución particular de la ecuación no homogénea (\ref{21}) en $\delta$. Entonces la combinación lineal

$$y(x) = c_{1}y_{1} + c_{2}y_{2} + \cdots + c_{n}y_{n} + y_{p} \label{22} \tag{22}$$

es solución a la ecuación no homogénea (\ref{21}), en $\delta$.

Demostración: Sea $y(x)$ la combinación lineal $y(x) = c_{1}y_{1} + c_{2}y_{2} + \cdots + c_{n}y_{n} + y_{p}$, si aplicamos el operador polinomial, tenemos

\begin{align*}
\mathcal{L} \{y(x)\} &= \mathcal{L} \{c_{1}y_{1} + c_{2}y_{2} + \cdots + c_{n}y_{n} + y_{p}\} \\
&= c_{1}\mathcal{L}\{y_{1}\} + c_{2}\mathcal{L}\{y_{2}\} + \cdots + c_{n}\mathcal{L}\{y_{n}\} + \mathcal{L}\{y_{p}\} \\
&= 0 + g(x) \\
&= g(x)
\end{align*}

Ya que $\mathcal{L}\{y_{i}\} = 0$ para cada $i = 1, 2, \cdots, n$ por ser cada $y_{i}$ solución a la ecuación homogénea, mientras que $\mathcal{L}\{y_{p}\} = g(x)$ por ser solución a la ecuación no homogénea. Entonces, como $\mathcal{L} \{y(x)\} = g(x)$, concluimos que la combinación lineal $y(x) = c_{1}y_{1} + c_{2}y_{2} + \cdots + c_{n}y_{n} + y_{p}$ es solución a la ecuación diferencial no homogénea.

$\square$

¿Y qué ocurre si las soluciones $y_{1}, y_{2}, \cdots, y_{n}$ forman un conjunto fundamental de soluciones?. La respuesta es que la combinación lineal $y(x) = c_{1}y_{1} + c_{2}y_{2} + \cdots + c_{n}y_{n} + y_{p}$ sería la solución general a la ecuación diferencial no homogénea (\ref{21}). Demostremos este resultado.

Teorema: Sea $y_{p}(x)$ cualquier solución particular de la ecuación diferencial lineal no homogénea de $n$-ésimo orden (\ref{21}) en el intervalo $\delta$, y sea $\{ y_{1}, y_{2}, \cdots, y_{n} \}$ un conjunto fundamental de soluciones de la ecuación diferencial homogénea asociada (\ref{1}) en $\delta$. Entonces la solución general de la ecuación no homogénea es

$$y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x) + y_{p}(x) \label{23} \tag{23}$$

en el intervalo $\delta$ y donde $c_{i}, i = 1, 2, \cdots, n$ son constantes arbitrarias.

Demostración: Sea $y(x)$ la solución general de la ecuación no homogénea (\ref{21}) y sea $y_{p}(x)$ una solución particular de la misma ecuación, ambas definidas en el intervalo $\delta$, de manera que $\mathcal{L} \{ y(x)\} = \mathcal{L} \{ y_{p}(x)\} = g(x)$, con $\mathcal{L}$ el operador polinomial. Nuestro objetivo es encontrar la forma explícita de $y(x)$.

Definamos la función $h(x) = y(x) -y_{p}(x)$ y notemos lo siguiente:

\begin{align*}
\mathcal{L} \{ h(x) \} &= \mathcal{L} \{ y(x) -y_{p}(x) \} \\
&= \mathcal{L} \{ y(x) \} -\mathcal{L} \{ y_{p}(x) \} \\
&= g(x) -g(x) \\
&= 0
\end{align*}

Esto es, $\mathcal{L} \{ h(x) \} = 0$, lo que significa que la función $h(x)$ es solución a la ecuación homogénea (\ref{1}) y por el teorema de la solución general de ecuaciones homogéneas podemos establecer que la función $h(x)$ tiene la siguiente forma:

$$h(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x) \label{24} \tag{24}$$

Con $\{ y_{1}, y_{2}, \cdots, y_{n} \}$ un conjunto fundamental de soluciones. Sustituyendo (\ref{24}) en $h(x) = y(x) -y_{p}(x)$ y despejando a la solución general $y(x)$ obtenemos finalmente que

$$y(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x) + y_{p}(x)$$

que es lo que queríamos demostrar.

$\square$

Recuerda, la diferencia entre las soluciones $(\ref{22})$ y $(\ref{23})$ es que en $(\ref{23})$ las $y_{i}, i = 1, 2, \cdots, n$ forman un conjunto fundamental de soluciones, es decir, son linealmente independientes entre sí, mientras que en (\ref{22}) no necesariamente forman una conjunto fundamental y sin embargo, también son solución a la ecuación (\ref{21}).

En el caso de las ecuaciones no homogéneas vemos que la solución general corresponde a la suma de la solución general de la ecuación homogénea asociada más una solución particular de la ecuación no homogénea, en este caso no homogéneo la solución general de la ecuación homogénea tiene un nombre especial.

Definición: La combinación lineal $y_{c}(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x)$, que es la solución general de la ecuación diferencial lineal homogénea de $n$-ésimo orden (\ref{1}), se denomina función complementaria para la ecuación no homogénea (\ref{21}).

Por tanto, resolver una ecuación lineal no homogénea implica resolver primero la ecuación homogénea asociada para obtener la función complementaria $y_{c}(x)$ y luego se encuentra una solución particular $y_{p}(x)$ a la ecuación no homogénea para finalmente sumarlas $y(x) = y_{c}(x) + y_{p}(x)$.

Realicemos un ejemplo.

Ejemplo: Probar que la función $y(x) = c_{1} e^{2x} + c_{2}x e^{2x} + x^{2} e^{2x} + x -2$, definida en el intervalo $\delta = (-\infty, \infty)$, es la solución general de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 2 e^{2x} + 4x -12$$

Solución: Primero probemos que las funciones $y_{1} = e^{2x}$ y $y_{2} = x e^{2x}$ forman un conjunto fundamental de soluciones de la ecuación homogénea asociada

$$\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 0$$

Para ello veamos que $y_{1}$ y $y_{2}$ son soluciones de la ecuación homogénea y que son linealmente independientes, es decir, que $W(y_{1}, y_{2}) \neq 0$. Calculemos las derivadas.

$$y_{1} = e^{2x} \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{dy_{1}}{dx} = 2 e^{2x} \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{d^{2}y_{1}}{dx^{2}} = 4 e^{2x}$$

$$y_{2} = xe^{2x} \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{dy_{2}}{dx} = e^{2x} + 2x e^{2x} \hspace{1cm} \Rightarrow \hspace{1cm} \dfrac{d^{2}y_{2}}{dx^{2}} = 4 e^{2x} + 4x e^{2x}$$

Si sustituyes los valores correspondientes en la ecuación homogénea notaras que se satisface la ecuación, es decir, se cumple que

$$\dfrac{d^{2}y_{1}}{dx^{2}} -4 \dfrac{dy_{1}}{dx} + 4y_{1} = 0 \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{2}}{dx^{2}} -4 \dfrac{dy_{2}}{dx} + 4y_{2} = 0$$

Por lo tanto, ambas son soluciones de la ecuación homogénea. Calculemos el Wronskiano.

\begin{align*}
W(y_{1}, y_{2}) = \begin{vmatrix}
e^{2x} & xe^{2x} \\
2 e^{2x} & e^{2x} + 2x e^{2x} \\
\end{vmatrix} = e^{2x}(e^{2x} + 2x e^{2x}) -xe^{2x}(2e^{2x}) = e^{4x} \neq 0
\end{align*}

Como $W(y_{1}, y_{2}) \neq 0$, $\forall x \in \delta$, por los teoremas vistos anteriormente concluimos que $\{y_{1} = e^{2x}, y_{2} = x e^{2x} \}$ forma un conjunto fundamental de soluciones de la ecuación homogénea asociada y que la solución general a dicha ecuación es $y_{c}(x) = c_{1} e^{2x} + c_{2}x e^{2x}$, donde el subíndice $c$ indica que es la función complementaria.

Ahora verifiquemos que la función $y_{p}(x) = x^{2} e^{2x} + x -2$ es una solución particular de la ecuación no homogénea. Calculemos la primera y segunda derivada.

$$\dfrac{dy_{p}}{dx} = 2x e^{2x} + 2x^{2} e^{2x} + 1$$

$$\dfrac{d^{2}y_{p}}{dx^{2}} = 2 e^{2x} + 8x e^{2x} + 4x^{2} e^{2x}$$

Sustituyendo en la ecuación diferencial tenemos

\begin{align*}
\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y &= (2 e^{2x} + 8x e^{2x} + 4x^{2} e^{2x}) -4(2x e^{2x} + 2x^{2} e^{2x} + 1) +4(x^{2} e^{2x} + x -2) \\
&= 2e^{2x} + (8x e^{2x} -8x e^{2x}) + (4x^{2} e^{2x} -8x^{2} e^{2x} + 4x^{2} e^{2x}) + 4x -12 \\
&= 2e^{2x} +4x -12
\end{align*}

Esto es, $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 2e^{2x} +4x -12$, que justo corresponde a la ecuación diferencial no homogénea, por lo tanto, efectivamente $y_{p}$ es una solución particular.

Como $\{y_{1} = e^{2x}, y_{2} = x e^{2x} \}$ es un conjunto fundamental de soluciones de la ecuación homogénea asociada y $y_{p}(x) = x^{2} e^{2x} + x -2$ es una solución particular de la ecuación no homogénea, por el teorema de la solución general de ecuaciones no homogéneas concluimos que la función $y(x) = c_{1} e^{2x} + c_{2}x e^{2x} + x^{2} e^{2x} + x -2$ es la solución general a la ecuación no homogénea.

$\square$

Hay algo muy curioso que ocurre en el ejemplo anterior. Mostramos que la función $y_{p}(x) = x^{2} e^{2x} + x -2$ es una solución particular de la ecuación no homogénea

$$\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 2e^{2x} +4x -12 = g(x)$$

Sin embargo, si haces los cálculos correspondientes notarás que la función $y_{p1}(x) = x^{2} e^{2x}$ es una solución particular de la ecuación

$$\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 2e^{2x} = g_{1}(x)$$

mientras que la función $y_{p2}(x) = x -2$ es una solución particular de la ecuación

$$\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 4y = 4x -12 = g_{2}(x)$$

Así, si superponemos las soluciones particulares $y_{p}(x) = y_{p1}(x) + y_{p2}(x)$ obtenemos en la ecuación diferencial la superposición de la funciones $g(x) = g_{1}(x) + g_{2}(x)$.

Lo anterior es efecto del principio de superposición para ecuaciones no homogéneas.

Teorema: Sean $y_{p1}, y_{p2}, \cdots, y_{pk}$, $k$ soluciones particulares de la ecuación diferencial lineal no homogénea de $n$-ésimo orden (\ref{21}) en un intervalo $\delta$, que corresponde, a su vez, a $k$ funciones diferentes $g_{1}, g_{2}, \cdots, g_{k}$. Es decir, se supone que $y_{pi}$ denota una solución particular de la ecuación diferencial correspondiente

$$a_{n}(x) \dfrac{d^{n}y}{dx^{n}} + a_{n-1}(x) \dfrac{d^{n-1}y}{dx^{n-1}} + \cdots + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = g_{i}(x) \label{25} \tag{25}$$

donde $i = 1, 2, , \cdots, k$. Entonces,

$$y_{p}(x) = y_{p1}(x) + y_{p2}(x) + \cdots + y_{pk}(x) \label{26} \tag{26}$$

es una solución particular de

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n-1} \dfrac{d^{n-1}y}{dx^{n-1}} + \cdots + a_{1} \dfrac{dy}{dx} + a_{0}y = g_{1}(x) + g_{2}(x) + \cdots + g_{k}(x) \label{27} \tag{27}$$

Demostración: Sea $\mathcal{L}$ el operador polinomial y sean $y_{pi}(x)$, $i = 1, 2, \cdots, k$, soluciones particulares de las ecuaciones no homogéneas $\mathcal{L} \{ y_{pi}(x) \} = g_{i}$, $i = 1, 2, \cdots, k$ respectivamente. Vamos a definir la función

$$y_{p}(x) = y_{p1}(x) + y_{p2}(x) + \cdots + y_{pk}(x)$$

Nuestro objetivo es demostrar que la función $y_{p}(x)$ es una solución particular de la ecuación (\ref{27}), es decir, que se cumple que

$$\mathcal{L} \{ y_{p}(x) \} = g_{1}(x) + g_{2}(x) + \cdots + g_{k}(x)$$

En efecto

\begin{align*}
\mathcal{L} \{ y_{p}(x)\} &= \mathcal{L} \{ y_{p1}(x) + y_{p2}(x) + \cdots + y_{pk}(x) \} \\
&= \mathcal{L} \{ y_{p1}(x) \} + \mathcal{L} \{ y_{p2}(x) \} + \cdots + \mathcal{L} \{ y_{pk}(x) \} \\
&= g_{1}(x) + g_{2}(x) + \cdots + g_{k}(x)
\end{align*}

Con esto queda probado que $y_{p}(x) = y_{p1}(x) + y_{p2}(x) + \cdots + y_{pk}(x)$ es solución de (\ref{27}).

$\square$

Corolario: Si la soluciones $y_{pi}(x)$ son soluciones particulares de (\ref{25}) para $i = 1, 2, \cdots, k$, entonces la combinación lineal

$$y_{p}(x) = c_{1}y_{p1}(x) + c_{2}y_{p2}(x) + \cdots + c_{k}y_{pk}(x) \label{28} \tag{28}$$

donde las $c_{i}$ son constantes, es también una solución particular de la ecuación diferencial

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n-1} \dfrac{d^{n-1}y}{dx^{n-1}} + \cdots + a_{1} \dfrac{dy}{dx} + a_{0}y = c_{1}g_{1} + c_{2}g_{2} + \cdots + c_{k}g_{k} \label{29} \tag{29}$$

Realicemos un último ejemplo.

Ejemplo: Probar que

  • $y_{p1}(x) = -4x^{2} \hspace{0.5cm}$ es solución particular de $\hspace{0.5cm} \dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = -16x^{2} + 24x -8$,
  • $y_{p2}(x) = e^{2x} \hspace{0.9cm}$ es solución particular de $\hspace{0.5cm} \dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = 2e^{2x}$,
  • $y_{p3}(x) = xe^{x} \hspace{0.9cm}$ es solución particular de $\hspace{0.5cm} \dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = 2x e^{x} -e^{x}$.

y probar que la superposición $y = y_{p1} + y_{p2} + y_{p3} = -4x^{2} + e^{2x} + xe^{x}$ es una solución de

$$\dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = -16x^{2} + 24x -8 + 2e^{2x} + 2xe^{x} -e^{x}$$

Solución: Sean $g_{1}(x) = -16x^{2} + 24x -8$, $g_{2}(x) = 2e^{2x}$ y $g_{3}(x) = 2x e^{x} -e^{x}$, de tarea moral muestra que efectivamente

$$\dfrac{d^{2}y_{p1}}{dx^{2}} -3\dfrac{dy_{p1}}{dx} + 4y_{p1} = g_{1}(x)$$

$$\dfrac{d^{2}y_{p2}}{dx^{2}} -3\dfrac{dy_{p2}}{dx} + 4y_{p2} = g_{2}(x)$$

$$\dfrac{d^{2}y_{p3}}{dx^{2}} -3\dfrac{dy_{p3}}{dx} + 4y_{p3} = g_{3}(x)$$

Por el principio de superposición para ecuaciones no homogéneas sabemos que la función $y(x) = y_{p1}(x) + y_{p2}(x) + y_{p3}(x)$ es solución de la ecuación

$$\dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = g_{1}(x) + g_{2}(x) + g_{3}(x)$$

Por lo tanto, la función $y(x) = -4x^{2} + e^{2x} + xe^{x}$ es solución a la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} + 4y = -16x^{2} + 24x -8 + 2e^{2x} + 2xe^{x} -e^{x}$$

Si gustas puedes calcular la primera y segunda derivada de $y(x)$ y verificar la ecuación anterior para asegurarte del resultado.

$\square$

Con esto concluimos nuestro estudio sobre algunas propiedades de las ecuaciones diferenciales lineales de orden superior. En la siguiente entrada conoceremos un primer método para resolver ecuaciones diferenciales de segundo orden.

Tarea moral

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

  1. Dadas las soluciones a las siguientes ecuaciones diferenciales lineales homogéneas de orden superior en el intervalo dado, calcular el Wronskiano para determinar si es un conjunto fundamental de soluciones y en caso de serlo dar la solución general.
  • $x^{3} \dfrac{d^{3}y}{dx^{3}} + 6x^{2} \dfrac{d^{2}y}{dx^{2}} + 4x \dfrac{dy}{dx} -4y = 0$, con soluciones

$\hspace{1cm} y_{1} = x, \hspace{0.6cm} y_{2} = \dfrac{1}{x^{2}}, \hspace{0.6cm} y_{3} = \dfrac{1}{x^{2}} \ln(x); \hspace{1cm} \delta = (0, \infty)$.

  • $\dfrac{d^{4}y}{dx^{4}} + \dfrac{d^{2}y}{dx^{2}} = 0$, con soluciones

$\hspace{1cm} y_{1} = 1, \hspace{0.6cm} y_{2} = x, \hspace{0.6cm} y_{3} = \cos(x), \hspace{0.6cm} y_{4} = \sin(x); \hspace{1cm} \delta = (\infty, \infty)$.

  1. Dadas las soluciones a las siguientes ecuaciones diferenciales lineales NO homogéneas de orden superior en el intervalo dado, probar que se trata de la solución general a la ecuación.
  • $\dfrac{d^{2}y}{dx^{2}} -7 \dfrac{dy}{dx} + 10y = 24 e^{x}$, con solución

$\hspace{1cm} y(x) = c_{1} e^{2x} + c_{2} e^{5x} + 6 e^{x}; \hspace{1cm} \delta = (\infty, \infty)$.

  • $2x^{2} \dfrac{d^{2}y}{dx^{2}} + 5x \dfrac{dy}{dx} + y = x^{2} -x$, con solución

$\hspace{1cm} y(x) = c_{1} \dfrac{1}{\sqrt{x}} + c_{2} \dfrac{1}{x} + \dfrac{1}{15}x^{2} -\dfrac{1}{6}x; \hspace{1cm} \delta = (0, \infty)$.

  1. Comprobar que las funciones $y_{p1}(x) = 3 e^{2x}$ y $y_{p2}(x) = x^{2} + 3x$ son, respectivamente, soluciones particulares de

$$\dfrac{d^{2}y}{dx^{2}} -6\dfrac{dy}{dx} + 5y = -9 e^{2x} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} -6\dfrac{dy}{dx} + 5y = 5x^{2} + 3x -16$$

  1. Usando el ejercicio anterior, encontrar la solución particular de la ecuación

$\hspace{1cm}\dfrac{d^{2}y}{dx^{2}} -6\dfrac{dy}{dx} + 5y = 5x^{2} + 3x -16 -9 e^{2x}$

$\hspace{1cm}$ y

$\hspace{1cm} \dfrac{d^{2}y}{dx^{2}} -6\dfrac{dy}{dx} + 5y = -10x^{2} -6x +32 + e^{2x}$.

Más adelante…

Ahora que ya conocemos algunas propiedades de las ecuaciones diferenciales de orden superior y sus soluciones, en particular de las lineales de segundo orden, es momento de comenzar a estudiar los distintos métodos de resolución de estas ecuaciones diferenciales.

En la siguiente entrada comenzaremos con un método que permite reducir una ecuación de segundo orden a una ecuación de primer orden de manera que podremos resolverla aplicando alguno de los métodos vistos en la unidad anterior. No es casualidad que dicho método se conozca como método de reducción de orden.

Entradas relacionadas