Aproximación de$\prod _{k=p+1}^{\infty } \cos \left(\frac{p \,\pi}{2 k}\right)$

Question

Después de este post, empecé a preguntarme acerca de las posibles aproximaciones de la infinita producto $$A_p=\prod _{k=p+1}^{\infty } \cos \left(\frac{p \,\pi}{2 k}\right)\tag 1$$ where $p$ es un número entero. Según lo que pude ver, no hay forma cerrada expresiones.

Así, como lo hice en el vinculado pregunta, he utilizado Bhaskara I aproximación (por $-\frac \pi 2 \leq x\leq\frac \pi 2$) $$\cos(x) \simeq\frac{\pi ^2-4x^2}{\pi ^2+x^2}\implies \cos\left(\frac{p\,\pi}{2k}\right)=\frac{4 \left(k^2-p^2\right)}{4 k^2+p^2}$$ y, a continuación, calculadas $$B_p=\prod _{k=p+1}^{\infty }\frac{4 \left(k^2-p^2\right)}{4 k^2+p^2}=\frac{\Gamma \left( p+1-i\frac p 2\right)\,\, \Gamma \left( p+1+i\frac p 2\right)}{(2p)!}\tag 2$$ que no parece ser muy malo (véase la tabla abajo).

Usando la aproximación de Stirling y la serie de Taylor para la supuesta grandes valores de $p$, que terminó con $$\log(B_p)=\frac 12\log \left(\frac{5\pi}{4}\right)-\left(\log \left(\frac{16}{5}\right)+\cot ^{-1}(2)\right)p+\frac 12 \log(p)+\frac{11}{120 p}+O\left(\frac{1}{p^3}\right)\tag 3$$ Lo que es interesante es notar que $$\log \left(\frac{16}{5}\right)+\cot ^{-1}(2)\approx 1.62680$$ which is not so far from $\frac \pi 2$ ($3.56$% diferencia relativa), y esto ya de llevar las preguntas (al menos, para mí).

Lo que es interesante es ver cómo están los números en una escala logarítmica. $$\left( \begin{array}{cccc} p & (2) & (3) & (1) \\ 1 & -0.85190 & -0.85120 & -0.84448 \\ 2 & -2.17732 & -2.17725 & -2.16333 \\ 3 & -3.61661 & -3.61660 & -3.59727 \\ 4 & -5.10720 & -5.10719 & -5.08299 \\ 5 & -6.62701 & -6.62700 & -6.59815 \\ 6 & -8.16570 & -8.16570 & -8.13233 \\ 7 & -9.71760 & -9.71760 & -9.67979 \\ 8 & -11.2793 & -11.2793 & -11.2371 \\ 9 & -12.8485 & -12.8485 & -12.8019 \\ 10 & -14.4236 & -14.4236 & -14.3727 \\ 20 & -30.3496 & -30.3496 & -30.2557 \\ 30 & -46.4164 & -46.4164 & -46.2798 \\ 40 & -62.5413 & -62.5413 & -62.3621 \\ 50 & -78.6981 & -78.6981 & -78.4764 \end{array} \right)$$

Sólo por curiosidad, para $1 \leq p \leq 50$, I ajuste de los parámetros para el modelo de $$\log(A_p)=a+b\,p+c\log(p)+\frac d p$$ y obtuvo un buen ajuste $$\begin{array}{clclclclc} \text{} & \text{Estimate} & \text{Standard Error} & \text{Confidence Interval} \\ a & +0.69204 & 0.00006 & \{+0.69193,+0.69216\} \\ b & -1.62255 & 0.00000 & \{-1.62256,-1.62255\} \\ c & +0.50040 & 0.00003 & \{+0.50035,+0.50045\} \\ d & +0.08598 & 0.00007 & \{+0.08583,+0.08613\} \\ \end{array}$$ mientras se $(3)$, los coeficientes se $(+0.68394,-1.62680,+0.50000,+0.09167)$ es decir muy muy cerca.

