...

テイラー展開 - OGIFEST

by user

on
Category: Documents
21

views

Report

Comments

Transcript

テイラー展開 - OGIFEST
テイラー展開(Taylor expansion)
例.(1.0007)15 を考える。
馬鹿正直に 1.0007 を 15 回かける
→具体的な値(=1.01055)
そんなもん大体 1 だろ!
→テイラー展開 0 次近似
大体「1+15×0.0007=1.0105」くらい?
→テイラー展開 1 次近似
x が 0.0007 みたいにすごく小さいとき
(1+x)15 の代わりに
1+15x を計算
sin(x)の代わりに
x を計算
cos(x)の代わりに
1-0.5x2 を計算
log(x)の代わりに
x-1 を計算
これがテイラー展開!
最初の例で 1 と 0.0007 を分離していることに注目する!
これはすごく小さな部分を分離して、まずは「大体 1」と考えている
これが「x=a の周りでテイラー展開をしなさい」ということ(今回は a=1)
x が a に近いんだったら f(x)は f(a)に近いはず
つまり f(x)~f(a)となる。
a が定数なので f(a)も定数、つまり f(x)を定数で近似していることになる
x=a の周りで考える、ってことは x が a に近いんだったら x-a は微小な量である
今回の場合だったら(x-a)=0.0007 となっている
もちろん微小な量である(x-a)を 2 乗 3 乗したものはもっと微小な数値となる
テイラー展開のポイント
微小な量である(x-a)の多項式 の形にするところ
仮の係数 c を用いて式にすると
f(x)~c0+c1(x-a)+c2(x-a)2+c3(x-a)3+c4(x-a)4+…
c の求め方は最後に記述するが、
c0=f(a)
…f(0)とは限らないことに注意
c1=f ’(a)/1!
…f ’(0)とは限らないことに注意
c2=f ’’(a)/2!
…f ’’(0)
〃
c3=f ’’’(a)/3!
…f ’’’(0)
〃
c4=f ’’’’(a)/4! …f ’’’’(0)
〃
…
cn=f(n)(a)/n! となる。
そして c が上記であることを踏まえた上で、
テイラー展開の完成型は以下のようになる
f(x)~c0+c1(x-a)+c2(x-a)2+c3(x-a)3+c4(x-a)4+…+cn(x-a)n+…
項が進むごとにどんどん微小な値を用いることになり、それは無限に続くので
どこで打ち切るかを決めてよい
0 次近似
f(x)~c0
微小な値(x-a)を無視すれば f(x)~f(a)となる
最初に説明した通り定数で近似したこととなり、今回の例で言うと 1 と近似される
1 次近似(線形近似)
f(x)~c0+c1(x-a)
ここで co=f(a)と c1=f ’(a)を代入すれば
f(x)~f(a)+f ’(a)(x-a) となる
これが 1 次近似、直線で近似しているので線形近似とも言う
今回の例の場合、a=1,x-a=0.0007 としているので
f(x)=x15
→ f ’(x)=15x14
→
f ’(a)=15a14=15
より
f(x)~f(a)+f ’(a)(x-a) =115+15×0.0007=1+15×0.0007 となり、最初の説明と一致する
2 次近似
cos(x)のような場合、x=0 の周りでテイラー展開すると 1 次近似しても結果が 1 となる
よって詳しく近似するなら 2 次の項まで取って近似すると cos(x)~1-0.5x2 となる
※c の求め方
f(x)~c0+c1(x-a)+c2(x-a)2+c3(x-a)3+c4(x-a)4+…
x に a を代入すると
f(x)~c0+c1(x-a)+c2(x-a)2+c3(x-a)3+c4(x-a)4+…
f(a)~ c0+c1(a-a)+c2(a-a)2+c3(a-a)3+c4(a-a)4+…
当然(a-a)は 0 となるので
c0=f(a) となる
一度 f(x)を微分してから a を代入する
f(x)~c0+c1(x-a)+c2(x-a)2+c3(x-a)3+c4(x-a)4+…
f ’(x)~0+c1+2 c2(x-a)+3 c3(x-a)2+4 c4(x-a)3+…
f ’(a)~c1+2 c2(a-a)+3 c3(a-a)2+4 c4(a-a)3+…
となって c1= f ’(a) となる
この展開を繰り返してゆくと
c0=f(a)
= f(a)
c1=f ’(a)
= f ’(a)
c2=f ’’(a)/2
= f ’’(a)/2!
c3=f ’’’(a)/(3・2)
= f ’’’(a)/3!
c4=f ’’’’(a)/(4・3・2)
= f ’’’’(a)/4!
よって
cn=f(n)(a)/n!
と書くことができる!
※ついでに f ’(a)などの求め方を一応
f(x)を微分したのが f ’(x)
それに x=a を代入したのが f ’(a)
f(x)を 2 回微分したのが f ’’(x)
それに x=a を代入したのが f ’’(a)
f(x)を 3 回微分したのが f ’’’(x)
それに x=a を代入したのが f ’’’(a)
f(x)を n 回微分したのが f(n)(x)
それに x=a を代入したのが f(n)(a)
これを n の階乗で割ると cn=f(n)(a)/n! を求めることができる
Fly UP