閉じる ホーム

Taylor展開の計算

≪Taylor展開の計算≫

二つの関数 g(x),h(x) の Taylor 多項式は計算できるものとしたときに、

g(x)+ h(x) 2g(x) g(x)h(x) g(x) f(x)
などの関数の Taylor 多項式をどのようにして求めることができるのか、すなわち「 Taylor 展開が計算できる関数たちの組み合わせ」として表わされる関数の Taylor展開が、それぞれの関数の Taylor展開からどのようにして計算できるのか、を考えてみよう。

まず、滑らかな関数g(x)
f(x)= f(0)+ f(0)'x + f''(0)2!x2 + f'''(0)3!x3 +⋯+ f(n)(0)n!xn + Rn(x)  ……(1)

Rn(x) = 1n! 0x (x-t)n f(n+1(t)dt  ……(2)
と書き表わすことができる。
いま、(1)式の右辺に現われる多項式の係数を
ak= f(k)(0) k!
として、Taylor多項式を単に

Pn(x)= a0+ a1x + a2x2 + a3x3 +⋯+ anxn  ……(3)
と書き表わすことにする。(以下、同様)
また(2)式の剰余項については、t=xs として、右辺に現われる積分の積分変数を t から s に変数変換すると (t(0x) は s(01),dtxds)

Rn(x) = 1n! 01 (x-xs)n f(n+1)(xs)xds

= xn+1n!  01 (1-s)n f(n+1)(xs)ds
ここで

f(n+1)(x) = 1n! 01 (1-s)n f(n+1)(xs)ds
とすると

Rn(x) = f(n+1)(x)·xn+1
したがって

f(x)= a0+ a1x + a2x2 + a3x3 +⋯+ anxn + f(n+1)(x)·xn+1  ……(4)

また「積分に関する平均値の定理」を用いると

f(n+1)(x)= f(n+1)(θ)(n+1)!    ……(5)
となるような実数 θ が 0 と x の間に存在する。
このとき、x0 は θ0 となるから

lim x0 f(n+1)(x)= lim θ0 f(n+1)(θ)(n+1)!  

 = f(n+1)(0)(n+1)!  = f(n+1)(0)  ……(6)
という有限な値になることが分かる。

以上、前準備ができたので、まず

f(x)= g(x)+ h(x)
の場合について考えてみよう。
いま、 g(x), h(x) の Taylor展開を、それぞれ

g(x)= b0+ b1x + b2x2 + b3x3 +⋯  ……(7)
h(x)= c0+ c1x + c2x2 + c3x3 +⋯  ……(8)
と表わすことにすると

f(x)= g(x)+ h(x)
 = (b0+ b1x + b2x2 + b3x3 +⋯)
  + (c0+ c1x + c2x2 + c3x3 +⋯)
  = (b0+c0)+ (b1+c1)x + (b2+c2)x2 + (b3+c3)x3 +⋯  ……(9)
と表わせることが分かるから、(9)式が Taylor展開を与えるのではないか、と思われる。

そこで、 g(x), h(x) の Taylor多項式を、それぞれ

Gn(x)= b0+ b1x + b2x2 +⋯+ bnxn  ……(10)
Hn(x)= c0+ c1x + c2x2 +⋯+ cnxn  ……(11)
と表わすことにすると、 g(x), h(x)は(4)式を適用して

g(x)= Gn(x)+ g(n+1)(x)·xn+1  ……(12)
h(x)= Hn(x)+ h(n+1)(x)·xn+1  ……(13)
と表わせるから

f(x)- (Gn(x) +Hn(x) )
 =( g(x)+ h(x))- (Gn(x) +Hn(x) )
 =( g(x)- Gn(x))+ (h(x)- Hn(x) )
 = g(n+1)(x)·xn+1 + h(n+1)(x)·xn+1
 ={ g(n+1)(x) + h(n+1)(x)}xn+1
ここで、k=0,1,2,3,…,n に対して

lim x0 f(x)- (Gn(x) +Hn(x) )       xk      

 = lim x0 ( { g(n+1)(x) + h(n+1)(x)}xn+1 )
 = { g(n+1)(0) + h(n+1)(0)}· 0
 =0
となることが分かる。
ところで、「Taylor多項式の考察」でみたように、k=0,1,2,3,…,n に対して

lim x0 f(x)- Qn(x)     xk     =0
を満たすような n 次の多項式 Qn(x) は、関数 f(x) の Taylor多項式 Pn(x) しか存在しないから

Pn(x) = Gn(x) + Hn(x)
となることが分かる。
したがって、f(x)=g(x)+h(x) の Taylor多項式は(9)式のように計算され、関数 f(x) の Taylor展開を与えていることが分かる。

(例) f(x)=ex+e-x の Taylor展開は

f(x) =( 1+x+x22! +x33! +⋯)
    +( 1-x+x22! -x33! +⋯)
=2( 1+x22! +x44! +⋯)

次に、スカラー倍 f(x)=Cg(x) をみてみよう。

いま(7)式の両辺に C を掛けてみると f(x)=Cg(x)
f(x)=Cg(x)
=C( b0+ b1x + b2x2 + b3x3 +⋯)
= Cb0+( Cb1)x +( Cb2)x2 +( Cb3)x3 +⋯

(12)式より
f(x)-CGn(x) = Cg(x)-CGn(x)
  = C(g(x)-Gn(x))
  = C(g(n+1)(x)·xn+1)
ここで、k=0,1,2,3,…,n に対して

lim x0 f(x)- CGn(x)       xk      

 = lim x0 ( C· g(n+1)(x) ·xn+1 )
 = C· g(n+1)(0) · 0
 =0
したがって、前と同様、関数 f(x) の Taylor多項式 Pn(x)
Pn(x) = CGn(x)
となることが分かる。

(例)
g(x) = 2sinx   の Taylor展開は
= 2( x-x33! +x55! -x77! +⋯)
= 2x-x33 +x560 -⋯)

