簡介這份資源面向STM32嵌入式開發(fā)者與數(shù)字信號處理學習者提供一套基于STM32H743、通過STM32CubeMX搭建的純凈FFT驗證例程解決全網(wǎng)缺少簡單直接FFT性能驗證工程的問題。壓縮包共1698個文件約49.05MB以C源碼與頭文件為主體輔以匯編啟動文件、鏈接腳本、庫文件及工程配置完整覆蓋FFT庫與相關支持庫并附STM32CubeMX源文件便于查看配置。代碼利用定時器精確統(tǒng)計實時FFT運算時間通過串口輸出耗時與自動計算結(jié)果開啟FPU進行硬件加速雙精度浮點運算1024點FFT僅需約300us仍有優(yōu)化空間。目前已有2409人學習下載適合希望快速驗證H743浮點性能、研究FFT庫集成與CubeMX配置的讀者參考借鑒。1. 從一次電機電流環(huán)的崩潰說起為什么我盯上了 STM32H743 的 FPU 加速 FFT去年調(diào)一臺三相電機的電流環(huán)采樣率拉到 20kHz每周期都要做一次 1024 點 FFT 看諧波。一開始跑在 M4 上軟件浮點硬算單次 FFT 干掉 1.8ms電流環(huán)直接崩了——采樣窗口還沒算完下一個周期就來了。后來換到 STM32H743打開 FPU掛上 CMSIS-DSP 的arm_rfft_fast_f32實測單次 1024 點 FFT 只要 300us 左右整整快了六倍。這份例程就是當時那套能跑通的工程STM32CubeMX 生成初始化代碼FPU 使能CMSIS-DSP 庫掛載FFT 1024 點實測耗時 300us帶完整的 FFT 庫和相關支持庫。如果你也在做電機控制、音頻分析、振動監(jiān)測這類需要實時頻譜的場景又不想在浮點運算上翻車這套東西可以直接拿去改。它適合兩類人一類是剛上手 H7 系列、想搞清楚 FPU 和 DSP 庫怎么配的新手另一類是被軟件浮點坑過、想看看硬件加速到底能壓到多少 us 的老手。2. 工程架構(gòu)與目錄拆解CubeMX 生成骨架CMSIS-DSP 填肉2.1 整體分層HAL 層、DSP 庫層、應用層這個工程不是那種把所有代碼塞進main.c的 Demo它按三層來組織。最底下是 STM32CubeMX 生成的 HAL 驅(qū)動層負責時鐘樹、GPIO、DMA、ADC 這些外設初始化中間是 CMSIS-DSP 庫層提供arm_rfft_fast_f32、arm_cmplx_mag_f32這些經(jīng)過匯編級優(yōu)化的浮點運算函數(shù)最上面是應用層放 FFT 的調(diào)用邏輯、數(shù)據(jù)搬運和耗時測量。這樣分層的好處是你換采樣源或者換 FFT 點數(shù)只動應用層就行底層不用碰。CubeMX 生成的部分集中在Core/Src和Core/IncDSP 庫放在Drivers/CMSIS/DSP下面。注意 H7 的 CMSIS-DSP 要用帶arm_math.h的那套別拿 F4 的舊庫直接替換指令集不一樣鏈接會報錯。2.2 目錄結(jié)構(gòu)一覽FFT_H743_Demo/ ├── Core/ │ ├── Inc/ │ │ ├── main.h │ │ ├── stm32h7xx_hal_conf.h │ │ └── stm32h7xx_it.h │ └── Src/ │ ├── main.c // 主循環(huán)FFT 調(diào)用入口 │ ├── stm32h7xx_hal_msp.c // HAL 底層初始化 │ ├── stm32h7xx_it.c // 中斷服務 │ └── system_stm32h7xx.c // 系統(tǒng)時鐘 ├── Drivers/ │ ├── CMSIS/ │ │ ├── DSP/ │ │ │ ├── Include/ // arm_math.h 等頭文件 │ │ │ └── Source/ // DSP 源碼可選編譯 │ │ └── Include/ // core_cm7.h │ └── STM32H7xx_HAL_Driver/ // HAL 庫 ├── MDK-ARM/ // Keil 工程文件 │ └── FFT_H743.uvprojx └── FFT_H743.ioc // CubeMX 配置文件這個結(jié)構(gòu)里FFT_H743.ioc是關鍵你后面要改時鐘、改外設直接雙擊它重新生成代碼不用手動改main.c里的初始化部分。MDK-ARM目錄下是 Keil 工程如果你用 STM32CubeIDE把源文件拖進去重建工程也行但注意 DSP 庫的路徑要重新指。2.3 關鍵配置FPU 和 DSP 庫怎么掛上去CubeMX 里使能 FPU 很簡單在Project Manager的Code Generator里勾上Copy only necessary library files然后在Pinout Configuration的System Core里找到CORTEX_M7把FPU設成Single Precision。但光這樣還不夠CMSIS-DSP 庫得手動加。我一般會在 Keil 的Options for Target里C/C選項卡的Define框加上ARM_MATH_CM7和__FPU_PRESENT1然后在Include Paths里把Drivers/CMSIS/DSP/Include加進去。庫文件用arm_cortexM7lfsp_math.lib這個 lfsp 后綴代表小端、單精度浮點別選錯。提示如果你在 CubeMX 里勾了Use MicroLIBDSP 庫的某些函數(shù)可能會和標準庫沖突建議先不勾等跑通了再試。3. 核心代碼逐段拆解從 ADC 采樣到 300us FFT 輸出3.1 時鐘樹配置讓 H743 跑到 480MHzH743 的 FFT 能壓到 300us一半功勞在時鐘。CubeMX 里我把 SYSCLK 設到 480MHzAHB 分頻 1APB1 分頻 4APB2 分頻 2。這樣 Cortex-M7 內(nèi)核跑滿FPU 的吞吐才上得來。具體配置在SystemClock_Config()里CubeMX 會自動生成但你要檢查PLL1的DIVN和DIVP對不對。常見做法是 HSE 用 25MHz 晶振DIVM設 5DIVN設 192DIVP設 2算下來 VCO 是 960MHz系統(tǒng)時鐘 480MHz。void SystemClock_Config(void) { RCC_OscInitTypeDef RCC_OscInitStruct {0}; RCC_ClkInitTypeDef RCC_ClkInitStruct {0}; // 使能 HSE25MHz 晶振 RCC_OscInitStruct.OscillatorType RCC_OSCILLATORTYPE_HSE; RCC_OscInitStruct.HSEState RCC_HSE_ON; RCC_OscInitStruct.PLL.PLLState RCC_PLL_ON; RCC_OscInitStruct.PLL.PLLSource RCC_PLLSOURCE_HSE; RCC_OscInitStruct.PLL.PLLM 5; // 25MHz / 5 5MHz RCC_OscInitStruct.PLL.PLLN 192; // 5MHz * 192 960MHz RCC_OscInitStruct.PLL.PLLP 2; // 960MHz / 2 480MHz RCC_OscInitStruct.PLL.PLLQ 4; RCC_OscInitStruct.PLL.PLLR 2; HAL_RCC_OscConfig(RCC_OscInitStruct); // 切換系統(tǒng)時鐘到 PLL1 RCC_ClkInitStruct.ClockType RCC_CLOCKTYPE_HCLK | RCC_CLOCKTYPE_SYSCLK | RCC_CLOCKTYPE_PCLK1 | RCC_CLOCKTYPE_PCLK2; RCC_ClkInitStruct.SYSCLKSource RCC_SYSCLKSOURCE_PLLCLK; RCC_ClkInitStruct.AHBCLKDivider RCC_SYSCLK_DIV1; // AHB 480MHz RCC_ClkInitStruct.APB1CLKDivider RCC_HCLK_DIV4; // APB1 120MHz RCC_ClkInitStruct.APB2CLKDivider RCC_HCLK_DIV2; // APB2 240MHz HAL_RCC_ClockConfig(RCC_ClkInitStruct, FLASH_LATENCY_4); }這段代碼里FLASH_LATENCY_4很關鍵480MHz 下 Flash 要插 4 個等待周期設少了會取指錯誤程序跑飛。PLLM、PLLN、PLLP這三個參數(shù)決定了最終頻率你換晶振的話要重新算。我一般會在main.c開頭加一句SystemCoreClockUpdate()確保SystemCoreClock變量和實際一致不然HAL_Delay會不準。3.2 FFT 初始化arm_rfft_fast_instance_f32 的坑CMSIS-DSP 的實數(shù) FFT 用arm_rfft_fast_instance_f32結(jié)構(gòu)體初始化的時候要指定點數(shù)。1024 點對應ARM_RFFT_FAST_1024這個枚舉值別寫錯寫錯了初始化會返回ARM_MATH_ARGUMENT_ERROR但很多人不看返回值后面算出來全是零。#include arm_math.h #define FFT_SIZE 1024 static arm_rfft_fast_instance_f32 fft_inst; static float32_t fft_input[FFT_SIZE]; static float32_t fft_output[FFT_SIZE]; static float32_t fft_mag[FFT_SIZE / 2]; void fft_init(void) { arm_status status; // 初始化 1024 點實數(shù) FFT 實例 status arm_rfft_fast_init_f32(fft_inst, FFT_SIZE); if (status ! ARM_MATH_SUCCESS) { // 初始化失敗通常是點數(shù)不支持或內(nèi)存不夠 while (1); } }arm_rfft_fast_init_f32的第二個參數(shù)是點數(shù)必須是 2 的冪32 到 4096 之間。H743 的 DTCM 有 128KB放這幾個數(shù)組綽綽有余但如果你把數(shù)組定義成局部變量棧會爆。我一般把它們放全局區(qū)或者用__attribute__((section(.DTCM)))指到 DTCM 里訪問速度比 AXI SRAM 快。3.3 執(zhí)行 FFT 并測耗時DWT 周期計數(shù)器測耗時別用HAL_GetTick()那個精度只有 1ms300us 根本測不出來。用 DWT 的周期計數(shù)器480MHz 下分辨率約 2ns足夠準。// 使能 DWT 周期計數(shù)器 void dwt_init(void) { CoreDebug-DEMCR | CoreDebug_DEMCR_TRCENA_Msk; DWT-CYCCNT 0; DWT-CTRL | DWT_CTRL_CYCCNTENA_Msk; } // 執(zhí)行一次 FFT 并返回耗時微秒 uint32_t fft_run_once(void) { uint32_t start, end; // 填充測試數(shù)據(jù)50kHz 正弦波假設采樣率 200kHz for (int i 0; i FFT_SIZE; i) { fft_input[i] arm_sin_f32(2.0f * PI * 50000.0f * i / 200000.0f); } start DWT-CYCCNT; // 實數(shù) FFT輸入 fft_input輸出 fft_output arm_rfft_fast_f32(fft_inst, fft_input, fft_output, 0); // 計算幅值輸出 512 個點 arm_cmplx_mag_f32(fft_output, fft_mag, FFT_SIZE / 2); end DWT-CYCCNT; // 480MHz 下周期數(shù)除以 480 就是微秒 return (end - start) / 480; }arm_rfft_fast_f32的第四個參數(shù)是ifftFlag0 表示正變換1 表示逆變換。輸出fft_output是交錯排列的復數(shù)arm_cmplx_mag_f32把它轉(zhuǎn)成幅值。實測下來FFT 本身約 280us加上幅值計算約 20us總共 300us 左右。如果你只做 FFT 不做幅值能壓到 280us 以內(nèi)。注意arm_cmplx_mag_f32的輸出長度是FFT_SIZE / 2因為實數(shù) FFT 的頻譜對稱只取前一半。別傳FFT_SIZE會越界。3.4 主循環(huán)把 FFT 掛到實際采樣流程里實際項目里不會在main里死循環(huán)跑 FFT一般是 ADC 采樣完一幀DMA 觸發(fā)中斷在中斷里置標志主循環(huán)檢測到標志再算。這樣不阻塞其他任務。volatile uint8_t fft_ready 0; int main(void) { HAL_Init(); SystemClock_Config(); dwt_init(); fft_init(); // 假設 ADC DMA 已經(jīng)配好每 1024 個點觸發(fā)一次中斷 // 中斷里把數(shù)據(jù)搬到 fft_input然后 fft_ready 1 while (1) { if (fft_ready) { fft_ready 0; uint32_t us fft_run_once(); // us 就是本次 FFT 耗時可以打到串口或者存起來 (void)us; } } }這段代碼里fft_ready用volatile修飾防止編譯器優(yōu)化掉。實際采樣中斷的配置在stm32h7xx_it.c里CubeMX 會生成 DMA 中斷入口你只需要在回調(diào)里搬數(shù)據(jù)。如果你用 ADC 的過采樣或者定時器觸發(fā)注意采樣率要和 FFT 點數(shù)匹配不然頻譜會混疊。4. 編譯、下載與實測Keil 和 CubeIDE 兩條路4.1 Keil MDK 下的編譯配置打開MDK-ARM/FFT_H743.uvprojx先檢查Options for Target里的Target選項卡ARM Compiler選Use default compiler version 6別用 AC5AC5 對 H7 的 FPU 支持有問題。C/C選項卡里Define加上ARM_MATH_CM7,__FPU_PRESENT1,__FPU_USED1Include Paths加上../Drivers/CMSIS/DSP/Include和../Drivers/CMSIS/Include。Linker選項卡里勾上Use Memory Layout from Target DialogH743 的 Flash 從0x08000000開始大小 2MBRAM 分 DTCM、AXI SRAM 幾塊默認配置就行。編譯的時候如果報arm_rfft_fast_f32未定義說明庫沒鏈上。在Manage Project Items里新建一個組把arm_cortexM7lfsp_math.lib加進去路徑在Drivers/CMSIS/Lib/ARM下面。這個 lib 是預編譯的不用把 DSP 源碼全加進工程編譯速度快很多。4.2 STM32CubeIDE 下的移植如果你用 CubeIDE新建一個 H743 的工程把Core/Src和Core/Inc里的文件替換成例程里的然后在Project Properties的C/C Build-Settings-MCU Settings里勾上Enable floating point unit。DSP 庫的添加方式不一樣CubeIDE 里要在Project Properties-C/C General-Paths and Symbols里加 Include 路徑然后在Source Location里把Drivers/CMSIS/DSP/Source整個目錄加進來讓 IDE 編譯源碼而不是鏈 lib。這樣編譯出來的代碼大一點但調(diào)試方便能單步進 DSP 函數(shù)。提示CubeIDE 編譯 DSP 源碼時注意把arm_rfft_fast_f32.c和arm_cmplx_mag_f32.c所在的目錄都加進去不然鏈接會缺函數(shù)。4.3 實測數(shù)據(jù)與驗證方法下載到板子上之后怎么確認 FFT 真的跑了 300us我一般用兩種方法交叉驗證。第一種是 DWT 計數(shù)器在fft_run_once里讀CYCCNT通過串口打出來。第二種是用 GPIO 翻轉(zhuǎn)在 FFT 前后拉高拉低一個 IO用示波器看脈寬。兩種方法測出來差不超過 5us說明數(shù)據(jù)可信。實測數(shù)據(jù)480MHz 主頻1024 點實數(shù) FFTarm_rfft_fast_f32耗時約 280usarm_cmplx_mag_f32耗時約 20us合計 300us。如果只做 FFT 不做幅值280us。這個數(shù)據(jù)和 ST 官方給的 benchmark 基本一致說明 FPU 和 DSP 庫都正常工作。如果你測出來是 1ms 以上先檢查 FPU 有沒有使能再看 DSP 庫是不是選錯了型號。5. 避坑與常見問題那些讓我熬夜的配置錯誤5.1 FPU 沒使能FFT 慢十倍現(xiàn)象FFT 耗時 3ms 以上和 M4 軟件浮點差不多。原因CubeMX 里CORTEX_M7的FPU沒設成Single Precision或者 Keil 的Define里漏了__FPU_PRESENT1。解決在 CubeMX 里重新生成代碼檢查stm32h7xx_hal_conf.h里__FPU_PRESENT是不是 1然后在 Keil 的Target選項卡里確認Floating Point Hardware選的是Single Precision。5.2 DSP 庫版本不匹配鏈接報錯現(xiàn)象編譯時報undefined symbol arm_rfft_fast_init_f32或者鏈接時提示arm_cortexM7lfsp_math.lib格式不對。原因用了 F4 的 DSP 庫或者 lib 文件選成了arm_cortexM7l_math.lib沒有 fsp 后綴。解決確認Drivers/CMSIS/DSP/Lib/ARM下的 lib 文件名帶lfspKeil 里Define加上ARM_MATH_CM7。如果還不行直接把 DSP 源碼加進工程編譯繞過 lib 兼容性問題。5.3 數(shù)組放棧上跑一次就 HardFault現(xiàn)象程序下載后一跑 FFT 就進 HardFault或者結(jié)果全是零。原因fft_input、fft_output定義在函數(shù)內(nèi)部1024 個 float 就是 4KB加上fft_output又是 4KB棧默認只有 1KB 左右直接溢出。解決把數(shù)組定義成全局變量或者用static修飾再或者把棧改大。我一般直接放全局區(qū)省事。5.4 采樣率和 FFT 點數(shù)不匹配頻譜混疊現(xiàn)象FFT 出來的峰值頻率和實際信號對不上或者出現(xiàn)多個假峰。原因采樣率低于信號最高頻率的 2 倍或者 FFT 點數(shù)對應的頻率分辨率沒算對。解決采樣率至少是信號最高頻率的 2.5 倍1024 點 FFT 的頻率分辨率是采樣率 / 1024。比如采樣率 200kHz分辨率就是 195Hz50kHz 的信號落在第 256 個 bin 上。5.5 DWT 計數(shù)器沒使能耗時測出來是零現(xiàn)象fft_run_once返回 0 或者一個固定的小值。原因DWT-CTRL的CYCCNTENA位沒置 1或者CoreDebug-DEMCR的TRCENA沒開。解決在dwt_init里先寫CoreDebug-DEMCR | CoreDebug_DEMCR_TRCENA_Msk再寫DWT-CTRL | DWT_CTRL_CYCCNTENA_Msk。注意這兩句的順序不能反先開調(diào)試單元再開計數(shù)器。6. 二次開發(fā)換點數(shù)、換窗函數(shù)、接實際傳感器6.1 改 FFT 點數(shù)從 1024 到 2048 或 512改點數(shù)不是把FFT_SIZE宏一改就完事。arm_rfft_fast_init_f32的第二個參數(shù)要跟著改ARM_RFFT_FAST_1024要換成ARM_RFFT_FAST_2048或ARM_RFFT_FAST_512。數(shù)組大小也要同步改fft_mag的長度是FFT_SIZE / 2。改完之后重新編譯如果初始化返回錯誤說明點數(shù)不支持CMSIS-DSP 的實數(shù) FFT 支持 32 到 4096 的 2 的冪。2048 點的耗時大約是 1024 點的 2.2 倍因為N log N的關系不是線性翻倍。#define FFT_SIZE 2048 // 改成 2048 static arm_rfft_fast_instance_f32 fft_inst; static float32_t fft_input[FFT_SIZE]; static float32_t fft_output[FFT_SIZE]; static float32_t fft_mag[FFT_SIZE / 2]; void fft_init(void) { // 點數(shù)變了枚舉值也要變 arm_rfft_fast_init_f32(fft_inst, FFT_SIZE); }6.2 加窗函數(shù)抑制頻譜泄漏實際傳感器信號不會正好整周期采樣不加窗的話頻譜泄漏很嚴重峰值旁邊全是裙邊。CMSIS-DSP 提供了arm_hanning_f32、arm_hamming_f32這些窗函數(shù)用法很簡單在 FFT 之前乘上去就行。static float32_t window[FFT_SIZE]; void apply_window(void) { // 生成漢寧窗 arm_hanning_f32(window, FFT_SIZE); // 逐點相乘 for (int i 0; i FFT_SIZE; i) { fft_input[i] * window[i]; } }arm_hanning_f32的第二個參數(shù)是點數(shù)生成的窗系數(shù)存在window數(shù)組里。加窗之后幅值會衰減需要乘一個補償系數(shù)漢寧窗的補償系數(shù)是 2.0。如果你做精確幅值測量這一步不能省。6.3 接實際傳感器ADC 配置和 DMA 搬運例程里用的是模擬正弦波實際項目要接 ADC。H743 的 ADC 支持 16 位最高采樣率取決于時鐘配置。我一般用定時器觸發(fā) ADCDMA 搬運到fft_input搬滿 1024 個點觸發(fā)中斷。CubeMX 里配置 ADC 的時候注意Continuous Conversion Mode要關掉用External Trigger觸發(fā)觸發(fā)源選定時器的TRGO事件。DMA 設成Circular模式這樣不用每次重新啟動。// ADC DMA 完成回調(diào)在 stm32h7xx_it.c 里 void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef *hadc) { // 一幀數(shù)據(jù)搬完了置標志讓主循環(huán)算 FFT fft_ready 1; }這個回調(diào)里不要做 FFT只置標志。FFT 耗時 300us在中斷里跑會阻塞其他中斷尤其是電流環(huán)這種對實時性要求高的場景。主循環(huán)檢測到fft_ready再算算完把結(jié)果存起來通過串口或者 CAN 發(fā)出去。6.4 驗證方法用已知頻率信號對一遍改完代碼怎么確認沒改錯我一般用信號發(fā)生器打一個已知頻率的正弦波比如 10kHz采樣率設 100kHz1024 點 FFT。理論上峰值應該落在第 102 個 bin 上10kHz / (100kHz / 1024) ≈ 102.4。如果落在別的地方先檢查采樣率對不對再看窗函數(shù)有沒有加錯。幅值也要對輸入 1Vpp 的正弦波加漢寧窗后峰值應該在 0.5 左右乘補償系數(shù)后回到 1.0。這套驗證流程我每次改 FFT 相關代碼都走一遍能擋掉大部分低級錯誤。從那以后我每次配 H7 的 DSP 庫都強制先跑一遍這個 1024 點 FFT 例程確認 FPU 和庫都正常再往上疊業(yè)務代碼。希望幫到你。本文還有配套的精品資源點擊獲取