Por cierto, todo esto muestra la alta calidad de la aproximación de $\cos(x)$ pero me pregunta sobre el valor de $A_p$. Si queremos truncar $(3)$ a $O\left(\frac{1}{p}\right)$,debemos tener $$A_p\sim \sqrt{\frac{5\pi p} 4 }\,e^{-\alpha \pi p}$$ with $\alpha \approx \frac 12$.

Me pregunto si podríamos encontrar una mejor aproximación de $A_p$ , incluso al precio de funciones más complejas. Cualquier idea será bienvenida.

Answer 1

Teniendo en cuenta que Bhaskara I aproximación $$\cos(x) \simeq\frac{\pi ^2-4x^2}{\pi ^2+x^2}\qquad \text{for} \qquad -\frac \pi 2 \leq x\leq\frac \pi 2$$ looks like a Padé approximant, I had the feeling that I could either pure Padé approximants built at $x=0$, es decir, $$f_n=\frac {1+\sum_{m=1}^n a_m x^{2m}}{1+\sum_{m=1}^n b_m x^{2m}}$$ o, puede ser mejor, $$g_n=\left(1-\frac{2x}{\pi }\right) \left(1+\frac{2x}{ \pi }\right)\frac {1+\sum_{m=1}^{n-1} c_m x^{2m}}{1+\sum_{m=1}^n d_m x^{2m}}$$ the second rational fraction being the Padé approximant of $\frac{\cos (x)}{\left(1-\frac{2x}{ \pi }\right) \left(1+\frac{2x}{ \pi }\right)}$; $g_n$ was chosen as potentially better than $f_n$ because the calculations will start closer and closer to $\frac \pi 2$. Para confirmar la validez de esta elección fueron calculadas $$I_n=\int_0^{\frac \pi 2} \big(\cos(x)-f_n\big)^2\,dx \qquad \text{and}\qquad J_n=\int_0^{\frac \pi 2} \big(\cos(x)-g_n\big)^2\,dx$$ La tabla a continuación muestra claramente la superioridad de la $g_n$. $$\left( \begin{array}{ccc} n & I_n & J_n \\ 1 & 7.086 \times 10^{-5} & 7.962 \times 10^{-6} \\ 2 & 6.754 \times 10^{-11} & 1.909 \times 10^{-12} \\ 3 & 2.375 \times 10^{-18} & 3.014 \times 10^{-20} \end{array} \right)$$ while $$\int_0^{\frac \pi 2} \Big(\cos(x)-\frac{\pi ^2-4x^2}{\pi ^2+x^2}\Big)^2\,dx =1.489\times 10^{-6}$$

Ahora, en sustitución de $x$ por $\frac {p \pi}{2k}$hace $$\cos \left(\frac{\pi p}{2 k}\right)=\frac {P_{m}(k^2)}{Q_{m}(k^2)}$$ where $P_m$ and $Q_m$ are homogeneous polynomials. $P_m$ shows only one real root and $Q_m$ only complex roots. By the end, limited to $n=2$, podemos escribir $$\cos \left(\frac{\pi p}{2 k}\right)=\frac{(k-r_1 p)\,(k+r_1p)\,(k-r_2p)\,(k+r_2p)} {(k-s_1 p)\,(k+s_1p)\,(k-s_2p)\,(k+s_2 p) }$$ y si $$A_p=\prod _{k=p+1}^{\infty } \cos \left(\frac{p \,\pi}{2 k}\right)$$, a continuación, $$A_p\sim\frac{\Gamma [1+(1-s_1) p]\, \Gamma [1+(1+s_1) p]\, \Gamma [1+(1-s_2) p]\, \Gamma [1+(1+s_2) p]}{\Gamma [1+(1-r_1) p]\, \Gamma [1+(1+r_1) p]\,\Gamma [1+(1-r_2) p]\, \Gamma [1+(1+r_2) p]}$$

