微積分 · テイラー級数

テイラー級数

テイラー級数は、ある点での導関数の値を使い、多項式近似を順に作ります。項を増やし、中心を動かして、近似がよくなる範囲と級数が収束しなくなる境界を調べましょう。

01

この可視化でのテイラー級数

関数ffの点aaにおけるテイラー級数(テイラー展開)は、その点での導関数の値から作るべき級数です。a=0a = 0の場合はマクローリン級数とも呼びます。中心はx0x_0とも書きます。nn次までの部分和がテイラー多項式PnP_nです。これは、aaでの関数値とnn階までの導関数の値がffと一致する、次数がnn以下のただ1つの多項式です。

テイラー多項式
Pn(x)=∑k=0nf(k)(a)k! (x−a)k.P_n(x) = \sum_{k=0}^{n} \frac{f^{(k)}(a)}{k!}\,(x-a)^k.

濃い曲線はff、黄色い曲線はPnP_n、青い点は中心aaです。xx軸上の黄色い棒は測定結果で、∣f(x)−Pn(x)∣<0.01|f(x) - P_n(x)| < 0.01を満たす最も広い区間∣x−a∣<r|x - a| < rを示します。aaの両側を0.0025間隔で調べ、初めて条件を外れる位置をさらに細かく求め、値を切り下げて表示します。青い括弧はa±Ra \pm Rにあり、Rは収束半径です。級数は∣x−a∣<R|x - a| < Rで収束し、∣x−a∣>R|x - a| > Rで発散します。ここで扱う7つの関数では、この区間内の和はffそのものです。またRRは、aaからffの最も近い特異点までの距離で、複素数の点も含めます。

あるxxでPnP_nがffからどれだけ離れているかは、剰余項で測れます。ffがaaとxxの間の区間全体でn+1n + 1階まで微分可能なら(その区間は、1/(1−x)1/(1 - x)の1での極などの特異点をまたげません)、次が成り立ちます。

剰余項(ラグランジュ形)
f(x)−Pn(x)=f(n+1)(ξ)(n+1)! (x−a)n+1ある ξ(a と x の間)に対して。f(x) - P_n(x) = \frac{f^{(n+1)}(\xi)}{(n+1)!}\,(x-a)^{n+1}\quad\text{ある } \xi \text{(} a \text{ と } x\text{ の間)に対して。}
02

注目するところ

  • ここでは収束半径の内側で部分和がffに収束するため、その中の各xxで、nnを大きくすると誤差は0に近づきます。ただし、毎回必ず小さくなるとは限りません。外側では各項が0に近づかず、Pn(x)P_n(x)も一定の値に近づきません。棒はRRに近づく目安であり、画面上で測ったものです。一方、RRは式から求まります。
  • 関数が実数全域で滑らかでも、収束半径が有限になることがあります。1/(1+x2)1/(1 + x^2)は複素数の点±i\pm iで定義されないため、aaを中心とする級数が収束するのは∣x−a∣<1+a2|x - a| < \sqrt{1 + a^2}だけです。
  • 区間の端での挙動はさまざまです。ln⁡(1+x)\ln(1 + x)の0での級数は、x=1x = 1ではln⁡2\ln 2に収束しますが、x=−1x = -1では発散します。1/(1−x)1/(1 - x)の級数は両端で発散し、1+x\sqrt{1 + x}の級数は両端で収束します。
  • a=0a = 0では、sin⁡x\sin xには奇数べきだけ、cos⁡x\cos xには偶数べきだけが現れます。そのため、正弦のP2mP_{2m}とP2m−1P_{2m-1}は同じ多項式です。n → ∞はこうした0の項を飛ばすので、各ステップで曲線が変わります。
  • 収束するだけでは、和がffに等しいとは限りません。x≠0x \ne 0でe−1/x2e^{-1/x^2}、0で0となる関数は、0でのすべての階数の導関数が0です。そのため、0でのテイラー級数は、すべてのxxに対して0になります。
03

どこで使うか

  • 物理では、ある点の近くで最初の数項を使います。振り子の小さな振動にはsin⁡θ≈θ\sin\theta \approx \thetaを使い、(1−v2/c2)−1/2≈1+v2/(2c2)(1 - v^2/c^2)^{-1/2} \approx 1 + v^2/(2c^2)を使うと、低速での相対論的運動エネルギーは12mv2\tfrac12 mv^2になります。
  • 剰余項は、計算に必要な項数を教えてくれます。1/k!1/k!をk=0k = 0から1010まで足すと、これはexe^xのP10(1)P_{10}(1)であり、すでにeeを誤差3×10−83 \times 10^{-8}未満で近似できます。
04

操作方法

  • 次数のつまみを回すか、n → ∞を押すと、次の0でない項を1つずつ、n=30n = 30まで追加します。最大まで進んでもう一度押すと、0から再生し直します。青い点または軸上の印をドラッグすると中心が動き、代表的な中心のキーでは正確な値に移動します。
  • グラフの下の帯は、各xでの∣f−Pn∣|f - P_n|を対数目盛りで示します。プレゼンテーションモードはプロジェクター向けの全画面表示です。スペースキーで項を追加し、← →で次数を変更、Hで収束半径を隠します。
05

関連項目