平均値の定理は x0 と x の間のある未知の c に対して f(x)=f(x0)+f′(c)(x−x0) と述べる。これは精度の粗い剰余項を持つ線形近似だ。テイラーの公式(Taylor’s Formula)はこの考えを全次数に拡張する:f を n 次多項式で近似し、残りの誤差を明示的に表す。
テイラー多項式
定義。 f の x0 における**n 次テイラー多項式(Taylor polynomial of degree n)**を
Pn(x):=k=0∑nk!f(k)(x0)(x−x0)k=f(x0)+f′(x0)(x−x0)+2!f′′(x0)(x−x0)2+⋯+n!f(n)(x0)(x−x0)n
と定める。
係数 k!f(k)(x0) は Pn(k)(x0)=f(k)(x0)(k=0,1,…,n)を成り立たせる唯一の値だ:この多項式は x0 において f とその n 次までの全導関数の値が一致する。
ラグランジュ剰余項によるテイラーの公式
定理(テイラーの公式)。 f が x0 と x を含む開区間上で n+1 回微分可能とする。このとき
f(x)=Pn(x)+Rn(x),(1)
が成り立つ。ここで**ラグランジュ剰余項(Lagrange remainder)**は
Rn(x)=(n+1)!f(n+1)(ξ)(x−x0)n+1(2)
であり、ξ は x0 と x の間のある点だ。
証明
x=x0 を固定し、J を x0 と x を端点とする閉区間とする。次の関数を定義する:
F(t):=f(x)−k=0∑nk!f(k)(t)(x−t)k,G(t):=(x−t)n+1.
F、G はともに J 上で連続、J の内部で微分可能だ。次のことに注目する:
- F(x)=0 かつ G(x)=0。
- F(x0)=f(x)−Pn(x) かつ G(x0)=(x−x0)n+1=0。
F の導関数の望遠鏡和。 積の微分法則を各項 k!f(k)(t)(x−t)k に適用すると:
dtd[k!f(k)(t)(x−t)k]=k!f(k+1)(t)(x−t)k−(k−1)!f(k)(t)(x−t)k−1.
k=0 から n まで和をとると、望遠鏡のように各項 (k−1)!f(k)(t)(x−t)k−1 が k↦k−1 に対応する第一項と消え合い、k=n の第一項だけが残る:
F′(t)=−n!f(n+1)(t)(x−t)n,G′(t)=−(n+1)(x−t)n.
ロールの定理の適用。 補助関数
ϕ(t):=F(t)−G(x0)F(x0)G(t)
を定める。ϕ(x0)=0 かつ ϕ(x)=F(x)−G(x0)F(x0)⋅0=0 だ。ロールの定理より、x0 と x の間に ϕ′(ξ)=0 を満たす ξ が存在する:
F′(ξ)−G(x0)F(x0)G′(ξ)=0.
ξ=x より (x−ξ)n=0 なので、F′、G′ を代入すると:
−n!f(n+1)(ξ)(x−ξ)n+G(x0)F(x0)(n+1)(x−ξ)n=0⟹G(x0)F(x0)=(n+1)!f(n+1)(ξ).
両辺に G(x0)=(x−x0)n+1 を掛けると F(x0)=Rn(x) が (2) の形で得られる。□
標準的なマクローリン展開
x0=0 の場合を**マクローリン展開(Maclaurin expansion)**という。n→∞ として Rn→0 を確認すると、全ての x∈R に対して以下が成り立つ:
ex=k=0∑∞k!xk=1+x+2!x2+3!x3+⋯
sinx=k=0∑∞(2k+1)!(−1)kx2k+1=x−6x3+120x5−⋯
cosx=k=0∑∞(2k)!(−1)kx2k=1−2x2+24x4−⋯
ln(1+x)=k=1∑∞k(−1)k−1xk=x−2x2+3x3−⋯(−1<x≤1)
(1+x)α=k=0∑∞(kα)xk=1+αx+2!α(α−1)x2+⋯(∣x∣<1)
ここで一般化二項係数は (kα):=k!α(α−1)⋯(α−k+1) だ。
剰余項による誤差評価
ラグランジュ形式 (2) は具体的な誤差上界を与える。x0 と x の間の全ての t に対して ∣f(n+1)(t)∣≤M ならば、
∣Rn(x)∣≤(n+1)!M∣x−x0∣n+1.
例。 x0=0 における P3 を使って e0.1 を近似する。
P3(0.1)=1+0.1+20.01+60.001≈1.105166.
剰余項は ∣R3(0.1)∣≤4!e0.1(0.1)4<243⋅10−4≈1.25×10−5 を満たす。この近似は小数点以下5桁まで正確だ。
極限の計算
テイラーの公式は不定形の極限を機械的に計算できるようにする。limx→0x3sinx−x を求めるには、sinx=x−6x3+O(x5) と展開すればよい:
x3sinx−x=x3−x3/6+O(x5)→−61.
まとめ
- x0 におけるテイラー多項式 Pn は x0 において f とその n 次までの全導関数の値が一致する。
- テイラーの公式:f(x)=Pn(x)+Rn(x) であり、ラグランジュ剰余項は Rn(x)=(n+1)!f(n+1)(ξ)(x−x0)n+1(ξ は x0 と x の間のある点)だ。
- 証明の方針:x0 での値が剰余項に等しい補助関数を定義し、コーシーの平均値定理を適用して最高次の導関数項以外を消去する。
- ラグランジュ剰余項は ∣f(n+1)∣≤M のとき上界 ∣Rn(x)∣≤(n+1)!M∣x−x0∣n+1 を与える。
- ex、sinx、cosx、ln(1+x)、(1+x)α の標準的な冪級数は n→∞ とすることで得られる。