本人所做的一些有关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)
|
QQ好友和群
QQ空间
腾讯微博
腾讯朋友
收藏
淘帖
顶
踩
