c语言中使用fft函数,DSP TMS320C6678中FFT函数的调用与具体使用方法（CCS中操作）...

首先明确一点，C6678中FFT的函数形式为：

N ： length of FFT in complex samples

ptr_x ： pointer to complex data input

ptr_w ： pointer to complex twiddle factor

ptr_y ： pointer to complex output data

brev ： pointer to bit reverse table containing 64 entries

n_min ： should be 4 if N can be represented as Power of 4 else, n_min should be 2

offset ： index in complex samples of sub-fft from start of main fft

n_max ： size of main fft in complex samples(可令n_max=N)

其中，64位bit反位表brev为：

unsigned char brev[64] =

{

0x00, 0x20, 0x10, 0x30, 0x08, 0x28, 0x18, 0x38,

0x04, 0x24, 0x14, 0x34, 0x0c, 0x2c, 0x1c, 0x3c,

0x02, 0x22, 0x12, 0x32, 0x0a, 0x2a, 0x1a, 0x3a,

0x06, 0x26, 0x16, 0x36, 0x0e, 0x2e, 0x1e, 0x3e,

0x01, 0x21, 0x11, 0x31, 0x09, 0x29, 0x19, 0x39,

0x05, 0x25, 0x15, 0x35, 0x0d, 0x2d, 0x1d, 0x3d,

0x03, 0x23, 0x13, 0x33, 0x0b, 0x2b, 0x1b, 0x3b,

0x07, 0x27, 0x17, 0x37, 0x0f, 0x2f, 0x1f, 0x3f

};

旋转矩阵(旋转因子)ptr_w为：

void tw_fft_gen (float *w, int n)

{

int i, j, k;

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;

}

}

}

DSPF_sp_fftSPxSP函数的调用方式为：

1.下载库

将下载好的库“dsplib_c66x_3_1_0_0”复制到CCS的安装路径下

2.添加头文件

右键点击工程–>Properties–>Build–>c6000 Compiler–>Include Options

3.添加库文件和库文件路径

右键点击工程–>Properties–>Build–>C6000 Linker–>File Search Path添加库文件和库文件路径

完成以上的准备工作后便可进行DSPF_sp_fftSPxSP函数的调用了。

参考博文：http://bbs.elecfans.com/jishu_1602988_1_1.html

文章中部分内容来自网络，如有侵权，请联系作者，立即删除！！！