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

¿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

Demostración ecuación diferencial lineal homogénea de segundo orden con Wronskiano

En esta ocasión, vengo a mostrarles la solución del problema 47 de ecuaciones diferenciales de segundo orden homogeneas del libro de  Ecuaciones diferenciales con aplicaciones 2da edición,  del autor  Dennis G. Zill. El problema es el siguiente: Sean $y_{1}$ y $y_{2}$ dos soluciones de: \[a_{2}(x)y''+a_{1}(x)y'+a_{0}y=0\] (a) Sí $W(y_{1},y_{2})$ es el wronskiano de $y_{1}$ y $y_{2}$, demuestre que \[a_{2}(x)\frac{dW}{dx}+a_{1}(x)W=0\] Como $y_{1}$ y $y_{2}$ son soluciones de la ecuación diferencial homogénea de segundo orden, luego también deben ser soluciones de las ecuaciones diferenciales: \[a_{2}(x)y_{1}''+a_{1}(x)y_{1}'+a_{0}y_{1}=0\] \[a_{2}(x)y_{2}''+a_{1}(x)y_{2}'+a_{0}y_{2}=0\] El wronskiano es de la forma: \[W(y_{1},y_{2})=\begin{vmatrix}y_{1} & y_{2}\\y_{1}' & y_{2}'\end{vmatrix}\] \[W(y_{1},y_{2})=y_{1}y_{2}'-y_{2}y_{1}'\] La derivada del Wronskiano es: \[\frac{dW}{dx}=(y_{1}y_{2}')'-(y_{2}y_{1}')'\...
Publicar un comentario
Leer más

¿Cómo hallar el centro de masa de una lámina con integrales dobles? 2

 A petición de un suscriptor y seguidor voy a mostrar dos ejercicios para el cálculo del centro de masa con integrales dobles, comencemos: El segundo ejercicio es: 2. Calcular el centro de masa de una lámina representada por la región $R$ que se encuentra por encima del eje $x$ y entre las líneas $y=x$; $y=-x$, $x^{2}+y^{2}=4y$; $x^{2}+y^{2}=6y$; $y>0$; donde la densidad es $\sqrt{x^{2}+y^{2}}$ Primero vamos a representar el área, al cuál vamos a hallar el centro de masa (voy a utilizar Geogebra). Este problema resulta más simple si lo resolvemos mediante coordenadas polares, así que la densidad $\rho(x,y)$ en coordenadas polares de acuerdo a las reglas de transformación $x=rcos\theta$ y $y=rsen\theta$ es: \[\rho(x,y)=\rho(rcos\theta,rsen\theta)=\sqrt{x^{2}+y^{2}}=\sqrt{(rcos\theta)^{2}+(rsen\theta)^{2}}=\sqrt{r^{2}}=r\] Como vimos en el primer ejercicio , las integrales que tiene cada punto del centro de masa se pueden representar en coordenadas polares, escogemos este tipo de ...
Publicar un comentario
Leer más