Processing math: 100%
スポンサーリンク

「ちゃんとしたテイラー展開の話」「代表例を計算してみよう!」【解析学の基礎シリーズ】1変数実数値関数の微分編 その11

微分法

本記事の内容

本記事はテイラー展開のイメージを説明し、代表的な例を実際に計算してみる記事です。

本記事を読むにあたり、テイラーの定理とベキ級数の収束半径を知っている必要があるため、以下の記事も合わせてご覧ください。

↓テイラー展開の記事

↓ベキ級数の収束半径の記事

テイラーの定理とテイラー展開のチャラい復習

テイラーの定理とテイラー展開をチャラく復習します。

テイラーの定理

テイラーの定理はなんだったっけ?というと、以下でした。

テイラーの定理 nNIRの区間、f:IRn階微分可能な関数、aIxIとするとき、 f(x)=f(a)+f(a)1!(xa)+f(a)2!(xa)2++f(n1)(a)(n1)!(xa)n1+Rn と書いてRnを定めれば、 Rn=f(n)(c)n!(xa)n を満たすようなcaxの間に、すなわちa<xならばc(a,x)に、x<aならばc(x,a)に、a=xならばc=aとして存在する。

この定理の証明は【解析学の基礎シリーズ】1変数実数値関数の微分編 その10を御覧ください。

このRnn次剰余項とよび、この形で表される剰余項をLagrange(ラグランジュ)剰余項と呼ぶのでした。

テイラー展開

テイラーの定理のRnに注目してみると、
Rn=f(k)(c)n!(xa)n
ですので、Rnxに依存します。
言い換えれば、変数をxとする関数です。
故に、以降はRnRk(x)と書き直します。

このRnがある点aの十分近くでlimnRn(x)=0が成り立てば、
f(x)=n=0f(n)(a)n!(xa)n
が成り立ちます。
これをfaを中心とするテイラー展開といったのでした。

しっかりとしたテイラー展開の話をしましょう。

先程、チャラくテイラー展開について復習しましたが、チャラさを置いておいても実は厳密ではありません。

「十分近くって何かネ?」やらと色々疑問が出てきます。
確かに、「何となく分かる」かなと思いますが、「なんとなく」というのは数学ではあまり意味を持ちません。
そこで厳密な話をしようと思います。

余談(読み飛ばしてOK) テイラー展開をググってみると、「こうやって計算するんだ!」「これは公式だ!」というように、ちょろっと計算してその結果を列挙している記事は多々見つかるのですが、収束半径まで踏み込んだ記事はあまりなさそうに見受けられます(筆者調べ)。
テイラー展開はある条件下の関数をベキ級数での表現なので、本来は収束半径やらを考察する必要があります。
本記事はそこまで踏み込んだ話をしようと思います

厳密なテイラー展開

さて、テイラー展開を厳密な定理として書くと、以下です。

定理1.(テイラー展開とその条件)
  1. テイラー展開
  2. 関数faを含む区間IにおいてC級(何度でも微分可能)で、Iの各点xlimkRk(x)=0 を満たすとする。このとき、fI上で f(x)=n=0f(n)(a)n!(xa)n の形の級数で表される。この形の級数をfaを中心とするテイラー展開(テイラー series expansion)という。
  3. テイラー展開可能な条件
  4. 定数C0M0が存在して、任意のnNおよびxIに対して |f(n)(x)|CMn を満たすとき、faを中心としてIでテイラー展開できる。

1.テイラー展開については、すでに前回証明しているので、【解析学の基礎シリーズ】1変数実数値関数の微分編 その10を御覧ください。

本記事の主たる目標の1つが2.テイラー展開可能な条件です。

これを示すためには、ratio test(証明の道中で明示します)が必要になります。

定理1.の証明

繰り返しになりますが、1.についてはすでに証明しているので、【解析学の基礎シリーズ】1変数実数値関数の微分編 その10を御覧ください。

