近似法

提供: testwiki
ナビゲーションに移動 検索に移動

近似法(きんじほう)とは関数の厳密値や方程式の厳密解を求めるときに、それが不可能または困難であるか、簡便のために近似値あるいは近似解を得る方法である。

無理数の近似

実数(特に無理数)を有理数により近似することはディオファントス近似として知られている[1]。例えば円周率 π=3.14159... の有理数による近似値として、古代エジプトでは 256/81 = 3.16049... が、古代バビロニアでは 25/8 = 3.125 が知られていた[2]アルキメデステンプレート:Indent を証明した[3](詳細は円周率の歴史を見よ)。ディオファントス近似は連分数と密接な関係がある[1]

テンプレート:詳細記事

テイラー展開

テイラー展開を用いる。

関数f (x ) のx = a の近傍における近似値を考える。f (x )をx = a においてテイラー展開すれば

f(x)=n=0f(n)(a)n!(xa)n

となる。x -a の値が十分小さければ、高次の項は無視することができる。とくに2次以上を無視すれば

f(x)f(a)+f(a)(xa)

となる。また、n 次の項まで考えたものをn 次近似と呼ぶ。すなわち上の例は1次近似である。

具体例

主要な関数のx0における2次近似を挙げておく。

  • ex1+x+x22
  • ln(1+x)xx22
  • (1+x)n1+nx+n(n1)2x2
  • sinxx
  • cosx1x22

漸近展開

特殊関数などの関数の振る舞いはしばしば漸近展開によって記述される。関数列 {φn}n=0 が点 x0 における漸近関数列であるとは、すべての nテンプレート:Indent が成り立つことをいう(点 x0 としては無限遠 x0= でもよい)[4]。そのうえで、ある関数 f の点 x0 における (広義の) 漸近級数とは、形式級数 n=0anφn(x) ですべての k に対して テンプレート:Indent を満足するもののことを言う[4]。このとき テンプレート:Indent と表記する[4]

漸近級数は一般には発散級数であってもよく、しかししばしばその有限項の和がもとの関数 f(x) をよく近似する[5]。例えば関数 テンプレート:Indentxスティルチェス級数を漸近展開として持つ[6]テンプレート:Indent f(10) の値はスティルチェス級数の第9項までの和により誤差 3.6×104 で評価できる[7]。ただしスティルチェス級数は発散級数であり、さらに多くの和を計算すると逆に誤差が増大する[7]

テンプレート:詳細記事

多項式近似

区間 [1,1] で定義された連続関数 f を多項式により近似することを考える。このような近似が可能であることはストーン=ワイエルシュトラスの定理によって保証されている。なお一般の区間 [a,b] の場合、変数変換 テンプレート:Indent により区間 [1,1] の問題に帰着できる[8]

ある多項式 p が問題の関数 f をどの程度「良く」近似できているのかは、区間 [1,1] における一様ノルム テンプレート:Indent によって表現される[9]。例えば区間 [a,b] を等間隔の n 区間 [xi,xi+1] (xi=a+(ba)i/n) に分割し、各小区間の端点 xi での関数値 fi(xi)ラグランジュ補間を近似多項式として採用することを考える。しばしばこの方式では区間 [1,1] の端付近で補間多項式の値が補間点以外ではもとの関数値から大きく外れるルンゲ現象が起き、その結果一様ノルムが大きな値を取るため好ましくない[10]

一様ノルムの観点で最良の近似多項式がミニマックス近似多項式である[9]n を高々 n 次の実係数多項式の全体とする。区間 [1,1] で定義された連続関数 f について、そのテンプレート:仮リンク ϕnPn とは、任意の pPn に対して テンプレート:Indent が成立するもののことを言う[9]。連続関数のミニマックス近似多項式は常に存在し一意であることが保証されている[9]。チェビシェフの等振動定理によると、ミニマックス多項式はもとの関数のある n+1 点でのラグランジュ補間に一致する[11]。ただし、どの点で補間を行えばよいのかは関数 f に依存する[11]。ミニマックス近似多項式を求めるRemezのアルゴリズムが知られている[12][13]。しかし、実際にそれを構成することはあまり容易ではない[9]

実用的には関数の多項式近似にはチェビシェフ補間を用いることが多い。n 次のチェビシェフ補間とは、Tn+1(x)n+1 個の零点における関数値をラグランジュ補間するものである。多項式近似におけるチェビシェフ多項式の有用性は次の定理によって示される[8]

n1 とする. 任意の実係数 n 次モニック多項式 p(x)テンプレート:Indent を満足する。等号成立は p(x)=21nTn(x) のときである.

従って n 次のモニック多項式のうち 0 を近似する最良のものがチェビシェフ補間であると言える[14]。この性質のため、チェビシェフ補間ではルンゲ現象は発生しない[11]。一般にはチェビシェフ補間がミニマックス近似多項式を与えるとは限らないものの、その良い近似を与えることが期待できる[15]

