17 votos

Integral de la $\int_0^\infty \log(1-e^{-a x})\cos (bx)\, dx=\frac{a}{2b^2}-\frac{\pi}{2b}\coth \frac{\pi b}{a}$

$$\mathcal{J}:=\int_0^\infty \log(1-e^{-x})\cos (bx)\, dx=\frac{a}{2b^2}-\frac{\pi}{2b}\coth \frac{\pi b}{a},\qquad \mathcal{Re}(a)>0, b>0. $$ Traté de escribir $$ \mathcal{J}=-\int_0^\infty \sum_{n=1}^\infty\frac{e^{-anx}}{n}\cos(bx)\,dx $$ pero el sastre de la serie, $\log (1-\xi)=-\sum_{n=1}^\infty \xi^n/n, \ |\xi|<1$, por lo que este no es tan útil para hacer la integral. Traté de escribir $$ \mathcal{J}=\frac{1}{b}\int_0^\infty \log(1-e^{-ax})d(\sen bx)=\frac{1}{b}\left(\log(1-e^{-ax})\sin (bx)\big|^\infty_0 - \int_0^\infty \frac{\sin (bx)}{{e^{ax}-1}}dx \right), $$ el término se desvanece así tenemos $$ \mathcal{J}=\frac{a}{b}\int_0^\infty \frac{\sin(bx)}{1-e^{ax}}dx=\frac{a}{b}\mathcal{Im}\bigg[\int_0^\infty \frac{e^{ibx}}{e^{ax}-1}dx\bigg] $$ lo que no estoy seguro de cómo resolver. Aviso hay singularidades en $x=2i\pi n/a, \ n\in \mathbb{Z}$.

Tenemos que calcular el residuo de todas las singularidades a lo largo del eje imaginario. El residuo de la contribución a la integral $$ 2\pi i\cdot \sum_{n= 0}^\infty \frac{ e^{-2\pi nb/a}}{e^{2 \pi n}}=2\pi i \sum_{n=0}^\infty e^{n( -2\pi b/a-2i\pi)}=\frac{2\pi i}{e^{-(2\pi b/a+2\pi i)}}$$ Tomando la parte imaginaria da y re-escribir la integral da un resultado diferente. ¿De dónde me salen mal? ¿Cómo podemos calcular esto? Gracias

19voto

Leucippus Puntos 11926

Considere la integral \begin{align} I = \int_{0}^{\infty} \ln(1-e^{-ax}) \ \cos(bx) \ dx. \end{align} Expanda el logaritmo para obtener \begin{align} I &= - \sum_{n=1}^{\infty} \frac{1}{n} \ \int_{0}^{\infty} e^{-a n x} \ \cos(bx) \ dx \\ &= - \sum_{n=1}^{\infty} \frac{1}{n} \ \frac{an}{ a^{2} n^{2} + b^{2} } \\ &= - \frac{1}{a} \ \sum_{n=1}^{\infty} \frac{1}{n^{2} + (b/a)^{2}}. \end{align} El uso de la expansión \begin{align} \coth(\pi x) = \frac{1}{\pi x} + \frac{2 x}{\pi} \sum_{n=1}^{\infty} \frac{1}{n^{2} + x^{2}} \end{align} a continuación, el valor de la integral se convierte en \begin{align} \int_{0}^{\infty} \ln(1-e^{-ax}) \ \cos(bx) \ dx = \frac{a}{2 b^{2}} - \frac{\pi}{2 b} \ \coth\left(\frac{b \pi}{a}\right). \end{align}

11voto

Lance Roberts Puntos 11164

Esto puede no ser el método más fácil, pero que parece estar interesado en un contorno manera de ir sobre él.

a partir de su $$-a/b\int_{0}^{\infty}\frac{\sin(bx)}{e^{ax}-1}dx$$

Considere la función$$f(z)=\frac{e^{ibz}}{e^{az}-1}$$

El uso de un rectángulo en el primer cuadrante con la altura de la $\displaystyle 2\pi i/a$ con cuarto de círculo sangrías alrededor de $2\pi i/a$$0$.

Será de 6 partes para armar:

$$I_{1}+I_{2}+I_{3}+I_{4}+I_{5}+I_{6}=0........(1)$$

