7 votos

Encontrar las soluciones de la ecuación diferencial $3x^2y''+5xy'+3xy=0$.

Encontrar todas las soluciones de la forma $y(x)= x^m \sum_{n=0}^{\infty} a_nx^n, \ x>0 (m \in \mathbb{R})$ de la ecuación diferencial $3x^2y''+5xy'+3xy=0$.

Eso es lo que he intentado:

Desde $x>0$ la ecuación diferencial puede escribirse como sigue.

$$y''+ \frac{5}{3x}y'+ \frac{1}{x}y=0$$

$$p(x)=\frac{5}{3x}, q(x)= \frac{1}{x}$$

El punto de $0$ es regular singular, es decir, las funciones $xp(x), x^2q(x)$ puede ser escrito como la potencia de la serie en una región de $0$.

Buscamos soluciones de la forma $y(x)=x^m \sum_{n=0}^{\infty} a_n x^n$ adecuado $m \in \mathbb{R}$ y para el adecuado $a_n \in \mathbb{R}$ $x \in (0,R)$ donde $R$ es un número positivo.

Entonces tenemos:

$$y'(x)= \sum_{n=0}^{\infty} a_n (n+m) x^{n+m-1} \\ \Rightarrow xy'=\sum_{n=0}^{\infty} a_n (n+m) x^{n+m} \\ \Rightarrow 5xy'=\sum_{n=0}^{\infty} 5a_n (n+m) x^{n+m}$$

y

$$y''(x)= \sum_{n=0}^{\infty} a_n (n+m)(n+m-1) x^{n+m-2} \\ \Rightarrow x^2y''=\sum_{n=0}^{\infty} a_n (n+m)(n+m-1) x^{n+m} \\ \Rightarrow 3x^2y''=\sum_{n=0}^{\infty} 3a_n (n+m)(n+m-1) x^{n+m}$$

$$3xy= \sum_{n=0}^{\infty} 3a_n x^{n+m+1}$$

Así que tiene que llevar a cabo las siguientes acciones:

$$\sum_{n=0}^{\infty} \left[ 3a_n (n+m)(n+m-1) x^{n+m}+5a_n (n+m) x^{n+m}+ 3a_n x^{n+m+1}\right]=0 \Rightarrow \sum_{n=0}^{\infty} \left[ 3a_n (n+m)(n+m-1) +5a_n (n+m) + 3a_n x\right]x^{n+m}=0 $$

Por lo que tiene que sostener que:

$$a_n=\frac{-3a_{n-1}}{3m+3n+2}, \forall n=1,2,3, \dots$$

EDIT: estoy mirando de nuevo al ejercicio. Para $m=0$ tengo las siguientes:

$$a_1=-\frac{3a_0}{5} \\ a_2=\frac{3^2 a_0}{5 \cdot 8} \\ a_3=-\frac{3^3 a_0}{5 \cdot 8 \cdot 11} \\ a-4= \frac{3^4 a_0}{5 \cdot 8 \cdot 11 \cdot 14}$$

Así que no es de $m=0$ la fórmula general para $a_n$ el siguiente? $$$$

$$a_n=(-1)^n \frac{a_0}{\prod_{i=0}^{n-1} (3i+5)}$$

Y para $m=-\frac{3}{2}$ no es la fórmula para $a_n$ el siguiente?

$$a_n=(-1)^n \frac{3^n a_0}{ \prod_{i=0}^{n-1} \frac{(6i+1)}{2}}$$

Si es así, entonces podríamos decir lo siguiente?

$$y_1(x)= x^0 \sum_{n=0}^{\infty} \frac{(-1)^n x^n}{\prod_{i=0}^{n-1}(3i+5)}$$

y

$$y_2(x)=x^{-\frac{3}{2}} \sum_{n=0}^{\infty} \frac{(-1)^n 3^n}{ \prod_{i=0}^{n-1} \frac{(6i+1)}{2}}x^n$$

son soluciones de la ecuación diferencial para $a_0=1$.

$$\left| \frac{\frac{(-1)^{n+1} x^{n+1}}{\prod_{i=0}^n (3i+5)}}{\frac{(-1)^{n} x^n}{\prod_{i=0}^{n-1} (3i+5)}}\right|=\left| \frac{x}{3n+5}\right| \to 0<1$$

Qué podemos deducir de esto último que el radio de convergencia es $+\infty$. Si es así, ¿podemos continuar de la manera sig?

De manera similar se demuestra que el radio de convergencia de $y_2(x)$$+\infty$.

$$$$