2.を証明します。
Rn(x)=f(n)(c)n!(xa)n
であり、定数C0M0が存在して、任意のnNおよびxIに対して
|f(n)(x)|CMn
を満たすとします。
xIは任意なので、テイラーの定理で存在が保証されているIの内点cIに対して、
|R(x)|=|f(n)(c)|n!|xa|nCMnn!|xa|n
です。
ここで、ratio testを使います。

定理2.(ratio test) 正項級数n=0anに対して、 limnan+1an=lR が存在して、l<1ならばn=0anは収束し、l>1ならばn=0anは発散する。 また、極限が存在しない場合、すなわちlimnan+1an=の場合もn=0anは発散する。

定理2.(ratio test)の証明は【解析学の基礎シリーズ】級数とベキ級数編 その1を御覧ください。

an=CMnn!|xa|n
として、定理2.(ratio test)を使うと、
limnan+1an=limnCMn+1(n+1)!|xa|n+1CMnn!|xa|n=limnMn|xa|=M|xa|limn1n+1=M|xa|0=0<1
です。
ただし、limn1n+1=0を使いました。
この事実の証明は【解析学の基礎シリーズ】実数の連続性編 その8を御覧ください。

従って、an+1an=0<1により、定理2.(ratio test)からn=0anは収束します。

さらに、次の事実を使います。

補題3. 級数n=0anが収束すれば、limnan=0である。

この証明は【解析学の基礎シリーズ】級数とべき級数編 その2の補題3.を御覧ください。

従って、
limnan=limnCMnn!|xa|n=0
となるので、1.の条件を満たすから、テイラー展開可能です。

定理1.の証明終わり

では、代表的な関数のテイラー展開とマクローリン展開の例を計算してみましょう!

テイラー展開とマクローリン展開の代表例

この節で紹介するテイラー展開とマクローリン展開の代表例は大きく分けて

  • 指数関数
  • 対数関数
  • 三角関数(tanxは飛ばします。理由は余談を参照してください。)

です。
テイラー展開を再掲しておきます。

定理1.(テイラー展開とその条件)
  1. テイラー展開
  2. 関数faを含む区間IにおいてC級(何度でも微分可能)で、Iの各点xlimkRk(x)=0 を満たすとする。このとき、fI上で f(x)=n=0f(n)(a)n!(xa)n の形の級数で表される。この形の級数をfaを中心とするテイラー展開(テイラー series expansion)という。
  3. テイラー展開可能な条件
  4. 定数C0M0が存在して、任意のnNおよびxIに対して |f(n)(x)|CMn を満たすとき、faを中心としてIでテイラー展開できる。

要は、上記の3つの種類の関数について、
f(x)=n=0f(n)(a)n!(xa)n
を収束半径とともに具体的に計算しよう、という話です。

さらに、収束半径は以下でした。

定理2. ベキ級数n=0an(xa)nに対して、次の1.および2.を満たす0R(Rの場合も有り得る)が存在する。
  1. |xa|<Rならば、n=0an(xa)nは絶対収束する。
  2. |xa|>Rならば、n=0an(xa)nは収束しない。
定理7. 極限 limn|anan+1|=R[0,] が存在すれば、Rn=0an(xa)nの収束半径である。

これらの証明は【解析学の基礎シリーズ】級数とベキ級数編 その2を御覧ください。

指数関数のテイラー展開とマクローリン展開

1. f(x)=kxx=aまわりのテイラー展開

f(x)=kxC級(何度でも微分可能)です。
そこで、f(n)(x)が計算してみます。

  • f(0)(x)=f(x)=kx
  • f(1)(x)=kx(logk)
  • f(2)(x)=kx(logk)2
  • f(3)(x)=kx(logk)3

というようになります。
従って、f(n)(x)=kx(logk)nです。
f(n)(x)は求まりましたが、テイラー展開可能なのか?ということを吟味します。

kRですので、k|k|<|k+1|です。
そして、logk<log(k+1)で、かつlog(k+1)>0です。
故に任意のt>0に対して、|xt|<1であれば、C=|kt|とし、M=log(k+1)とすれば、定理1.の2.により
|f(n)(x)|=kx(logk)n|kt|(log(k+1))n
となります。
また、tは任意の実数ですので、テイラー展開可能です。

故に、x=aまわりのテイラー展開は、
f(x)=n=0f(n)(a)n!(xa)n=n=0ka(logk)nn!(xa)n
です。

特に、a=0の周りのテイラー展開、すなわちマクローリン展開は、
f(x)=n=0f(n)(0)n!xn=n=0(logk)nn!xn

次に収束半径です。
an=ka(logk)nn!とします。
このとき、
limn|anan+1|=limn|ka(logk)nn!ka(logk)n+1(n+1)!|=limn|1n!logk(n+1)!|=limn|n+1logk|=
ですので、R=、すなわち収束半径はです。
すなわち、任意のxRに対して、級数
kx=n=0ka(logk)nn!(xa)n
の右辺は収束するため、xにどんな実数を代入してもその値が実数値で求まる、ということです。

2. f(x)=ex(eはNapier数)のテイラー展開とマクローリン展開(k=eの場合)

上記のテイラー展開において、k=eとすればOKです。
すなわち、
ex=n=0ean!(xa)n
です。
特にa=0のときのテイラー展開、すなわちマクローリン展開は
f(x)=n=0f(n)(0)n!xn=n=01n!xn
です。
収束半径は、an=ean!としたとき、
limn|anan+1|=limn|n+1logk|=limn|n+1|=
となり、収束半径R=です。
すなわち、任意のxRに対して、級数
ex=n=0ean!(xa)n
の右辺は収束するため、xにどんな実数を代入してもその値が実数値で求まる、ということです。

対数関数のテイラー展開とマクローリン展開

f(x)=log(1+x)C級(何度でも微分可能)です。
そこで、f(n)(x)が計算してみます。
どうしてlogxではなくlog(1+x)なのか、というのは後述します。

  • f(0)(x)=f(x)=log(1+x)
  • f(1)(x)=11+x
  • f(2)(x)=1(1+x)2
  • f(3)(x)=2(1+x)3
  • f(4)(x)=6(1+x)4

というようになります。
従って、f(n)(x)=(1)n1(n1)!(1+x)nです。
f(n)(x)は求まりましたが、テイラー展開可能なのか?ということを吟味します。

Rn(x)=f(n)(c)n!(xa)n=(1)n1(n1)!(1+c)n1n!(xa)n
により、
|Rn(x)|=1n|xa1+c|n
です。
これは、|xa|<1であれば、nで収束します。
実際、|xa|<1により、1<xa<1だからa1<x<a+1です。
caxの間に存在するため、a<c<xかまたはx<c<aです。
仮にa<c<xであれば、a<c<x<a+1<c+1ですので、xa1+c<1です。
また、x<c<aであれば、xa<x<c<c+1により、xa1+c<1です。
従って収束します。

以上のことから、log(1+x)はテイラー展開可能です。

故に、x=aまわりのテイラー展開は、
f(x)=n=0f(n)(a)n!(xa)n=n=0(1)n11n(1+a)n(xa)n
です。

特に、a=0の周りのテイラー展開、すなわちマクローリン展開は、
f(x)=n=0f(n)(0)n!xn=n=0(1)n1nxn
です。
ちなみに、ここでlogxx>0でなければなりません。
従って、logxx=0まわりでテイラー展開できません。
つまり、マクローリン展開ができません。
従って、x=0でテイラー展開ができるようにするためにlog(1+x)を考えるのです。

次に収束半径です。
an=(1)n11n(1+a)nとします。
このとき、
limn|anan+1|=limn|1n(1+a)n1(n+1)(1+a)n|=limn|n+1n(1+a)|=limn|(1+1n)(1+a)|=1+a
ですので、R=1+a、すなわち収束半径は1+aです。
すなわち、|xa|<1+aを満たすxRでないと
log(1+x)=n=0(1)n11n(1+a)n(xa)n
の右辺は収束しません。

