标题: stm32f103官方DSP库测试程序 可做128点、256点的fft运算,时间很短 [打印本页]
作者: dingxan    时间: 2019-12-24 20:23
标题: stm32f103官方DSP库测试程序 可做128点、256点的fft运算,时间很短
基于stm32f103测试，103存储较少，自己写的fft程序不适合在103上运行，使用STM32官方的库，可以做128点、256点的fft，运算时间很短；运算所需要的官方库文件在user目录下。


单片机源程序如下:
  1. #include "led.h"
  2. #include "delay.h"
  3. #include "key.h"
  4. #include "sys.h"
  5. #include "lcd.h"
  6. #include "usart.h"         
  7. #include "adc.h"
  8. #include "timer.h"
  9. #include "oled.h"
  10. #include "bmp.h"
  11. #include "stm32_dsp.h"
  12. #include "table_fft.h"
  13. #include "GetADValueBuffer.h"

  15. /*感谢大壮提供的频率采集部分程序*/

  17. /***************************************************************************************************  很重要
  18. //#define Fs 400 //采样频率，在函数中没定义，但是出的波形和它有关，具体体现在AD采样频率上，调节AD的采样周期，
  19.                                                                                 和采样函数GetADValueBuffer()的delay时间，可以设置其大小。数据装入MATLAB，横轴为点数，乘以频域的分辨率即可。
  20.                                                                                 Fs=1/T;Tp=NT;Fs=NF；
  21.                                                                                 其中：
  22.                                                                                 Fs--时域采样频率；
  23.                                                                                 Tp--数据时间长度；
  24.                                                                                 T--时域采样间隔；
  25.                                                                                 N--采样点数；
  26.                                                                                 F--频域采样间隔；！！！！！！分辨率
  27. PS:记得#include<math.h>，不然做sqrt时会出错！并且编译器不报错！！！
  28. ***************************************************************************************************/

  30. #define NPT 256 //fft点数，即N

  32. extern u16        count;
  33. extern u16        pinlv;

  35. extern long lBufInArray[NPT]; //输入缓冲区
  36. extern long lBufOutArray[NPT/2]; //输出缓冲区
  37. extern long lBufMagArray[NPT/2];  //变换后缓冲区


  40. int main(void)
  41. {
  42.         int i;
  43.   u16 adcx;
  44.         int gewei,shifenwei,baifenwei,qianfenwei;
  45.         float temp;
  46.         delay_init();                     //延时函数初始化         
  47.         NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);//设置中断优先级分组为组2：2位抢占优先级，2位响应优先级
  48.         uart_init(115200);                 //串口初始化为115200
  49.                         
  50.          Adc_Init();                                  //ADC初始化
  51.          
  52.         TIM3_Int_Init(9999,7199);//改成1s
  53.          TIM5_Cap_Init(0XFFFF,72-1);        //以1Mhz的频率计数
  54.          
  55.                  OLED_Init();                        //初始化OLED  
  56.                 OLED_Clear()          ;



  60.         while(1)
  61.         {
  62. //                adcx=Get_Adc_Average(ADC_Channel_1,10);

  64.                 GetADValueBuffer();

  66.                
  67. //                for(i=0;i<NPT;i++)
  68. //                {
  69. //                        printf("%ld\r\n",lBufInArray[i]);
  70. //                }
  71.                
  72.                
  73.                 cr4_fft_256_stm32(lBufOutArray, lBufInArray, NPT);
  74.                 GetPowerMag();
  75.                
  76. //                for(i=0;i<NPT/2;i++)
  77. //                {
  78. //                        printf("%ld\r\n",lBufOutArray[i]);
  79. //                }
  80.                
  81.                         for(i=0;i<NPT/2;i++)
  82.                 {
  83.                         printf("%ld\r\n",lBufMagArray[i]);
  84.                 }
  85.                
  86.                 printf("helloworld");
  87.                
  88. //                temp=(float)adcx*(3.3/4096);
  89. //                temp=temp/10;
  90.                 //printf("v=%.3f A\r\n",temp);
  91.                 //printf("f:%d Hz\r\n",pinlv);//打印频率
  92.                
  93.                 gewei=(int)temp;
  94.                 shifenwei=(int)(temp*10);
  95.                 baifenwei=(int)(temp*100);
  96.                 qianfenwei=(int)(temp*1000);
  97. //               
  98. //                OLED_ShowCHinese(0,0,0);
  99. //                OLED_ShowCHinese(18,0,1);
  100. //                OLED_ShowCHinese(36,0,2);
  101. //                OLED_ShowCHinese(54,0,3);
  102. //                OLED_ShowCHinese(72,0,4);
  103. //                OLED_ShowCHinese(90,0,5);
  104. //                OLED_ShowCHinese(108,0,6);//中文，不得出现校名，注释了
  105.                
  106.                 OLED_ShowString(0,3,"Ipp:");
  107.                 OLED_ShowNum(46,3,gewei,1,16);
  108.                 OLED_ShowString(54,3,".");
  109.                 OLED_ShowNum(62,3,shifenwei,1,16);
  110.                 OLED_ShowNum(70,3,baifenwei,1,16);
  111.                 OLED_ShowNum(78,3,qianfenwei,1,16);
  112.                 OLED_ShowString(90,3,"A");//电流
  113.                
  114.                 OLED_ShowString(0,6,"Fre:");
  115.                 OLED_ShowNum(46,6,pinlv,5,16);//最多5位数字
  116.                 OLED_ShowString(90,6,"Hz");

  118.                
  119.                 LED0=!LED0;
  120.                 delay_ms(1000);
  121.         }
  122. }
复制代码

所有资料51hei提供下载:
STM官方DSP库测试.7z (238.55 KB, 下载次数: 145)

作者: ar_fa    时间: 2019-12-28 08:51
收下了，谢楼主分享
作者: 蒜毫炒肉    时间: 2020-6-1 11:17
官方库不是 4基的嘛？还能测128点？
作者: lwgold    时间: 2023-3-23 18:33
收下了，谢楼主分享，很NB的
作者: Lumos6fly    时间: 2023-7-20 16:56
感谢楼主分享，比官方的点猛




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