级数求和与定积分计算报告

给定级数

要求完成以下五件事：

1. 计算 $S(1000)$ 的数值，并说明保留的精度；
2. 求 $N\to\infty$ 时 $S(N)$ 的极限，给出精确表达式和数值；
3. 计算定积分 $\displaystyle\int_{0}^{1}\frac{\ln(1+x)}{x}\,\mathrm{d}x$ 的精确值；
4. 判断第 3 题的定积分与第 2 题的极限之间是什么关系；
5. 对所有结果执行三重验证（代回反算 / 量级检查 / 边界试算），给出验证结论；
6. 用专业版和通俗版解释：该级数为什么收敛，以及第 3 题积分为什么与第 2 题极限有关系。

已知条件清单：求和指标 $k\in\{1,\dots,N\}$；积分区间 $[0,1]$；被积函数在 $x=0$ 处有可去奇点（$\lim_{x\to0^+}\ln(1+x)/x=1$）。

隐含约束：报告数值须来源可追溯；中间结果不中途舍入；精度特征（有效数字 / 机器精度 / 误差界）须标注。

本报告采用三通道并行 + 交叉验证策略：

验证方案设计：对每个结果执行三类验证——① 代回反算（用独立方法/恒等式重算并比对残差）；② 量级检查（与上下界、常识基准比对）；③ 边界试算（特殊值、极端参数、尾部估计）。

计算过程

按定义直接求和。为审计精度，用三种独立方式同时计算并交叉比对：

1. 精确有理数：$S(1000)=\sum_{k=1}^{1000}\frac{1}{k^2}$ 逐项化为 $\frac{a}{b}$（分子、分母各 866 位十进制数），代数精确、零舍入；
2. mpmath 60 位十进制浮点累加（mp.fsum，高精度补偿求和）；
3. IEEE 754 双精度math.fsum（对比工程精度）。

两种高精度结果逐位一致：

精确有理数   S(1000) = a/b   （分子/分母各 866 位）
mpmath 60dps S(1000) = 1.64393456668155980313905802382221558965210344649368531671724
精确有理数在 60 位下  = 1.64393456668155980313905802382221558965210344649368531671724
  二者之差              = 0.0   （逐位一致，0 误差）
float64 fsum S(1000) = 1.6439345666815597

结果（保留 30 位有效数字）

精度说明

结论：报告中给出的 30 位有效数字均经过"精确有理数 + 60 位高精度"双通道逐位核对，所有 30 位数字可靠。浮点双精度最多只能保证约 16 位。

精确表达式

这是经典巴塞尔问题（Basel problem）：1650 年前后由 Pietro Mengoli 提出，1735 年由欧拉（L. Euler, De summis serierum reciprocarum, E41）首次证明。

证明要点（Fourier 级数法，一个现代标准证法）

$x^2$ 在 $[-\pi,\pi]$ 上的 Fourier 级数为

令 $x=\pi$，因 $\cos(n\pi)=(-1)^{n}$，故 $(-1)^{n}\cos(n\pi)=1$，得

（欧拉的原证基于 $\sin x$ 的无穷乘积展开与系数比较，结论相同。）

数值（60 位）

pi^2/6  (mpmath 60dps) = 1.64493406684822643647241516664602518921894990120679843773556
zeta(2) (mpmath 60dps) = 1.64493406684822643647241516664602518921894990120679843773556   ← 一致
pi^2/6  (SymPy 50dps)  = 1.6449340668482264364724151666460251892189499012068
pi^2/6  (float64)      = 1.6449340668482264

精确表达式

推导（幂级数展开 + 逐项积分）

对 $x\in(-1,1]$ 有 $ \ln(1+x)=\sum_{n=1}^{\infty}(-1)^{n+1}\dfrac{x^{n}}{n} $，故

级数 $\sum (-1)^{n+1}x^{n-1}/n$ 在 $[0,1]$ 上一致收敛（Abel 判别法；$x=0$ 处为可去奇点、和为 $1$），可逐项积分：

