欧美极品高清xxxxhd,国产日产欧美最新,无码AV国产东京热AV无码,国产精品人与动性XXX,国产传媒亚洲综合一区二区,四库影院永久国产精品,毛片免费免费高清视频,福利所导航夜趣136
標題:
本人所做的一些有關FFT的實驗STM32源碼
[打印本頁]
作者:
zhuyun
時間:
2018-3-31 11:02
標題:
本人所做的一些有關FFT的實驗STM32源碼
本人所做的一些有關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 上傳
點擊文件名下載附件
下載積分: 黑幣 -5
歡迎光臨 (http://m.raoushi.com/bbs/)
Powered by Discuz! X3.1