21 votos

Constante de Euler-Mascheroni en la integral de la función de Bessel

Actualmente estoy haciendo malabares con algunas integrales. En un libro de texto de física, Chaikin-Lubensky [1] En el capítulo 6, (6.1.26), encontré una integral que dice \begin{equation} \int_0^{1} \textrm{d} y\, \frac{1 - J_0(y)}{y} - \int_{1}^{\infty} \textrm{d} y\, \frac{J_0(y)}{y} = -.116. \end{equation} Sólo dan el resultado como valor en coma flotante sin nombrar las fuentes. El valor se parece sospechosamente a $\gamma - \ln(2)$ a mí ( $\gamma$ siendo la constante de Euler-Mascheroni), lo que resolvería un problema que tengo en otra parte. No estoy familiarizado con las manipulaciones típicas que se utilizan en este tipo de integrales ni con las distintas definiciones de la constante de Euler-Mascheroni. He trasteado un poco con integrales de coseno $\textrm{Ci}(y)$ pero no llegó muy lejos con ella. Así que me alegro de recibir sugerencias.

26voto

metamorphy Puntos 186

Una forma relativamente elemental es empezar por lo conocido $$\gamma=\int_0^1\frac{1-\cos t}{t}\,dt-\int_1^\infty\frac{\cos t}{t}\,dt.$$

Ponga $t=ax$ para $a>0$ y hacer algunos reordenamientos, para obtener $$\int_0^1\frac{1-\cos ax}{x}\,dx-\int_1^\infty\frac{\cos ax}{x}\,dx=\gamma+\log a.$$

Ahora la representación integral $J_0(y)=\frac2\pi\int_0^{\pi/2}\cos(y\cos x)\,dx$ produce $$\int_0^1\frac{1-J_0(y)}{y}\,dy-\int_1^\infty\frac{J_0(y)}{y}\,dy=\frac2\pi\int_0^{\pi/2}(\gamma+\log\cos x)\,dx$$ después de intercambiar las integraciones (lo que no es difícil de justificar).

El resultado se deduce ahora de $\int_0^{\pi/2}\log\cos x\,dx\color{gray}{=\int_0^{\pi/2}\log\sin x\,dx}=-(\pi/2)\log2$ .

18voto

psychotik Puntos 171

Para una amplia clase de funciones $f$ se puede justificar rigurosamente el siguiente cálculo heurístico: 1)

\begin{align*} &\int_{0}^{\infty} \frac{f(x) - f(0)\mathbf{1}_{[0,1]}(x)}{x} \, \mathrm{d}x \\ &= \int_{0}^{\infty} \left( f(x) - f(0)\mathbf{1}_{[0,1]}(x) \right) \left( \int_{0}^{\infty} e^{-xs} \, \mathrm{d}s \right) \, \mathrm{d}x \\ &= \int_{0}^{\infty} \left[ \int_{0}^{\infty} \left( f(x) - f(0)\mathbf{1}_{[0,1]}(x) \right) e^{-xs} \, \mathrm{d}x \right] \, \mathrm{d}s \tag{“Fubini”} \\ &= \int_{0}^{\infty} \left[ \mathcal{L}\{f\}(s) - f(0)\mathcal{L}\{\mathbf{1}_{[0,1]}\}(s) \right] \, \mathrm{d}s \\ &= \lim_{R \to \infty} \int_{0}^{R} \left[ \mathcal{L}\{f\}(s) - f(0)\mathcal{L}\{\mathbf{1}_{[0,1]}\}(s) \right] \, \mathrm{d}s \\ &= \lim_{R\to\infty} \left( \int_{0}^{R} \mathcal{L}\{f\}(s) \, \mathrm{d}s - f(0)(\gamma + \log R) \right) \end{align*}

donde $\mathcal{L}\{f\}$ denota la transformada de Laplace de $f$ y utilizamos el conocido resultado de que

$$ \int_{0}^{R} \mathcal{L}\{\mathbf{1}_{[0,1]}\}(s) \, \mathrm{d}s = \int_{0}^{R} \frac{1 - e^{-s}}{s} \, \mathrm{d}s = \gamma + \log R + o(1) $$

como $ R \to \infty $ . (Sugerencia: Realice la integración por partes y utilice $\int_{0}^{\infty} e^{-s}\log s \, \mathrm{d}s = -\gamma$ .) Así, cuando $f = J_0$ es la función de Bessel del 1er tipo y orden $0$ ,