其中 $\eta(s)=\sum_{n\ge1}(-1)^{n+1}/n^{s}$ 为 Dirichlet eta 函数。由恒等式 $\eta(s)=(1-2^{1-s})\,\zeta(s)$，取 $s=2$：

符号计算直接验证

SymPy 精确积分 : integrate(log(1+x)/x, (x,0,1))  =  pi**2/12      ← 符号级精确
eta(2) - zeta(2)/2 = 0                                        ← 恒等式精确成立
数值           : 0.822467033424113218236207583323012594609474950603399218867779
mpmath quad    : 0.822467033424113218236207583323012594609474950603399218867779
scipy.quad     : 0.8224670334241132    （估计误差 ≈ 9.1e-15）

关系结论

即：第 3 题的定积分恰好等于第 2 题极限（巴塞尔级数和）的一半，比例严格为 $1:2$。

为什么是"恰好一半"（奇偶项分离）

记 $\zeta(2)=\sum_{n=1}^{\infty}1/n^{2}$，将其拆成奇数项与偶数项：

直觉：在 $\eta(2)$ 中，所有偶数项被交替符号翻成负号；而偶数项之和恰为整体的 $1/4$，因此 $\zeta-\eta$ 恰等于 $2\cdot(1/4)\zeta=\zeta/2$，即 $\eta=\zeta/2$。等价地，从 $\eta(2)=\frac12\zeta(2)$ 看，系数 $1-2^{1-2}=1/2$ 直接给出比值。

数值核验

ratio  (pi^2/12) / (pi^2/6) = 0.5
|2·I − L| = 0.0          （|2×0.8224670334241132… − 1.6449340668482264…| = 0）
pi^2/12 × 2 = 1.64493406684822643647241516665  =  pi^2/6  ✔

辅助：SymPy Sum 直接求 $S(1000)=1.643934566681559803139058023822215589652\cdots$，与主结果逐位一致。

验证总结论

三个数值对象（$S(1000)$、$\pi^{2}/6$、$\pi^{2}/12$）与一个关系式（$I=\frac12\lim S$）的全部 15 项检验均通过（PASS），无任何异常，无需回炉。所有数值均来自真实工具执行结果，可复现、可追溯。

容差说明：代回反算以机器精度 $\varepsilon\approx2.2\times10^{-16}$ 为基准放大 10~100 倍；60 位通道内部以 $10^{-47}$ 为误差界；所有比对残差均远小于容差。

专业版 为什么 $\sum 1/k^{2}$ 收敛

由单调收敛定理：$S(N)$ 单调递增。对 $k\ge2$，利用分项不等式

于是

部分和有上界 2。单调递增且有上界 ⇒ 收敛（Cauchy/单调收敛定理）。更精确地，积分判别给出收敛半径相同、且 $ \sum_{k=N+1}^{\infty}\frac1{k^{2}}\sim\frac1N $。该极限不是任意常数——由 Fourier 级数或 $\sin x$ 的无穷乘积可知其精确值 $=\zeta(2)=\frac{\pi^{2}}{6}$（Euler, 1735）。

专业版 为什么积分 $=\frac12\lim S$

被积函数 $ \frac{\ln(1+x)}{x} $ 的幂级数在 $[0,1]$ 一致收敛（Abel 判别），故逐项积分合法，得到交错平方倒数级数：

而奇偶分离表明 $ \eta(2)=\zeta(2)-2\sum_{m}\frac1{(2m)^{2}}=\zeta(2)-\frac12\zeta(2)=\frac12\zeta(2) $，即积分恰为巴塞尔和的一半。本质上：积分以"面积"的形式实现了交错求和，而交错求和仅"翻转"偶数项，由于偶数项占全体的 $1/4$，两次翻转后净剩 $1/2$。

通俗版 为什么这个级数会"加得完"（收敛）

想象往盒子里倒越来越小的沙粒：第 1 颗 = 1，第 2 颗 = 0.25，第 3 颗 ≈ 0.111，第 4 颗 = 0.0625……每颗都是"上一颗的几分之一"，缩得飞快。