Esto conduce a los siguientes resultados (para los logaritmos) $$\left( \begin{array}{cccc} p & \text{using } f_n & \text{using } g_n& \text{exact} \\ 1 & -0.844481 & -0.844481 & -0.844481 \\ 2 & -2.163326 & -2.163329 & -2.163327 \\ 3 & -3.597265 & -3.597275 & -3.597270 \\ 4 & -5.082976 & -5.082997 & -5.082989 \\ 5 & -6.598132 & -6.598168 & -6.598155 \\ 6 & -8.132292 & -8.132346 & -8.132328 \\ 7 & -9.679735 & -9.679810 & -9.679787 \\ 8 & -11.23699 & -11.23708 & -11.23706 \\ 9 & -12.80178 & -12.80190 & -12.80187 \\ 10 & -14.37255 & -14.37270 & -14.37266 \\ 20 & -30.25528 & -30.25578 & -30.25570 \\ 30 & -46.27897 & -46.27993 & -46.27979 \\ 40 & -62.36080 & -62.36228 & -62.36208 \\ 50 & -78.47457 & -78.47661 & -78.47636 \end{array} \right)$$ que se ve mucho mejor que en mi post.

Answer 2

Como por los comentarios, pero con referencias y enlaces. El uso de esta desigualdad y debido a $0<\frac{p \pi}{2k}<\frac{\pi}{2}$ para $k\geq p+1$, tenemos $$\color{red}{\prod\limits_{k=p+1}\cos{\frac{p\pi}{2k}}\leq} \frac{1}{e^{\frac{\pi^2 \cdot p^2}{8}\sum\limits_{k=p+1}\frac{1}{k^2}}}= \frac{1}{e^{\frac{\pi^2 \cdot p^2}{8}\left(\zeta(2)-\sum\limits_{k=1}^p\frac{1}{k^2}\right)}}=... \etiqueta{1}$$ y (por ejemplo, Apostol "Introducción a la Teoría Analítica de números", página 55) $$\sum\limits_{n\leq x}\frac{1}{n^s}=\frac{x^{1-s}}{1-s}+\zeta(s)+O(x^{-s}) \tag{2}$$ para $s>0, s\ne 1$, llevando a $$\sum\limits_{n\leq p}\frac{1}{n^2}=-\frac{1}{p}+\zeta(2)+O\left(\frac{1}{p^2}\right) \tag{3}$$ y, sustituyendo en la $(1)$ $$...=\frac{1}{e^{\frac{\pi^2 \cdot p^2}{8}\left(\frac{1}{p}-O\left(\frac{1}{p^2}\right)\right)}}= \frac{1}{e^{\frac{\pi^2 \cdot p}{8}-\frac{\pi^2 \cdot O(1)}{8}}}= \color{red}{\frac{e^{\frac{\pi^2}{8}\cdot O(1)}}{e^{\frac{\pi^2}{8}\cdot p}}}= \frac{O(1)}{e^{\frac{\pi^2}{8}\cdot p}}$$

Lo que conduce a $\sim \frac{1}{e^{\alpha \pi p}}$ en lugar de $\sim \frac{\sqrt{p}}{e^{\alpha \pi p}}$.

---

Sin embargo, este no parece ser muy importante para un gran $p$, simplemente porque (yo uso $-b$ a ha $b=1.62255>0$): $$A_p=e^{\ln{A_p}}= e^{a-bp+c\ln{p}+\frac{d}{p}}= e^{a}\cdot e^{\frac{d}{p}} \cdot \frac{p^c}{e^{pb}}= O(1)\cdot \frac{1}{e^{(b-c)p}}\cdot \frac{p^c}{e^{cp}}\leq ...$$ debido a $\frac{p^c}{e^{cp}}\rightarrow 0, p\rightarrow\infty$ $$...\leq \frac{O(1)}{e^{(b-c)p}}$$ Pero $\frac{\pi^2}{8}> b-c$, por lo tanto, en una forma $$\frac{O(1)}{e^{\frac{\pi^2}{8}\cdot p}}<\frac{O(1)}{e^{(b-c)p}}$$ desde algunos lo suficientemente grande como $p$ en adelante.