$y_1, y_2$ son linealmente independientes en $(0,+\infty)$. Porque si $c_1, c_2 \in \mathbb{R}$ $c_1y_1(x)+c_2y_2(x)=0 \forall x \in (0,+\infty)$ desde $c_1 y_1(x)+ c_2y_2(x)$ es una potencia de la serie con radio de convergencia $+\infty$ tenemos $0= c_1 y_1(x)+c_2y_2(x)= \sum_{n=0}^{\infty} d_n x^n$ algunos $d_n \in \mathbb{R}$ e lo $d_n=0 \forall n=0,1,2, \dots$

Sin embargo $d_0=c_1=0$$d_1=-\frac{3}{5} c_2=0 \Rightarrow c_2=0$.

Por lo tanto, la solución general de la ecuación de sincronización es:

$$y(x)=c_1 \sum_{n=0}^{\infty} \frac{(-1)^n x^n}{\prod_{i=0}^{n-1}(3i+5)}+ c_2 x^{-\frac{3}{2}} \sum_{n=0}^{\infty} \frac{(-1)^n 3^n}{ \prod_{i=0}^{n-1} \frac{(6i+1)}{2}}x^n, c_1, c_2 \in \mathbb{R}$$

EDITAR: Me volvió a hacer los cálculos de $m=0$ y ahora tengo las siguientes: $$$$ Para $n=1$: $a_1=-\frac{3a_0}{1 \cdot 5}$

$$$$ Para $n=2$: $a_2=\frac{3^2 a_0}{2 \cdot 5 \cdot 8}$ $$$$ Para $n=3$: $a_3=-\frac{3^3 a_0}{2 \cdot 3 \cdot 5 \cdot 8 \cdot 11}$ $$$$ Para $n=4$: $a_4=\frac{3^4 a_0}{2 \cdot 3 \cdot 4 \cdot 5 \cdot 8 \cdot 11 \cdot 14}$ $$$$

Son mi actual cálculos derecho o fueron a la anterior, ¿correcto?

$$$$ Si están en lo correcto ¿cómo podemos escribir la fórmula de $a_n$ sin el uso de la función Gamma?

EDIT:me vuelve a intentar de nuevo. No podríamos escribir la fórmula general para $a_n$ al $m=0$ como sigue? $$$$ $$a_n=\frac{(-1)^n a_0}{n! \prod_{i=1}^n (3i+2)}$$ O estoy equivocado?

También debería ser $m_2=-\frac{2}{3}$. O estoy equivocado?

Si es que como que no es la fórmula general para $a_n$ en este caso el siguiente?

$$a_n=(-1)^n \frac{a_0}{n! \prod_{i=0}^{n-2} (2 \cdot 2+3 \cdot i)}$$

O estoy equivocado?

10voto

Roger Hoover Puntos 56

Sugerencia: si se reemplaza $y(x)$$x^{-1/3}g(2\sqrt{x})$, $g$ satisface una ecuación diferencial de Bessel.

4voto

andy.holmes Puntos 518

Usted necesita llevar a cabo un índice de cambio a fin de que todos los poderes de $x$ tienen el mismo exponente en cada término. $$ 3xy= \sum_{n=0}^{\infty} 3a_n x^{n+m+1}= \sum_{n=1}^{\infty} 3a_{n-1} x^{n+m}= \sum_{n=0}^{\infty} 3a_{n-1} x^{n+m} $$ La última tiene que definir un auxiliar del coeficiente de $a_{-1}=0$. Si no te gusta coeficientes negativos índice, usted tiene que tratar el plazo para $n=0$ por separado.

4voto

Markus Scheuer Puntos 16133

Creo que su elaboración de su pregunta, junto con algunos de sus comentarios es muy cerca a una solución correcta.

Vamos a calcular el $y_1(x)$.

Empezamos con $a_n$ que se ha especificado correctamente en uno de sus comentarios:

\begin{align*} a_n=\frac{-3a_{n-1}}{(m+n)(3m+3n+2)}\qquad\forall n=1,2,3,\ldots \end{align*} Establecimiento $m=0$ obtenemos \begin{align*} a_n&=\frac{-3a_{n-1}}{n(3n+2)}\\ &=\frac{-a_{n-1}}{n\left(n+\frac{2}{3}\right)}\\ &=\frac{a_{n-2}}{n(n-1)\left(n+\frac{2}{3}\right)\left(n-\frac{1}{3}\right)}\\ &=\ldots\\ &=\frac{(-1)^na_0}{n!\left(n+\frac{2}{3}\right)\left(n+\frac{2}{3}-1\right)\cdots\left(n+\frac{2}{3}-(n-1)\right)}\\ &=\frac{(-1)^na_0}{n!{\left(n+\frac{2}{3}\right)}_n}\tag{1}\\ &=(-1)^n\frac{a_0}{n!}\frac{\Gamma\left(\frac{5}{3}\right)}{\Gamma\left(n+\frac{5}{3}\right)}\tag{2} \end{align*}