• 不超线：把从第 2 颗往后的每一颗都"凑"成两块之差（$\frac1{k^2}<\frac1{k-1}-\frac1k$），这些差像多米诺骨牌一样逐项抵消（裂项），加完只剩"第 1 块 − 最后 1 块"，所以总和永远 < 1 + 1 = 2。
• 只增不减：每加一项和都在变大，却又永远超不过 2。
• 结论：一个只增不减、又封顶的数字，必然"停下来"。它停在哪？——π²/6 ≈ 1.6449，这叫巴塞尔问题，欧拉 1735 年算出来的，数学史上最有名的求和之一。

通俗版 为什么那个积分和这个极限"有关系"

把对数展开成"无限长多项式"：$\ln(1+x)=x-\frac{x^2}{2}+\frac{x^3}{3}-\frac{x^4}{4}+\cdots$。除以 $x$ 后是 $1-\frac{x}{2}+\frac{x^2}{3}-\frac{x^3}{4}+\cdots$，再逐块求面积（积分），每一项 $\frac{x^{n-1}}{n}$ 都贡献 $1/n^2$，于是

而原来的级数是 $1+\frac14+\frac19+\frac1{16}+\cdots$（全加号）。差别只在一个地方：偶数项 $1/4,1/16,1/36,\dots$ 的符号被"翻"成了减号。

• 偶数项占整个和的 $1/4$（因为 $\frac1{(2m)^2}=\frac14\cdot\frac1{m^2}$）。
• 把偶数项从"加"改成"减"，相当于把 $1/4$ 的贡献减去两次：原来算一份，现在先去掉、再倒扣，净少了一半。
• 所以交错和恰好 = 非交错和的一半：$\pi^2/12 = \frac12\times\pi^2/6$。

一句话类比：原始级数是一张 $π²/6$ 的总账单；那个积分是同一张账单的"半价券"——它把带偶数编号的条目全部翻成负的，正好抵消掉一半。

外部依据核验（联网检索）：巴塞尔问题的历史与证明（Mengoli 提出、Euler 1735 证明、$\sum1/n^2=\pi^2/6$）[cite:d3f75d5f-57][cite:d3f75d5f-11][cite:d3f75d5f-83]；Dirichlet eta 函数定义与 $\eta(s)=(1-2^{1-s})\zeta(s)$、$\eta(2)=\pi^2/12$[cite:a834dec9-2][cite:a834dec9-5][cite:a834dec9-6]。以上属 A 级标准结果，多来源一致。其余公式（幂级数展开、Abel 逐项积分、Fourier 级数证法）为训练知识中的标准教材内容，均以本报告数值计算独立复核。

补充说明

• 精度限制：双精度（IEEE 754, GB/T 17966）只能保证约 16 位有效数字；如需更多位数应使用 mpmath/任意精度通道（本报告已如此执行）。
• 特殊值处理：$x=0$ 处被积函数为 $0/0$ 型，按可去奇点处理，极限值 $1$；本报告所有数值工具均据此正确计算，未产生 NaN/Inf/除零异常。
• 误差提示：$S(1000)$ 距极限的差约为 $9.995\times10^{-4}$（尾差 $\sim1/N$），若要用部分和逼近 $\pi^2/6$ 需按 $1/N$ 量级估计误差；用 Euler–Maclaurin 尾修正可显著改善。
• 扩展方向：$\zeta(2k)$ 的一般值均为 $\pi^{2k}$ 与 Bernoulli 数的有理组合；类似技巧可算 $\int_0^1\ln(1-x)/x\,\mathrm dx=-\pi^2/6$、$\int_0^1\frac{\ln(1+x)}{x}\mathrm dx$ 的推广等。

本报告公式在支持 LaTeX 渲染的环境（Typora / Obsidian / Jupyter / 浏览器 + MathJax）中显示；若离线打开，公式将以原始 LaTeX 文本显示，全部关键数值已在正文中以纯文本形式给出。