テイラー展開とは、関数を無限個の多項式の和で表す手法のことである。
一次関数や二次関数、三次関数といった多項式関数は、xの値を与えれば、関数の値を計算するのは難しくない。
しかし、指数関数や三角関数といった関数については、特殊なxの値以外については、関数の値を計算するのは難しい。
そこで、あらゆる関数を多項式で表すことができれば値が計算できるのに…となるが、そのような場合に特に役に立つのがテイラー展開である。
以下、f(x)を、定義域内の任意の点で任意の回数だけ微分可能なxの関数とする。
ここで、ある一点x=aを決め、x=aの近傍で、f(x)を(x-a)のk乗の和(kは0以上の整数)で近似することを考える。(x-a)kの係数をak(a)とし(係数がaに依存することから)、n次の項まで考える。
当然、n次の項までの和とf(x)には誤差が生じるので、この誤差である剰余項をRn+1と表現すると(n+1乗以上の項をまとめたようなものと考えると分かりやすい)、
f(x)=Σ(k=0→n){ak(a)・(x-a)k}+Rn+1 …①
となる。
剰余項には多種多様な表現があるので、詳細は解析学の教科書等を参照のこと。
ここで、n→∞としたとき、剰余項Rn+1が0に収束するならば、f(x)を無限個の(x-a)の冪で表現することが可能であり、式①は
f(x)=Σ(n=0→∞){an(a)・(x-a)n} …②
と表現される。
肝心の係数an(a)であるが、結論を言えば、f(n)(x)をf(x)のn次導関数として、
an(a)= f(n)(a)/n!
となるので、式②は
f(x)=Σ(n=0→∞)[{f(n)(a)/n!}・(x-a)n] …③
と表現される。これを「f(x)のx=aのまわりでのテイラー展開」などと表現する。(なお、0次微分f(0)(x)はf(x)そのものとし、0!=1、(x-a)0=1にも注意。)
厳密な説明ではないが、③の両辺をn回微分してx=aを代入すれば両辺が一致し、係数an(a)の表式が妥当であることが確かめられる。
n→∞のとき剰余項Rn+1が収束し、f(x)のx=aのまわりでのテイラー展開が可能となるxの範囲であるが、
場合に分けられる。2.の場合におけるρを収束半径という。なお1.の場合は、収束半径は∞という。
テイラー展開の中でも、x=0の周りのテイラー展開を考える場合が多く、これをマクローリン展開という。③にa=0を代入すると、関数f(x)のマクローリン展開は
f(x)=Σ(n=0→∞)[{f(n)(0)/n!}・xn] …④
と表現される。
実際にテイラー展開を用いて関数の値を計算するときは、当然無限項の和は計算できないので、有限個の項で打ち切って計算することとなる。項数は求める精度によるが、一般的にはxがaに十分近い場合や、係数の分母に階乗がある場合は少ない項数でも高精度で計算することができる。
なお、ここでは変数xを実数として考えたが、変数を複素数に拡張してテイラー展開を考えることもできる。
以下に、いくつかの関数のx=0のまわりのテイラー展開(マクローリン展開)を掲載する。主に、表式が簡単となる関数について掲載する。
なお、三角関数は弧度法とし、対数関数は自然対数とする。
xの変域がない式は収束半径が∞であり、任意のxについて成立する。
マクローリン展開の式でも、微分でsinx→cosx→-sinx→-cosx→sinx→…と変化することが分かる。
sinxは奇関数なので奇数次のみ、cosxは偶関数なので偶数次のみとなる。
収束半径は1。収束半径上の点x=±1においても成立。奇関数なので奇数次のみ。
logxがx=0において定義されないので、log(1+x)としてx=0のまわりの展開を考える。
収束半径は1。収束半径上の点x=-1では発散するが、x=1においては成立する。
sinhx=(ex-e-x)/2、coshx=(ex+e-x)/2という表式通りの形となる。
三角関数の式の符号を全てプラスにした形となる。
収束半径は1。収束半径上の点x=±1では発散する。
logを使った表式のため、logの値の計算に利用でき、値によっては収束が速くて有用。
掲示板
掲示板に書き込みがありません。
急上昇ワード改
最終更新:2024/04/20(土) 04:00
最終更新:2024/04/20(土) 04:00
ウォッチリストに追加しました!
すでにウォッチリストに
入っています。
追加に失敗しました。
ほめた!
ほめるを取消しました。
ほめるに失敗しました。
ほめるの取消しに失敗しました。