8.6: ODE inhomogéneas

Ahora consideramos la ecuación de oda lineal no homogénea general de segundo orden\ ref {8.1}:

\[\ddot{x}+p(t) \dot{x}+q(t) x=g(t), \nonumber \]

con condiciones iniciales\(x\left(t_{0}\right)=x_{0}\) y\(\dot{x}\left(t_{0}\right)=u_{0}\). Hay un método de solución de tres pasos cuando el término no homogéneo\(g(t) \neq 0\). (i) Encontrar la solución general de la ecuación homogénea

\[\ddot{x}+p(t) \dot{x}+q(t) x=0 . \nonumber \]

Denotemos la solución homogénea por

\[x_{h}(t)=c_{1} X_{1}(t)+c_{2} X_{2}(t), \nonumber \]

donde\(X_{1}\) y\(X_{2}\) son soluciones linealmente independientes de\(Equation \ref{8.22}\),\(c_{1}\) y\(c_{2}\) son todavía constantes indeterminadas. (ii) Encontrar una solución particular\(x_{p}\) de la ecuación no homogénea Ecuación\ ref {8.21}. Una solución particular se encuentra fácilmente cuando\(p(t)\) y\(q(t)\) son constantes, y cuando\(g(t)\) es una combinación de polinomios, exponenciales, senos y cosenos. iii) Escribir la solución general de la Ecuación\ ref {8.21} como la suma de las soluciones homogéneas y particulares,

\[x(t)=x_{h}(t)+x_{p}(t), \nonumber \]

y aplicar las condiciones iniciales para determinar las constantes\(c_{1}\) y\(c_{2}\). Obsérvese que debido a la linealidad de\(Equation \ref{8.21}\),

\[\begin{aligned} \ddot{x}+p \dot{x}+q x &=\frac{d^{2}}{d t^{2}}\left(x_{h}+x_{p}\right)+p \frac{d}{d t}\left(x_{h}+x_{p}\right)+q\left(x_{h}+x_{p}\right) \\ &=\left(\ddot{x}_{h}+p \dot{x}_{h}+q x_{h}\right)+\left(\ddot{x}_{p}+p \dot{x}_{p}+q x_{p}\right) \\ &=0+g \\ &=g \end{aligned} \nonumber \]

para que la Ecuación\ ref {8.23} resuelva\(Equation \ref{8.21}\), y las dos constantes libres en se\(x_{h}\) puedan usar para satisfacer las condiciones iniciales.

Consideraremos aquí solo el caso del coeficiente constante. Ahora ilustramos el método de solución con un ejemplo.

Ejemplo\(\PageIndex{1}\)

Resolver\(\ddot{x}-3 \dot{x}-4 x=3 e^{2 t}\) con\(x(0)=1\) y\(\dot{x}(0)=0\).

Solución

Ver tutorial en YouTube

Primero, resolvemos la ecuación homogénea. La ecuación característica es

\[\begin{aligned} r^{2}-3 r-4 &=(r-4)(r+1) \\ &=0, \end{aligned} \nonumber \]

para que

\[x_{h}(t)=c_{1} e^{4 t}+c_{2} e^{-t} . \nonumber \]

Segundo, encontramos una solución particular de la ecuación no homogénea. La forma de la solución particular se elige de tal manera que lo exponencial cancelará de ambos lados de la oda. El ansatz que elegimos es

\[x(t)=A e^{2 t} \nonumber \]

donde\(A\) es un coeficiente aún indeterminado. Al sustituir\(x\) en la oda, diferenciando usando la regla de la cadena, y cancelando la exponencial, obtenemos

\[4 A-6 A-4 A=3 \nonumber \]

a partir de la cual determinamos\(A=-1 / 2\). Obtener una solución para\(A\) independiente de\(t\) justifica la Ecuación ansatz\ ref {8.24}. Tercero, escribimos la solución general a la oda como la suma de las soluciones homogéneas y particulares, y determinamos\(c_{1}\) y\(c_{2}\) que satisfacen las condiciones iniciales. Tenemos

\[x(t)=c_{1} e^{4 t}+c_{2} e^{-t}-\frac{1}{2} e^{2 t} ; \nonumber \]

y tomando el derivado,

\[\dot{x}(t)=4 c_{1} e^{4 t}-c_{2} e^{-t}-e^{2 t} . \nonumber \]

Aplicando las condiciones iniciales,

\[\begin{aligned} &c_{1}+c_{2}-\frac{1}{2}=1 \\ &4 c_{1}-c_{2}-1=0 \end{aligned} \nonumber \]

\[\begin{aligned} c_{1}+c_{2} &=\frac{3}{2} \\ 4 c_{1}-c_{2} &=1 \end{aligned} \nonumber \]

Este sistema de ecuaciones lineales puede resolverse sumando\(c_{1}\) las ecuaciones a obtener\(c_{1}=1 / 2\), después de lo cual se\(c_{2}=1\) pueden determinar a partir de la primera ecuación. Por lo tanto, la solución para\(x(t)\) eso satisface tanto la oda como las condiciones iniciales viene dada por

\[\begin{aligned} x(t) &=\frac{1}{2} e^{4 t}-\frac{1}{2} e^{2 t}+e^{-t} \\ &=\frac{1}{2} e^{4 t}\left(1-e^{-2 t}+2 e^{-5 t}\right) \end{aligned} \nonumber \]

donde hemos agrupado los términos en la solución para mostrar mejor el comportamiento asintótico para grandes\(t\).

