[论文解读] Fast exact digital differential analyzer for circle generation

Jan L. Cieśliński, Leonid Moroz|arXiv (Cornell University)|Apr 17, 2013

Advanced Electrical Measurement Techniques被引用 1

一句话总结

本文提出一种新颖的两步数字微分分析器（DDA）算法，利用显式中点法则实现精确的圆生成。通过利用移位和加法操作，该方法在保持计算效率的同时实现了最大精度——精确到舍入误差范围内，其精度和成本均优于传统的单步DDA方法。

ABSTRACT

In the first part of the paper we present a short review of applications of digital differential analyzers (DDA) to generation of circles showing that they can be treated as one-step numerical schemes. In the second part we present and discuss a novel fast algorithm based on a two-step numerical scheme (explicit midpoint rule). Although our algorithm is as cheap as the simplest one-step DDA algoritm (and can be represented in terms of shifts and additions), it generates circles with maximal accuracy, i.e., it is exact up to round-off errors.

研究动机与目标

• 开发一种计算高效且高精度的数字圆生成方法，适用于计算机图形学和数值控制领域。
• 克服传统单步DDA方法固有的不准确性，这些方法由于离散化误差会产生对数螺旋或椭圆螺旋。
• 设计一种两步数值格式，使轨迹在舍入误差范围内完全保持圆形，同时适用于实时应用。
• 通过引入周期保持参数δ(h)的多项式逼近，探索计算成本与精度之间的权衡。

提出的方法

• 将圆生成建模为连续系统 dx/dϑ = -y, dy/dϑ = x 的两步显式中点法则离散化。
• 推导出递推关系：xn+2 = xn - 2δyn+1, yn+2 = yn + 2δxn+1，其中δ是步长h的函数。
• 证明当δ(h) = sin(h)时，该格式在圆x² + y² = r²上精确生成点，同时保持轨迹和周期。
• 引入sin(h)的多项式逼近，如δ = h - h³/6，以用位移和加法替代昂贵的乘法运算。
• 通过证明所有特征值位于单位圆上，分析稳定性，确保扰动有界且长期数值稳定。
• 通过符号展开和渐近展开验证方法，包括当h = 2⁻ᵐ时√(1 - h²)的级数展开，从而支持高效的定点算术。

实验结果

研究问题

• RQ1 能否设计一种两步数值格式，在保持低计算成本的同时实现圆的精确生成（舍入误差范围内）？
• RQ2 在离散DDA格式中，如何在不使用昂贵三角函数运算的情况下保持圆周运动的周期？
• RQ3 在定点算术中，何种多项式逼近sin(h)能在精度与效率之间实现最佳平衡？
• RQ4 为何基于Matsushiro的三阶DDA（a = 1 - h²/2, c = h - h³/8）在所有三阶多项式DDA方案中精度最高？
• RQ5 在圆生成算法中，使用位移和加法在多大程度上可替代乘法而不损失精度？

主要发现

• 所提出的基于显式中点法则的两步DDA算法，当δ(h) = sin(h)时，如定理3.6所证明，可精确生成圆（舍入误差范围内）。
• 当δ(h) = sin(h)时，公式(48)保持了圆形轨迹和周期T = 2πh / arcsin(h)，确保了正确的角速度。
• 通过用δ = h - h³/6逼近sin(h)，该方法实现了三阶精度，并如公式(51)所示，用位移替代了昂贵的乘法运算。
• 采用δ = h - h³/6的算法计算效率高，仅需移位和加法即可实现，适用于嵌入式和实时系统。
• 系统矩阵的所有特征值位于单位圆上，表明初始扰动不会增长，从而保证了长期稳定性。
• 该方法在精度上优于所有单步DDA方案，包括精度最高的三阶DDA（Matsushiro的方案），后者仍会产生k ≈ h⁵/128的螺旋偏差。