STM32F1 FFT库的使用
ST公司为了方便客户使用FFT,自己做了一个库,不过这个库是有限制的。点数必须是4的次方,分别是64、256和1024个点。速度完全满足客户的要求。
1、第一步必须添加使用FFT的库文件到inc和src中。附上百度云网盘链接(http://pan.baidu.com/s/1gfHkS0b)
2、导入相应的头文件啊。。。
3、导入相应的文件之后嘛,就是直接调用库函数这东东了。不过前面还得先定义一下参数。
参数定义呢就像上面那张图片。lBUFIN1呢就是你需要进行FFT的数组哦,后面的N就是你需要变换的点数啦。
lBUFOUT1呢就是FFT变换 之后的数据啦。在这个参数数组呢高16位呢就是你的实部,低16位呢就是你的虚部。哈哈,直接贴代码可能会直观点。
long lBUFIN1[N]; /* Complex input vector */long lBUFOUT1[N]; /* Complex output vector */#define N 512void dsp_asm_powerMag(void) { s16 X1,Y1; u32 i; Xmax = Ymax = 0; cr4_fft_256_stm32(lBUFOUT1, lBUFIN1, N); for(i=0;i < N/2;i++) { X1 = (lBUFOUT1[i] << 16) >> 16; Y1 = (lBUFOUT1[i] >> 16); { float X = N * ((float)X1) /32768; float Y = N * ((float)Y1) /32768; float Mag = sqrt(X*X + Y*Y)/N; lBUFMAG1[i] = (u32)(Mag * 65536); } } }
因为FFT之后只有一半的数据是可以使用的,而且他们是对称关系的。所以就N/2。