La integral se puede establecer en 0 porque no hay polos en el interior del contorno.

A lo largo de la parte inferior horizontal en el eje x: $$I_{1}=\int_{\epsilon}^{R}\frac{e^{ibx}}{e^{ax}-1}dx$$

arriba a la derecha lado vertical:

$$\left|\frac{e^{ibR}}{e^{aR}-1}\right|\to 0, \;\ as \;\ R\to \infty$$

$$I_{2}=0$$

a lo largo de la parte superior horizontal: $$I_{3}=-\int_{\epsilon}^{R}\frac{e^{ib(x+2\pi i/a)}}{e^{a(x+2\pi i/a)}-1}dx=-e^{-2\pi b/a}\int_{\epsilon}^{r}\frac{e^{ibx}}{e^{ax}-1}dx$$

top cuarto de círculo alrededor de la sangría en $2\pi i/a$,

donde x varía de

$(\epsilon, \epsilon+\frac{2\pi i}{a})$ $(0,\frac{2\pi i}{a}-\frac{2\pi i}{a}\epsilon)$

$$I_{4}=\frac{-\pi i}{2}Res\left(f(z), \frac{2\pi a}{b}\right)=\frac{-\pi i}{2}\cdot \frac{e^{ib(2\pi i/a)}}{ae^{a(2\pi i/a)}}=\frac{-\pi i}{2a}e^{-2\pi b/a}$$

Abajo a la izquierda lado vertical. parametrizar con $\displaystyle z=iy, \;\ dz=idy$

$$I_{5}=-i\int_{\epsilon}^{2\pi/a}\frac{e^{-by}}{e^{ayi}-1}dy$$

Cuarto de círculo de sangría alrededor del origen con x variando de $\displaystyle (0,i\epsilon)$ $\displaystyle (\epsilon, 0)$.

$$I_{6}=\frac{-\pi i}{2}Res(f,0)=\frac{-\pi i}{2}\cdot \frac{e^{ib(0)}}{ae^{a(0)}}=\frac{-\pi i}{2a}$$

Ahora, ensamblar todas las partes conectando a todos en (1):, y deje $\displaystyle \epsilon\to 0, \;\ R\to \infty$

$$\int_{C}\frac{e^{ibz}}{e^{az}-1}dz=\int_{0}^{\infty}\frac{e^{ibx}}{e^{ax}-1}dx+I_{2}-e^{-2\pi b/a}\int_{0}^{\infty}\frac{e^{ibx}}{e^{ax}-1}dx$$ $$-\frac{\pi i}{2a}e^{-2\pi b/a}-\frac{\pi i}{2a}-i\int_{0}^{2\pi /a}\frac{e^{-by}}{e^{ayi}-1}dy=0$$

$$\rightarrow (1-e^{-2\pi b/a})\int_{0}^{\infty}\frac{\sin(bx)}{e^{ax}-1}dx+\int_{0}^{2\pi/a}\frac{(-i)e^{-by}}{e^{ayi}-1}dy=\frac{\pi i}{2a}(1+e^{-2\pi b/a})$$

Tomando imaginaria, la última integral(el que va hacia abajo de la izquierda lado vertical) puede ser demostrado ser igual a

$$\int_{0}^{2\pi/a}\frac{e^{-by}}{2}dy=\frac{1-e^{-2\pi b/a}}{2b}$$

resolviendo la integral en cuestión, finalmente tenemos:

$$\int_{0}^{\infty}\frac{\sin(bx)}{e^{ax}-1}dx=\frac{\frac{\pi}{2a}(1+e^{-2\pi b/a})-\frac{1-e^{-2\pi b/a}}{2b}}{1-e^{-2\pi b/a}}$$

$$=\frac{\pi}{2a}\coth(\frac{\pi b}{a})-\frac{1}{2b}$$

multiplicando esto por el $-a/b$ desde el principio, se reduce a una forma que puede ser escrito en términos de funciones trigonométricas hiperbólicas como la solución sugiere.

y, al final obtenemos:

$$\frac{a}{2b^{2}}-\frac{\pi}{2b}\coth(\frac{\pi b}{a})$$

7voto

Felix Marin Puntos 32763

