【安富莱——DSP教程】第31章复数FFT的逆变换实现

ti

第31章复数FFT的逆变换实现

本章主要讲解复数FFT的逆变换。

本章节使用的复数FFT函数来自ARM官方库的TransformFunctions部分

31.1 复数FFT的逆变换实现

31.2 总结

硬汉Eric2013 · 2015-7-4 14:16:59

31.1 复数FFT的逆变换实现

本小节主要讲解复数FFT的逆变换实现，通过函数arm_cfft_f32实现浮点数的逆变换。

31.1.1 arm_cfft_f32逆变换

函数定义如下：

void arm_cfft_f32(

constarm_cfft_instance_f32 * S,

float32_t * p1,

uint8_tifftFlag,

uint8_tbitReverseFlag)

参数定义：

[in] *S points to an instance of the floating-point CFFT structure.

[in, out] *p1 points to the complex data buffer of size2*fftLen. Processing

occurs in-place.

[in]ifftFlag flag that selects forward (ifftFlag=0) orinverse (ifftFlag=1) transform.

[in]bitReverseFlag flag that enables(bitReverseFlag=1) or disables (bitReverseFlag=0) bit

reversal of output.

注意事项：

结构const arm_cfft_instance_f32的定义如下（在文件arm_math.h文件）：

typedefstruct

{

uint16_t fftLen;

constfloat32_t *pTwiddle;

constuint16_t *pBitRevTable;

uint16_t bitRevLength;

}arm_cfft_instance_f32;

下面通过函数arm_cfft_f32计算一个正弦波的FFT，然后再使用函数arm_cfft_f32做FFT逆变换，并使用 Matlab计算变换前后的结果对比。

/*
*********************************************************************************************************
* 函数名: arm_cfft_f32_app
* 功能说明: 调用函数arm_cfft_f32_app计算CFFT正变换和你变换。
* 形参：无
* 返回值: 无
*********************************************************************************************************
*/
static void arm_cfft_f32_app(void)
{
uint16_t i;
fftSize = 1024;
ifftFlag = 0;
doBitReverse = 1;
/* 按照实部，虚部，实部，虚部..... 的顺序存储数据 */
for(i=0; i<1024; i++)
{
testInput_f32_10khz[i*2+1] = 0;
/* 50Hz正弦波，采样率1KHz */
testInput_f32_10khz[i*2] = arm_sin_f32(2*3.1415926f*50*i/1000);
printf("%frn", testInput_f32_10khz[i*2]);
}
/* CFFT变换：FFT正变换 */
arm_cfft_f32(&arm_cfft_sR_f32_len1024, testInput_f32_10khz, ifftFlag, doBitReverse);
ifftFlag = 1;
/* CFFT变换：FFT逆变换 */
arm_cfft_f32(&arm_cfft_sR_f32_len1024, testInput_f32_10khz, ifftFlag, doBitReverse);
printf("****************************分割线***************************************rn");
/* 串口打印求解的模值 */
for(i=0; i<1024; i++)
{
printf("%frn", testInput_f32_10khz[i*2]);
}
}