LINEスタンプ制作代行サービス・LINEスタンプの作り方!

お電話でのお問い合わせ:03-6869-8600

stampfactory大百科事典

差分法

数値解析における有限差分法(ゆうげんさぶんほう、)あるいは単に差分法は、微分方程式を解くために微分を有限差分近似(差分商)で置き換えて得られる差分方程式で近似するという離散化手法を用いる数値解法である。18世紀にオイラーが考案したと言われる。今日ではFDMは偏微分方程式の数値解法として支配的な手法である解の誤差とは、真の解析解と近似解との間の差として定義される。有限差分法における誤差の原因は丸め誤差および打ち切り誤差または離散化誤差である。問題に対する解の近似に有限差分法を用いるためには、まず初めに問題の領域を離散化しなければならない。これは普通は、その領域を一様な格子に分ければよい。これは有限差分法がしばしば「時間刻み」な仕方で微分に対する離散的な数値近似の集合を提供することを意味することに注意。一般に注目すべきはで、典型的にはこれをO-記法で表す。局所打ち切り誤差は、各点における誤差について言うもので、真値 と近似値 との差として書ける。この誤差の評価には、テイラー展開の剰余項を見るのが簡便である。式 に対するテイラー展開のラグランジュ型剰余項から、局所打ち切り誤差の支配項が求められる。例えば、一階差分近似 () を考えればであり、適当な代数的操作でと書きなおすと、この左辺は有限差分法で得られる近似値であり、右辺が真値と剰余項の和である。明らかにこの剰余項が局所打ち切り誤差であり、そのオーダーがとなることが分かる。この場合、局所打ち切り誤差は刻み幅に比例するということになる。有限差分法の近似解の精度と計算量は方程式の離散化の仕方や刻み幅の取り方に依存する。これらは刻み幅を小さくするにつれ著しく増加するから、実用上は必要な精度と計算時間を天秤にかけて十分合理的な条件で近似を行う。時間の刻み幅が大きければ多くの場合に計算速度は早くなるが、大きくしすぎると不安定性を生じ、データの精度に問題がでる。数値モデルの安定性を決定するために、フォン・ノイマンの安定性解析を用いるのが普通である。最も簡単な例として、次の1階常微分方程式を考える:これを解くには、差分商を用いてと近似する。この方法をオイラー法という。この最後の方程式のように、微分方程式の微分を差分商に置き換えたものを、差分方程式(さぶんほうていしき、)と呼ぶ。偏微分方程式の例として、一様ディリクレ境界条件に従う1次元規格化熱伝導方程式を考える:左辺は時間"t" による微分、右辺は座標"x" による2階微分である。また、境界条件および初期条件は以下とする:これを数値的に解く1つの方法は、すべての微分を差分で近似することである。空間の領域をメッシュformula_13で、時間の領域をメッシュformula_14で分割しよう。どちらの分割も等間隔とし、空間点の間隔を"h"、時刻の間隔を"k" とする。formula_15の数値的近似をformula_16で表す。時刻formula_17には前進差分を用い、空間点formula_18で2次微分に対して2次中央差分を用いれば、次の漸化式:が得られる。これを陽解法という。formula_20の値は次のように得られる:ただしここでformula_22 (拡散数と呼ばれる)である。ゆえに、時刻"n" での値がわかれば、対応する時刻"n" + 1 での値も漸化式を用いて求められる。formula_23とformula_24には境界条件(この例ではどちらも0)を適用する。この陽解法は、"r" ≤ 1/2 であれば数値的に安定で収束することが知られている。誤差は時間ステップ数と空間ステップ数の2乗とに比例する:時刻formula_26に後退差分を用い、空間点formula_18 で2階中央差分を用いれば、漸化式:が得られる。これを陰解法という。線形方程式系:を解けば、formula_20が得られる。この方法は常に数値的に安定で収束するが、時刻ごとに方程式系を解く必要があるため、陽解法よりも繁雑である。誤差は時間ステップ数と空間ステップ数の4乗とに比例する。さいごに、時刻formula_31で中央差分を、空間点formula_18での空間微分に2階中央差分を用いれば、漸化式:が得られる。これをクランク・ニコルソン法()という。線形方程式系:を解けば、formula_20が得られる。この方法は常に数値的に安定で収束するが、各時刻で方程式系を解く必要があるので繁雑なことが多い。誤差は時間ステップ数の4乗と空間ステップ数の2乗とに比例する:しかし、境界付近では誤差はO("h" ) でなくO("h" ) となることが多い。クランク・ニコルソン法は時間ステップ数が少なければたいてい最も正確な方法である。陽解法はそれより正確でなく不安定でもあるが、最も実行しやすく、繁雑さも最も少ない。陰解法は時間ステップ数が多い場合に最も優れている。

出典:wikipedia

LINEスタンプ制作に興味がある場合は、
下記よりスタンプファクトリーのホームページをご覧ください。