$\newcommand{\+}{^{\daga}} \newcommand{\ángulos}[1]{\left\langle\, nº 1 \,\right\rangle} \newcommand{\llaves}[1]{\left\lbrace\, nº 1 \,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\, nº 1 \,\right\rbrack} \newcommand{\ceil}[1]{\,\left\lceil\, nº 1 \,\right\rceil\,} \newcommand{\dd}{{\rm d}} \newcommand{\down}{\downarrow} \newcommand{\ds}[1]{\displaystyle{#1}} \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 (\, nº 1 \,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\parcial #3^{#1}}} \newcommand{\pp}{{\cal P}} \newcommand{\raíz}[2][]{\,\sqrt[#1]{\vphantom{\large Un}\,#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} \newcommand{\wt}[1]{\widetilde{#1}}$ $\ds{{\cal J}\equiv\int_{0}^{\infty}\ln\pars{1 - \expo{-ax}}\cos\pars{bx}\, \dd x={a \sobre 2b^{2}} - {\pi \sobre 2b}\,\coth\pars{\pi b \a través de una}:\ {\large ?}.\quad a > 0,\ b > 0}$.

\begin{align} {\cal J}&\equiv \Re\int_{0}^{\infty}\ln\pars{1 - \expo{-ax}}\expo{\ic bx}\,\dd x =\Re\int_{0}^{\infty}\sum_{n = 1}^{\infty}\pars{-\,{\expo{-nax} \over n}} \expo{\ic bx}\,\dd x \\[3mm]&=-\,\Re\sum_{n = 1}^{\infty}{1 \over n}\int_{0}^{\infty} \expo{-\pars{na - \ic b}x}\,\dd x =-\,\Re\sum_{n = 1}^{\infty}{1 \over n}\,{1 \over na - \ic b} \\[3mm]&=-\,{1 \over a}\Re\sum_{n = 0}^{\infty} {1 \over \pars{n + 1}\pars{n + 1 - \ic b/a}} =-\,{1 \over a}\, \Re\bracks{\Psi\pars{1} - \Psi\pars{1 - \ic b/a} \over 1 - \pars{1 - \ic b/a}} \end{align} donde $\ds{\Psi\pars{z}}$ es la La Función Digamma ${\bf\mbox{6.3.1}}$.

\begin{align} {\cal J}&=-\,{1 \over b}\,\Im\bracks{\Psi\pars{1} - \Psi\pars{1 - {b \over a}\,\ic}} ={1 \over b}\,\Im\Psi\pars{1 - {b \over a}\,\ic} \end{align}

Con la identidad ${\bf\mbox{6.3.13}}$: $$ {\cal J}={1 \over b}\llaves{% -\,{1 \over 2\pars{-b/a}} + \media\,\pi\coth\pars{\pi\bracks{-\,{b \a través de una}}}} $$

$$\color{#00f}{\large% {\cal J}\equiv\int_{0}^{\infty}\ln\pars{1 - \expo{-ax}}\cos\pars{bx}\,\dd x ={a \sobre 2b^{2}} - {\pi \sobre 2b}\,\coth\pars{\pi b \a través de una}} $$

2voto

schooner Puntos 1602

Vamos $$ I(a)=\int_{0}^{\infty}\ln(1-e^{-ax})\cos(bx)dx. $$ Entonces \begin{eqnarray} I'(a)=&=&\int_{0}^{\infty}\frac{xe^{-ax}}{1-e^{-ax}}\cos(bx)dx\\ &=&\int_{0}^{\infty}\sum_{n=0}^\infty xe^{-a(n+1)x}\cos(bx)dx\\ &=&\sum_{n=0}^\infty \frac{(a(n+1)-b)(a(n+1)+b)}{(a^2(n+1)^2+b^2)^2}\\ &=&\sum_{n=1}^\infty \left(\frac{1}{a^2n^2+b^2}-\frac{2b^2}{(a^2n^2+b^2)^2}\right)\\ &=&\frac{1}{2}\left(\frac{1}{b^2}-\frac{\pi^2}{a^2\sinh^2(\frac{b\pi}{a})}\right). \end{eqnarray} Así $$ I(a)=\frac{a}{2b^2}-\frac{\pi\cot(\frac{b\pi}{a})}{2b}. $$

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