本人所做的一些有關(guān)FFT的實驗STM32源碼

ID:299814 · 發(fā)表于 2018-3-31 11:02

本人所做的一些有關(guān)FFT的實驗

單片機源程序如下:

#include "fft.h"
#include "math.h"
u32 lBufInArray[NPT];
u32 lBufOutArray[NPT];
u32 lBufMagArray[NPT];
float PI2=6.28318530717959;
u32 Fs=192000;
/******************************************************************
函數(shù)名稱:InitBufInArray()
函數(shù)功能:模擬采樣數(shù)據(jù)，采樣數(shù)據(jù)中包含3種頻率正弦波
參數(shù)說明:
備注:在lBufInArray數(shù)組中，每個數(shù)據(jù)的高16位存儲采樣數(shù)據(jù)的實部，
低16位存儲采樣數(shù)據(jù)的虛部(總是為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;
}
}
/******************************************************************
函數(shù)名稱:GetPowerMag()
函數(shù)功能:計算各次諧波幅值
參數(shù)說明:
備　　注:先將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);
}
}

復(fù)制代碼

所有資料51hei提供下載:

FFT.rar (414.22 KB, 下載次數(shù): 47)

帳號		自動登錄	找回密碼
密碼			立即注冊

久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

本人所做的一些有關(guān)FFT的實驗STM32源碼

評分