常微分方程式

本節では常微分方程式の近似解法について述べる。摂動論によって近似解を得ることができる状況がしばしばあるが、最高階微分項の係数が微小パラメータ ϵ であるような方程式の場合、ϵ=0 とした方程式ともとの方程式が質的に異なる特異摂動問題であり、通常の摂動論によって解くことができない[16]

境界層理論

流体力学ナビエ–ストークス方程式)における境界層のように、解析対象となる領域のうち、異なる領域によって有効な近似方程式が異なる場合がある。このような場合、各領域で境界条件を満足する適切な近似解を構成し、それを漸近接続することによって大域的な解を構成することができる[17]。この手法は境界層理論として知られている[18]

WKB近似

ϵ を微小パラメータとする2階常微分方程式 テンプレート:Indent は、A の変化がゆっくりとみなせるならば、指数関数型の解 テンプレート:Indent により近似できる[19]。この考察に基づく常微分方程式の漸近級数解の理論はWKB近似として知られる[20]。これはもともとジョゼフ・リウヴィルらによって19世紀から用いられていた[19]が、量子力学においてシュレディンガー方程式の近似解法としてWentzelらによって用いられたことからWKB近似として知られるようになった[21]。また光学幾何光学物理光学)における幾何光学近似あるいは物理光学近似とも対応している[22]

複スケール解析

複数の「時間」スケール(ここでは常微分方程式の独立変数を時間とみなす)を持つ問題の場合、独立変数として各時間スケールに対応する変数を導入することによって特異摂動問題を解く複スケール解析が知られている[23]。この手法は例えば擬調和振動子の摂動における永年項問題に適用できる[24]だけでなく、境界層問題等の他の特異摂動問題にも適用できる[25]

脚注

テンプレート:Reflist

参考文献

和書
  • テンプレート:Cite book
  • テンプレート:Cite book
  • 山内二郎、森口繁一、一松信 (編):「電子計算機のための数値計算法 II」、培風館(1967年3月5日)の第I編 "関数計算" の中の第1章、第2章、第3章。
  • 竹ノ内脩、西白保敏彦:「近似理論:関数の近似」、培風館(現代数学レクチャーズB-11)、ISBN 4-563-00431-6 (1986年6月10日).
  • 浜田穂積:「近似式のプログラミング」、培風館、ISBN 4-563-01382-X (1995年9月30日).
  • テンプレート:Cite book
  • テンプレート:Cite book
  • E. W. Cheney:「近似理論入門」、共立出版(1977年10月20日). ※ 原著は Introduction to Approximation Theory, McGraw-Hill (1966).
  • Theodore J. Rivlin: An Introductoion to the Approximation of Functions, Dover pub., ISBN 0-486-64069-8 (1981). ※初版は Blaisdel publishing (1969年)
  • 有本卓:「数値解析(1)」、コロナ社、ISBN 4-339-00124-4 (1981年4月30日)。 第5章「関数近似と補間法」。
洋書
  • E. W. Cheney: Multivariate Approximation Theory: Selected Topics, SIAM (CBMS 51), ISBN 0-89871-207-6 (1986).
  • Charles K. Chui: Multivariate Splines, SIAM (CBMS 54), ISBN 0-89871-226-2 (1988).

関連項目

テンプレート:Analysis-stub

  1. 1.0 1.1 Hensley, p. 13.
  2. 中村, pp. 51-52.
  3. 中村, pp. 53-54.
  4. 4.0 4.1 4.2 柴田, p. 33.
  5. 柴田, p. 22.
  6. 柴田, p. 23.
  7. 7.0 7.1 柴田, pp. 24-25.
  8. 8.0 8.1 Gil, Segura & Temme, p. 55.
  9. 9.0 9.1 9.2 9.3 9.4 Gil, Segura & Temme, p. 51.
  10. Gil, Segura & Temme, p. 54.
  11. 11.0 11.1 11.2 Gil, Segura & Temme, p. 63.
  12. テンプレート:Cite web
  13. 浜田望:「コンピューティングの玉手箱(80) ルメの最良近似アルゴリズム」、共立出版「bit」、vol.23,no.9, pp.1286-1287 (1991年8月号)
  14. Gil, Segura & Temme, p. 62.
  15. テンプレート:Cite book
  16. 柴田, p. 86.
  17. 柴田, p. 89.
  18. 柴田, p. 90.
  19. 19.0 19.1 柴田, p. 125.
  20. 柴田, pp. 124-125.
  21. 柴田, p. 120.
  22. 柴田, p. 126.
  23. 柴田, p. 151.
  24. 柴田, pp. 154-158.
  25. 柴田, pp. 151, 167-174.