STM32 官方DSP的 FFT库使用

最早我实现FFT是移植的网上搜罗来的FFT 的 C程序代码，奈何运算效率太低，只好放弃

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include "FFT.h"
#include "sample.h"
#define N DEPTH_ADC_MAX



complex x[N + 10];// *W; /*输入序列,变换核*/
//complex W[N];
int size_x=0;      /*输入序列的大小，在本程序中仅限2的次幂*/
double PI = 3.1415926535;         /*圆周率*/
void fft(void);     /*快速傅里叶变换*/

void change(void); /*变址*/
void add(complex,complex,complex *);   /*复数加法*/
void mul(complex,complex,complex *);   /*复数乘法*/
void sub(complex,complex,complex *);   /*复数减法*/



complex fft_result(uint16_t index )
{
   
   
    if(index >= N)
        index = 0;
    return x[index];
}

int fft_analyze(uint16_t * sample_data,uint16_t sample_cnt)
{
   
   
    int i;                             /*输出结果*/

    //PI=atan(1)*4;
    PI = 3.1415926535;

    if(sample_cnt > N)
        return -1;

    size_x = sample_cnt;



    for(i=0; i<size_x; i++)
    {
   
   
        x[i].real = (double)(*( sample_data + i));
        x[i].img = 0;
    }


//	initW();//调用变换核
    fft();//调用快速傅里叶变换

    return 0;
}



/*快速傅里叶变换*/

void fft()
{
   
   
    int i=0,j=0,k=0,l=0;
    complex W1;
    int W_Index = 0;
    complex up,down,product;


    change();  //调用变址函数
    for(i=0; i< log(size_x)/log(2) ; i++)      /*一级蝶形运算 stage */
    {
   
   
        l=1<<i;
        for(j=0; j<size_x; j+= 2<