ST公司為了方便客戶使用FFT,自己做了一個庫,不過這個庫是有限制的。點數必須是4的次方,分別是64、256和1024個點。速度完全滿足客戶的要求。 1、第一步必須添加使用FFT的庫文件到inc和src中。附上百度雲網盤鏈接(http://pan.baidu.com/s/1gfHkS0b) 2、導入 ...
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 512
void 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。