[论文解读] Formulae of numerical differentiation

Maxim Dvornikov|ArXiv.org|Jun 5, 2003

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

一句话总结

本文推导出使用任意数量等距离散点进行一阶和二阶导数中心数值微分的显式闭式公式，无需构建插值多项式。该方法通过权重系数的谱分析，证明在广泛的空间频率范围内具有高精度，且在点数趋于无穷时收敛至精确导数，已在 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.

研究动机与目标

• 开发一种通用方法，用于在不构建插值多项式的情况下，计算给定离散等距点上函数的一阶和二阶导数。
• 解决现有数值微分公式的局限性，即通常仅适用于固定且较少的插值点数量。
• 分析微分公式中权重系数的谱特性，以评估在不同空间频率下的精度。
• 通过计算 sin(x) 的导数并考察当点数趋于无穷时的收敛性，验证该方法。
• 为任意数量的插值点提供显式、闭式表达的权重系数公式，从而在数值模拟中实现高效实现。

提出的方法

• 该方法利用函数值在对径点上的对称差分，以消除偶次幂多项式系数，从而分离出奇次幂项用于导数近似。
• 通过差分 $ d_m = f_m - f_{-m} $ 构造奇次幂 $ m $ 的和，建立关于奇次幂多项式系数的线性方程组。
• 利用克莱姆法则求解该方程组，得到权重系数 $ c_{2k+1} $ 的显式表达式，其以 $ \alpha_m^{(2k+1)}(n) $ 表示。
• 通过范德蒙德型矩阵的行列式计算一阶和二阶导数的权重系数 $ \alpha_m^{(1)}(n) $ 和 $ \alpha_m^{(2)}(n) $。
• 通过计算离散傅里叶变换分析权重系数序列的谱特性，并与理想导数谱进行比较。
• 在 MATLAB 中实现该方法，采用一种通过递归乘积迭代计算系数的算法，避免完整矩阵求逆。

实验结果

研究问题

• RQ1 能否在不构建插值多项式的情况下，为任意数量的插值点推导出数值微分的显式闭式公式？
• RQ2 权重系数序列的谱特性如何影响数值微分在不同空间频率下的精度？
• RQ3 当插值点数增加时，数值导数的收敛行为如何，特别是在 $ n \to \infty $ 的极限下？
• RQ4 对于光滑函数如 $ \sin x $，所推导公式的精度与解析导数相比如何？
• RQ5 该方法能否在高精度格点模拟中实现高效且可扩展的应用，例如在量子场论的模拟中？

主要发现

• 所推导的一阶和二阶导数公式对任意数量 $ n $ 的插值点均以闭式表达，无需插值多项式。
• 谱分析表明，一阶导数的权重系数序列在广泛的空间频率范围内近似于理想导数谱 $ (2\pi r/N)^2 $，尤其在 $ n $ 增大时表现更优。
• 当 $ n=11 $ 时，一阶导数近似在空间频率 $ r \lesssim 750 $ 内保持高精度；当 $ n=21 $ 时，精度可扩展至 $ r \lesssim 800 $，表明具有良好的高频精度。
• 该方法在 $ n \to \infty $ 的极限下收敛至 $ \sin x $ 的精确导数，确认了理论一致性。
• 二阶导数公式同样表现出高精度，其谱偏离理想抛物响应的范围在宽频带内保持较小。
• MATLAB 实现通过递归乘积高效计算权重系数，使该方法在大规模模拟中具备实际应用价值。