stm32f103 fft变换有老哥懂吗，求助求助！！！

道桦 · 发表于 2018-9-25 10:11:31

在网上看到了@依旧淡然的博客，有关用dsp库进行fft的，写得很好，但是不知道为什么我运算出来的结果跟博主的完全不一样的，还很奇怪
下面附上我的程序以及运算结果，希望有大神可以帮我解答一下
#include "led.h"
#include "delay.h"
#include "key.h"
#include "sys.h"
#include "lcd.h"
#include "usart.h"
#include "stm32_dsp.h"
#include "table_fft.h"
#include "stm32f10x.h"
#include "math.h"
unsigned long lBufOutArray[300];
unsigned long lBufInArray[300];
unsigned short int lBufMagArray[300];
float X[300],Y[300],Mag[300];
unsigned int NPT=256;
/******************************************************************
函数名称:InitBufInArray()
函数功能:模拟采样数据，采样数据中包含3种频率正弦波(350Hz，8400Hz，18725Hz)
参数说明:
备注:在lBufInArray数组中，每个数据的高16位存储采样数据的实部，
      低16位存储采样数据的虚部(总是为0)
*******************************************************************/
void InitBufInArray()
{
  unsigned short i;
float fx;
  int Fs;
  Fs=44800;
for(i=0; i<NPT; i++)
{
      fx = 1500 * sin(6.28318530717959* i * 350.0 / Fs) +2700 * sin(6.28318530717959* i * 8400.0 / Fs) +4000 * sin(6.28318530717959* i * 18725.0 / Fs);
lBufInArray[i] = ((signed short)fx) << 16;
}
}
/******************************************************************
函数名称:dsp_fft()
函数功能:调用DSP库函数进行FFT
参数说明:其中，参数pssOUT表示FFT输出数组指针，参数pssIN表示要进行FFT运算的输入数组指针，参数Nbin表示了点数。
备注:lBufOutArray中每个元素的数据格式为；高16位存储虚部，低16位存储实部。
*******************************************************************/
void dsp_fft()
{
cr4_fft_256_stm32(&lBufOutArray, &lBufInArray, NPT);
}
/******************************************************************
函数名称:GetPowerMag()
函数功能:计算各次谐波幅值
参数说明:
备　　注:先将lBufOutArray分解成实部(X)和虚部(Y)，然后计算幅值(sqrt(X*X+Y*Y)
*******************************************************************/
void GetPowerMag()
{
  signed short lX,lY;
unsigned short i;
for(i=0; i<NPT/2; i++)
{
      lX  = (lBufOutArray[i] << 16) >> 16;
      lY  = (lBufOutArray[i] >> 16);
      X[i]= NPT * ((float)lX) / 32768;
      Y[i]= NPT * ((float)lY) / 32768;
  Mag[i]=sqrt(X[i]* X[i]+ Y[i]* Y[i]) / NPT;
      if(i == 0)
         lBufMagArray[i] = (unsigned long)(Mag[i] * 32768);
      else
         lBufMagArray[i] = (unsigned long)(Mag[i] * 65536);
}
}
/******************************************************************
函数名称:int main()
函数功能:将lBufMagArray数组中各次谐波的幅值（即各个频率分量的幅值）输出打印出来
参数说明:
备　　注:
*******************************************************************/
int main(void)
{
int a;
int i[130],p[130];
delay_init();    //延时函数初始化
NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);  //设置NVIC中断分组2:2位抢占优先级，2位响应优先级
uart_init(115200); //串口初始化为115200

InitBufInArray();
cr4_fft_256_stm32(&lBufOutArray[0], &lBufInArray[0], NPT);
GetPowerMag();



  for(a=0;a<=127;a=a+1)
{
  i[a]=a;
}
for(a=0;a<=127;a=a+1)
{
  p[a]=a*175;
}

printf("\ti,\t\tp,\tmag,\t\tx,\t\ty\r\n");
for(a=0;a<=127;a++)
{
  printf("%12d,%12d,%12.7d,%12.7d,%12.7d\r\n",i[a],p[a],Mag[a],X[a],Y[a]);
}

while(1)
{
// printf("i\tp\tmag\tx\ty\n");
// printf("\r\n");
// for(a=0;a<=127;a++)
// {
//  printf("%12d,%12d,%12d,%12d,%12d",i[a],p[a],Mag[a],X[a],Y[a]);
//  printf("\r\n");
// }
}
}