CORDIC 算法中的每一步都添加cos(arctan(2^-i))
(或1/sqrt(1+2^-2i)
) 的缩放,因此对于 4 步 CORDIC,总缩放为:
cos(arctan(2^-0))*cos(arctan(2^-1))*cos(arctan(2^-2))*cos(arctan(2^-3)) = 0.60883
如果添加更多迭代,它会达到 0.607252935 和一些。至于如何处理该因素,取决于您和您实现的功能。您可能想要删除它并给出未缩放的结果,或者您的消费者可能能够更轻松地处理它。
在 的情况下tanh
,我认为存在具有相同缩放比例的两个值的除法,因此在这种情况下会抵消。我真的不太确定那个,因为我没有在这种模式下使用 CORDIC,所以你必须自己弄清楚。
更新
CORDIC 的优点是每次迭代都是一个 shift+add onX
和Y
。在您的示例中,您合并了比例因子,它给出了正确的答案,但不是必需的,并且不允许简单的 shift + add 操作。因此,我们将缩放保留为后期操作(或初始化),或者如果后续模块可以处理缩放答案,则完全退出。
获得 20 的 sin/cos 的方法是,从X = 1; Y = 0; Z = 20
:
Z0 = 20 - 45.00 = -25.00; X0 = 1 - 0 = 1; Y0 = 0 + 2^-0 = 1
Z1 = Z0 + 26.57 = 1.57; X1 = X0 + Y0*2^-1 = 1.5000; Y1 = Y0 - X0*2^-1 = 0.5000
Z2 = Z1 - 14.04 = -12.47; X2 = X1 - Y1*2^-2 = 1.3750; Y2 = Y1 + X1*2^-2 = 0.8750
Z3 = Z2 + 7.13 = - 5.35; X3 = X2 + Y2*2^-3 = 1.4844; Y3 = Y2 - X2*2^-3 = 0.7031
Z4 = Z3 + 3.58 = - 1.77; X4 = X3 + Y3*2^-4 = 1.5283; Y4 = Y3 - X3*2^-4 = 0.6104
Z5 = Z4 + 1.79 = 0.02; X5 = X4 + Y4*2^-5 = 1.5474; Y5 = Y4 - X4*2^-5 = 0.5626
Z6 = Z5 - 0.90 = - 0.88; X6 = X5 - Y5*2^-6 = 1.5386; Y6 = Y5 + X5*2^-6 = 0.5868
最终答案按比例缩放:cos(20) = X6*0.607 = 0.9339
和sin(20) = Y6*0.607 = 0.3562
。忽略缩放因子的好处应该是很明显的,每次迭代都是Xi = Xi-1 +- Yi-1 * 2^-i
,这是一个移位和加法。您也不需要将比例因子存储在 ROM 或其他地方;您需要的唯一内存是 arctan(2^-i)。
请注意,您可以通过以X = 0.607
代替开头来获得未缩放的结果X = 1
。