[论文解读] A Dual Number Approach for Numerical Calculation of derivatives and its use in the Spherical 4R Mechanism

F. Pe, R. Pe+2 more|arXiv (Cornell University)|Jan 1, 2013

Robotic Mechanisms and Dynamics参考文献 23被引用 3

一句话总结

本文提出了一种基于双重数的高效方法，可在单步计算中同时求解向量函数的一阶和二阶导数，利用嵌套双重数算术。该方法在 Fortran 中实现，可在运动学分析中实现精确且自动化的微分，通过在球面 4R 机构上的应用得到验证，显著提升了计算效率和导数精度。

ABSTRACT

This paper proposes a methodology to calculate both the rst and second derivatives of a vector function of one variable in a single computation step. The method is based on the nested application of the dual number approach for rst order derivatives. It has been implemented in Fortran language, a module which contains the dual version of elementary functions as well as more complex functions, which are common in the

研究动机与目标

• 开发一种统一的计算方法，可在一次操作中计算向量函数的一阶和二阶导数。
• 将常用于一阶导数的双重数方法扩展至高阶导数。
• 在 Fortran 中实现该方法，并提供基本函数和复杂数学函数的双重数版本库。
• 将该方法应用于球面 4R 机构的运动学分析，该机构为复杂的空间机构，需要高阶导数。
• 提升机制设计与分析中数值导数计算的计算效率和精度。

提出的方法

• 该方法采用嵌套双重数，每个双重数同时携带函数值和最高至二阶的导数信息。
• 一阶导数通过标准双重数算术计算；二阶导数通过将双重数算术递归应用于一阶导数分量获得。
• 通过自定义 Fortran 模块实现基本函数（如 sin、cos、exp）和运动学中使用的复合函数的双重数版本。
• 使用双重数算术对感兴趣的向量函数进行求值，从而同时获得一阶和二阶导数。
• 该方法避免了符号微分，并相比有限差分法或分步执行的算法微分，显著降低了计算开销。
• 通过在球面 4R 机构上的应用对该方法进行验证，高阶导数在动态分析和敏感性分析中至关重要。

实验结果

研究问题

• RQ1 双重数算术能否扩展至在单次计算遍历中计算向量函数的二阶导数？
• RQ2 与传统的有限差分法或分步计算一阶/二阶导数相比，嵌套双重数方法在精度和效率方面表现如何？
• RQ3 该方法对球面 4R 等复杂机构的运动学分析有何实际影响？
• RQ4 能否有效实现一个稳健的 Fortran 双重数函数库以用于工程应用？
• RQ5 该方法在机制设计中基于导数的计算可靠性方面提升了多少？

主要发现

• 该方法成功通过嵌套双重数算术在单次求值中同时计算出向量函数的一阶和二阶导数。
• 在 Fortran 中的实现展示了计算效率，避免了多次函数调用或符号运算的需要。
• 双重数库能够对运动学中常用的初等函数和复合函数实现精确的导数计算。
• 该方法在球面 4R 机构上得到验证，高阶导数在动态分析和敏感性分析中至关重要。
• 与有限差分近似相比，该方法具有更高的精度，并且相比分步导数计算，显著降低了计算开销。
• 结果证实了双重数方法在机械系统分析中用于高阶导数的可行性与有效性。