Ahora encontramos soluciones particulares para algunos términos relativamente simples no homogéneos utilizando este método de coeficientes indeterminados.

Ejemplo\(\PageIndex{2}\)

Encuentre una solución particular de\(\ddot{x}-3 \dot{x}-4 x=2 \sin t\).

Solución

Ver tutorial en YouTube

Se muestran dos métodos para encontrar una solución particular. El primer método más directo prueba el ansatz

\[x(t)=A \cos t+B \sin t \nonumber \]

donde el argumento de coseno y seno debe estar de acuerdo con el argumento de seno en el término no homogéneo. Se requiere el término coseno porque la derivada del seno es coseno. Tras la sustitución en la ecuación diferencial, obtenemos

\[(-A \cos t-B \sin t)-3(-A \sin t+B \cos t)-4(A \cos t+B \sin t)=2 \sin t \nonumber \]

o términos de reagrupación,

\[-(5 A+3 B) \cos t+(3 A-5 B) \sin t=2 \sin t \nonumber \]

Esta ecuación es válida para todos\(t\), y en particular para\(t=0\) y\(\pi / 2\), para los cuales desaparecen las funciones seno y coseno. Para estos dos valores de\(t\), encontramos

\[5 A+3 B=0, \quad 3 A-5 B=2 ; \nonumber \]

y resolviendo para\(A\) y\(B\), obtenemos

\[A=\frac{3}{17}, \quad B=-\frac{5}{17} \nonumber \]

Por lo tanto, la solución particular viene dada por

\[x_{p}=\frac{1}{17}(3 \cos t-5 \sin t) . \nonumber \]

El segundo método de solución hace uso de la relación\(e^{i t}=\cos t+i \sin t\) para convertir el término seno no homogéneo en una función exponencial. Introducimos la función compleja\(z(t)\) dejando

\[z(t)=x(t)+i y(t) \nonumber \]

y reescribir la ecuación diferencial en forma compleja. Podemos reescribir la ecuación de una de dos maneras. Por un lado, si usamos\(\sin t=\operatorname{Re}\left\{-i e^{i t}\right\}\), entonces la ecuación diferencial se escribe como

\[\ddot{z}-3 \dot{z}-4 z=-2 i e^{i t} ; \nonumber \]

y al equiparar las partes real e imaginaria, esta ecuación se convierte en las dos ecuaciones diferenciales reales

\[\ddot{x}-3 \dot{x}-4 x=2 \sin t, \quad \ddot{y}-3 \dot{y}-4 y=-2 \cos t . \nonumber \]

La solución que estamos buscando, entonces, es\(x_{p}(t)=\operatorname{Re}\left\{z_{p}(t)\right\}\).

Por otro lado, si escribimos\(\sin t=\operatorname{Im}\left\{e^{i t}\right\}\), entonces la compleja ecuación diferencial se convierte en

\[\ddot{z}-3 \dot{z}-4 z=2 e^{i t} \nonumber \]

que se convierte en las dos ecuaciones diferenciales reales

\[\ddot{x}-3 \dot{x}-4 x=2 \cos t, \quad \ddot{y}-3 \dot{y}-4 y=2 \sin t . \nonumber \]

Aquí, la solución que estamos buscando es\(x_{p}(t)=\operatorname{Im}\left\{z_{p}(t)\right\}\).

Procederemos aquí resolviendo la Ecuación\ ref {8.26}. Como ahora tenemos una función exponencial como el término no homogéneo, podemos hacer que el ansatz

\[z(t)=C e^{i t}, \nonumber \]

donde ahora esperamos\(C\) que sea una constante compleja. Tras la sustitución en la oda Ecuación\ ref {8.26} y usando\(i^{2}=-1\):

\[-C-3 i C-4 C=2 \nonumber \]

o resolviendo para\(C\):

\[\begin{aligned} C &=\frac{-2}{5+3 i} \\ &=\frac{-2(5-3 i)}{(5+3 i)(5-3 i)} \\ &=\frac{-10+6 i}{34} \\ &=\frac{-5+3 i}{17} \end{aligned} \nonumber \]

Por lo tanto,

\[\begin{aligned} x_{p} &=\operatorname{Im}\left\{z_{p}\right\} \\ &=\operatorname{Im}\left\{\frac{1}{17}(-5+3 i)(\cos t+i \sin t)\right\} \\ &=\frac{1}{17}(3 \cos t-5 \sin t) \end{aligned} \nonumber \]

Ejemplo\(\PageIndex{3}\)

Encuentre una solución particular de\(\ddot{x}+\dot{x}-2 x=t^{2}\).

Solución

Ver tutorial en YouTube

El correcto ansatz aquí es el polinomio

\[x(t)=A t^{2}+B t+C . \nonumber \]

Al sustituirnos en la oda, tenemos

\[2 A+2 A t+B-2 A t^{2}-2 B t-2 C=t^{2} \nonumber \]

\[-2 A t^{2}+2(A-B) t+(2 A+B-2 C) t^{0}=t^{2} . \nonumber \]

Equiparar poderes de\(t\)

\[-2 A=1, \quad 2(A-B)=0, \quad 2 A+B-2 C=0 ; \nonumber \]

y resolviendo,

\[A=-\frac{1}{2}, \quad B=-\frac{1}{2}, \quad C=-\frac{3}{4} \nonumber \]

Por lo tanto, la solución particular es

\[x_{p}(t)=-\frac{1}{2} t^{2}-\frac{1}{2} t-\frac{3}{4} . \nonumber \]