En (1) se aplica el símbolo de Pochhammer $(x)_n=x(x-1)\cdots (x-n+1)$ y en (2) utilizamos la identidad $$(x)_n=\frac{\Gamma(x+1)}{\Gamma(x-n+1)}$$

Así, encontramos la solución

$$y_1(x)=a_0\Gamma\left(\frac{5}{3}\right)\sum_{n=0}^{\infty}\frac{(-1)^n}{\Gamma\left(n+\frac{5}{3}\right)}\frac{x^n}{n!}$$

De acuerdo a un comentario de @JackDAurizio vamos a echar un vistazo a la función de Bessel de primera especie

\begin{align*} J_\alpha(x)=\sum_{n=0}^{\infty}\frac{(-1)^n}{n!\Gamma(n+\alpha+1)}\left(\frac{x}{2}\right)^{2n+\alpha} \end{align*}

Observamos

\begin{align*} x^{-\frac{1}{3}}J_{\frac{2}{3}}(2\sqrt{x})=x^{-\frac{1}{3}}\sum_{n=0}^{\infty}\frac{(-1)^n}{n!\Gamma\left(n+\frac{5}{3}\right)}x^{n+\frac{1}{3}} =\sum_{n=0}^{\infty}\frac{(-1)^n}{n!\Gamma\left(n+\frac{5}{3}\right)}x^{n} \end{align*}

Llegamos a la conclusión de

\begin{align*} y_1(x)=a_0\Gamma\left(\frac{5}{3}\right)x^{-\frac{1}{3}}J_{\frac{2}{3}}(2\sqrt{x}) \end{align*}

Del mismo modo podemos encontrar $$y_2(x)=a_0\Gamma\left(\frac{1}{3}\right)x^{-\frac{1}{3}}J_{-\frac{2}{3}}(2\sqrt{x})$$

y la solución general es, por tanto, $$y(x)=c_1\Gamma\left(\frac{5}{3}\right)x^{-\frac{1}{3}}J_{\frac{2}{3}}(2\sqrt{x})+c_2\Gamma\left(\frac{1}{3}\right)x^{-\frac{1}{3}}J_{-\frac{2}{3}}(2\sqrt{x})$$


[2015-06-04]: Nota de acuerdo a la OPs comentario

Evitar la función Gamma $\Gamma(x)$ una representación a través del símbolo de Pochhammer (1) es conveniente:

$$a_n=\frac{(-1)^na_0}{n!{\left(n+\frac{2}{3}\right)}_n}$$

Observamos para $n=3$:

\begin{align*} a_3&=-\frac{3^3a_0}{2\cdot3\cdot5\cdot8\cdot11}\\ &=-\frac{3^3a_0}{3!\cdot5\cdot8\cdot11}\\ &=-\frac{a_0}{3!\cdot\frac{5}{3}\cdot\frac{8}{3}\cdot\frac{11}{3}}\\ &=-\frac{a_0}{3!\cdot\left(3-\frac{4}{3}\right)\left(3-\frac{1}{3}\right)\left(3+\frac{2}{3}\right)}\\ &=-\frac{a_0}{3!\left(3+\frac{2}{3}\right)_3} \end{align*} De manera similar el caso de $n=4$: \begin{align*} a_4&=\frac{3^4a_0}{2\cdot3\cdot4\cdot5\cdot8\cdot11\cdot14}\\ &=\frac{3^4a_0}{4!\cdot5\cdot8\cdot11\cdot14}\\ &=\frac{a_0}{4!\cdot\frac{5}{3}\cdot\frac{8}{3}\cdot\frac{11}{3}\cdot\frac{14}{3}}\\ &=\frac{a_0}{4!\left(4-\frac{7}{3}\right)\left(4-\frac{4}{3}\right)\left(4-\frac{1}{3}\right)\left(4+\frac{2}{3}\right)}\\ &=\frac{a_0}{4!\left(4+\frac{2}{3}\right)_4} \end{align*}

Obtenemos la siguiente representación para $y_1(x)$: \begin{align*} y_1(x)&=x^m\sum_{n=0}^{\infty}a_nx^n\\ &=x^0\sum_{n=0}^{\infty}\frac{(-1)^na_0}{n!{\left(n+\frac{2}{3}\right)}_n}x^n\\ &=a_0\sum_{n=0}^{\infty}\frac{(-1)^n}{{\left(n+\frac{2}{3}\right)}_n}\frac{x^n}{n!} \end{align*}


[2015-06-04]: nota Adicional de acuerdo a la OPs comentarios

Desde

