找回密碼
 立即注冊

QQ登錄

只需一步，快速開始

搜索
查看: 5870|回復(fù): 4
收起左側(cè)

stm32f103官方DSP庫測試程序 可做128點、256點的fft運算,時間很短

  [復(fù)制鏈接]
ID:672840 發(fā)表于 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 //采樣頻率，在函數(shù)中沒定義，但是出的波形和它有關(guān)，具體體現(xiàn)在AD采樣頻率上，調(diào)節(jié)AD的采樣周期，
  19.                                                                                 和采樣函數(shù)GetADValueBuffer()的delay時間，可以設(shè)置其大小。數(shù)據(jù)裝入MATLAB，橫軸為點數(shù)，乘以頻域的分辨率即可。
  20.                                                                                 Fs=1/T;Tp=NT;Fs=NF；
  21.                                                                                 其中：
  22.                                                                                 Fs--時域采樣頻率；
  23.                                                                                 Tp--數(shù)據(jù)時間長度；
  24.                                                                                 T--時域采樣間隔；
  25.                                                                                 N--采樣點數(shù)；
  26.                                                                                 F--頻域采樣間隔；�。。。。�！分辨率
  27. PS:記得#include<math.h>，不然做sqrt時會出錯！并且編譯器不報錯�。。�
  28. ***************************************************************************************************/

  30. #define NPT 256 //fft點數(shù)，即N

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

  35. extern long lBufInArray[NPT]; //輸入緩沖區(qū)
  36. extern long lBufOutArray[NPT/2]; //輸出緩沖區(qū)
  37. extern long lBufMagArray[NPT/2];  //變換后緩沖區(qū)


  40. int main(void)
  41. {
  42.         int i;
  43.   u16 adcx;
  44.         int gewei,shifenwei,baifenwei,qianfenwei;
  45.         float temp;
  46.         delay_init();                     //延時函數(shù)初始化         
  47.         NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);//設(shè)置中斷優(yōu)先級分組為組2：2位搶占優(yōu)先級，2位響應(yīng)優(yōu)先級
  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的頻率計數(shù)
  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);//中文，不得出現(xiàn)校名，注釋了
  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位數(shù)字
  116.                 OLED_ShowString(90,6,"Hz");

  118.                
  119.                 LED0=!LED0;
  120.                 delay_ms(1000);
  121.         }
  122. }
復(fù)制代碼
TS%R9BS{V7__UUBXY2[J{KO.png
所有資料51hei提供下載:
STM官方DSP庫測試.7z (238.55 KB, 下載次數(shù): 145)

評分

參與人數(shù) 1黑幣 +50 收起 理由
admin + 50 共享資料的黑幣獎勵！

查看全部評分

回復(fù)

使用道具 舉報

ID:652775 發(fā)表于 2019-12-28 08:51 | 顯示全部樓層
收下了，謝樓主分享
回復(fù)

使用道具 舉報

ID:334838 發(fā)表于 2020-6-1 11:17 | 顯示全部樓層
官方庫不是 4基的嘛？還能測128點？
回復(fù)

使用道具 舉報

ID:778733 發(fā)表于 2023-3-23 18:33 | 顯示全部樓層
收下了，謝樓主分享，很NB的
回復(fù)

使用道具 舉報

ID:1089115 發(fā)表于 2023-7-20 16:56 | 顯示全部樓層
感謝樓主分享，比官方的點猛
回復(fù)

使用道具 舉報

您需要登錄后才可以回帖 登錄 | 立即注冊

本版積分規(guī)則

小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術(shù)交流QQ群281945664

Powered by 單片機教程網(wǎng)

快速回復(fù) 返回頂部 返回列表