標題: 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輸入數(shù)組
  13. float fft_outputbuf[FFT_LENGTH];        //FFT輸出數(shù)組

  14. u8 timeout;//定時器溢出次數(shù)
  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);//設置系統(tǒng)中斷優(yōu)先級分組2
  26.         delay_init(168);  //初始化延時函數(shù)
  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計數(shù)頻率,最大計時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執(zhí)行時間
  44.         POINT_COLOR=BLUE;        //設置字體為藍色   
  45.         arm_cfft_radix4_init_f32(&scfft,FFT_LENGTH,0,1);//初始化scfft結構體,設定FFT相關參數(shù)
  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定時器的計數(shù)器值
  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);        //把運算結果復數(shù)求模得幅值
  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, 下載次數(shù): 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出來的都是一些隨機數(shù),keil仿真里 看就知道了。
作者: Accelerat    時間: 2020-8-13 22:27
問一下各位,這個是干什么的
作者: chenzm905    時間: 2021-8-1 14:53
有沒有頻率




歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1