2 votos

Una solución aproximada de una ecuación diferencial de segundo orden

En el libro, Métodos matemáticos para estudiantes de física y campos afines, segunda edición de Sadri Hassani, página 667, el autor ha afirmado que, para la siguiente ecuación diferencial

$y''(x) - x^2 y(x) \approx 0$ ,

donde $x \to \infty$ se puede fácilmente obtener una solución aproximada de la forma $e^{\pm x^2/2}$ .

¿Existe algún enfoque para obtener esta solución? además de resolver la ecuación diferencial exacta $y''(x) - x^2 y(x) = 0$ por el método de Frobenius, y luego tomando el límite de las soluciones (polinomios de Hermite) como $x \to \infty$ ?

2 votos

Una idea para pensar en un enfoque heurístico: $(D+x)(D-x)=D^2-Dx+xD-x^2$ Así que a la inversa $D^2-x^2=(D+x)(D-x)+Dx-xD$ . Ahora $Dx-xD=1$ (el operador de identidad), por lo que su ecuación es $((D+x)(D-x)+1)y=0$ . ¿Puede argumentar que el $1$ término es "sólo una perturbación" para grandes $x$ ? Si es así, el problema de "orden principal" es mucho más fácil ya que se puede resolver por el método de integración de factores.

0 votos

@Ian Gracias por tu perspicaz comentario.

1 votos

¿El método de Frobenius consiste en suponer que la solución de la ecuación de segundo orden se puede expandir a series infinitas y así intentar obtener la solución de esta manera?

1voto

andy.holmes Puntos 518

En general, la ecuación $y''-qy=0$ con $0<q(x)\to\infty$ para $x\to\infty$ tiene la aproximación WKB (ver wikipedia, este es un ejemplo estándar) de soluciones de base $$ y(x)=q(x)^{-\frac14}\exp\left(\pm\int\sqrt{q(x)}dx\right) $$ donde el primer factor es el término de segundo orden en la expansión. Así que, efectivamente, en primer orden se obtendría $y(x)=\exp(\pm\frac{x^2}2)$ mientras que en el segundo orden, habría además un factor $\frac1{\sqrt{x}}$ .


El enfoque consiste en establecer $y=\exp(S)$ con una expansión $S=S_0+S_1+S_2+...$ donde $S_0\gg S_1\gg S_2\gg...$ para $x\to\infty$ . Esta escala también se traduce en las derivadas. Juego para la brevedad $S'=s$ entonces aislando los componentes de igual escala en $s'+s^2-q=0$ da $$ s_0^2=q\\ s_0'+2s_0s_1=0\\\vdots $$ lo que implica $s_0=\pm\sqrt q$ y $s_1=-\frac{s_0'}{2s_0}\implies S_1=-\frac12\ln|s_0|=-\frac14\ln(q)$ .

0 votos

Gracias por escribir una respuesta. Si queremos comprobar las soluciones, es decir, $e^{\pm x^2/2}$ e insertarlos en la ecuación diferencial: $y''(x) - x^2 y(x) \approx 0$ ; para la solución de signo negativo, obtenemos: $- e^{-x^2/2} \approx 0$ donde como $x$ va al infinito, la exponencial desaparece, por lo que es una solución correcta. Mi pregunta es: ¿Por qué la solución de signo positivo también es aceptable? (Porque no satisface la ecuación diferencial).

1 votos

El resto $e^{x^2/2}$ es pequeño frente a los otros términos de tamaño $x^2e^{x^2/2}$ para $x\approx\infty$ . Esto no es mucho, pero también es sólo la primera aproximación asintótica.

1 votos

@Novice De hecho, en términos de error relativo, la situación acaba siendo la misma: el error relativo decae como una potencia negativa de $x$ como $|x| \to \infty$ . Esto es así porque el "carácter exponencial" de la solución ya está correctamente capturado por la aproximación WKB de orden principal.

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