テイラーの公式

Basis
最終更新: タグ: 微積分, 微分

平均値の定理x0x_0xx の間のある未知の cc に対して f(x)=f(x0)+f(c)(xx0)f(x) = f(x_0) + f'(c)(x - x_0) と述べる。これは精度の粗い剰余項を持つ線形近似だ。テイラーの公式(Taylor’s Formula)はこの考えを全次数に拡張する:ffnn 次多項式で近似し、残りの誤差を明示的に表す。

テイラー多項式

定義。 ffx0x_0 における**nn 次テイラー多項式(Taylor polynomial of degree nn)**を

Pn(x)    k=0nf(k)(x0)k!(xx0)k=f(x0)+f(x0)(xx0)+f(x0)2!(xx0)2++f(n)(x0)n!(xx0)nP_n(x) \;\coloneqq\; \sum_{k=0}^{n} \frac{f^{(k)}(x_0)}{k!}\,(x - x_0)^k = f(x_0) + f'(x_0)(x-x_0) + \frac{f''(x_0)}{2!}(x-x_0)^2 + \cdots + \frac{f^{(n)}(x_0)}{n!}(x-x_0)^n

と定める。

係数 f(k)(x0)k!\dfrac{f^{(k)}(x_0)}{k!}Pn(k)(x0)=f(k)(x0)P_n^{(k)}(x_0) = f^{(k)}(x_0)k=0,1,,nk = 0, 1, \ldots, n)を成り立たせる唯一の値だ:この多項式は x0x_0 において ff とその nn 次までの全導関数の値が一致する。

ラグランジュ剰余項によるテイラーの公式

定理(テイラーの公式)。 ffx0x_0xx を含む開区間上で n+1n+1 回微分可能とする。このとき

f(x)  =  Pn(x)  +  Rn(x),(1)f(x) \;=\; P_n(x) \;+\; R_n(x), \tag{1}

が成り立つ。ここで**ラグランジュ剰余項(Lagrange remainder)**は

Rn(x)  =  f(n+1)(ξ)(n+1)!(xx0)n+1(2)R_n(x) \;=\; \frac{f^{(n+1)}(\xi)}{(n+1)!}\,(x - x_0)^{n+1} \tag{2}

であり、ξ\xix0x_0xx の間のある点だ。

証明

xx0x \neq x_0 を固定し、JJx0x_0xx を端点とする閉区間とする。次の関数を定義する:

F(t)    f(x)k=0nf(k)(t)k!(xt)k,G(t)    (xt)n+1.F(t) \;\coloneqq\; f(x) - \sum_{k=0}^{n} \frac{f^{(k)}(t)}{k!}(x - t)^k, \qquad G(t) \;\coloneqq\; (x - t)^{n+1}.

FFGG はともに JJ 上で連続、JJ の内部で微分可能だ。次のことに注目する:

  • F(x)=0F(x) = 0 かつ G(x)=0G(x) = 0
  • F(x0)=f(x)Pn(x)F(x_0) = f(x) - P_n(x) かつ G(x0)=(xx0)n+10G(x_0) = (x - x_0)^{n+1} \neq 0

FF の導関数の望遠鏡和。 積の微分法則を各項 f(k)(t)k!(xt)k\tfrac{f^{(k)}(t)}{k!}(x-t)^k に適用すると:

ddt[f(k)(t)k!(xt)k]=f(k+1)(t)k!(xt)k    f(k)(t)(k1)!(xt)k1.\frac{d}{dt}\left[\frac{f^{(k)}(t)}{k!}(x-t)^k\right] = \frac{f^{(k+1)}(t)}{k!}(x-t)^k \;-\; \frac{f^{(k)}(t)}{(k-1)!}(x-t)^{k-1}.

k=0k = 0 から nn まで和をとると、望遠鏡のように各項 f(k)(t)(k1)!(xt)k1\tfrac{f^{(k)}(t)}{(k-1)!}(x-t)^{k-1}kk1k \mapsto k-1 に対応する第一項と消え合い、k=nk = n の第一項だけが残る:

F(t)  =  f(n+1)(t)n!(xt)n,G(t)  =  (n+1)(xt)n.F'(t) \;=\; -\frac{f^{(n+1)}(t)}{n!}(x - t)^n, \qquad G'(t) \;=\; -(n+1)(x-t)^n.

ロールの定理の適用。 補助関数

ϕ(t)    F(t)F(x0)G(x0)G(t)\phi(t) \;\coloneqq\; F(t) - \frac{F(x_0)}{G(x_0)}\,G(t)

を定める。ϕ(x0)=0\phi(x_0) = 0 かつ ϕ(x)=F(x)F(x0)G(x0)0=0\phi(x) = F(x) - \tfrac{F(x_0)}{G(x_0)} \cdot 0 = 0 だ。ロールの定理より、x0x_0xx の間に ϕ(ξ)=0\phi'(\xi) = 0 を満たす ξ\xi が存在する:

F(ξ)F(x0)G(x0)G(ξ)=0.F'(\xi) - \frac{F(x_0)}{G(x_0)}\,G'(\xi) = 0.

ξx\xi \neq x より (xξ)n0(x - \xi)^n \neq 0 なので、FF'GG' を代入すると:

f(n+1)(ξ)n!(xξ)n+F(x0)G(x0)(n+1)(xξ)n=0        F(x0)G(x0)=f(n+1)(ξ)(n+1)!.-\frac{f^{(n+1)}(\xi)}{n!}(x-\xi)^n + \frac{F(x_0)}{G(x_0)}(n+1)(x-\xi)^n = 0 \;\;\Longrightarrow\;\; \frac{F(x_0)}{G(x_0)} = \frac{f^{(n+1)}(\xi)}{(n+1)!}.

両辺に G(x0)=(xx0)n+1G(x_0) = (x-x_0)^{n+1} を掛けると F(x0)=Rn(x)F(x_0) = R_n(x)(2)(2) の形で得られる。\square

標準的なマクローリン展開

x0=0x_0 = 0 の場合を**マクローリン展開(Maclaurin expansion)**という。nn \to \infty として Rn0R_n \to 0 を確認すると、全ての xRx \in \mathbb{R} に対して以下が成り立つ:

ex  =  k=0xkk!  =  1+x+x22!+x33!+e^x \;=\; \sum_{k=0}^{\infty} \frac{x^k}{k!} \;=\; 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots sinx  =  k=0(1)k(2k+1)!x2k+1  =  xx36+x5120\sin x \;=\; \sum_{k=0}^{\infty} \frac{(-1)^k}{(2k+1)!} x^{2k+1} \;=\; x - \frac{x^3}{6} + \frac{x^5}{120} - \cdots cosx  =  k=0(1)k(2k)!x2k  =  1x22+x424\cos x \;=\; \sum_{k=0}^{\infty} \frac{(-1)^k}{(2k)!} x^{2k} \;=\; 1 - \frac{x^2}{2} + \frac{x^4}{24} - \cdots ln(1+x)  =  k=1(1)k1kxk  =  xx22+x33(1<x1)\ln(1+x) \;=\; \sum_{k=1}^{\infty} \frac{(-1)^{k-1}}{k} x^k \;=\; x - \frac{x^2}{2} + \frac{x^3}{3} - \cdots \quad (-1 < x \leq 1) (1+x)α  =  k=0(αk)xk  =  1+αx+α(α1)2!x2+(x<1)(1+x)^\alpha \;=\; \sum_{k=0}^{\infty} \binom{\alpha}{k} x^k \;=\; 1 + \alpha x + \frac{\alpha(\alpha-1)}{2!}x^2 + \cdots \quad (|x| < 1)

ここで一般化二項係数は (αk)α(α1)(αk+1)k!\binom{\alpha}{k} \coloneqq \dfrac{\alpha(\alpha-1)\cdots(\alpha-k+1)}{k!} だ。

剰余項による誤差評価

ラグランジュ形式 (2)(2) は具体的な誤差上界を与える。x0x_0xx の間の全ての tt に対して f(n+1)(t)M|f^{(n+1)}(t)| \leq M ならば、

Rn(x)    M(n+1)!xx0n+1.|R_n(x)| \;\leq\; \frac{M}{(n+1)!}\,|x - x_0|^{n+1}.

例。 x0=0x_0 = 0 における P3P_3 を使って e0.1e^{0.1} を近似する。

P3(0.1)  =  1+0.1+0.012+0.0016    1.105166.P_3(0.1) \;=\; 1 + 0.1 + \frac{0.01}{2} + \frac{0.001}{6} \;\approx\; 1.10516\overline{6}.

剰余項は R3(0.1)e0.14!(0.1)4<3241041.25×105|R_3(0.1)| \leq \dfrac{e^{0.1}}{4!}(0.1)^4 < \dfrac{3}{24} \cdot 10^{-4} \approx 1.25 \times 10^{-5} を満たす。この近似は小数点以下5桁まで正確だ。

極限の計算

テイラーの公式は不定形の極限を機械的に計算できるようにする。limx0sinxxx3\lim_{x \to 0} \dfrac{\sin x - x}{x^3} を求めるには、sinx=xx36+O(x5)\sin x = x - \dfrac{x^3}{6} + O(x^5) と展開すればよい:

sinxxx3  =  x3/6+O(x5)x3    16.\frac{\sin x - x}{x^3} \;=\; \frac{-x^3/6 + O(x^5)}{x^3} \;\to\; -\frac{1}{6}.

まとめ

  • x0x_0 におけるテイラー多項式 PnP_nx0x_0 において ff とその nn 次までの全導関数の値が一致する。
  • テイラーの公式f(x)=Pn(x)+Rn(x)f(x) = P_n(x) + R_n(x) であり、ラグランジュ剰余項Rn(x)=f(n+1)(ξ)(n+1)!(xx0)n+1R_n(x) = \dfrac{f^{(n+1)}(\xi)}{(n+1)!}(x-x_0)^{n+1}ξ\xix0x_0xx の間のある点)だ。
  • 証明の方針x0x_0 での値が剰余項に等しい補助関数を定義し、コーシーの平均値定理を適用して最高次の導関数項以外を消去する。
  • ラグランジュ剰余項は f(n+1)M|f^{(n+1)}| \leq M のとき上界 Rn(x)M(n+1)!xx0n+1|R_n(x)| \leq \dfrac{M}{(n+1)!}|x-x_0|^{n+1} を与える。
  • exe^xsinx\sin xcosx\cos xln(1+x)\ln(1+x)(1+x)α(1+x)^\alpha の標準的な冪級数は nn \to \infty とすることで得られる。