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

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動先: 案内検索

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

具体例[編集]

以下、いろいろな積分と、それに対応する二重指数関数型の変換を示す。

台形公式への適用例[編集]

積分

の場合、変数変換

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

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

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

ここで、は被積分関数の関数値を評価する回数である。は、離散化誤差()と打ち切り誤差()がほぼ等しくなるように決める。

関連項目[編集]

参考文献[編集]

  • 森正武 『数値解析 (第2版)』 共立出版、2002年
  • 高橋秀俊:「数値積分法の迷信」、数学セミナー、1971年3月号、日本評論社。また高橋秀俊:「数理と現象」,岩波書店(1975年1月10日)、頁157-168にも再録。
  • Takahasi, H. and Mori, M. (1974). “Double exponential formulas for numerical integration”. Publications of the Research Institute for Mathematical Sciences (京都大学) 9 (3): 721--741. 
  • P.Rabinowitz, 森 正武 (翻訳):「計算機による数値積分法」、日本コンピュータ協会(コンピュータ・サイエンス研究書シリーズ 37)、 (1981年2月)。※原著の内容に、訳者による二重指数関数型数値積分法の解説を追加。
  • Ooura, T. and Mori, M. (1991). “The double exponential formula for oscillatory functions over the half infinite interval”. Journal of Computational and Applied Mathematics (Elsevier) 38 (1-3): 353--360. 
  • 森正武:二重指数関数型変換のすすめ (PDF) ,京都大学数理解析研究所講究録、第1040巻(1998),pp.143-153.
  • 森正武:数値解析における二重指数関数型変換の最適性, 数学(日本数学会),Vol.50, No.3 (1998), pp.248-264. doi:10.11429/sugaku1947.50.248
  • 渡辺二太(研究論文):「二重指数関数型数値積分公式について」 (PDF) (Comment on Double Exponential Formula), (1990年3月).
  • 大浦拓哉:DE公式と同じ漸近性能を持つ IMT型積分公式 (PDF)