标题: 本人所做的一些有关FFT的实验STM32源码 [打印本页]
作者: zhuyun    时间: 2018-3-31 11:02
标题: 本人所做的一些有关FFT的实验STM32源码
本人所做的一些有关FFT的实验

单片机源程序如下:
  1. #include "fft.h"
  2. #include "math.h"


  5. u32 lBufInArray[NPT];
  6. u32 lBufOutArray[NPT];
  7. u32 lBufMagArray[NPT];
  8. float PI2=6.28318530717959;
  9. u32 Fs=192000;
  10. /******************************************************************
  11. 函数名称:InitBufInArray()
  12. 函数功能:模拟采样数据，采样数据中包含3种频率正弦波
  13. 参数说明:
  14. 备    注:在lBufInArray数组中，每个数据的高16位存储采样数据的实部，
  15.           低16位存储采样数据的虚部(总是为0)
  16. *******************************************************************/
  17. void InitBufInArray(void)
  18. {
  19.     unsigned short i;
  20.     float fx;
  21.     for(i=0; i<NPT; i++)
  22.     {
  23.         fx = 1000 * sin(PI2 * i * 1000.0 / Fs) +
  24.              3000 * sin(PI2 * i * 3000.0 / Fs) +
  25.              4000 * sin(PI2 * i * 5000.0 / Fs);
  26.         lBufInArray[i] = ((signed short)fx) << 16;
  27.     }
  28. }
  29. /******************************************************************
  30. 函数名称:GetPowerMag()
  31. 函数功能:计算各次谐波幅值
  32. 参数说明:
  33. 备　　注:先将lBufOutArray分解成实部(X)和虚部(Y)，然后计算幅值(sqrt(X*X+Y*Y)
  34. *******************************************************************/
  35. void GetPowerMag(void)
  36. {
  37.     signed short lX,lY;
  38.     float X,Y,Mag;
  39.     unsigned short i;
  40.     for(i=0; i<NPT/2; i++)
  41.     {
  42.         lX  = (lBufOutArray[i] << 16) >> 16;
  43.         lY  = (lBufOutArray[i] >> 16);
  44.         X = NPT * ((float)lX) / 32768;
  45.         Y = NPT * ((float)lY) / 32768;
  46.         Mag = sqrt(X * X + Y * Y) / NPT;
  47.         if(i == 0)
  48.             lBufMagArray[i] = (unsigned long)(Mag * 32768);
  49.         else
  50.             lBufMagArray[i] = (unsigned long)(Mag * 65536);
  51.     }
  52. }

复制代码

所有资料51hei提供下载:
FFT.rar (414.22 KB, 下载次数: 47)






欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1