10 votos

Informática $\int_0^\infty\mathrm{d} x\frac{x}{e^x+1}$ con la integración de contorno

Vamos a configurar: $$ \int_0^\infty\mathrm{d}x\frac{x}{e^x+1}=I. $$ me gustaría calcular utilizando, probablemente, los métodos de análisis complejo y la integración de contorno.

8voto

Felix Marin Puntos 32763

$\newcommand{\+}{^{\daga}} \newcommand{\ángulos}[1]{\left\langle #1 \right\rangle} \newcommand{\llaves}[1]{\left\lbrace #1 \right\rbrace} \newcommand{\bracks}[1]{\left\lbrack #1 \right\rbrack} \newcommand{\ceil}[1]{\,\left\lceil #1 \right\rceil\,} \newcommand{\dd}{{\rm d}} \newcommand{\down}{\downarrow} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\equalby}[1]{{#1 \cima {= \cima \vphantom{\enorme}}}} \newcommand{\expo}[1]{\,{\rm e}^{#1}\,} \newcommand{\fermi}{\,{\rm f}} \newcommand{\piso}[1]{\,\left\lfloor #1 \right\rfloor\,} \newcommand{\mitad}{{1 \over 2}} \newcommand{\ic}{{\rm i}} \newcommand{\iff}{\Longleftrightarrow} \newcommand{\imp}{\Longrightarrow} \newcommand{\isdiv}{\,\left.\a la derecha\vert\,} \newcommand{\cy}[1]{\left\vert #1\right\rangle} \newcommand{\ol}[1]{\overline{#1}} \newcommand{\pars}[1]{\left( #1 \right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\parcial #3^{#1}}} \newcommand{\pp}{{\cal P}} \newcommand{\raíz}[2][]{\,\sqrt[#1]{\,#2\,}\,} \newcommand{\sech}{\,{\rm sech}} \newcommand{\sgn}{\,{\rm sgn}} \newcommand{\totald}[3][]{\frac{{\rm d}^{#1} #2}{{\rm d} #3^{#1}}} \newcommand{\ul}[1]{\underline{#1}} \newcommand{\verts}[1]{\left\vert\, nº 1 \,\right\vert}$ $\large\tt\mbox{Just another one !!!}$

En el primer paso integramos por partes: \begin{align} &\color{#00f}{\large\int_{0}^{\infty}{x \over \expo{x} + 1}\,\dd x} =\int_{0}^{\infty}\ln\pars{1 + \expo{-x}}\,\dd x= \int_{0}^{\infty}\sum_{\ell = 1}^{\infty}{\pars{-1}^{\ell + 1} \over \ell}\, \expo{-\ell x}\,\dd x \\[3mm]&=\sum_{\ell = 1}^{\infty}{\pars{-1}^{\ell + 1} \over \ell}\ \overbrace{\int_{0}^{\infty}\expo{-\ell x}\,\dd x}^{\ds{1 \over \ell}} =\sum_{\ell = 1}^{\infty}{\pars{-1}^{\ell + 1} \over \ell^{2}} =\sum_{\ell = 0}^{\infty} \bracks{{1 \over \pars{2\ell + 1}^{2}} - {1 \over \pars{2\ell + 2}^{2}}} \\[3mm]&=\bracks{\sum_{\ell = 1}^{\infty}{1 \over \ell^{2}} -\sum_{\ell = 1}^{\infty}{1 \over \pars{2\ell}^{2}}} -{1 \over 4}\sum_{\ell = 1}^{\infty}{1 \over \ell^{2}} =\half\ \underbrace{\sum_{\ell = 1}^{\infty}{1 \over \ell^{2}}} _{\ds{\zeta\pars{2} = {\pi^{2} \over 6}}} =\ \color{#00f}{\large{\pi^{2} \over 12}} \end{align}

$\zeta\pars{z}$ es la Riemann Zeta Función .

7voto

Ron Gordon Puntos 96158

Considere la integral de contorno

$$\oint_C dz \frac{z^2}{e^z+1}$$

donde $C$ es el rectángulo con vértices en $0$, $R$, $R+i 2 \pi$, y $i 2 \pi$, en ese orden, con una pequeña hendidura semicircular de radio $\epsilon$ $z=i \pi$ en el rectángulo. Por lo tanto, el contorno de la integral es igual a

$$\int_0^R dx \frac{x^2-(x+i 2 \pi)^2}{e^x+1} + i \int_0^{2 \pi} dy \frac{(R+i y)^2}{e^{R+i y}+1} \\+i PV \int_0^{2 \pi} dy \frac{y^2}{e^{i y}+1} + i \epsilon \int_{\pi/2}^{-\pi/2} d\phi \, e^{i \phi} \frac{(i \pi+\epsilon e^{i \phi})^2}{-e^{\epsilon e^{i \phi}}+1}$$

donde $PV$ denota el valor principal de Cauchy. Como $R \to \infty$, la segunda integral se desvanece. Como $\epsilon \to 0$, la cuarta parte integral de los enfoques $-i \pi^3$. Por Cauchy teorema, el contorno de la integral es cero. Por lo tanto, hemos

$$-i 4 \pi \int_0^{\infty} dx \frac{x}{e^x+1} + 4 \pi^2 \int_0^{\infty} \frac{dx}{e^x+1}\\+i PV \int_0^{2 \pi} dy \frac{y^2}{e^{i y}+1} - i \pi^3 = 0$$

Ahora,

$$i PV \int_0^{2 \pi} dy \frac{y^2}{e^{i y}+1} = \frac12 PV \int_0^{2 \pi} dy \, y^2 \, \tan{\frac{y}{2}} + i \frac12 \int_0^{2 \pi} dy \, y^2 $$

La equiparación de la imaginaria, tenemos que

$$-4 \pi \int_0^{\infty} dx \frac{x}{e^x+1} = \pi^3 - \frac{4 \pi^3}{3}$$

o

$$\int_0^{\infty} dx \frac{x}{e^x+1} = \frac{\pi^2}{12}$$

5voto

Marko Riedel Puntos 19255

Por medio de la diversidad y desde que admitir variable compleja métodos en cuenta que lo que tenemos aquí es un valor especial de Mellin de transformación: $$\mathfrak{M}\left(\frac{1}{e^x+1}; s\right) = \int_0^\infty \frac{1}{e^x+1} x^{m-1} dx.$$ Para calcular esta transformación podemos proceder en un punto poco ortodoxa manera por la expansión de la función que se transforma en una serie: $$\int_0^\infty \frac{1}{e^x+1} x^{m-1} dx = \int_0^\infty \frac{e^{-x}}{1+e^{-x}} x^{m-1} dx = \int_0^\infty \sum_{q\ge 1} (-1)^{q+1} e^{-qx} x^{m-1} dx \\ = \sum_{q\ge 1} (-1)^{q+1} \int_0^\infty e^{-qx} x^{m-1} dx = \Gamma(s) \sum_{q\ge 1} \frac{(-1)^{q+1}}{q^s} = \Gamma(s) \left(1-\frac{2}{2^s}\right) \zeta(s).$$ Poner $s=2$ obtener $$\Gamma(2) \times \frac{1}{2} \times \zeta(2) = \frac{\pi^2}{12}.$$ Aquí hemos utilizado la representación integral de la función gamma, que supongo que es admisible en un complejo de la variable tipo de método de respuesta.

4voto

Aaron Puntos 1

Partimos de la integral de contorno evaluado a través del teorema de Cauchy: $$ \oint_{\Gamma_{\epsilon,R}}\mathrm{d}z\frac{z^2}{e^z+1} = 0 $$ donde $\Gamma_{\epsilon, R}$ es el rectángulo de vértices $0,\ R,\ R+i2\pi,\ i2\pi$, con una sangría en la singularidad $i\pi$ con un semicírculo de radio $\epsilon$.

División en sus ramas naturales de los rendimientos: $$ \int_0^R\mathrm{d}x\frac{x^2}{e^x+1}+ \int_0^{2\pi}i\mathrm{d}y \frac{(iy+R)^2}{e^{iy+R}+1}+ \int_R^0\mathrm{d}x\frac{(x+i2\pi)^2}{e^{x+i2\pi}+1}+ \left(\int_{2\pi}^{\pi+\epsilon}+ \int_{\pi\epsilon}^0\right)i\mathrm{d}y\frac{(iy)^2}{e^{i y}+1}+ \int_{\pi/2}^{-\pi/2}i\epsilon e^{i\theta}\mathrm{d}\theta\frac{(\epsilon e^{i\theta}+i\pi)^2}{e^{i\pi + \epsilon e^{i\theta}}+1}=0.$$

Ahora observamos que la primera integral se cancela con el primer término de la expansión de la tercera integral, además el segundo $\to0$$R\to \infty$. Nuestra deseada integral aparece como la rectangular término de dicha expansión.

Ahora estamos a la izquierda con los siguientes cálculos: $$ -\int_{-\pi/2}^{\pi/2}i\epsilon e^{i\theta}\mathrm{d}\theta\frac{(\epsilon e^{i\theta}+i\pi)^2}{e^{i\pi + \epsilon e^{i\theta}}+1}=i\pi^2\int_{-\pi/2}^{\pi/2}\mathrm{d}\theta(-1+O(\epsilon)) \to_{\epsilon\to0}-i\pi^3 $$ y, utilizando la fórmula de Euler y la división en real y la parte imaginaria: $$ i\int_0^{2\pi}\mathrm{d}y\frac{y^2}{\cos y+1+i\pecado y}=\int_0^{2\pi}\mathrm{d}y\ \frac{y^2}{2}- \int_0^{2\pi}\mathrm{d}y\frac{y^2\pecado y}{2(\cos y +1}. $$ Por lo tanto, el uso de la fácil-a-compruebe $\int\mathrm{d}x\frac{1}{e^x+1}=\ln 2$: $$ 0-i4\pi I + 4\pi^2\ln2+i\int_0^{2\pi}\mathrm{d}y\ \frac{y^2}{2}-\int_0^{2\pi}\mathrm{d}y\frac{y^2\pecado y}{2(\cos y +1}-i\pi^3=0. $$ Por comparación de lo real y parte imaginaria: $$ I=\frac{\pi^2}{12} $$ $$ \int_0^{2\pi}\mathrm{d}y\frac{y^2\pecado y}{2(\cos y +1}=4\pi^2\ln 2 $$

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