标题: stm32STM32F407_FFTDSP库_ADC+DMA源码 [打印本页]

作者: 鸡米花    时间: 2018-7-21 20:52
标题: stm32STM32F407_FFTDSP库_ADC+DMA源码
stm32f407dma和adc

单片机源程序如下:
  1. #include "sys.h"
  2. #include "delay.h"  
  3. #include "usart.h"  
  4. #include "led.h"
  5. #include "key.h"
  6. #include "lcd.h"
  7. #include "timer.h"
  8. #include "math.h"
  9. #include "arm_math.h"
  10. #include "adc.h"

  11. #define FFT_LENGTH                         1024        //FFT长度,默认是1024点FFT

  12. float fft_inputbuf[FFT_LENGTH*2];        //FFT输入数组
  13. float fft_outputbuf[FFT_LENGTH];        //FFT输出数组

  14. u8 timeout;//定时器溢出次数
  15. extern u16 AD_Value[1024];
  16. extern u16 After_filter[3];
  17. u16 AD_FFT[1024];
  18. int main(void)
  19. {
  20.   arm_cfft_radix4_instance_f32 scfft;
  21.         u8 key,t=0;
  22.         float time;
  23.         u8 buf[50];
  24.         u16 i,tt;

  25.         NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);//设置系统中断优先级分组2
  26.         delay_init(168);  //初始化延时函数
  27.         uart_init(115200);                //初始化串口波特率为115200
  28.        
  29.         LED_Init();                                        //初始化LED
  30.         KEY_Init();                                        //初始化按键
  31.         LCD_Init();                                        //初始化LCD
  32.        
  33.         DMA_Configuration();
  34.         ADC1_Configuration();
  35.        
  36.         TIM3_Int_Init(65535,84-1);        //1Mhz计数频率,最大计时65ms左右超出
  37.           POINT_COLOR=RED;
  38.         LCD_ShowString(30,50,200,16,16,"Explorer STM32F4");       
  39.         LCD_ShowString(30,70,200,16,16,"DSP FFT TEST");       
  40.         LCD_ShowString(30,90,200,16,16,"ATOM@ALIENTEK");
  41.         LCD_ShowString(30,110,200,16,16,"2014/7/2");       
  42.         LCD_ShowString(30,130,200,16,16,"KEY0:Run FFT");//显示提示信息
  43.         LCD_ShowString(30,160,200,16,16,"FFT runtime:");//显示FFT执行时间
  44.         POINT_COLOR=BLUE;        //设置字体为蓝色   
  45.         arm_cfft_radix4_init_f32(&scfft,FFT_LENGTH,0,1);//初始化scfft结构体,设定FFT相关参数
  46.         while(1)
  47.         {
  48.                 key=KEY_Scan(0);
  49.                 if(key==KEY0_PRES)
  50.                 {
  51.                         for(i=0;i<FFT_LENGTH;i++)
  52.                                 AD_FFT[i]=AD_Value[i];
  53.                         for(i=0,tt=0;i<FFT_LENGTH;i++)//来自ADC信号序列
  54.                         {
  55.                                  fft_inputbuf[2*i] = AD_FFT[tt];
  56.                                  fft_inputbuf[2*i+1]=0;//虚部全部为0
  57.                                  tt++;
  58.                         }
  59. //                        for(i=0;i<FFT_LENGTH;i++)//手动生成信号序列
  60. //                        {
  61. //                                 fft_inputbuf[2*i]=100+
  62. //                                                   10*arm_sin_f32(2*PI*i/FFT_LENGTH)+50*arm_cos_f32(2*PI*i*128/FFT_LENGTH)+
  63. //                                                                   30*arm_sin_f32(2*PI*i*4/FFT_LENGTH)+
  64. //                                                   50*arm_cos_f32(2*PI*i*8/FFT_LENGTH);        //生成输入信号实部
  65. //                                 fft_inputbuf[2*i+1]=0;//虚部全部为0
  66. //                        }                       
  67.                         TIM_SetCounter(TIM3,0);//重设TIM3定时器的计数器值
  68.                         timeout=0;
  69.                        
  70.                         arm_cfft_radix4_f32(&scfft,fft_inputbuf);        //FFT计算(基4)
  71.                        
  72.                         time=TIM_GetCounter(TIM3)+(u32)timeout*65536;                         //计算所用时间
  73.                         sprintf((char*)buf,"%0.3fms\r\n",time/1000);       
  74.                         LCD_ShowString(30+12*8,160,100,16,16,buf);        //显示运行时间               
  75.                        
  76.                         arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);        //把运算结果复数求模得幅值
  77.                        
  78.                         printf("\r\n%d point FFT runtime:%0.3fms\r\n",FFT_LENGTH,time/1000);
  79.                         printf("FFT Result:\r\n");
  80.                         for(i=0;i<FFT_LENGTH;i++)
  81.                         {
  82.                                 printf("fft_outputbuf[%d]:%f\r\n",i,fft_outputbuf[i]);
  83.                                 LCD_DrawLine(i, 0,i,fft_outputbuf[i]/10000);
  84.                                 LCD_ShowxNum(30,5+i*20,fft_outputbuf[i],8,16,0);
  85.                         }
  86.                 }else delay_ms(10);
  87.                 t++;
  88.                 if((t%10)==0)LED0=!LED0;                  
  89.         }
  90. }

复制代码

所有资料51hei提供下载:
DSP-FFT-ADC我的(好的).rar (1.09 MB, 下载次数: 221)



作者: z1026563274    时间: 2018-7-21 23:45
请问这个可以测得频率吗 咳咳 老哥
作者: 1245455936    时间: 2018-7-23 00:31
谢谢分享
作者: Jimnice    时间: 2018-7-27 10:28
谢谢分享
作者: hzpj01    时间: 2018-10-15 19:52
1024个点做完大概用时多少?
作者: bouna    时间: 2020-3-19 21:11
where is please the schematic ?
作者: 我是你爸爸丶    时间: 2020-6-3 22:41
这代码有点问题啊,adc压根没用到,fftoutput出来的都是一些随机数,keil仿真里 看就知道了。
作者: Accelerat    时间: 2020-8-13 22:27
问一下各位,这个是干什么的
作者: chenzm905    时间: 2021-8-1 14:53
有没有频率




欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1