C++ 如何为 CORDIC 生成 Kvalues 序列

我在CORDIC上进行了网络搜索，并查看了不同的页面，try 了许多变体[没有显示出来].

从这里:https://nonagon.org/ExLibris/cordic-c-javascript我得到了公式:

prev = prev / sqrt(1 + (1.0 / pow(2, 2.0 * i)))

归根结底，对于您的代码来说，它是更改:

pow(2, -i)

进入:

pow(2, -2.0 * i)

以下是修改后的有效代码:

#include <stdio.h>
#include <math.h>

#define PRECISION   6

#define CALL(_fnc) \
    do { \
        printf("\n" #_fnc "\n"); \
        _fnc(); \
    } while (0)

void
orig(void)
{
    double prev = 1.0;

    for (int i = 0; i < PRECISION; i++)
        printf("%f\n", prev = prev / sqrt(1 + pow(2, -i)));
}

void
fix1(void)
{
    float prev = 1;

    for (int i = 0; i < PRECISION; i++)
        printf("%f\n",prev = prev / sqrt(1 + pow(2, -2.0 * i)));
}

int
main(void)
{

    CALL(orig);
    CALL(fix1);

    return 0;
}

以下是程序输出:

orig
0.707107
0.577350
0.516398
0.486864
0.472328
0.465116

fix1
0.707107
0.632456
0.613572
0.608834
0.607648
0.607352

我建议使用double而不是float，并提高printf格式的精度(例如%.14f):

#include <stdio.h>
#include <math.h>

#define PRECISION   6

#define CALL(_fnc) \
    do { \
        printf("\n" #_fnc "\n"); \
        _fnc(); \
    } while (0)

void
orig(void)
{
    double prev = 1.0;

    for (int i = 0; i < PRECISION; i++)
        printf("%.14f\n", prev = prev / sqrt(1 + pow(2, -i)));
}

void
fix1(void)
{
    double prev = 1;

    for (int i = 0; i < PRECISION; i++)
        printf("%.14f\n",prev = prev / sqrt(1 + pow(2, -2.0 * i)));
}

int
main(void)
{

    CALL(orig);
    CALL(fix1);

    return 0;
}

以下是增强后的输出:

orig
0.70710678118655
0.57735026918963
0.51639777949432
0.48686449556015
0.47232793895472
0.46511640189705

fix1
0.70710678118655
0.63245553203368
0.61357199107790
0.60883391251775
0.60764825625617
0.60735177014130