sin x のテイラー展開(マクローリン展開)
x を弧度法で測るとき、すべての実数 x について
sinx=x−3!x3+5!x5−7!x7+⋯n! は 1⋅2⋅3⋯n(階乗)。
この記事の要点
- limxsinx=1 は「x が小さいとき sinx≒x」という意味。その誤差が 6x3 くらいであることが、微分を重ねた不等式で分かる
- 同じ方法を続けると、sinx を x の多項式で表す規則が見える。これがテイラー展開
- 項を増やすほど、0 から遠いところまで sinx に重なる
limxsinx=1 の証明は、x が 0 に近いとき sinx が x とほぼ同じ値になる、ということを示していました。このコラムでは「ほぼ同じ」のずれを調べ、その先にある大学数学のテイラー展開へ進みます。
高校の問題:sin x と x のずれはどれくらいか
問題
x→0limx3x−sinx を求めよ。
x→0 で分母も分子も 0 に近づきます。x−sinx が x3 と同じ速さで 0 に近づくなら、極限は 0 でない値になるはずです。
微分で不等式をはさんで解く
x>0 とします。次の2つの事実を使います。
- sinx<x(面積のはさみうちで得た不等式)
- f′(x)>0 かつ f(0)=0 なら、x>0 で f(x)>0(f は増加関数だから)
1. cosx を下から押さえる
f1(x)=cosx−1+2x2 とおくと f1(0)=0 です。微分すると
f1′(x)=−sinx+x=x−sinx>0したがって f1(x)>0、つまり cosx>1−2x2 です。
2. sinx を下から押さえる
f2(x)=sinx−x+6x3 とおくと f2(0)=0 で、
f2′(x)=cosx−1+2x2=f1(x)>0なので f2(x)>0、つまり sinx>x−6x3 です。
3. cosx を上から押さえる
f3(x)=1−2x2+24x4−cosx とおくと f3(0)=0 で、
f3′(x)=−x+6x3+sinx=f2(x)>0なので cosx<1−2x2+24x4 です。
4. sinx を上から押さえる
f4(x)=x−6x3+120x5−sinx とおくと f4(0)=0 で、
f4′(x)=1−2x2+24x4−cosx=f3(x)>0なので sinx<x−6x3+120x5 です。
5. はさみうちの原理を使う
2 と 4 をあわせて、x>0 で
x−6x3<sinx<x−6x3+120x5各辺を −1 倍して x を足し、x3 で割ると
61−120x2<x3x−sinx<61x→+0 で左辺は 61 に近づくので、x→0limx3x−sinx=61
です。x3x−sinx は x→−0 でも同じ値(x を −x にしても変わらない)なので、極限は 61 です。
つまり x−sinx≒6x3、すなわち sinx≒x−6x3 という、x だけの近似より一段よい近似が得られました。
規則が見える:sin x のテイラー展開
上の流れで出てきた式を並べます。
sinx≒x,sinx≒x−6x3,sinx≒x−6x3+120x5分母は 6=3!、120=5! で、符号は +,−,+ と交互に変わります。このまま項を増やしたものが、冒頭の式です。
項を増やすと、原点から遠いところまで sin x と重なる数値でも確かめます。x=1 のとき sin1=0.841471⋯ で、
| 近似式 | 値 | sin1 とのずれ |
|---|
| x | 1 | 0.1585 |
| x−3!x3 | 0.8333⋯ | 0.0081 |
| +5!x5 | 0.84167⋯ | 0.00020 |
| −7!x7 | 0.841468⋯ | 0.0000027 |
項を1つ足すごとに、ずれが大きく減ります。
係数はどう決まるのか
高校の不等式では、式を「当たりをつけて」書きました。大学では、係数を機械的に決めます。f(x)=c0+c1x+c2x2+c3x3+⋯ と表せたとして、x=0 を代入すると f(0)=c0。微分してから x=0 を代入すると f′(0)=c1。もう一度微分すると f′′(0)=2c2、さらに微分して f′′′(0)=3⋅2⋅c3 となり、一般に
cn=n!f(n)(0)です。f(n) は n 回微分した関数。f(x)=sinx なら、(sinx)′=cosx を繰り返して、x=0 での値は 0, 1, 0, −1, 0, 1, ⋯ の順になります。偶数番目が 0、奇数番目が ±1 で、3!x3 の符号が − になる理由も、これで分かります。
その先:大学数学へ
x=0 の近くだけでなく、x=a のまわりで同じことをした
f(x)=f(a)+f′(a)(x−a)+2!f′′(a)(x−a)2+⋯をテイラー展開、a=0 の場合をマクローリン展開と呼びます。1715年にブルック・テイラー(イギリス)が一般の式を発表し、1742年にコリン・マクローリン(スコットランド)が著書でその特別な場合を詳しく扱いました。
- べき級数:x の累乗を無限に足した式。ex=1+x+2!x2+3!x3+⋯ や cosx=1−2!x2+4!x4−⋯ も展開できます。eiθ の項を実部と虚部に分けると cosθ+isinθ になり、オイラーの公式も確かめられます。
- 収束:無限に足した結果が決まった値に近づくこと。sinx は全ての x で収束しますが、1−x1=1+x+x2+⋯ は ∣x∣<1 のときだけです。
- ロピタルの定理・ランダウの記号:問題の極限は、sinx=x−6x3+(x5 以上の項) を代入すれば一行で出ます。大学の極限計算の基本道具です。
- 物理:単振り子の運動方程式で sinθ≒θ とするのは、テイラー展開の第1項だけを使うことです。θ=0.2rad(約 11∘)でも誤差は約 0.7% です(単振動)。
注意
sinx≒x−6x3 が使えるのは、x が 0 に近いときだけです。図のように、x が大きくなると多項式は sinx から離れていきます。また、高校では「無限に足す」ことの厳密な扱いを学びません。このコラムの展開は、考え方の紹介です。
この記事のまとめ
- x→0limx3x−sinx=61 は、f′>0 と f(0)=0 を繰り返す不等式のはさみうちで求まる。
- その流れを続けると sinx=x−3!x3+5!x5−⋯ が見える。係数は n!f(n)(0)。
- テイラー展開は、関数を多項式で表す大学数学の基本で、ex・オイラーの公式・物理の近似につながる。