$$ J_0(0) = 1 \qquad \text{and}\qquad \mathcal{L}\{J_0\}(s) = \frac{1}{\sqrt{s^2 + 1}}, $$

y así,

\begin{align*} \int_{0}^{\infty} \frac{J_0(x) - \mathbf{1}_{(0,1)}(x)}{x} \, \mathrm{d}x &=\lim_{R\to\infty} \left[ \int_{0}^{R} \frac{\mathrm{d}s}{\sqrt{s^2+1}} - (\gamma + \log R) \right] \\ &=\lim_{R\to\infty} \left[ \log(R + \sqrt{R^2 + 1}) - (\gamma + \log R) \right] \\ &= \log 2 - \gamma. \end{align*}

Reordenando esto se obtiene

$$ \int_{0}^{\infty} \frac{\mathbf{1}_{(0,1)}(x) - J_0(x)}{x} \, \mathrm{d}x = \gamma - \log 2. $$

14voto

Claude Leibovici Puntos 54392

Utilización de la función G de Meijer $$\int \frac{J_0(y)}{y}\,dy=-\frac{1}{2} G_{1,3}^{2,0}\left(\frac{y^2}{4}| \begin{array}{c} 1 \\ 0,0,0 \end{array} \right)$$ $$\int_0^1 \frac{1-J_0(y)}{y}\,dy=\frac{1}{8} \, _2F_3\left(1,1;2,2,2;-\frac{1}{4}\right)$$ $$\int_1^\infty \frac{J_0(y)}{y}\,dy=\frac{1}{8} \, _2F_3\left(1,1;2,2,2;-\frac{1}{4}\right)-\gamma +\log (2)$$ Buena anulación de las condiciones.

Editar

Hacerlo más general $$\int_0^a \frac{1-J_0(y)}{y}\,dy=\frac{1}{8} a^2 \, _2F_3\left(1,1;2,2,2;-\frac{a^2}{4}\right)$$ $$\int_a^\infty \frac{J_0(y)}{y}\,dy=\frac{1}{8} a^2 \, _2F_3\left(1,1;2,2,2;-\frac{a^2}{4}\right)-\log (a)-\gamma +\log (2)$$ $$\int_0^a \frac{1-J_0(y)}{y}\,dy-\int_a^\infty \frac{J_0(y)}{y}\,dy=\gamma+\log \left(\frac{a}{2}\right)$$

11voto

Yuriy S Puntos 179

Intentaré una prueba más directa.

$$ I=\int_0^{1} \textrm{d} y\, \frac{1 - J_0(y)}{y} - \int_{1}^{\infty} \textrm{d} y\, \frac{J_0(y)}{y} $$

Tenemos:

$$J_0(y)=\frac{2}{\pi}\int_0^1 \frac{\cos (x y) dx}{\sqrt{1-x^2}}$$

Intentemos la expansión en serie (formal) de la función coseno, entonces obtenemos:

$$J_0(y)=\frac{2}{\pi} \sum_{n=0}^\infty \frac{(-1)^n y^{2n}}{(2n)!} \int_0^1 \frac{x^{2n} dx}{\sqrt{1-x^2}}$$

Ahora podemos subsumir (formalmente) esta serie en la integral original, por ahora consideraré la primera parte:

$$I_1=\frac{2}{\pi} \sum_{n=1}^\infty \frac{(-1)^{n-1} }{(2n)!} \int_0^1 \frac{x^{2n} dx}{\sqrt{1-x^2}} \int_0^1 y^{2n-1} dy$$

$$I_1=\frac{1}{\pi} \sum_{n=1}^\infty \frac{(-1)^{n-1} }{(2n)! n} \int_0^1 \frac{x^{2n} dx}{\sqrt{1-x^2}}$$

La última integral es sólo la función Beta:

$$\int_0^1 \frac{x^{2n} dx}{\sqrt{1-x^2}}=\frac{1}{2} \int_0^1 \frac{u^{n-1/2} du}{\sqrt{1-u}}=\frac{1}{2} B\left(n+\frac{1}{2},\frac{1}{2}\right)=\frac{1}{2} \frac{\Gamma(n+1/2)\Gamma(1/2)}{\Gamma(n+1)}=\\=\frac{\sqrt{\pi}}{2} \frac{\Gamma(n+1/2)}{n!}$$

