|
FFT官方库的使用准备工作 其中，stm32_dsp.h和table_fft.h两个文件是必须添加的。stm32_dsp.h是STM32的DSP库的头文件。 进行256点的FFT，只需要调用STM32 DSP库函数中的cr4_fft_256_stm32()函数即可。该函数的原型为： void cr4_fft_256_stm32(void *pssOUT, void *pssIN, uint16_t Nbin);
需要说明的是：按照FFT官方库的说明，pssOUT和pssIN都必须是32位的数据类型，其中高16位存储实部，低16位存储虚部。对于pssIN来说，低16位存储的虚部总是为0。 代码实例调用FFT官方库的代码为： //填充数组
for(i=0;i<NPT;i++)
//这里因为单片机的ADC只能测正的电压 所以需要前级加直流偏执
//加入直流偏执后，需要在软件上减去2048即一半，达到负半周期测量的目的（需要根据具体情况来进行配置）
lBufInArray[i] = ((signed short)(adc_buf[i]-2048)) << 16;
cr4_fft_256_stm32(lBufOutArray, lBufInArray, NPT);
同时，计算各次谐波幅值的函数为：
void GetPowerMag()
{
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);
//除以32768再乘65536是为了符合浮点数计算规律
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);
}
}
| 
楼主
标题置顶
标题高亮
