通过MCP连接 →

输入计算

数学公式

广告

结果

积分值
0.7468241328
定积分在 [a, b] 上的近似值
方法 Tanh-Sinh(双指数)
细分数 n 100
函数求值次数 201

什么是 Tanh-Sinh 求积?

Tanh-Sinh 求积,又称双指数(DE)公式,是一种擅长计算有限区间 \([a, b]\) 上定积分的数值积分方法——尤其是当被积函数在端点处存在奇异性时表现突出。它通过变量替换 \(u = \tanh\!\left(\tfrac{\pi}{2}\sinh(t)\right)\),把区间端点映射到 \(t = \pm\infty\)。在这些端点附近,被积函数的贡献会以双指数速度衰减,因此即便是在 \(a\) 或 \(b\) 处发散的函数(例如 \(\tfrac{1}{\sqrt{1-x^2}}\)),也能被精确地积分出来。本工具属于通用数学方法,世界各地都适用。

S 形的 tanh-sinh 变换曲线,向负一和正一处的渐近线趋于平缓
双指数变量替换将实数轴映射到 \((-1, 1)\),使节点向两端聚集。

使用方法

用标准记法输入你的函数 \(f(x)\)(运算符 + - * / ^、括号,以及 sin、cos、exp、log、sqrt、abs 等函数,还可使用常数 pi 和 e)。再填写积分下限 \(a\)、上限 \(b\),以及控制节点密度的细分数 \(n\)。\(n\) 越大精度越高,但计算量也随之增加;常用的实际取值范围为 50~400。被积函数需在开区间内解析(端点处有奇异性没关系),且不能是周期函数。

公式解析

首先用 \(x = \tfrac{b-a}{2}u + \tfrac{a+b}{2}\)、\(dx = \tfrac{b-a}{2}\,du\) 把区间归一化到 \([-1, 1]\)。随后 DE 公式采用节点 \(t_k = k\,h\)、横坐标 \(u_k = \tanh\!\left(\tfrac{\pi}{2}\sinh(t_k)\right)\) 以及权重 \(w_k = \dfrac{\tfrac{\pi}{2}\cosh(t_k)}{\cosh^2\!\left(\tfrac{\pi}{2}\sinh(t_k)\right)}\)。积分近似为 $$\int_{a}^{b} f(x)\,dx \approx \frac{b-a}{2}\, h \sum_{k=-N}^{N} w_k\, f\!\left(\frac{a+b}{2} + \frac{b-a}{2}\,x_k\right)$$ 当某个节点的权重下溢为零(即饱和端点)时会被跳过,从而避免在可能奇异的边界处求值。

从 a 到 b 的区间,阴影区域下方两端附近密集聚集的求积节点
节点在端点 \(a\) 和 \(b\) 附近聚集,因此能精确处理端点奇异性。

实例演算

计算 \(f(x) = \exp(-x^2)\) 在 \([0, 1]\) 上的积分。精确值为 \(\tfrac{\sqrt{\pi}}{2}\cdot\operatorname{erf}(1) \approx 0.7468241\)。即使步长较粗(\(h = 0.5\),\(N = 4\)),该公式也能给出约 \(0.7467\);而采用默认值 \(n = 100\) 时,结果可精确到大约十二位有效数字。

常见问题

可以对端点有奇异性的函数积分吗? 可以——这正是该方法的核心优势。\(a\) 或 \(b\) 处的可积奇异性都能被妥善处理。

为什么周期性会有影响? 双指数公式是为非周期被积函数量身设计的;对于周期函数,梯形公式收敛得更快,而 DE 公式可能并不准确。

如果 a 等于 b 会怎样? 积分为零。如果 \(a > b\),由于因子 \(\tfrac{b-a}{2}\) 携带符号,结果会被正确地取负。

最后更新: