§D1.16テイラーの定理と剰余項の評価

最終更新

高校では「関数を多項式で近似する」テイラー展開を形式的に使いました。本記事では、テイラーの定理を証明し、多項式による近似の誤差を正確に評価します。テイラーの定理は、平均値の定理の一般化です。

1 コーシーの平均値定理

テイラーの定理の証明に用いるために、平均値の定理を2つの関数に拡張した定理を用意します。

定理 1.1 (コーシーの平均値定理).F,GF, Gが[a,b][a,b]で連続、(a,b)(a,b)で微分可能とする。このときあるc∈(a,b)c \in (a,b)で(F(b)−F(a))G′(c)=(G(b)−G(a))F′(c)\big(F(b)-F(a)\big)G'(c) = \big(G(b)-G(a)\big)F'(c)。

証明.h(x)=(F(b)−F(a))G(x)−(G(b)−G(a))F(x)h(x) = \big(F(b)-F(a)\big)G(x) - \big(G(b)-G(a)\big)F(x)とおくとhhは[a,b][a,b]で連続・(a,b)(a,b)で微分可能で、h(a)=F(b)G(a)−G(b)F(a)=h(b)h(a) = F(b)G(a)-G(b)F(a) = h(b)。§D1.14 定理 2.2よりあるccでh′(c)=0h'(c)=0、これが主張の等式。▨

2 テイラーの定理

定理 2.1 (テイラーの定理).ffがaaとxxを含む閉区間でf(n−1)f^{(n-1)}まで連続、その開区間でf(n)f^{(n)}が存在するとする。このときaaとxxの間のある点ccで

Rn(x)=f(n)(c)n!(x−a)n.R_n(x) = \frac{f^{(n)}(c)}{n!}(x-a)^n.

証明.xxを固定し、ttの関数

F(t)=f(x)−∑k=0n−1f(k)(t)k!(x−t)k,G(t)=(x−t)nF(t) = f(x) - \sum_{k=0}^{n-1} \frac{f^{(k)}(t)}{k!}(x-t)^k, \qquad G(t) = (x-t)^n

を考える。F(x)=0F(x) = 0、F(a)=Rn(x)F(a) = R_n(x)、G(x)=0G(x)=0、G(a)=(x−a)nG(a)=(x-a)^n。FFを微分すると、和が望遠鏡状に打ち消し合って

F′(t)=−f(n)(t)(n−1)!(x−t)n−1F'(t) = -\frac{f^{(n)}(t)}{(n-1)!}(x-t)^{n-1}

だけが残る(ddt[f(k)(t)k!(x−t)k]=f(k+1)(t)k!(x−t)k−f(k)(t)(k−1)!(x−t)k−1\frac{d}{dt}\big[\frac{f^{(k)}(t)}{k!}(x-t)^k\big] = \frac{f^{(k+1)}(t)}{k!}(x-t)^k - \frac{f^{(k)}(t)}{(k-1)!}(x-t)^{k-1}が隣接項で相殺)。またG′(t)=−n(x−t)n−1G'(t) = -n(x-t)^{n-1}。コーシーの平均値定理(定理 1.1)を[a,x][a,x]上のF,GF, Gに適用すると、あるccで

=−f(n)(c)(n−1)!(x−c)n−1−n(x−c)n−1=f(n)(c)n!.= \frac{-\frac{f^{(n)}(c)}{(n-1)!}(x-c)^{n-1}}{-n(x-c)^{n-1}} = \frac{f^{(n)}(c)}{n!}.

左辺はRn(x)−0(x−a)n−0=Rn(x)(x−a)n\dfrac{R_n(x) - 0}{(x-a)^n - 0} = \dfrac{R_n(x)}{(x-a)^n}。よってRn(x)=f(n)(c)n!(x−a)nR_n(x) = \dfrac{f^{(n)}(c)}{n!}(x-a)^n。▨

n=1n=1とすると、テイラーの定理は平均値の定理そのものになります。テイラーの定理では、「1点の傾き」ではなく「高次の導関数まで」を使って関数ffを近似し、その近似の誤差を最高次の導関数で押さえます。

3 剰余項の評価とテイラー級数の収束

剰余項の形が分かると、近似の誤差を具体的に評価できます。∣f(n)(c)∣|f^{(n)}(c)|を区間上でMMで抑えれば、∣Rn(x)∣≤Mn!∣x−a∣n|R_n(x)| \le \dfrac{M}{n!}|x-a|^nが成り立ちます。この不等式は、次の二つに用いることができます。

  • テイラー級数が本当に関数に収束するか(Rn→0R_n \to 0か)の判定
  • sin⁡x,ex\sin x, e^xなどの値の、必要な精度での計算(数値計算の基礎)

例 3.1.f(x)=exf(x) = e^xはf(n)(t)=etf^{(n)}(t) = e^tゆえ、a=0a = 0とすると00とxxの間のccで∣Rn(x)∣=ecn!∣x∣n≤e∣x∣n!∣x∣n|R_n(x)| = \dfrac{e^c}{n!}|x|^n \le \dfrac{e^{|x|}}{n!}|x|^n。固定したxxに対し∣x∣nn!→0\dfrac{|x|^n}{n!} \to 0(階乗が指数より速く増える)だからRn(x)→0R_n(x) \to 0、すなわちex=∑k=0∞xkk!e^x = \sum_{k=0}^\infty \frac{x^k}{k!}がすべてのxxで成り立つ。とくにx=1x = 1でe=∑k=0∞1k!e = \sum_{k=0}^\infty \frac{1}{k!}、nn項で打ち切った誤差はe1n!<3n!\dfrac{e^1}{n!} < \dfrac{3}{n!}以下である。

4 滑らかな関数とテイラー級数

補題 4.1. 関数θ:R→R\theta:\R\to\Rを

θ(t)={e−1/t,t>0,0,t≤0\theta(t)= \begin{cases} e^{-1/t},&t>0,\\ 0,&t\leq0 \end{cases}

で定める。このときθ\thetaはR\R上で何回でも微分可能であり、各階の導関数は連続である。すなわち、θ\thetaはC∞C^\infty級(滑らか)である。さらに、すべての非負整数nnについてθ(n)(0)=0\theta^{(n)}(0)=0である。ただしθ(0)=θ\theta^{(0)}=\thetaとする。

証明.NNを非負整数、s>0s>0とする。指数関数の各階の導関数は指数関数自身であるから、定理 2.1をa=0a=0、x=sx=s、n=N+2n=N+2として適用すると、あるc∈(0,s)c\in(0,s)に対して

es=∑k=0N+1skk!+ec(N+2)!sN+2≥sN+1(N+1)!e^s=\sum_{k=0}^{N+1}\frac{s^k}{k!} +\frac{e^c}{(N+2)!}s^{N+2} \geq\frac{s^{N+1}}{(N+1)!}

である。したがって

0≤sNe−s≤(N+1)!s⟶0(s→∞).0\leq s^Ne^{-s}\leq\frac{(N+1)!}{s}\longrightarrow0 \qquad(s\to\infty).

任意の実係数多項式PPは非負整数冪の有限和であるから、P(s)e−s→0P(s)e^{-s}\to0も成り立つ。

多項式PnP_nをP0(s)=1P_0(s)=1と

Pn+1(s)=s2(Pn(s)−Pn′(s))(n≥0)P_{n+1}(s)=s^2\bigl(P_n(s)-P_n'(s)\bigr) \qquad(n\geq0)

によって定める。t>0t>0では、積と合成の微分により

ddt(Pn(1/t)e−1/t)=t−2(Pn(1/t)−Pn′(1/t))e−1/t=Pn+1(1/t)e−1/t\frac{d}{dt}\bigl(P_n(1/t)e^{-1/t}\bigr) =t^{-2}\bigl(P_n(1/t)-P_n'(1/t)\bigr)e^{-1/t} =P_{n+1}(1/t)e^{-1/t}

である。したがってθ\thetaは(0,∞)(0,\infty)上で何回でも微分可能であり、θ(n)(t)=Pn(1/t)e−1/t\theta^{(n)}(t)=P_n(1/t)e^{-1/t}である。一方、t<0t<0では各階の導関数は00である。

N=0N=0の減衰評価からθ(t)→0\theta(t)\to0がt↓0t\downarrow0で成り立つため、θ\thetaは原点でも連続である。n≥0n\geq0とし、θ(n)\theta^{(n)}がR\R上で連続で、θ(n)(0)=0\theta^{(n)}(0)=0であると仮定する。h>0h>0では

θ(n)(h)−θ(n)(0)h=h−1Pn(1/h)e−1/h⟶0(h↓0)\frac{\theta^{(n)}(h)-\theta^{(n)}(0)}{h} =h^{-1}P_n(1/h)e^{-1/h}\longrightarrow0 \qquad(h\downarrow0)

である。ここでは多項式sPn(s)sP_n(s)に減衰評価を適用した。一方、h<0h<0では差商は00である。したがって両側からの差商の極限は00であり、θ(n+1)(0)\theta^{(n+1)}(0)が存在して00に等しい。また、多項式Pn+1P_{n+1}に対する減衰評価からθ(n+1)(h)→0\theta^{(n+1)}(h)\to0がh↓0h\downarrow0で成り立つ。一方、h<0h<0ではθ(n+1)(h)=0\theta^{(n+1)}(h)=0である。したがってθ(n+1)\theta^{(n+1)}は原点でも連続である。帰納法により、すべての非負整数nnについて、θ(n)\theta^{(n)}はR\R上で連続であり、θ(n)(0)=0\theta^{(n)}(0)=0である。▨

例 4.2.補題 4.1の関数θ\thetaと多項式x↦x2x\mapsto x^2の合成をf:R→Rf:\R\to\Rとする。すなわち

f(x)=θ(x2)={e−1/x2,x≠0,0,x=0f(x)=\theta(x^2)= \begin{cases} e^{-1/x^2},&x\ne0,\\ 0,&x=0 \end{cases}

である。このとき、すべての非負整数nnに対して、ある多項式qn,0,…,qn,nq_{n,0},\ldots,q_{n,n}が存在して

f(n)(x)=∑k=0nqn,k(x)θ(k)(x2)f^{(n)}(x)=\sum_{k=0}^n q_{n,k}(x)\theta^{(k)}(x^2)

となる。実際、n=0n=0ではq0,0=1q_{0,0}=1とすればよく、各項の導関数は

qn,k′(x)θ(k)(x2)+2xqn,k(x)θ(k+1)(x2)q_{n,k}'(x)\theta^{(k)}(x^2) +2xq_{n,k}(x)\theta^{(k+1)}(x^2)

であるから、帰納的に同じ形の有限和が得られる。補題 4.1により各項は連続で、x=0x=0では00となるため、ffはC∞C^\infty級であり、すべての非負整数nnについてf(n)(0)=0f^{(n)}(0)=0である。したがって、00を中心とするffのテイラー級数は恒等的に00であり、x≠0x\ne0におけるf(x)>0f(x)>0と一致しない。固定したx≠0x\ne0に対して、剰余項はすべてのn≥1n\geq1でRn(x)=f(x)>0R_n(x)=f(x)>0であるため、00へ収束しない。このように、C∞C^\infty級であることから、テイラー級数が関数と局所的に一致するという解析性は従わない。なお、複素解析では正則関数は解析的である。

前提記事