用STM32的DSP库做的FFT测试，测试的频率没毛病，但是幅值就GG了

只要在路上就行

用STM32的DSP库做的FFT测试，测试的频率没毛病，但是幅值就GG了，有没有做过的啊




#include "stm32f10x.h"
#include "led.h"
#include "delay.h"
#include "USART1.h"
#include "adc.h"
#include "TIME2.h"
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "stm32_dsp.h"
#include "table_fft.h"

#define PI2 6.28318530717959
#define Fs 1024*2                   //采样频率 Hz
#define NPT 1024                    //采样点数

u16 ADC_input[NPT];

u32 input[NPT];
u32 output[NPT],Mag[NPT];

float Fn[NPT/2];                        //各个点的频率
float Mag_max;
u16 ADC_DataNum = 0;
s16 Fn_Num;
u16 ADC_flag = 0;


extern __IO uint16_t ADC_ConvertedValue[];    // ADC1转换的电压值通过MDA方式传到SRAM   
/**********************原程序*********************************************/

/***************计算幅值*****************/
void PowerMag(u16 FFT_NPT)  
{

    s16 lx, ly;
    u32 i;
    float mag, x, y;
   
    for (i = 0; i < FFT_NPT/2; i++ )
    {

        lx = (output[i]<<16)>>16;
        ly = (output[i]>>16);
        
         x = FFT_NPT*((float)lx)/32768;
         y = FFT_NPT*((float)ly)/32768;
        mag = sqrt(x*x+y*y)/FFT_NPT;        
        Mag[i] = (u32)(mag*65536);//交流信号        
    }
    Mag[0] = Mag[0]/2;//直流信号   
}

/***********取最大值************/
void Compare_Max()
{
    u16 i;
    Mag_max = Mag[1];
    for(i=1; i < NPT/2; i++)
    {
        if(Mag_max < Mag[i])
        {
            Fn_Num = i;
            Mag_max = Mag[i];
            
        }
    }

}

/********FFT变换**********/
void DSP_FFT1024()
{
    u16 i;   
   
    for (i = 0; i <NPT ; i++)
    {
         input[i] = ADC_input[i]<<16;
    }
   
    cr4_fft_1024_stm32(output,input,NPT);   //FFT变换
    PowerMag(NPT);

}

int main(void)
{  
    u16 i ;
    delay_init(72);   
    TIME2_Init();
    ADC1_Init();     
    USART1_Configuration() ;
    while (1)
    {  
        if(ADC_flag == 1)
        {
            TIM_Cmd(TIM2,DISABLE);
            DMA_Cmd(DMA1_Channel1,DISABLE);
            ADC_flag = 0;
            DSP_FFT1024();
            for(i = 0; i < NPT/2; i++)
            {
                Fn[i] = ((float)i*Fs/NPT);            
            }

            Compare_Max();
            printf("F = %f,%d,  %f\r\n",Fn[Fn_Num],Mag[Fn_Num],Mag_max);   
            
            DMA_Cmd(DMA1_Channel1,ENABLE);
            TIM_Cmd(TIM2,ENABLE);    //开启定时器2
        }
    }
}
/***********************************************************************************/




void DMA1_Channel1_IRQHandler()  
{  
    int  i;
    if(DMA_GetITStatus(DMA_IT_TC))                      //判断DMA传输完成中断  
    {      
            for(i=0;i<NPT;i++)
            {
                ADC_input[i] = ADC_ConvertedValue[i];
            }
            ADC_flag = 1;
    }
    DMA_ClearITPendingBit(DMA_IT_TC);                   //清除DMA中断标志位  
}  

星空pp

你好,请问可以分享一下dsp库吗，谢谢

十一呀

？？？？