Ir al contenido principal

¿Cómo resolver una ecuación diferencial no homogénea por el método de variación de parámetros? 1

Tenemos la ecuación diferencial no homogénea:
\[y''+9y=sec(x)\]
Primero hallamos la solución complementaria $y_{c}$, que se obtiene de resolver la ecuación diferencial homogénea:
\[y''+9y=0\]
Proponemos una solución de la forma $y=e^{mx}$ y sus derivadas $y'=me^{mx}$, $y''=m^{2}e^{mx}$, reemplazamos en nuestra ecuación diferencial:
\[m^{2}e^{mx}+9e^{mx}=0\]
Factorizamos $e^{mx}$ y como este ultimo no puede ser cero, la ecuación cuadrática a resolver para determinar $m$ es de la forma:
\[e^{mx}(m^{2}+9)=0\]
\[m^{2}+9=0\]
Las soluciones a esa ecuación cuadrática son de la forma:
\[m_1=i3 \quad m_2=-i3\]
La solución particular es de la forma:
\[y_c=C_1e^{i3x}+C_2e^{-i3x}\]
Aplicamos la identidad de Euler y obtenemos:
\[y_c=c_1cos(3x)+c_2sen(3x)\]
Para hallar la solución particular $y_p$ utilizamos el método de variación de parámetros, que emplea el Wronskiano:
\[W(f(x),g(x))=\begin{vmatrix}f(x) & g(x)\\f'(x) & g'(x)\end{vmatrix}\]
Identificamos las funciones $f(x)=cos(3x)$ y $g(x)=sen(3x)$, hallamos las derivadas $f'(x)=-3sen(3x)$ y $g'(x)=3cos(3x)$, reemplazamos en nuestro Wronskiano:
\[W(cos(3x),sen(3x))=\begin{vmatrix}cos(3x) & sen(3x)\\-3sen(3x) & 3cos(3x)\end{vmatrix}\]
Si el determinante es diferente de cero, luego las soluciones son linealmente independientes:
\[W(cos(3x),sen(3x))=\begin{vmatrix}cos(3x) & sen(3x)\\ -3sen(3x) & 3cos(3x)\end{vmatrix}\]
\[=(cos(3x))(3cos(3x))-(sen(3x))(-3sen(3x))\]
\[=3cos^{2}(x)+3sen^{2}(x)=3(cos^{2}(3x)+sen^{2}(3x))\]
\[=3\]
Luego el determinante nos da diferente de cero.
La variación de parámetros nos dice que para hallar la solución $y_p$ debemos integrar los términos:
\[u'_1=-\frac{y_2h(x)}{W} \quad u'_2=-\frac{y_1h(x)}{W}\]
Donde $y_1=cos(3x)$ y $y_2=sen(3x)$ y $h(x)=sec(x)$ (El término que hace no homogénea a nuestra ecuación diferencial $y''+9y=sec(x)$ y $W$ es el resultado del determinante (Wronskiano)
Calculamos $u'_1$ y $u'_2$ :
\[u'_1=-\frac{sen(3x)sec(x)}{3} \quad u'_2=\frac{cos(3x)sec(x)}{3}\]
Integramos cada una de las $u_n$ anteriores respecto a $x$, para eso vamos a tener la forma:
\[u'_1=-\frac{sen(3x)}{3cos(x)} \quad u'_2=\frac{cos(3x)}{3cos(x)}\]
Integramos:
\[u_1=\int -\frac{sen(3x)}{3cos(x)}dx \quad u_2=\int \frac{cos(3x)}{3cos(x)}dx\]
Y obtenemos las soluciones de las integrales:
\[u_1=-\frac{1}{3}ln(|cos(x)|) -\frac{2}{3}cos^{2}(x)\]
\[u_2=-\frac{1}{3}sen(2x)-\frac{1}{3}x\]
La solución particular es de la forma:
\[y_p=u_1y_1+u_2y_2\]
\[y_p=\left(-\frac{1}{3}ln(|cos(x)|) -\frac{2}{3}cos^{2}(x)\right)cos(3x)+\left(-\frac{1}{3}sen(2x)-\frac{1}{3}x\right)sen(3x)\]
\[y_p=-\left(\frac{1}{3}ln(|cos(x)|) cos(3x)+\frac{2}{3}cos^{2}(x)cos(3x)+\frac{1}{3}sen(2x)sen(3x)+\frac{1}{3}xsen(3x)\right)\]
Y la solución de la ecuación diferencial no homogénea es:
\[y=y_c+y_p\]
\[y=c_1cos(3x)+c_2sen(3x)-\left(\frac{1}{3}ln(|cos(x)|) cos(3x)+\frac{2}{3}cos^{2}(x)cos(3x)+\frac{1}{3}sen(2x)sen(3x)+\frac{1}{3}xsen(3x)\right)\]

Etiquetas:

Comentarios

Publicar un comentario

Entradas populares de este blog

Problema de probabilidad 1

A un mono se le dan 9 bloques: 3 en forma de cuadrados ,3 como rectángulos y 3 como triángulos. si saca tres de cada clase en orden, es decir, tres triángulos, luego la misma cantidad de cuadrados y así sucesivamente. ¿ sospecharía usted que el mono asoció figuras que tengan forma idéntica?. calcule la probabilidad de este evento. Como tenemos 3 figuras, debemos tener en cuenta todas las posibles formas de organizarlas, Luego una forma de representar este orden es de la siguiente manera: Tomamos $C$ como Cuadrado, $R$ como Rectángulo, y $T$ como Triángulo \[CRT \quad RTC \quad TCR\] \[CTR \quad RCT \quad TRC\] Que corresponden a $6$ formas de organizar las tres figuras, este procedimiento se puede representar fácilmente con la operación factorial, para nuestro caso queda definido como: \[3!=3 \times 2 \times 1 = 6\] Luego por cada vez que el mono saca tres bloques, estos vendrán organizados de $3!$ formas, y como son 3 tandas de veces que el mono va a sacar cantidades
Publicar un comentario
Leer más

¿Cómo resolver la ecuación diferencial de un oscilador armónico simple con gravedad?

En este caso vamos a solucionar la ecuación diferencial para un oscilador armónico unidimensional que solo tiene el movimiento en el eje y, y además está en presencia de la gravedad, su ecuación diferencial es la siguiente: \[\ddot{y}+\omega^{2}y=g\] Esta es una ecuación diferencial no homogénea y la vamos a resolver por el método de variación de parámetros. Para este método necesitamos resolver la ecuación diferencial homogenea: \[\ddot{y}+\omega^{2}y=0\] Que ya la hemos resuelto en este enlace  aunque resuelta con la coordenada $x$, pero eso no va a importar, porque va a ser la misma solución, solo que cambiamos de $x$ a $y$: \[y(t)=c_{1}cos(\omega t)+c_{2}sen(\omega t)\] Para resolver la ecuación diferencial debemos hallar el wronskiano: \[W(f(t),g(t))=\begin{vmatrix}f(t) & g(t)\\f'(t) & g'(t)\end{vmatrix}\] Identificamos las funciones $f(t)=c_{1}cos(\omega t)$ y $g(t)=c_{2}sen(\omega t)$, hallamos las derivadas $f'(t)=-c_{1}\omega sen(\omega t)$ y $g'(t)=c_{
Publicar un comentario
Leer más

¿Cómo realizar una sustitución para integrar una función en términos de arcoseno?

La integral que queremos resolver es la siguiente: \[\sqrt{\frac{m}{2E}}\int \frac{dx}{\sqrt{1-\frac{1}{2}\frac{k}{E}x^{2}}}\] Realizamos un cambio de variable de la siguiente forma: \[\frac{k}{2E}x^{2}=y^{2}\] Para poder hallar la derivada sacamos raíz cuadrada a ambos términos: \[\sqrt{\frac{k}{2E}}x=y\] Sacamos las derivadas y obtenemos: \[\sqrt{\frac{k}{2E}}dx=dy\] Luego pasando el factor que multiplica el $dx$ al lado del dy: \[dx=\sqrt{\frac{2E}{k}}dy\] Por tanto la integral queda de la forma: \[\sqrt{\frac{m}{2E}}\sqrt{\frac{2E}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}\] \[\sqrt{\frac{m}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}\] Que corresponde a una integral de arcoseno, por lo tanto la respuesta nos queda como: \[\sqrt{\frac{m}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}=\sqrt{\frac{m}{k}}arcsen(y)+C\] Deshaciendo el cambio de variable, obtenemos finalmente el resultado de nuestra integral: \[\sqrt{\frac{m}{2E}}\int \frac{dx}{\sqrt{1-\frac{1}{2}\frac{k}{E}x^{2}}}=\sqrt{\frac{m}{k}}arcse
Publicar un comentario
Leer más