さらに、f(x)=g(x)h(x)    の場合、
(7),(8)式より
f(x) = g(x)h(x)
 =( b0+ b1x + b2x2 + b3x3 +⋯)
  ·( c0+ c1x + c2x2 + c3x3 +⋯)
 =( b0+c0)+( b1c0+ b0c1) x
  +( b0c2+ b1c1+ b2c0) x2 +⋯  ……(14)
さて、Taylor多項式の積 Gn(x)·Hn(x) の多項式において( l,m=0,1,2,⋯,n ) n次式以下の部分と (n+1)次式以上の部分をそれぞれ
Qn(x) = 0l+mn blcm xl+m
Q_n(x) = n+1l+m2n blcm xl+m
とすると
Gn(x)·Hn(x) = Qn(x)+ Q_n(x)  ……(15)
また(12),(13)式より
g(x)h(x) = (Gn(x)+ g(n+1)(x)xn+1)
  · (Hn(x)+ h(n+1)(x)xn+1)
 = Gn(x)Hn(x)
  + Gn(x)· h(n+1)(x)xn+1
  + Hn(x)· g(n+1)(x)xn+1
  + g(n+1)(x)xn+1· h(n+1)(x)xn+1
したがって
f(x) - Qn(x) = g(x)h(x) - Qn(x)
 = Gn(x)Hn(x) - Qn(x)  ⋯ (= Q_n(x))
  + Gn(x)· h(n+1)(x)xn+1
  + Hn(x)· g(n+1)(x)xn+1
  + g(n+1)(x)xn+1· h(n+1)(x)xn+1  ……(16)
ところで、(16)式の右辺に現われるすべての項が xn+1 で括れるから
k=0,1,2,3,…,n に対して

lim x0 f(x)- Qn(x)      xk     =0
となることが分かる。したがって

Pn(x) = Qn(x)
となることが分かる。

(例)  g(x) = ex h(x) = sinx  とすると (14)式より  ex sinx  の Taylor展開は

ex sinx =( 1+x+x22!+x33!+⋯)
  ·( x-x33!+x55!-x77!+⋯)
 = x+x2+(12!-13!)x3+⋯
 = x+x2+13x3+⋯
というように計算される。
最後に、h(0)0 として
f(x) = g(x) h(x)
の Taylor展開を考えてみよう。

いま、 f(x) = g(x) · 1 h(x)  とすると積の計算であるから、問題は  1h(x) を求めることにある、と考えることができる。
そこで、 h^(x)= 1h(x)  として関数 h^(x) の Taylor展開を
h^(x) = d0+ d1x + d2x2 + d3x3 +⋯  ……(17)
h(x) = c0+ c1x + c2x2 + c3x3 +⋯  ……(18)
とすると、 いま、h(x) h^(x) =1 より

1= h(x) h^(x)
 =( c0+ c1x + c2x2 + c3x3 +⋯)
  ·( d0+ d1x + d2x2 + d3x3 +⋯)
 = c0d0+ ( c0d1+ c1d0 )x
  +( c0d2+ c1d1+ c2d0 )x2 +⋯
 = 1+0·x+0·x2 +⋯
係数を比較して
c0d0 =1
c0d1+ c1d0 =0
c0d2+ c1d1+ c2d0 =0
 ⋮
いま、h(0)0 と仮定しているから c00 に注意すると上記の連立方程式を(18)式の係数を適用して上から順番に解くことができる。
例えば h(x)=cosx として h^(x)=1cosx とすると

1= h(x) h^(x)
 =( 1-x22! +x44! -x66! +x88! +⋯)
  ·( d0+ d1x + d2x2 + d3x3 +⋯)
 = d0+ d1x+ (d2-d02!)x2+ (d3-d12!)x3+ (d4+d04!-d22!)x4 +⋯
 = 1+0·x+0·x2+0·x3+0·x4 +⋯
より
d0=1
d1=0
d2=d02=12
d3=d12!=0
d4=d22!-d04! =524
 ⋮
したがって
h^(x)=1cosx
 = 1+12x2+524x4 +⋯  ……(19)
ところで、

tanx =sinx ·1cosx   より

 =( x-x33! +x55!+⋯)

  ·( 1+x22! +524x4+⋯)

 = x+(12!-13!)x3+(15!-13!·2!+524)x5+⋯

 = x+13x3 +215x5+⋯
というような計算もできる。

<「牛腸作 数学IB演習」・独習ノートより>

(1)テーラー展開とはなにか
(2)テーラー展開の注意点
(3)部分積分とテーラーの定理
(4)テーラーの定理・剰余項の考察
(5)テーラー多項式の考察
(7)合成関数のテーラー展開
(8)近似式としてのテーラー展開
(9) a のまわりでのテーラー展開
(10)テーラーの定理・極限
(11)多変数のテーラー展開
補論・積分に関する「平均値の定理」
補論・発散のスピード
補論・級数の収束判定
補論・ロルの定理・考察
補論・偏微分(1)偏導関数
補論・偏微分(2)ヤングの定理
余録・バーゼル問題とテーラー展開

  © nagaizumicyou.net