设为首页收藏本站
开启辅助访问 切换到宽版

 找回密码
 立即注册

QQ登录

只需一步，快速开始

热门： 51单片机 | 24小时必答区 | 单片机教程 | 单片机DIY制作 | STM32 | Cortex M3 | 模数电子 | 电子DIY制作 | 音响/功放 | 拆机乐园 | Arduino | 嵌入式OS | 程序设计

搜索
»论坛 › 嵌入式/单片机论坛 › STM32/8 › stm32f103官方DSP库测试程序 可做128点、256点的fft运算 ...
返回列表 发新帖
查看: 4532|回复: 4
收起左侧

stm32f103官方DSP库测试程序 可做128点、256点的fft运算,时间很短

  [复制链接]
ID:672840 发表于 2019-12-24 20:23 | 显示全部楼层 |阅读模式
基于stm32f103测试，103存储较少，自己写的fft程序不适合在103上运行，使用STM32官方的库，可以做128点、256点的fft，运算时间很短；运算所需要的官方库文件在user目录下。
QFL%){$R0{$T8Y`0IV~NPH2.png

单片机源程序如下:
  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. }
复制代码
TS%R9BS{V7__UUBXY2[J{KO.png
所有资料51hei提供下载:
STM官方DSP库测试.7z (238.55 KB, 下载次数: 135)

评分

参与人数 1黑币 +50 收起 理由
admin + 50 共享资料的黑币奖励！

查看全部评分

回复

使用道具 举报

ID:652775 发表于 2019-12-28 08:51 | 显示全部楼层
收下了，谢楼主分享
回复

使用道具 举报

ID:334838 发表于 2020-6-1 11:17 | 显示全部楼层
官方库不是 4基的嘛？还能测128点？
回复

使用道具 举报

ID:778733 发表于 2023-3-23 18:33 | 显示全部楼层
收下了，谢楼主分享，很NB的
回复

使用道具 举报

ID:1089115 发表于 2023-7-20 16:56 | 显示全部楼层
感谢楼主分享，比官方的点猛
回复

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表