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