Esto nos da:

$$I_1=\frac{1}{2\sqrt{\pi}} \sum_{n=1}^\infty \frac{(-1)^{n-1} \Gamma(n+1/2)}{(2n)! n! n}= \frac{1}{2\sqrt{\pi}} \sum_{n=0}^\infty \frac{(-1)^{n} \Gamma(n+3/2)}{(2n+2)! (n+1)! (n+1)}$$

Se trata claramente de series hipergeométricas generalizadas. Vamos a encontrar el coeficiente de orden 0:

$$c_0=\frac{\Gamma(3/2)}{4\sqrt{\pi}}=\frac18$$

Ahora la relación de coeficientes:

$$\frac{c_{n+1}}{c_n}=\frac{(n+3/2)(n+1)(n+1)}{(2n+3)(2n+4)(n+2)(n+2)} \frac{-1}{n+1}=$$

$$=\frac{(n+1)(n+1)}{(n+2)(n+2)(n+2)} \frac{-1}{4(n+1)}$$

De esto concluimos:

$$I_1=\frac18 {_2F_3} \left(1,1;2,2,2; -\frac14 \right)$$

Creo que la segunda integral puede hacerse de forma similar, pero habrá que poner algunos límites, que es donde entran las constantes.

Continuaré cuando tenga tiempo.

9voto

Thierry Lam Puntos 1079

Teorema maestro de Ramanujan muestra que $$\int_{0}^{\infty} J_{0}(t) t^{s-1} \, \mathrm dt= \frac{2^{s-1}\Gamma \left(\frac{s}{2} \right)}{\Gamma \left(1-\frac{s}{2} \right)}, \quad 0 < \Re(s) < 3/2. \tag{1}$$

La función $ J_{0}(t)-1$ tiene la misma transformada de Mellin, pero la región de convergencia es $-2 < \Re(s) < 0$ .

(Véase el teorema 8.1 del documento Teorema maestro de Ramanujan por T. Amdeberhan et al.)

Por lo tanto, tenemos $$ \begin{align} \int_{0}^{\infty} \frac{\mathbf{1}_{[0,1]}(t)-J_{0}(t)}{t} dt &=\lim_{s \to 0^{-}} \int_{0}^{\infty}\left( \mathbf{1}_{[0,1]}(t)-J_{0}(t) \right)t^{s-1} \, \mathrm dt \\&= \lim_{s \to 0^{-}} \left( \int_{0}^{\infty}\left(1-J_{0}(t) \right) t^{s-1} \, \mathrm dt - \int_{1}^{\infty} t^{s-1} \, \mathrm dt\right) \\ &= \lim_{s \to 0^{-}} \left( -\frac{2^{s-1}\Gamma \left(\frac{s}{2} \right)}{\Gamma \left(1-\frac{s}{2} \right)} + \frac{1}{s}\right) \\ &= \lim_{s \to 0^{-}} \left(\frac{-2^{s} \Gamma \left(1+\frac{s}{2}\right)+\Gamma \left(1- \frac{s}{2} \right)}{s \, \Gamma \left(1- \frac{s}{2} \right)}\right) \\ &\overset{(1)}{=} \lim_{s \to 0^{-}} \frac{- 2^{s} \ln(2) \Gamma\left(1+ \frac{s}{2} \right) - 2^{s-1} \Gamma'\left(1+ \frac{s}{2}\right) - \frac{1}{2} \Gamma'\left(1- \frac{s}{2} \right)}{\Gamma \left(1- \frac{s}{2} \right) - \frac{s}{2}\Gamma'\left(1- \frac{s}{2} \right)} \\ & = \frac{-\ln(2) - \frac{1}{2} \Gamma'(1)- \frac{1}{2} \Gamma'(1)}{1}\\ & = \gamma - \ln(2) . \end{align}$$

$(1)$ La regla de L'Hôpital


El mismo planteamiento muestra también, por ejemplo, que $$\int_{0}^{\infty} \frac{\frac{t}{2}\mathbf{1}_{[0,1]}(t)- J_{1}(t)}{t^{2}} \, \mathrm dt = \frac{2 \gamma - 2 \ln(2) -1}{4}. $$

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