Ir al contenido principal

¿Cómo resolver una ecuación diferencial que involucra circuitos por el método de series de potencias, hallando sus constantes?

La ecuación diferencial a resolver es de la forma:
\[Lq''(t)+Rq'(t)+\frac{1}{C}q(t)=E(t)\]
Los valores de las constantes son los siguientes $L=\frac{1}{4}H$, $R(t)=\left(1+\frac{t}{8}\right)\Omega$, $C=4F$, $E(t)=0$, reemplazando dichos valores la ecuación diferencial a resolver es de la forma:
\[\frac{1}{4}q''(t)+\left(1+\frac{t}{8}\right)q'(t)+\frac{1}{4}q(t)=0\]
Se propone una solución en series de potencias y sus correspondientes derivadas de hasta segundo orden:
\[q(t)=\sum_{n=0}^{\infty}c_{n}t^{n}\]
\[q(t)'=\sum_{n=1}^{\infty}c_{n}nt^{n-1}\]
\[q(t)''=\sum_{n=2}^{\infty}c_{n}n(n-1)t^{n-2}\]
Reemplazando en nuestra ecuación diferencial tenemos:
\[\frac{1}{4}\sum_{n=2}^{\infty}c_{n}n(n-1)t^{n-2}+\left(1+\frac{t}{8}\right)\sum_{n=1}^{\infty}c_{n}nt^{n-1}+\frac{1}{4}\sum_{n=0}^{\infty}c_{n}t^{n}=0\]
Distribuyendo para la segunda sumatoria tenemos:
\[\frac{1}{4}\sum_{n=2}^{\infty}c_{n}n(n-1)t^{n-2}+\sum_{n=1}^{\infty}c_{n}nt^{n-1}+\frac{t}{8}\sum_{n=1}^{\infty}c_{n}nt^{n-1}+\frac{1}{4}\sum_{n=0}^{\infty}c_{n}t^{n}=0\]
Multiplicando el término $t$ que le aporta a la tercera sumatoria tendremos:
\[\frac{1}{4}\sum_{n=2}^{\infty}c_{n}n(n-1)t^{n-2}+\sum_{n=1}^{\infty}c_{n}nt^{n-1}+\frac{1}{8}\sum_{n=0}^{\infty}c_{n}nt^{n}+\frac{1}{4}\sum_{n=0}^{\infty}c_{n}t^{n}=0\]
Hacemos $k=n-2$ y $n=k+2$ para la primera sumatoria, $k=n-1$ y $n=k+1$ para la segunda sumatoria y $k=n$ para la tercera y cuarta sumatoria, reemplazamos y tenemos:
\[\frac{1}{4}\sum_{k=0}^{\infty}c_{k+2}(k+2)(k+1)t^{k}+\sum_{k=0}^{\infty}c_{k+1}(k+1)t^{k}+\frac{1}{8}\sum_{k=0}^{\infty}c_{k}kt^{k}+\frac{1}{4}\sum_{k=0}^{\infty}c_{k}t^{k}=0\]
Que se puede escribir de la forma:
\[\sum_{k=0}^{\infty}\left[\frac{1}{4}c_{k+2}(k+2)(k+1)+c_{k+1}(k+1)+\frac{1}{8}c_{k}k+\frac{1}{4}c_{k}\right]t^{k}=0\]
Luego:
\[\sum_{k=0}^{\infty}\underbrace{\left[\frac{1}{4}c_{k+2}(k+2)(k+1)+c_{k+1}(k+1)+\frac{1}{8}c_{k}k+\frac{1}{4}c_{k}\right]}_{=0}\underbrace{t^{k}}_{\neq0}=0\]
Al despejar la constante de mayor valor en términos de las constantes de valor tendremos:
\[\frac{1}{4}c_{k+2}(k+2)(k+1)+c_{k+1}(k+1)+\frac{1}{8}c_{k}k+\frac{1}{4}c_{k}=0\]
\[\frac{1}{4}c_{k+2}(k+2)(k+1)+c_{k+1}(k+1)+\frac{1}{4}c_{k}\left[\frac{k}{2}+1\right]=0\]
\[\frac{1}{4}c_{k+2}(k+2)(k+1)=-c_{k+1}(k+1)-\frac{1}{4}c_{k}\left[\frac{k}{2}+1\right]\]
\[c_{k+2}=-\frac{4(k+1)c_{k+1}}{(k+2)(k+1)}-\frac{4}{4}\frac{c_{k}}{(k+2)(k+1)}\left[\frac{k}{2}+1\right]\]
\[c_{k+2}=-\frac{4(k+1)c_{k+1}}{(k+2)(k+1)}-\frac{4}{4}\frac{c_{k}}{(k+2)(k+1)}\left[\frac{k+2}{2}\right]\]
\[c_{k+2}=-\frac{4c_{k+1}}{(k+2)}-\frac{c_{k}}{2(k+1)}\]
Le damos valores desde $k=0$ hasta el valor $k=3$ dado que se piden calcular los primeros $5$ términos de cada solución, por lo tanto:
Para $k=0$:
\[c_{0+2}=-\frac{4c_{0+1}}{(0+2)}-\frac{c_{0}}{2(0+1)}\]
\[c_{2}=-\frac{4c_{1}}{2}-\frac{c_{0}}{2(1)}\]
\[c_{2}=-\frac{4c_{1}}{2}-\frac{c_{0}}{2}\]
\[c_{2}=-2c_{1}-\frac{c_{0}}{2}\]
Para $k=1$:
\[c_{1+2}=-\frac{4c_{1+1}}{(1+2)}-\frac{c_{1}}{2(1+1)}\]
\[c_{3}=-\frac{4c_{2}}{3}-\frac{c_{1}}{2(2)}\]
\[c_{3}=-\frac{4c_{2}}{3}-\frac{c_{1}}{4}\]
\[c_{3}=-\frac{4}{3}\left(-2c_{1}-\frac{c_{0}}{2}\right)-\frac{c_{1}}{4}\]
\[c_{3}=\frac{8c_{1}}{3}+\frac{4c_{0}}{6}-\frac{c_{1}}{4}\]
\[c_{3}=\frac{8c_{1}}{3}-\frac{c_{1}}{4}+\frac{4c_{0}}{6}\]
\[c_{3}=\frac{32c_{1}}{12}-\frac{3c_{1}}{12}+\frac{4c_{0}}{6}\]
\[c_{3}=\frac{29c_{1}}{12}+\frac{2c_{0}}{3}\]
Para $k=2$:
\[c_{2+2}=-\frac{4c_{2+1}}{(2+2)}-\frac{c_{2}}{2(2+1)}\]
\[c_{4}=-\frac{4c_{3}}{4}-\frac{c_{2}}{2(3)}\]
\[c_{4}=-\frac{4c_{3}}{4}-\frac{c_{2}}{6}\]
\[c_{4}=-\frac{4}{4}\left(\frac{29c_{1}}{12}+\frac{2c_{0}}{3}\right)-\frac{1}{6}\left(-2c_{1}-\frac{c_{0}}{2}\right)\]
\[c_{4}=-\frac{29c_{1}}{12}-\frac{2c_{0}}{3}+\frac{4c_{1}}{12}+\frac{c_{0}}{12}\]
\[c_{4}=-\frac{29c_{1}}{12}-\frac{4c_{1}}{12}+\frac{2c_{0}}{3}+\frac{c_{0}}{12}\]
\[c_{4}=-\frac{29c_{1}}{12}-\frac{4c_{1}}{12}+\frac{8c_{0}}{12}+\frac{c_{0}}{12}\]
\[c_{4}=-\frac{33c_{1}}{12}+\frac{9c_{0}}{12}\]
\[c_{4}=-\frac{11c_{1}}{4}+\frac{3c_{0}}{4}\]
Para $k=3$:
\[c_{3+2}=-\frac{4c_{3+1}}{(3+2)}-\frac{c_{3}}{2(3+1)}\]
\[c_{5}=-\frac{4c_{4}}{5}-\frac{c_{3}}{2(4)}\]
\[c_{5}=-\frac{4c_{4}}{5}-\frac{c_{3}}{8}\]
\[c_{5}=-\frac{4}{5}\left(-\frac{11c_{1}}{4}+\frac{3c_{0}}{4}\right)-\frac{1}{6}\left(\frac{29c_{1}}{12}+\frac{2c_{0}}{3}\right)\]
\[c_{5}=\frac{100c_{1}}{60}-\frac{12c_{0}}{20}-\frac{29c_{1}}{96}-\frac{4c_{0}}{42}\]
\[c_{5}=\frac{10c_{1}}{6}-\frac{6c_{0}}{10}-\frac{29c_{1}}{96}-\frac{2c_{0}}{21}\]
\[c_{5}=\frac{5c_{1}}{3}-\frac{3c_{0}}{5}-\frac{29c_{1}}{96}-\frac{2c_{0}}{21}\]
\[c_{5}=\frac{5c_{1}}{3}-\frac{29c_{1}}{96}-\frac{3c_{0}}{5}-\frac{2c_{0}}{21}\]
\[c_{5}=\frac{160c_{1}}{96}-\frac{29c_{1}}{96}-\frac{63c_{0}}{105}-\frac{10c_{0}}{105}\]
\[c_{5}=\frac{131c_{1}}{96}-\frac{39c_{0}}{105}\]
\[c_{5}=\frac{131c_{1}}{96}-\frac{13c_{0}}{35}\]
Luego la solución en serie de potencias corresponde a:
\[q(t)=c_{0}+c_{1}t+c_{2}t^{2}+c_{3}t^{3}+c_{4}t^{4}+c_{5}t^{5}+...\]
Reemplazamos las constantes y tendremos nuestra solución expresada en términos de $c_{0}$ y $c_{1}$:
\[q(t)=c_{0}+c_{1}t+\left(-2c_{1}-\frac{c_{0}}{2}\right)t^{2}+\left(\frac{29c_{1}}{12}+\frac{2c_{0}}{3}\right)t^{3}+\left(-\frac{11c_{1}}{4}+\frac{3c_{0}}{4}\right)t^{4}+\left(\frac{131c_{1}}{96}-\frac{13c_{0}}{35}\right)t^{5}+...\]
Factorizamos en terminos de $c_{0}$ y $c_{1}$ y tendremos:
\[q(t)=c_{0}\left(1-\frac{t^{2}}{2}+\frac{2t^{3}}{3}+\frac{3t^{4}}{4}-\frac{13t^{5}}{35}+...\right)+c_{1}\left(t-2t^{2}+\frac{29t^{3}}{12}-\frac{11t^{4}}{4}+\frac{131t^{5}}{96}+...\right)\]
Aplicando condiciones iniciales en $q(0)=2$ tendremos:
\[q(0)=c_{0}\left(1-\frac{(0)^{2}}{2}+\frac{2(0)^{3}}{3}+\frac{3(0)^{4}}{4}-\frac{13(0)^{5}}{35}+...\right)+c_{1}\left((0)-2(0)^{2}+\frac{29(0)^{3}}{12}-\frac{11(0)^{4}}{4}+\frac{131(0)^{5}}{96}+...\right)\]
\[q(0)=2=c_{0}(1)+c_{1}(0)\]
Ahora cuando $\frac{dq(0)}{dt}=0$ lo que equivale a derivar la expresión anterior:
\[\frac{dq(t)}{dt}=2\left(-t-\frac{3t^{2}}{3}+\frac{12t^{3}}{4}+\frac{215t^{4}}{105}+...\right)+c_{1}\left(1-4t+\frac{87t^{2}}{12}-\frac{44t^{3}}{4}+\frac{655t^{4}}{96}+...\right)\]
\[\frac{dq(t)}{dt}=0=2\left(-(0)-\frac{3(0)^{2}}{3}+\frac{12(0)^{3}}{4}+\frac{215(0)^{4}}{105}+...\right)+c_{1}\left(1-4(0)+\frac{87(0)^{2}}{12}-\frac{44(0)^{3}}{4}+\frac{655(0)^{4}}{96}+...\right)\]
\[0=2c_{1}(1)\]
\[0=c_{1}\]
Luego la correspondiente solución a mi ecuación diferencial es:
\[q(t)=2\left(1-\frac{t^{2}}{2}+\frac{2t^{3}}{3}+\frac{3t^{4}}{4}-\frac{13t^{5}}{35}+...\right)\]
De acuerdo a las condiciones iniciales dadas.

Comentarios

Entradas populares de este blog

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}')'\...

¿Cómo son las integrales del tipo arcoseno o arcocoseno?

Del curso de Cálculo Integral aprendemos que las integrales de tipo arcoseno u arcocoseno son de la forma: \[\int \frac{du}{\sqrt{1-u^{2}}}=arcsen(u)\] ó \[\int -\frac{du}{\sqrt{1-u^{2}}}=arccos(u)\] Hasta es posible ver dicha integral de un poco más complicada sin importar que sea positiva (arcoseno) u negativa (arcocoseno): \[\int \pm \frac{f'(u)du}{\sqrt{1-[f(u)]^{2}}}\] Sólo en esos casos podemos conocer algunas integrales que podemos resolver, pero existen otras integrales que con un cambio de variable, u organización de términos especifico, puede darnos en términos de arcoseno u arcocoseno, un ejemplo puede ser el siguiente: \[\int \frac{du}{\sqrt{a^{2}-u^{2}}}\] Acá simplemente dentro de la raíz cuadrada debemos dejarlo de la forma: 1-término al cuadrado para que nos quede fácil de identificar, y así se pueda hacer fácil la integración: \[\int \frac{du}{\sqrt{a^{2}\left(1-\frac{u^{2}}{a^{2}}\right)}}\] El término $a^{2}$, sale de la raíz cuadrada como $a$,  \[\frac{1}{a}\int...

¿Cuál es la integral de arcotangente?

Nuestro problema ahora es hallar la integral a la función arcotangente: \[\int arctan(w)dw\] Debemos aplicar integración por partes: Uniforme de Vaca es para Una Vaca Vestida de Uniforme : \[\int UdV=UV-\int VdU\] Acá $U=arctan(w)$, su derivada es: $dU=\frac{1}{1+w^{2}}$, y $dV=dw$, su integral es: $V=w$, aplicamos de acuerdo a la integración por partes: \[\int arctan(w)dw=w \cdot arctan(w)-\int \frac{w}{1+w^{2}}dw\] Para la última integral utilizamos una sustitución, sí $p=1+w^{2}$, $dp=2wdw$, entonces $\frac{dp}{2}=wdw$, y la respectiva integral queda como: \[\frac{1}{2}\int \frac{dp}{p}=\frac{1}{2}ln(p)+C\] Deshacemos la sustitución: \[\frac{1}{2}ln(1+w^{2})+C\] Reemplazamos este resultado para obtener el resultado de la integral arcotangente: \[\int arctan(w)dw=w \cdot arctan(w)-\frac{1}{2}ln(1+w^{2})+C\] Finalmente obteniendo el resultado deseado.