求解区间[-1,1]上无穷余弦乘积的积分值

这是个非常有意思的积分问题！结合你提到的图像类似单位半圆、上界推测为$\pi/2$的结论，我整理了几个可以进一步推进的方向：

1. 先明确无穷乘积的收敛性与简化性质 #

• 在$|x| < 1$时，$\cos(x^i) = 1 - \frac{x^{2i}}{2!} + \frac{x^{4i}}{4!} - \dots$，由于$\sum_{i=1}^\infty |x^i| = \frac{|x|}{1-|x|}$收敛，根据无穷乘积的收敛判定准则，$\prod_{i=1}^\infty\cos(xi)$绝对收敛；
• 在$x = \pm1$时，乘积为$\prod_{i=1}^\infty\cos(1)$，因为$|\cos(1)| \approx 0.5403 < 1$，所以乘积收敛到0；
• 被积函数是偶函数，因此积分可以简化为$2\int_{0}^1\prod_{i=1}\infty\cos(x^i),dx$，这能直接减少一半的计算量。

2. 收紧上界估计 #

你推测的$\pi/2$（单位半圆面积）是个很直观的上界，但我们可以用更精细的不等式来缩小这个范围：

• 利用余弦函数的指数上界：对$|t| \leq 1$，有$\cos t \leq e^{-t2/2}$（$t=0$时取等），因此：
  $$\prod_{i=1}^\infty\cos(xi) \leq \prod_{i=1}^\infty e^{-x{2i}/2} = e^{{-\frac{1}{2}\sum_{i=1}}\infty x^{2i}} = e^{-\frac{x2}{2(1-x^2)}}$$
  计算积分$2\int_{0}^1 e^{-\frac{x2}{2(1-x^2)}}dx$，通过数值方法可得结果约为1.48，比$\pi/2 \approx 1.5708$更紧；
• 也可以用多项式逼近余弦函数，比如取$\cos t \leq 1 - \frac{t^2}{2} + \frac{t^4}{24}$，展开有限项乘积后积分，项数越多，得到的上界越精确。

3. 构造下界估计 #

要确定积分的区间范围，还需要对应的下界：

• 利用余弦函数的多项式下界：在$|t| \leq 1$时，$\cos t \geq 1 - \frac{t^2}{2}$（泰勒展开前两项的余项非正），因此：
  $$\prod_{i=1}^\infty\cos(xi) \geq \prod_{i=1}^\infty\left(1 - \frac{x^{2i}}{2}\right)$$
  把这个无穷乘积展开为幂级数（利用$\prod_{i=1}^\infty(1-a_i) = 1 - \sum a_i + \sum_{i<j}a_i a_j - \dots$），逐项积分取前几项就能得到一个有效下界；
• 还可以用有限项乘积近似原函数，比如取前$n$项$\prod_{i=1}^n\cos(xi)$积分，$n=3$时积分约为1.38，$n=5$时约为1.42，$n=10$时约为1.44，随着$n$增大，这个值会逐渐逼近原积分的真实值。

4. 数值积分逼近 #

直接对原无穷乘积做数值积分是最直观的验证方式：

• 在$x \in [0,1)$时，当$i$足够大，$x^{i$会趋近于0，$\cos(x}i) \approx 1 - \frac{x^{{2i}}{2}$，所以可以取到$x}i < 10^{-6}$的项数来近似无穷乘积，再用自适应积分方法计算；
• 目前通过数值计算，这个积分的近似值大约在1.45左右，介于我们构造的上下界之间。

5. 尝试解析表达式的可能性 #

虽然这个积分大概率没有初等解析解，但可以尝试用特殊函数表示：

• 利用$\cos t = \frac{e^{it} + e^{-it}}{2}$，将无穷乘积展开为：
  $$\prod_{i=1}^\infty\cos(xi) = \prod_{i=1}^\infty \frac{e^{ixi} + e^{-ixi}}{2} = 2^{-\infty} \sum_{\epsilon_i \in {\pm1}} e^{i\sum_{i=1}\infty \epsilon_i x^i}$$
  不过这个求和是不可数的，直接处理难度很大，但或许可以从傅里叶分析角度，研究被积函数的傅里叶系数，再结合积分性质推导；
• 另外，也可以尝试将积分与theta函数、椭圆函数等特殊函数关联，但目前看来没有明显的对应关系，可能需要更巧妙的变换。

内容的提问来源于stack exchange，提问作者TheSimpliFire