二重指数関数型数値積分公式

二重指数関数型数値積分公式（にじゅうしすうかんすうがたすうちせきぶんこうしき、テンプレート:Lang-en-short, 略してDE公式）とは変数変換に基づく数値積分の公式の一つである。この公式は森正武、高橋秀俊によって提案された。変換後の被積分関数が端点で二重指数関数的に減衰することが特徴である。数値積分の効率性の観点で、この公式がいろいろな点で使いやすく、非常に応用が利くと言われている。また、この公式は変換前の被積分関数が端点で特異性を持つときにも有効である。ただし、被積分関数によって適用できない場合があるので注意が必要である。

具体例

以下、いろいろな積分と、それに対応する二重指数関数型の変換を示す(テンプレート:Harvtxt)。

$\underset{-1}{\overset{1}{\int }}f(x)dx\Rightarrow x=\tanh (\frac{\pi }{2}\sinh t)$

$\underset{0}{\overset{\mathrm{\infty }}{\int }}f(x)dx\Rightarrow x=\exp (\frac{\pi }{2}\sinh t)$

$\underset{0}{\overset{\mathrm{\infty }}{\int }}{f}_1(x)\exp (-x)dx\Rightarrow x=\exp (t-\exp (-t))$

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(x)dx\Rightarrow x=\sinh (\frac{\pi }{2}\sinh t)$

台形公式への適用例

積分

$I=\underset{-1}{\overset{1}{\int }}f(x)dx$

の場合、変数変換

$x=\varphi (t)=\tanh (\frac{\pi }{2}\sinh t)$

によって積分は次のような形になる。

$I=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(\varphi (t)){\varphi }^{\prime }(t)dt$

これに、きざみ幅が等間隔$h$である台形公式を適用すると、

$I_h=h\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}f(\varphi (kh)){\varphi }^{\prime }(kh)$

を得る。さらに、この和を有限項までで打ち切ると、以下の数値積分公式が得られる：

$I_h^{(N)}=\frac{\pi }{2}h\sum _{k=-N_{-}}^{N_{+}}f(\tanh (\frac{\pi }{2}\sinh (kh)))\frac{\cosh (kh)}{{\cosh }^2(\frac{\pi }{2}\sinh (kh))},N=N_{-}+N_{+}+1$。

ここで、$N$は被積分関数$f(x)$の関数値を評価する回数である。$N_{-}$と$N_{+}$は、離散化誤差($\mathrm{\Delta }I=I-I_h$)と打ち切り誤差($\epsilon =I_h-I_h^{(N)}$)がほぼ等しくなるように決める。

特殊関数への応用

二重指数関数型積分公式は、ガンマ関数や変形ベッセル関数、行列値関数などの特殊関数の高精度計算・精度保証付き数値計算に応用されている。

関連項目

• 数値積分
• 台形公式
• 精度保証付き数値計算
• 森正武
• 高橋秀俊

参考文献

• テンプレート:Cite book
• Hidetosi Takahasi and Masatake Mori: Error Estimation in the Numerical Integration of Analytic Functions, Report of the Computer Centre, University of Tokyo, Vol.3 (1970), pp.41-108.
• H. Takahasi and M. Mori: Estimation of Errors in the Numerical Quadrature of Analytic Functions, Applicable Analysis, Vol.1 (1971) 201-229.
• Masatake Mori and Makoto Natori: Error Estimation in the Linear Approximation of Analytic Functions, Rep. Compt. Centre, Univ. Tokyo, Vol.4 (1971-1972), pp.1--17.
• 高橋秀俊「数値積分法の迷信」、数学セミナー、1971年3月号、日本評論社。また高橋秀俊：「数理と現象」，岩波書店（1975年1月10日）、頁157–168にも再録。
• テンプレート:Cite journal
• テンプレート:Cite book
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal

外部リンク

• DE-Sinc数値計算法