c - 在不使用标准功能的情况下，在C中生成正弦信号

我想在C中生成正弦信号而不使用标准函数sin()来触发LED亮度的正弦形变化。我的基本想法是使用40点和插值的查找表。

这是我的第一种方法:

const int sine_table[40] = {0, 5125, 10125, 14876, 19260, 23170, 26509, 29196,
31163, 32364, 32767,  32364, 31163, 29196, 26509, 23170, 19260, 14876, 10125,
5125, 0, -5126, -10126,-14877, -19261, -23171, -26510, -29197, -31164, -32365,
-32768, -32365, -31164, -29197, -26510, -23171, -19261, -14877, -10126, -5126};

int i = 0;
int x1 = 0;
int x2 = 0;
float y = 0;

float sin1(float phase)
{
    x1 = (int) phase % 41;
    x2 = x1 + 1;
    y = (sine_table[x2] - sine_table[x1])*((float) ((int) (40*0.001*i*100) % 4100)/100 - x1) + sine_table[x1];
    return y;
}

int main()
{
    while(1)
    {
    printf("%f      ", sin1(40*0.001*i)/32768);
    i = i + 1;
    }
}

不幸的是，该函数有时返回的值远大于1。此外，插值似乎不太好(我用它来创建LED的正弦形亮度变化，但这些变化非常不平滑)。

有谁有更好的主意在C中实现正弦发生器？

最佳答案

OP的主要问题是为表查找生成索引。

OP的代码尝试访问导致undefined behavior的外部数组sine_table[40]。至少解决该问题。

const int sine_table[40] = {0, 5125, 10125, ...
    ...
    x1 = (int) phase % 41;                     // -40 <= x1 <= 40
    x2 = x1 + 1;                               // -39 <= x2 <= 41
    y = (sine_table[x2] - sine_table[x1])*...  // bad code, consider x1 = 40 or x2 = 40,41

建议变更

    x1 = (int) phase % 40;   // mod 40, not 41
    if (x1 < 0) x1 += 40;    // Handle negative values
    x2 = (x1 + 1) % 40;      // Handle wrap-around
    y = (sine_table[x2] - sine_table[x1])*...

存在更好的方法，但下面将重点介绍OP的方法。

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

const int sine_table[40] = { 0, 5125, 10125, 14876, 19260, 23170, 26509, 29196,
31163, 32364, 32767, 32364, 31163, 29196, 26509, 23170, 19260, 14876, 10125,
5125, 0, -5126, -10126, -14877, -19261, -23171, -26510, -29197, -31164, -32365,
-32768, -32365, -31164, -29197, -26510, -23171, -19261, -14877, -10126, -5126 };

int i = 0;
int x1 = 0;
int x2 = 0;
float y = 0;

float sin1(float phase) {
  x1 = (int) phase % 40;
  if (x1 < 0) x1 += 40;
  x2 = (x1 + 1) % 40;
  y = (sine_table[x2] - sine_table[x1])
      * ((float) ((int) (40 * 0.001 * i * 100) % 4100) / 100 - x1)
      + sine_table[x1];
  return y;
}

int main(void) {
  double pi = 3.1415926535897932384626433832795;
  for (int j = 0; j < 1000; j++) {
    float x = 40 * 0.001 * i;
    float radians = x * 2 * pi / 40;
    printf("%f %f %f\n", x, sin1(x) / 32768, sin(radians));
    i = i + 1;
  }
}

输出量

         OP's     Reference sin()
0.000000 0.000000 0.000000
0.040000 0.006256 0.006283
0.080000 0.012512 0.012566
...
1.960000 0.301361 0.303035
2.000000 0.308990 0.309017
2.040000 0.314790 0.314987
...
39.880001 -0.020336 -0.018848
39.919998 -0.014079 -0.012567
39.959999 -0.006257 -0.006283

更好的代码不会将i, x1, x2, y值作为全局变量传递，而是作为函数参数或函数变量传递。也许那是OP调试的产物。

这是相当广泛的。在速度，精度，代码空间，可移植性或可维护性方面更好？ sine()函数很容易制作。高质量的需要更多的努力。

尽管模糊，但OP使用较小的查找表是一个不错的开始-尽管我认为无需任何浮点运算就可以完成操作。我建议OP构建一个经过测试的有效解决方案，并将其发布在Code Review中以提出改进意见。

关于c - 在不使用标准功能的情况下，在C中生成正弦信号，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/47906401/