本人所做的一些有關FFT的實驗
單片機源程序如下:
- #include "fft.h"
- #include "math.h"
- u32 lBufInArray[NPT];
- u32 lBufOutArray[NPT];
- u32 lBufMagArray[NPT];
- float PI2=6.28318530717959;
- u32 Fs=192000;
- /******************************************************************
- 函數名稱:InitBufInArray()
- 函數功能:模擬采樣數據，采樣數據中包含3種頻率正弦波
- 參數說明:
- 備 注:在lBufInArray數組中，每個數據的高16位存儲采樣數據的實部，
- 低16位存儲采樣數據的虛部(總是為0)
- *******************************************************************/
- void InitBufInArray(void)
- {
- unsigned short i;
- float fx;
- for(i=0; i<NPT; i++)
- {
- fx = 1000 * sin(PI2 * i * 1000.0 / Fs) +
- 3000 * sin(PI2 * i * 3000.0 / Fs) +
- 4000 * sin(PI2 * i * 5000.0 / Fs);
- lBufInArray[i] = ((signed short)fx) << 16;
- }
- }
- /******************************************************************
- 函數名稱:GetPowerMag()
- 函數功能:計算各次諧波幅值
- 參數說明:
- 備 注:先將lBufOutArray分解成實部(X)和虛部(Y)，然后計算幅值(sqrt(X*X+Y*Y)
- *******************************************************************/
- void GetPowerMag(void)
- {
- signed short lX,lY;
- float X,Y,Mag;
- unsigned short i;
- for(i=0; i<NPT/2; i++)
- {
- lX = (lBufOutArray[i] << 16) >> 16;
- lY = (lBufOutArray[i] >> 16);
- X = NPT * ((float)lX) / 32768;
- Y = NPT * ((float)lY) / 32768;
- Mag = sqrt(X * X + Y * Y) / NPT;
- if(i == 0)
- lBufMagArray[i] = (unsigned long)(Mag * 32768);
- else
- lBufMagArray[i] = (unsigned long)(Mag * 65536);
- }
- }
所有資料51hei提供下載:
FFT.rar
(414.22 KB, 下載次數: 47)
2018-4-1 22:57 上傳
點擊文件名下載附件
|
