嵌入式开发者社区

标题: 6678 DSPLIB 库函数的使用 [打印本页]

作者: 戴亮军 时间: 2016-11-8 15:02
标题: 6678 DSPLIB 库函数的使用
现在想利用库里面提供的DSPF_sp_fftSPxSP函数做FFT，再使用这个函数前需要做两个准备，第一个是定义位反转比特表：
unsigned char brev[64] = {
0x0, 0x20, 0x10, 0x30, 0x8, 0x28, 0x18, 0x38,
0x4, 0x24, 0x14, 0x34, 0xc, 0x2c, 0x1c, 0x3c,
0x2, 0x22, 0x12, 0x32, 0xa, 0x2a, 0x1a, 0x3a,
0x6, 0x26, 0x16, 0x36, 0xe, 0x2e, 0x1e, 0x3e,
0x1, 0x21, 0x11, 0x31, 0x9, 0x29, 0x19, 0x39,
0x5, 0x25, 0x15, 0x35, 0xd, 0x2d, 0x1d, 0x3d,
0x3, 0x23, 0x13, 0x33, 0xb, 0x2b, 0x1b, 0x3b,
0x7, 0x27, 0x17, 0x37, 0xf, 0x2f, 0x1f, 0x3f
};

第二个是声明旋转因子函数：
void tw_gen (float *w, int n)
{
int i, j, k;
// const double PI = 3.141592654;

for (j = 1, k = 0; j <= n >> 2; j = j << 2)
{
      for (i = 0; i < n >> 2; i += j)
      {
#ifdef _LITTLE_ENDIAN
         w[k]    = (float) sin (2 * PI * i / n);
         w[k + 1] = (float) cos (2 * PI * i / n);
         w[k + 2] = (float) sin (4 * PI * i / n);
         w[k + 3] = (float) cos (4 * PI * i / n);
         w[k + 4] = (float) sin (6 * PI * i / n);
         w[k + 5] = (float) cos (6 * PI * i / n);

#else
         w[k]    = (float)  cos (2 * PI * i / n);
         w[k + 1] = (float) -sin (2 * PI * i / n);
         w[k + 2] = (float)  cos (4 * PI * i / n);
         w[k + 3] = (float) -sin (4 * PI * i / n);
         w[k + 4] = (float)  cos (6 * PI * i / n);
         w[k + 5] = (float) -sin (6 * PI * i / n);
#endif
         k += 6;
      }
}
}

现在发现一个问题：在调用w_gen (float *w, int n)这个函数时发现，参数n的取值最大只能是1024，否则这个函数调用就不成功，请问有人知道为什么吗？

作者: IRONMAN 时间: 2016-11-8 18:28
不应该随着点数的不同而不同吗？

作者: 希望缄默 时间: 2016-11-9 08:59
这是 TI DSPLIB 的 Bug
实数 FFT 是没有问题的

作者: 小灰灰 时间: 2016-12-10 22:05
可以呀！我做的复数FFT正确呀！

欢迎光临嵌入式开发者社区 (https://51ele.net/)