\begin{align*} \prod_{j=1}^{n}(3j+2)&=5\cdot8\cdot11\cdot\ldots\cdot(3n+2)\\ &=(3\cdot1+2)(3\cdot2+2)(3\cdot3+2)\cdot\ldots\cdot(3n+2)\\ &=3^n\left(1+\frac{2}{3}\right)\left(2+\frac{2}{3}\right)\left(3+\frac{2}{3}\right)\cdot\ldots\cdot\left(n+\frac{2}{3}\right)\\ &=3^n\left(n+\frac{2}{3}\right)_n \end{align*} es perfectamente válido, para escribir \begin{align*} a_n=\frac{(-1)^n3^na_0}{n!\prod_{j=1}^{n}(3j+2)}=\frac{(-1)^na_0}{n!\left(n+\frac{2}{3}\right)_n} \end{align*}

Nota: por lo general, es más conveniente utilizar el símbolo de Pochhammer, desde consecutivos de los factores que tienen el bonito diferencia $1$ y no tenemos necesidad de escribir el producto símbolo junto con el factor de $3^n$. Es similar a la situación que preferimos escribir $n!$ en lugar de $\prod_{j=1}^nj$.

Cuando se mira en el coeficiente de \begin{align*} a_n=\frac{-3a_{n-1}}{(n+m)(3n+3m+2)}\qquad n\geq 1,m\in \mathbb{R} \end{align*} calculado por la OP en uno de sus comentarios, se observa que el $m=-\frac{2}{3}$ es una buena elección, además de a $m=0$ a fin de obtener una segunda solución linealmente independiente, ya que en este caso tenemos la simple relación de recurrencia \begin{align*} a_n&=\frac{-3a_{n-1}}{(n-\frac{2}{3})3n}=\frac{-a_{n-1}}{(n-\frac{2}{3})n}\\ &=\frac{a_{n-2}}{\left(n-\frac{2}{3}\right)\left(n-\frac{5}{3}\right)n(n-1)}\\ &=\ldots\\ &=\frac{(-1)^na_0}{\left(n-\frac{2}{3}\right)\left((n-1)-\frac{2}{3}\right)\cdot\ldots\cdot\left((n-(n-1))-\frac{2}{3}\right)n!}\\ &=\frac{(-1)^na_0}{\left(n-\frac{2}{3}\right)\left(n-\frac{5}{3}\right)\cdot\ldots\cdot\left(\frac{1}{3}\right)n!}\\ &=\frac{(-1)^na_0}{n!\left(n-\frac{2}{3}\right)_n}\tag{3} \end{align*} Por supuesto, también podemos utilizar el producto símbolo en (3) en lugar del símbolo de Pochhammer \begin{align*} \left(n-\frac{2}{3}\right)_n&=\prod_{j=0}^{n-1}\left(n-j-\frac{2}{3}\right)=\prod_{j=0}^{n-1}\left(n-(n-1-j)-\frac{2}{3}\right)\\ &=\prod_{j=0}^{n-1}\left(j+\frac{1}{3}\right)=\frac{1}{3^n}\prod_{j=0}^{n-1}\left(3j+1\right) =\frac{1}{3^n}\prod_{j=1}^{n}\left(3j-2\right)\\ \end{align*}

Así, obtenemos $y_2(x)$ \begin{align*} y_2(x)&=a_0x^{-\frac{2}{3}}\sum_{n=0}^{\infty}\frac{(-1)^n}{\left(n-\frac{2}{3}\right)_n}\frac{x^n}{n!}\\ &=a_0x^{-\frac{2}{3}}\sum_{n=0}^{\infty}\frac{(-1)^n3^n}{\prod_{j=1}^{n}\left(3j-2\right)}\frac{x^n}{n!} \end{align*}

Conclusión: Alternativamente a la representación de $y(x)$ a través de funciones de Bessel de primera especie obtenemos como solución general para $y(x)$:

\begin{align*} y(x)&= c_0\sum_{n=0}^{\infty}\frac{(-1)^n}{{\left(n+\frac{2}{3}\right)}_n}\frac{x^n}{n!}+c_1x^{-\frac{2}{3}}\sum_{n=0}^{\infty}\frac{(-1)^n}{\left(n-\frac{2}{3}\right)_n}\frac{x^n}{n!}\\ &=c_0\sum_{n=0}^{\infty}\frac{(-1)^n3^n}{\prod_{j=1}^{n}(3j+2)}\frac{x^n}{n!}+c_1x^{-\frac{2}{3}}\sum_{n=0}^{\infty}\frac{(-1)^n3^n}{\prod_{j=1}^{n}\left(3j-2\right)}\frac{x^n}{n!}\qquad c_0,c_1\in\mathbb{R}\\ \end{align*}

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X