透過 MCP 連接 →

輸入計算

數學公式

廣告

結果

回傳的 Tanh-Sinh 節點
10
node/weight pairs (order n = 20)
實際採用的 t_a 4.2
步長 h 0.442105
所有權重總和 1.9999993431
i t_i 節點 x_i 權重 w_i
1 0.2211 0.3364317911573048 0.6309622363150247
2 0.6632 0.8074765118645584 0.296772693493876
3 1.1053 0.9711342024624363 0.0662076633937352
4 1.5474 0.9982615398799233 0.0059249094153592
5 1.9895 0.9999745093540499 0.0001318520493753
6 2.4316 0.9999999602027466 0.0000003168563807
7 2.8737 0.999999999998166 0.0000000000226176
8 3.3158 1 0
9 3.7579 1 0
10 4.2 1 0
11

這個計算器的功能

Tanh-Sinh 求積節點與權重計算器會在標準區間 [-1, 1] 上,產生 Tanh-Sinh(又稱雙指數,DE)積分法所需的橫座標(節點)\(x_i\) 與對應的權重 \(w_i\)。取得這些數對之後,你就能把任意定積分化為一個簡單的加權和:f(x) 在 [-1, 1] 上的積分,近似等於所有 \(w_i\) 乘以 \(f(x_i)\) 的總和。

從 -1 到 1 的數線,求積節點向端點聚集
tanh-sinh 節點向 [-1, 1] 的端點聚集,能很好地處理奇異點。

方法與公式

Tanh-Sinh 求積法採用變數變換 \(x = \tanh\!\left(\tfrac{\pi}{2}\sinh t\right)\),把整條實數線 \(t\) 映射到開區間 (-1, 1)。變換後的被積函數會以雙指數速度衰減,因此最普通的梯形法則也能以驚人的速度收斂。把 \(t\) 截斷到 \([-t_a, t_a]\),再取 \(n\) 個等距取樣點、步長 \(h = 2 t_a / (n - 1)\),則每個點的座標為 \(t_i = -t_a + (i - 1) h\)、節點 \(x_i = \tanh\!\left(\tfrac{\pi}{2}\sinh t_i\right)\),權重則為 \(w_i = h\,\tfrac{\pi}{2}\cosh t_i\) 除以 cosh 的平方(\(\tfrac{\pi}{2}\sinh t_i\))。

$$ x_i = \tanh\!\left(\tfrac{\pi}{2}\sinh t_i\right), \qquad w_i = \frac{h\,\tfrac{\pi}{2}\cosh t_i}{\cosh^{2}\!\left(\tfrac{\pi}{2}\sinh t_i\right)} $$
Advertisement
tanh-sinh 變換 x 對 t 的圖形,以及對 t 的鐘形權重曲線
雙指數映射 x(t) 與快速衰減的權重 w(t) 作為 t 的函數。

使用方式

先選定階數 \(n\)(梯形取樣點的數量),再決定 \(t_a\) 要依你設定的精度自動推算,還是手動輸入,最後選擇要顯示幾位有效數字。在自動模式下,半寬度為 \(t_a = \mathrm{round}\!\left[(\text{digits} + 1)^{0.46},\,1\right]\);以 22 位為例,會得到文件中記載的預設值 \(t_a = 4.2\)。「僅半側」選項會利用對稱性 \(x_{-i} = -x_i\)、\(w_{-i} = w_i\),只回傳非負的一側;「全部」則會列出從接近 -1 到接近 +1 的每一個節點。

範例演算

當 \(n = 3\)、手動設 \(t_a = 4\),並選擇「全部」時:\(h = 8 / 2 = 4\)。三個 \(t\) 值分別為 -4、0、4。在 \(t = 0\) 時,\(x = \tanh(0) = 0\),$$ w = \tfrac{\pi}{2}\,h = 1.5707963 \times 4 = 6.2831853. $$在 \(t = \pm 4\) 時,\(\tfrac{\pi}{2}\sinh(4)\) 這個引數極大,因此 \(x\) 飽和趨近 \(\pm 1\),權重則因下溢而幾乎變成 0。當 \(n\) 較大且搭配適當的 \(t_a\) 時,權重總和會逼近 2,正好是 \(f = 1\) 在 [-1, 1] 上的精確積分值。

常見問題

為什麼兩端的權重幾乎為零?雙指數衰減會使 cosh 的平方(\(\tfrac{\pi}{2}\sinh t\))在邊緣附近趨於溢位,因此那些權重會消失——而這正是這套方法如此精確的原因。

這裡的「階數 n」是什麼意思?它指的是在 \([-t_a, t_a]\) 上等距取的梯形取樣點數量;點數越多並搭配合適的 \(t_a\),精度就越高。

我能在一般區間 [a, b] 上積分嗎?可以——只要重新縮放:代入 \(x = \tfrac{b - a}{2} x_i + \tfrac{a + b}{2}\),並把每個權重乘上 \(\tfrac{b - a}{2}\) 即可。

最後更新: