logo

テイラー級数の応用 📂微分積分学

テイラー級数の応用

説明

テイラー級数（マクローリン級数）とは、与えられた関数を冪級数で近似したもので、関数ffのテイラー級数は次のようになる。

∑n=0∞f(n)(a)n!(x−a)n=f(a)+f′(a)(x−a)+f′′(a)2!(x−a)2+⋯ \sum\limits_{n=0}^{\infty} \dfrac{f^{(n)}(a)}{n!} (x-a)^{n} = f(a) + f^{\prime}(a)(x-a) + \dfrac{f^{\prime \prime}(a)}{2!} (x-a)^{2} + \cdots

良好な条件下では、ffとそのテイラー級数は実際に一致する。

f(x)=∑n=0∞f(n)(a)n!(x−a)n=f(a)+f′(a)(x−a)+f′′(a)2!(x−a)2+⋯ f(x) = \sum\limits_{n=0}^{\infty} \dfrac{f^{(n)}(a)}{n!} (x-a)^{n} = f(a) + f^{\prime}(a)(x-a) + \dfrac{f^{\prime \prime}(a)}{2!} (x-a)^{2} + \cdots

テイラー級数は極限を含んでいて一見するとより難しく見えるが、実際にはそうではない。関数ffが複雑な場合、テイラー級数で扱うほうが便利なこともある。基本的にテイラー級数は（無限次の）多項式なので、微分や積分、そして関数値計算が容易だ。

様々な状況で関数を冪級数として表現することがどのように有用かを見てみよう。

関数値の近似

例として、サイン関数sin⁡\sinを考えてみよう。こうした超越関数の関数値を計算するのは容易ではない。sin⁡x\sin xの値が簡単に分かるのはxx、00、π6\dfrac{\pi}{6}などの特殊角の場合だけだ。直ちにsin⁡0.672\sin 0.672の値を計算しろと言われても難しいだろう。このような場合、テイラー級数でsin⁡\sinの近似値を計算すればよい。sin⁡\sinのテイラー級数は次のようになる。

sin⁡x=x−x33!+x55!−x77!+⋯=∑n=0∞(−1)nx2n+1(2n+1)! \sin x = x - \dfrac{x^{3}}{3!} + \dfrac{x^{5}}{5!} - \dfrac{x^{7}}{7!} + \cdots = \sum\limits_{n=0}^{\infty} (-1)^{n} \dfrac{x^{2n+1}}{(2n+1)!}

つまり次のような近似式を使えば、適切な掛け算と足し算だけでsin⁡0.672\sin 0.672の値を簡単に計算することができる。

sin⁡x≈x−x33!+x55!−x77! \sin x \approx x - \dfrac{x^{3}}{3!} + \dfrac{x^{5}}{5!} - \dfrac{x^{7}}{7!}

実際、プログラミング言語では超越関数をテイラー級数で近似して実装することもある。ジュリアのコードを見てみれば、サイン関数がテイラー級数の7番目の項までで近似されていることが確認できる。

定数を用いて関数を直接定義して計算してみると、ジュリアで実装されたサイン関数とほぼ同じ値が得られる。

julia> function mysin(x)
           S = [-1.66666666666666324348e-01,
                8.33333333332248946124e-03,
                -1.98412698298579493134e-04,
                2.75573137070700676789e-06,
                -2.50507602534068634195e-08,
                1.58969099521155010221e-10]
           sinx = x +  S[1]*x^3 + S[2]*x^5 + S[3]*x^7 + S[4]*x^9 + S[5]*x^11 + S[6]*x^13
           return sinx
       end
mysin (generic function with 1 method)

julia> sin(0.672)
0.622552387351666

julia> mysin(0.672)
0.6225523873516658

区間[0,π/4][0, \pi/4]でsin⁡x\sin xとs6=∑n=06(−1)nx2n+1(2n+1)!s_{6} = \displaystyle \sum\limits_{n=0}^{6} (-1)^{n} \dfrac{x^{2n+1}}{(2n+1)!}を比較してみると、ほとんど違いがないことが分かる。

S = [-1.66666666666666324348e-01,
         8.33333333332248946124e-03,
         -1.98412698298579493134e-04,
         2.75573137070700676789e-06,
         -2.50507602534068634195e-08,
         1.58969099521155010221e-10]

x = 0:0.01:π/4
y = sin.(x)
s₆ = x .+ S[1]*x.^3 .+ S[2]*x.^5 .+ S[3]*x.^7 .+ S[4]*x.^9 .+ S[5]*x.^11 .+ S[6]*x.^13

using Plots
plot(x, y, label="sin(x)", linewidth=8, dpi=300, size=(728,300))
plot!(x, s₆, label="s_6", linewidth=3, ls=:dash, lc=:red)

無理数の近似

三角関数の計算も難しいのに、逆三角関数の計算はさらに難しくないか？アークタンジェントの関数値はアークタンジェントのテイラー級数を使って近似できるだろう。

tan⁡−1x=∑n=0∞(−1)nx2n+12n+1 \tan^{-1} x = \sum\limits_{n=0}^{\infty} (-1)^{n} \dfrac{x^{2n+1}}{2n+1}

上の式から面白い公式を得ることができ、x=1x = 1を代入すると次のようになる。

π4=tan⁡−11=1−13+15−17+⋯ \dfrac{\pi}{4} = \tan^{-1} 1 = 1 - \dfrac{1}{3} + \dfrac{1}{5} - \dfrac{1}{7} + \cdots

  ⟹  π=4(1−13+15−17+⋯ ) \implies \pi = 4 \left( 1 - \dfrac{1}{3} + \dfrac{1}{5} - \dfrac{1}{7} + \cdots \right)

つまり、無理数π\piを近似できる級数を得ることができる。

積分

関数をテイラー級数として表現することは、関数値の計算だけでなく積分をする際にも有利だ。それもそのはず、テイラー級数は多項式だからだ。多項式の積分は非常に簡単なので、複雑な形の関数はテイラー級数として表して積分すると簡単だ。例えばf(x)=sin⁡xx2f(x) = \dfrac{\sin x}{x^{2}}のような関数の不定積分を求めるのは難しい。（部分積分法が思い浮かぶかもしれないが、実際にやってみると行き詰まってしまう。）このような場合、sin⁡\sinのテイラー級数を積分するほうが簡単だ。

∫sin⁡xx2dx=∫1x2∑n=0∞(−1)nx2n+1(2n+1)!dx=∫1x2(x−x33!+x55!−x77!+⋯ )dx=∫(1x−x3!+x35!−x57!+⋯ )dx=C+ln⁡∣x∣−x22⋅3!+x44⋅5!−x66⋅7!+⋯=C+ln⁡∣x∣+∑n=1∞(−1)nx2n2n(2n+1)! \begin{align*} \int \dfrac{\sin x}{x^{2}} dx &= \int \dfrac{1}{x^{2}} \sum\limits_{n=0}^{\infty} (-1)^{n} \dfrac{x^{2n+1}}{(2n+1)!} dx \\ &= \int \dfrac{1}{x^{2}} \left( x - \dfrac{x^{3}}{3!} + \dfrac{x^{5}}{5!} - \dfrac{x^{7}}{7!} + \cdots \right) dx \\ &= \int \left( \dfrac{1}{x} - \dfrac{x}{3!} + \dfrac{x^{3}}{5!} - \dfrac{x^{5}}{7!} + \cdots \right) dx \\ &= C + \ln |x| - \dfrac{x^{2}}{2 \cdot 3!} + \dfrac{x^{4}}{4 \cdot 5!} - \dfrac{x^{6}}{6 \cdot 7!} + \cdots \\ &= C + \ln |x| + \sum\limits_{n=1}^{\infty} (-1)^{n} \dfrac{x^{2n}}{2n(2n+1)!} \end{align*}

別の例としてex2e^{x^{2}}のような関数を見てみよう。ex2e^{x^{2}}の積分というのは、微分積分学の基本定理により、微分してex2e^{x^{2}}になる関数を見つけることだ。このような関数を見つけようとすると、指数にx2x^{2}が含まれていて容易ではない。このような場合、ex2e^{x^{2}}を冪級数で表して不定積分するほうが簡単だ。

∫ex2dx=∫(1+x2+x42!+⋯+x2nn!+⋯ )dx=C+x+x33⋅1!+x55⋅2!+x77⋅3!+⋯+x2n+1(2n+1)n!+⋯=C+∑n=0∞x2n+1(2n+1)n! \begin{align*} \int e^{x^{2}} dx &= \int \left( 1+x^2+\dfrac{x^4}{2!}+\cdots +\dfrac{x^{2n}}{n!}+\cdots \right) dx \\ &= C + x + \dfrac{x^3}{3 \cdot 1!} + \dfrac{x^5}{5 \cdot 2!}+\dfrac{x^7}{7 \cdot 3!}+ \cdots + \dfrac{x^{2n+1}}{(2n+1)n!} + \cdots \\ &= C + \sum \limits_{n=0}^\infty \dfrac{x^{2n+1}}{(2n+1)n!} \end{align*}

調和振動子

任意の関数を冪級数で表現する方法はいろいろな分野で応用される。物理学では、単純調和振動子のポテンシャルエネルギーをV(x)=a0+a1x+a2x2+⋯V(x) = a_{0} + a_{1}x + a_{2}x^{2} + \cdotsのような級数で与え、フックの法則F(x)=−kxF(x) = -kxと単純調和振動子の運動方程式x¨+kxx=0\ddot{x} + \dfrac{k}{x}x = 0を導くことができる。

微分方程式の解

いくつかの微分方程式は解が冪級数で表現される。ベッセル微分方程式とその解（ベッセル関数）は次の通り。

x2y′′+xy′+(x2−ν2)y=0 x^2 y^{\prime \prime} +xy^{\prime} +(x^2-\nu^2)y = 0

Jν(x)=∑n=0∞(−1)nΓ(n+1)Γ(n+ν+1)(x2)2n+ν J_{\nu}(x)=\sum \limits_{n=0}^{\infty} \frac{(-1)^{n} }{\Gamma (n+1) \Gamma (n+\nu+1)} \left(\frac{x}{2} \right)^{2n+\nu}