例えば、x=1だったとすると、
n=0(1)n11n(1+a)n(1a)n=n=0(1)n11n(1+a)n(1)n(1+a)n=n=01n=
となって、発散してしまいます。

三角関数のテイラー展開とマクローリン展開

まずは、sinxについてです。

sinxのテイラー展開とマクローリン展開

f(x)=sinxC級(何度でも微分可能)です。
そこで、f(n)(x)が計算してみます。

  • f(0)(x)=f(x)=sinx
  • f(1)(x)=cosx
  • f(2)(x)=sinx
  • f(3)(x)=cosx
  • f(4)(x)=sinx
  • f(5)(x)=cosx

というようになります。
従って、f(n)(x)は、mN{0}として、

  • f(4m)(x)=sinx
  • f(4m+1)(x)=cosx
  • f(4m+2)(x)=sinx
  • f(4m+3)(x)=cosx

です。
すなわち、f(n)(x)=sin(x+nπ2)です。

f(n)(x)は求まりましたが、テイラー展開可能なのか?ということを吟味します。

Rn(x)=f(n)(c)n!(xa)n=1n!sin(c+nπ2)(xa)n
です。
任意のxRに対して、|sinx|1です。
従って、
|Rn(x)||(xa)nn!|
です。
ここで、次の事実を使います。

命題8. 任意の正の実数αに対して、limnαnn!=0が成り立つ。
命題8.の証明

なんてことありません。
アルキメデスの原理から、α<Nという自然数Nが存在します。
アルキメデスの原理は何だったか、というと、以下でした。

アルキメデスの原理 (aR>0)(bR>0)(N0N) s.t. b<N0a ただし、R>0={xRx>0}、すなわち正の実数の集合である。

アルキメデスの原理の証明は【解析学の基礎シリーズ】実数の連続性編 その12を御覧ください。

このとき、n>Nであれば、
αnn!=αnαn1αNαN1α2α
と表すことができます。
後半のαN1α2αは定数です。
従って、
limnαn!=limnαnαn1αNαN1α2α=αN1α2αlimnαnαn1αN=αN1α2αlimnαnlimnαn1limnαN=0
従って、成り立ちます。

命題8.の証明終わり

命題8.により
0limn|Rn(x)|limn|(xa)nn!|=0
となって、はさみうちの原理から、limnRn(x)=0です。

ここで、はさみうちの原理とは、

はさみうちの原理
{an}nN{bn}nN{cn}nNを実数の数列とする。 このとき、任意のnNに対してancnbnが成り立ち、かつある実数Aに対してlimnan=limnbn=Aが成り立つならば、limncn=Aである。

でした。
はさみうちの原理の証明は【解析学の基礎シリーズ】実数の連続性編 その15を御覧ください。

以上のことから、sinxはテイラー展開可能です。

故に、x=aまわりのテイラー展開は、
f(x)=n=0f(n)(a)n!(xa)n=n=01n!sin(a+nπ2)(xa)n
です。

特に、a=0の周りのテイラー展開、すなわちマクローリン展開は、
f(x)=n=0f(n)(0)n!xn=n=01n!sin(nπ2)xn
です。
ここで、sinxの性質から、

  • f(0)(0)0!=10!sin0=0
  • f(1)(0)1!=11!cos0=1
  • f(2)(0)2!=12!sin0=0
  • f(3)(0)3!=13!cos0=13!
  • f(4)(0)4!=14!sin0=0
  • f(5)(0)5!=12!cos0=15!

です。
従って、
f(x)=n=0f(n)(0)n!xn=n=0(1)n+1(2n1)!xn
です。

次に収束半径です。
an=1n!sin(a+nπ2)とします。
このとき、
limn|anan+1|=limn|sin(a+nπ2)n!sin(a+(n+1)π2)(n+1)!|=limn|sin(a+nπ2)sin(a+(n+1)π2)n|
で、かつ、任意のxNに対して、|sinx|1ですので、R=、すなわち収束半径はです。
すなわち、任意のxRに対して、級数
sinx=n=0f(n)(a)n!(xa)n=n=01n!sin(a+nπ2)(xa)n
の右辺は収束するため、xにどんな実数を代入してもその値が実数値で求まる、ということです。

cosxのテイラー展開とマクローリン展開

sinxとほぼ同じです。

f(x)=cosxC級(何度でも微分可能)です。
そこで、f(n)(x)が計算してみます。

  • f(0)(x)=f(x)=cosx
  • f(1)(x)=sinx
  • f(2)(x)=cosx
  • f(3)(x)=sinx
  • f(4)(x)=cosx
  • f(5)(x)=sinx

というようになります。
従って、f(n)(x)は、mN{0}として、

  • f(4m)(x)=cosx
  • f(4m+1)(x)=sinx
  • f(4m+2)(x)=cosx
  • f(4m+3)(x)=sinx

です。
すなわち、f(n)(x)=cos(x+nπ2)です。

f(n)(x)は求まりましたが、テイラー展開可能なのか?ということを吟味します。

Rn(x)=f(n)(c)n!(xa)n=1n!cos(c+nπ2)(xa)n
です。
任意のxRに対して、|cosx|1です。
従って、
|Rn(x)||(xa)nn!|
です。

sinxのときと同様にlimnRn(x)=0です。
従ってsinxはテイラー展開可能です。

故に、x=aまわりのテイラー展開は、
f(x)=n=0f(n)(a)n!(xa)n=n=01n!cos(a+nπ2)(xa)n
です。

特に、a=0の周りのテイラー展開、すなわちマクローリン展開は、
f(x)=n=0f(n)(0)n!xn=n=01n!cos(nπ2)xn
です。

ここで、cosxの性質から、

  • f(0)(0)0!=10!cos0=1
  • f(1)(0)1!=11!sin0=0
  • f(2)(0)2!=12!cos0=12!
  • f(3)(0)3!=13!sin0=0
  • f(4)(0)4!=14!cos0=14!
  • f(5)(0)5!=12!sin0=0

です。
従って、
f(x)=n=0f(n)(0)n!xn=n=0(1)n+12(n1)!xn
です。

次に収束半径です。
an=1n!cos(a+nπ2)とします。
このとき、
limn|anan+1|=limn|cos(a+nπ2)n!cos(a+(n+1)π2)(n+1)!|=limn|cos(a+nπ2)cos(a+(n+1)π2)n|
で、かつ、任意のxNに対して、|cosx|1ですので、R=、すなわち収束半径はです。
すなわち、任意のxRに対して、級数
sinx=n=0f(n)(a)n!(xa)n=n=01n!cos(a+nπ2)(xa)n
の右辺は収束するため、xにどんな実数を代入してもその値が実数値で求まる、ということです。

余談 tanxのテイラー展開、マクローリン展開にはベルヌーイ数という数が登場します。つまり、sinxcosxのようにスッキリとした形では表現できません。ベルヌーイ数まで踏み込むと量が膨大になるだけでなく、さらなる予備知識が要求されるので、今回は割愛します。

本当に近似できてるのかネ?

できてます。
それを視覚的に理解するため、筆者が作ったGeogebraのノートを貼ります。
このノートは変数をカーソルで動かせるようになっているため、バーを動かして色々試してみてください。

下のノートにおいて、f(x)=sinxgfのマクローリン展開、hfx=2π周りのテイラー展開です。
次数というのは、何次の項まで計算するか(何階の微分係数まで使うか)を表す数です。

今回はテイラー展開の代表例を真面目に計算してみました。
私調べですが、殆どの記事が「公式だ!」というように表示していて、「本当にテイラー展開可能なのか?」というところまで踏み込んでいなかったと思います。
そこで、真面目に展開が可能なのか、というところから計算をしてみました。

次回はコーシーの第2平均値定理です。

乞うご期待!質問、コメントなどお待ちしております!

コメントをする

タイトルとURLをコピーしました