[論文レビュー] Formulae of numerical differentiation

Maxim Dvornikov|ArXiv.org|Jun 5, 2003

Iterative Methods for Nonlinear Equations参考文献 7被引用数 22

ひとこと要約

本稿では、補間多項式を構築する必要なく、任意の数の等間隔離散点を用いた1階および2階微分の明示的閉形式式を導出する。この手法は重み係数のスペクトル解析を用いて、広範な空間周波数範囲で高い精度を示し、無限に多くの点の極限において正確な微分に収束することを検証している。検証は sin(x) を用いて実施された。

ABSTRACT

We derived the formulae of central differentiation for the finding of the first and second derivatives of functions given in discrete points, with the number of points being arbitrary. The obtained formulae for the derivative calculation do not require direct construction of the interpolating polynomial. As an example of the use of the developed method we calculated the first derivative of the function having known analytical value of the derivative. The result was examined in the limiting case of infinite number of points. We studied the spectral characteristics of the weight coefficients sequence of the numerical differentiation formulae. The performed investigation enabled one to analyze the accuracy of the numerical differentiation carried out with the use of the developed technique.

研究の動機と目的

• 与えられた離散的かつ等間隔の点における関数の1階および2階微分を、補間多項式を構築することなく計算する一般化された手法を開発すること。
• 従来の数値微分公式が通常、固定で少数の補間点に制限されているという限界を克服すること。
• 微分公式における重み係数のスペクトル的性質を分析し、異なる空間周波数における精度を評価すること。
• sin(x) の微分を計算し、点数が無限に近づく際の収束を検証することで、手法を検証すること。
• 任意の数の補間点に対して有効な、明示的かつ閉形式の重み係数の式を提供することで、数値シミュレーションにおける効率的実装を可能にすること。

提案手法

• 関数値の対称的差分を用いて、偶数次多項式係数を消去し、微分近似に寄与する奇数次項に集中する。
• 差分 $ d_m = f_m - f_{-m} $ を用いて、奇数次多項式係数のための一次方程式系を定式化し、これは $ m $ の奇数乗の和として表現される。
• この方程式系の解はクラメルの公式を用いて導出され、重み係数 $ \alpha_m^{(2k+1)}(n) $ を用いた $ c_{2k+1} $ の明示的表現が得られる。
• 1階および2階微分の重み係数 $ \alpha_m^{(1)}(n) $ および $ \alpha_m^{(2)}(n) $ は、バテルマン行列に類似た行列の行列式を用いて計算される。
• 重み係数系列のスペクトル特性は、離散フーリエ変換を計算することで分析され、理想の微分スペクトルと比較される。
• MATLAB において、フル行列の逆行列計算を回避する再帰的積分による反復的アルゴリズムを用いて、重み係数が実装された。

実験結果

リサーチクエスチョン

• RQ1 補間多項式を構築することなく、任意の数の補間点に対して数値微分の明示的閉形式式を導出可能か？
• RQ2 重み係数系列のスペクトル的性質は、異なる空間周波数における数値微分の精度にどのように影響するか？
• RQ3 補間点数を増加させた際の数値微分の収束挙動はいかなるものか、特に $ n \to \infty $ の極限においては？
• RQ4 滑らかな関数（例：$ \sin x $）に対して、導出された公式の精度は解析的微分と比べてどの程度か？
• RQ5 この手法は、量子場の理論の高精度格子シミュレーションへの応用を想定した、効率的かつスケーラブルに実装可能か？

主な発見

• 1階および2階微分の導出式は、任意の数 $ n $ の補間点に対して閉形式で表され、多項式補間の必要がなくなる。
• スペクトル解析により、1階微分の重み係数系列が、空間周波数の広い範囲で理想の微分スペクトル $ (2\pi r/N)^2 $ をよく近似していることが示された。特に $ n $ が増加するにつれてその近似が向上する。
• $ n=11 $ の場合、空間周波数 $ r \lesssim 750 $ まで精度が保たれ、$ n=21 $ の場合、$ r \lesssim 800 $ まで精度が維持される。これは高周波数領域でも高い精度を示している。
• $ n \to \infty $ の極限において、$ \sin x $ の微分が正確に収束することが確認され、理論的整合性が裏付けられた。
• 2階微分の公式に対しても、広い周波数帯域にわたり、理想の放物型応答からのスペクトルずれが小さく、高い精度を示した。
• MATLAB 実装では、再帰的積分を用いて重み係数を効率的に計算でき、大規模シミュレーションにおける実用的応用が可能となった。