Integración de $\int_0^{\pi/2} \frac{\sin^2 ax}{\sin^2 x}\,dx$

Question

Estoy buscando formas de evaluar:

$$\int_0^{\pi/2} \frac{\sin^2 ax}{\sin^2 x}\,dx$$

donde, $a$ es cualquier número real positivo. (Técnicas de integración real o compleja, cualquiera sirve).

Nota al margen de los votos negativos/cerrados: Se discutió en la sala de chat 36 de SE ayer (en caso de que no hayas revisado el comentario de robjohn) y decidí publicarlo en main para que robjohn pueda compartir su impresionante solución.

Answer 1

Utilizaremos $$ H(x)=\sum_{k=1}^\infty\left(\frac1k-\frac1{k+x}\right)\tag{1} $$ que es la extensión analítica estándar de los Números Armónicos. $$ \begin{align} &\int_0^{\pi/2}\frac{\sin^2(ax)}{\sin^2(x)}\,\mathrm{d}x\tag{2}\\ &=\frac12\int_{-\pi/2}^{\pi/2}\frac{\left(e^{iax}-e^{-iax}\right)^2}{\left(e^{ix}-e^{-ix}\right)^2}\,\mathrm{d}x\tag{3}\\ &=\frac1{4i}\int_{-\pi}^\pi\frac{\left(e^{iax}-1\right)^2}{\left(e^{ix}-1\right)^2}e^{-iax}\,\mathrm{d}e^{ix}\tag{4}\\ &=\frac1{4i}\int_\gamma\frac{z^a-2+z^{-a}}{\left(z-1\right)^2}\,\mathrm{d}z\tag{5}\\ &=\frac1{4i}\int_0^1\frac{t^a\left(e^{-\pi ia}-e^{\pi ia}\right)+t^{-a}\left(e^{\pi ia}-e^{-\pi ia}\right)}{\left(t+1\right)^2}\,\mathrm{d}t\tag{6}\\ &=\frac{\sin(\pi a)}2\int_0^1\frac{t^{-a}-t^a}{\left(t+1\right)^2}\,\mathrm{d}t\tag{7}\\ &=\frac{\sin(\pi a)}2\sum_{k=0}^\infty(-1)^k(k+1)\int_0^1\left(t^{k-a}-t^{k+a}\right)\,\mathrm{d}t\tag{8}\\ &=\frac{\sin(\pi a)}2\sum_{k=0}^\infty(-1)^k(k+1)\left(\frac1{k+1-a}-\frac1{k+1+a}\right)\tag{9}\\ &=\frac{\sin(\pi a)}2\sum_{k=0}^\infty(-1)^ka\left(\frac1{k+1-a}+\frac1{k+1+a}\right)\tag{10}\\ &=\frac{a\sin(\pi a)}2\sum_{k=1}^\infty\left(\frac1{2k-1-a}-\frac1{2k-a}+\frac1{2k-1+a}-\frac1{2k+a}\right)\tag{11}\\[3pt] &=\frac{a\sin(\pi a)}4\left(-H\left(\frac{-1-a}2\right)+H\left(\frac{-a}2\right)-H\left(\frac{-1+a}2\right)+H\left(\frac{a}2\right)\right)\tag{12}\\[6pt] &=\frac{a\sin(\pi a)}4\left(-\psi\left(\frac{1-a}2\right)+\psi\left(\frac{-a}2\right)-\psi\left(\frac{1+a}2\right)+\psi\left(\frac{a}2\right)\right)\tag{13} \end{align} $$ Explicación:
$\phantom{1}(3)$ : $\sin(x)=\frac{e^{ix}-e^{-ix}}{2i}$
$\phantom{1}(4)$ : sustituto $x\mapsto x/2$
$\phantom{1}(5)$ : $z=e^{ix}$ y $\gamma$ es una circunferencia unitaria ccw y el eje real negativo es un corte de rama
$\phantom{1}(6)$ : colapsar el contorno en el eje real negativo
$\phantom{(10)\text{:}}$ set $\color{#C00000}{z=te^{-i\pi}}$ abajo y $\color{#00A000}{z=te^{i\pi}}$ por encima del eje
$\phantom{1}(7)$ : $\sin(x)=\frac{e^{ix}-e^{-ix}}{2i}$
$\phantom{1}(8)$ : aplicar la serie de potencias para $\frac1{(1+x)^2}$
$\phantom{1}(9)$ : integrar
$(10)$ : álgebra
$(11)$ : separa los términos positivos y negativos
$(12)$ : aplicar $(1)$
$(13)$ : $H(x)=\gamma+\psi(1+x)=\gamma+\psi(x)+\frac1x$

Answer 2

Creo que he encontrado una manera de abordar esta integral que es un poco diferente de la bonita respuesta de @robjohn. Para empezar anotamos la siguiente indentidad algebraica simple

$$ \frac{1}{\sin^2(x)}=\frac{-4}{(e^{i x}-e^{-i x})^2}=\frac{-4 e^{2 ix}}{(e^{i x}-1)^2}=-2 i\partial_x\frac{1}{e^{2 ix}-1} $$

por lo que podemos reescribir la integral como

$$ I(a)=-2i\int_0^{\pi/2}\sin^2(a x)\left(\partial_x\frac{1}{e^{2 ix}-1}\right) $$

esto es ahora ideal para una integración por partes que nos deja

$$ I(a)=-2 i\color{\red}{\sin^2(a \pi/2)}+2ia\underbrace{\int_0^{\pi/2}\frac{\sin(2a x)}{e^{2 ix}-1}}_{J(a)} \quad \color{\red}{(*)} $$

expandiendo el denominador en una serie geométrica e intercambiando integración y suma obtenemos

$$ J(a)=-\sum_{n=0}^{\infty}\int_0^{\pi/2}\sin(2a x)e^{2 i xn} $$

las integrales resultantes se realizan fácilmente

$$ J(a)=\frac{1}{2}\left[a\underbrace{\sum_{n=0}^{\infty}\frac{1}{ (a^2-n^2)}}_{\color{\orange}{S_1(a)}}-a\cos (\pi a)\underbrace{\sum_{n=0}^{\infty}\frac{(-1)^n }{(a^2-n^2)}}_{\color{\green}{S_2(a)}}+i\sin (\pi a)\underbrace{\sum_{n=0}^{\infty}\frac{(-1)^n n }{(a^2-n^2)}}_{\color{\violet}{S_3(a)}}\right] $$

La suma $\color{\violet}{S_3(a)}$ ya fue calculado por Rob, y no repetimos los pasos correspondientes aquí.

Lo que queda por demostrar es esencialmente que $\frac{a}{2}\left(\color{\orange}{S_1(a)}-\cos(\pi a)\color{\green}{S_2(a)}\right)=\color{\red}{\sin^2(a \pi/2)}$ para que las partes imaginarias en $\color{\red}{(*)}$ anular. Empezamos con $\color{\green}{S_2(a)}$ :

$$ \color{\green}{S_2(a)}=\frac{1}{2a}\left(\sum_{n=0}^{\infty}\frac{(-1)^n}{n+a}-\frac{(-1)^n}{n-a}\right)=\frac{1}{2a}\left(\sum_{n\in\mathbb{Z}}\frac{(-1)^n}{a-n}-\frac{1}{a}\right) $$

Ahora podemos emplear, por ejemplo la expansión del polo de $\frac{1}{\sin(z)}$ y obtener $$ \color{\green}{S_2(a)}=\frac{1}{2a^2}+\frac{\pi}{2\sin(\pi a)a} $$

Para calcular $\color{\orange}{S_1(a)}$ utilizamos la expansión de polos de $\cot(z)$ y obtener

$$ \color{\orange}{S_1(a)}=\frac{1}{2a^2}+\frac{\pi \cot(\pi a) }{ 2a} $$

Utilizando el hecho de que $1-\cos(x)=2\sin^2(x/2)$ y $\cot(x)=\frac{\cos(x)}{\sin(x)}$ podemos llegar a la conclusión de que

$$ \frac{a}{2}\left(\color{\orange}{S_1(a)}-\cos(\pi a)\color{\green}{S_2(a)}\right)=\frac{a}{2}\left(2 \frac{\color{\red}{\sin^2(a \pi/2)}}{a}\right)=\color{\red}{\sin^2(a \pi/2)} $$

por lo que la parte imaginaria de $\color{red}{(*)}$ es cero y por lo tanto

$$ I(a)=-a\sin (\pi a)\color{\violet}{S_3(a)} $$

que coincide con la otra respuesta dada

Answer 3

Parece que hay muchos métodos para lanzar en este caso.

Utilice las piezas dejando

$u=\sin^{2}(ax), \;\ dv=\frac{1}{\sin(x)}dx, \;\ du=2a\sin(ax)\cos(ax)dx, \;\ v=-cot(x)dx$

$$-\sin^{2}(ax)\cot(x)|_{0}^{\frac{\pi}{2}}+2a\int_{0}^{\frac{\pi}{2}}\sin(ax)\cos(ax)\cot(x)dx$$

Como la parte izquierda es $0$ y utilizando la identidad $2\sin(u)\cos(u)=\sin(2u)$ entonces tenemos..:

$$a\int_{0}^{\frac{\pi}{2}}\sin(2ax)\cot(x)dx\tag{1}$$

Parece que hay muchas maneras de hacerlo. Voy a utilizar la famosa relación:

$$\int_{a}^{b}p(x)\cot(x)dx=2\sum_{n=1}^{\infty}\int_{a}^{b}p(x)\sin(2nx)dx$$

Está bien en este caso dejando $p(x)=\sin(2ax), \;\ p(0)=0$

Así,

$$2a\sum_{n=1}^{\infty}\int_{0}^{\frac{\pi}{2}}\sin(2ax)\sin(2nx)dx$$

$$\frac{a}{2}\underbrace{\sum_{n=1}^{\infty}\frac{\sin(\pi(a-n))}{a-n}}_{[1]}-\frac{a}{2}\underbrace{\sum_{n=1}^{\infty}\frac{\sin(\pi(a+n))}{a+n}}_{[2]}$$

He cortado por lo sano usando la tecnología, pero estas sumas son algo famosas y sumando [1]-[2] resulta:

$$[1]=\frac{a}{2}\left(\frac{-1}{2}\sin(\pi a)\left(\psi\left(1/2-a/2\right)-\psi\left(1-a/2\right)\right)\right)$$

$$[2]=\frac{a}{2}\left(\frac{-1}{2}\sin(\pi a)\left(\psi\left(a/2+1\right)-\psi\left(a/2+1/2\right)\right)\right)$$

$$[1]-[2]=\text{required result}$$

Otros pensamientos:

volver a utilizar las piezas en (1)

o

Considere $1/2\int_{0}^{\pi}\frac{1-\cos(ax)}{1-\cos(x)}dx$$

Answer 4

$\newcommand{\angles}[1]{\left\langle\, #1 \,\right\rangle} \newcommand{\braces}[1]{\left\lbrace\, #1 \,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\, #1 \,\right\rbrack} \newcommand{\dd}{\mathrm{d}} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,} \newcommand{\half}{{1 \over 2}} \newcommand{\ic}{\mathrm{i}} \newcommand{\iff}{\Leftrightarrow} \newcommand{\imp}{\Longrightarrow} \newcommand{\ol}[1]{\overline{#1}} \newcommand{\pars}[1]{\left(\, #1 \,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\root}[2][]{\,\sqrt[#1]{\, #2 \,}\,} \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}} \newcommand{\verts}[1]{\left\vert\, #1 \,\right\vert}$ Pregunta : $\ds{\int_{0}^{\pi/2}{\sin^{2}\pars{ax} \over \sin^{2}\pars{x}}\,\dd x\,,\quad a \in \mathbb{R}}$

---

\begin{align} &\color{#f00}{\int_{0}^{\pi/2}{\sin^{2}\pars{ax} \over \sin^{2}\pars{x}}\,\dd x} = \int_{0}^{\pi/2}{1 - \cos\pars{2ax} \over 1 - \cos\pars{2x}}\,\dd x = \half\int_{0}^{\pi}{1 - \cos\pars{ax} \over 1 - \cos\pars{x}}\,\dd x \\[3mm] = &\ {1 \over 4}\int_{-\pi}^{\pi}{1 - \cos\pars{ax} \over 1 - \cos\pars{x}}\,\dd x = {1 \over 4}\,\Re\int_{-\pi}^{\pi}{1 + \verts{a}x\ic - \expo{\verts{a}x\ic} \over 1 - \cos\pars{x}}\,\dd x \\[3mm] = &\ {1 \over 4}\,\Re\oint_{\verts{z} = 1}{1 + \verts{a}\ln\pars{z} - z^{\verts{a}} \over 1 - \pars{z + z^{-1}}/2}\,{\dd z \over \ic z} = \half\,\Im\oint_{\verts{z} = 1} {z^{\verts{a}} - \verts{a}\ln\pars{z} - 1 \over \pars{z - 1}^{2}}\,\dd z \end{align}

Elegiremos el $z^{\verts{a}}$ y $\ln\pars{z}$ cortes de rama a lo largo del negativo $x$ -eje con $z = 0$ como su puntos de ramificación . Entonces, el contorno de integración es un ojo de la cerradura que no encierra ningún polo, de manera que toda la contribución a la integral proviene de las integrales por encima y por debajo de los cortes de rama. A saber, \begin{align} &\color{#f00}{\int_{0}^{\pi/2}{\sin^{2}\pars{ax} \over \sin^{2}\pars{x}}\,\dd x} \\[5mm] = &\ -\,\half\,\Im\int_{-1}^{0}{\pars{-x}^{\verts{a}}\exp\pars{\pi\verts{a}\ic} - \verts{a}\bracks{\ln\pars{-x} + \ic\pi} - 1 \over \pars{x - 1}^{2}}\,\dd x \\[3mm] & -\ \half\,\Im\int_{0}^{-1}{\pars{-x}^{\verts{a}}\exp\pars{-\pi\verts{a}\ic} - \verts{a}\bracks{\ln\pars{-x} - \ic\pi} - 1 \over \pars{x - 1}^{2}}\,\dd x \\[5mm] = &\ -\,\half\,\Im\int_{0}^{1}{x^{\verts{a}}\exp\pars{\pi\verts{a}\ic} - \verts{a}\bracks{\ln\pars{x} + \ic\pi} - 1 \over \pars{1 + x}^{2}}\,\dd x \\[3mm] & +\ \half\,\Im\int_{0}^{1}{x^{\verts{a}}\exp\pars{-\pi\verts{a}\ic} - \verts{a}\bracks{\ln\pars{x} - \ic\pi} - 1 \over \pars{1 + x}^{2}}\,\dd x \\[5mm] = &\ -\sin\pars{\pi\verts{a}}\int_{0}^{1}{x^{\verts{a}} \over \pars{1 + x}^{2}} \,\dd x + \pi\verts{a}\int_{0}^{1}{\dd x \over \pars{1 + x}^{2}} \\[3mm] & = \half\,\pi\verts{a} - \sin\pars{\pi\verts{a}}\int_{0}^{1}{x^{\verts{a}} \over \pars{1 + x}^{2}}\,\dd x \end{align}

Integrar por partes: \begin{align} &\color{#f00}{\int_{0}^{\pi/2}{\sin^{2}\pars{ax} \over \sin^{2}\pars{x}}\,\dd x} \\[3mm] = &\ \half\,\pi\verts{a} + \half\,\sin\pars{\pi\verts{a}} - \verts{a}\sin\pars{\pi\verts{a}}\int_{0}^{1}{x^{\verts{a} - 1} \over 1 + x} \,\dd x \\[3mm] = &\ \half\,\pi\verts{a} + \half\,\sin\pars{\pi\verts{a}} + \verts{a}\sin\pars{\pi\verts{a}}\int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 + x} \,\dd x - \verts{a}\sin\pars{\pi\verts{a}}\int_{0}^{1}{\dd x \over 1 + x} \\[3mm] = &\ \half\,\pi\verts{a} + \bracks{\half - \verts{a}\ln\pars{2}} \sin\pars{\pi\verts{a}} + a\sin\pars{\pi a}\color{#00f}{\int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 + x} \,\dd x}\tag{1} \end{align}

El " $\color{#00f}{blue\ integral}$ " se evalúa de la siguiente manera: \begin{align} &\color{#00f}{% \int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 + x}\,\dd x} = \int_{0}^{1}\pars{1 - x^{\verts{a} - 1}} \pars{{1 \over 1 + x} + {1 \over 1 - x}}\,\dd x - \int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 - x}\,\dd x \\[3mm] = &\ 2\int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 - x^{2}}\,\dd x - \int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 - x}\,\dd x = \int_{0}^{1}{x^{-1/2} - x^{\verts{a}/2 - 1} \over 1 - x}\,\dd x - \int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 - x}\,\dd x \\[3mm] & = \int_{0}^{1}{1 - x^{\verts{a}/2 - 1} \over 1 - x}\,\dd x - \int_{0}^{1}{1 - x^{-1/2} \over 1 - x}\,\dd x - \int_{0}^{1}{1 - x^{\verts{a} - 1} \over 1 - x}\,\dd x \\[3mm] = &\ \bracks{\Psi\pars{{\verts{a} \over 2}} + \gamma}\ -\ \overbrace{\bracks{\Psi\pars{\half} + \gamma}}^{\ds{-2\ln\pars{2}}}\ -\ \bracks{\Psi\pars{\verts{a}} + \gamma} \\[3mm] = &\ \color{#00f}{% 2\ln\pars{2} + \bracks{\Psi\pars{{\verts{a} \over 2}} - \Psi\pars{\verts{a}}}} \end{align}

En la sustitución de este resultado en $\pars{1}$ encontramos: \begin{align} &\color{#f00}{\int_{0}^{\pi/2}{\sin^{2}\pars{ax} \over \sin^{2}\pars{x}}\,\dd x} \\[3mm] = &\ \color{#f00}{\half\,\pi\verts{a} + \bracks{\verts{a}\ln\pars{2} + \half}\sin\pars{\pi\verts{a}} + a\sin\pars{\pi a} \bracks{\Psi\pars{{\verts{a} \over 2}} - \Psi\pars{\verts{a}}}} \end{align}