簡介面向 STM32F407 學(xué)習(xí)與開發(fā)者以定時(shí)器觸發(fā) ADC 采樣用 DMA 搬運(yùn)數(shù)據(jù)后執(zhí)行 FFT并通過串口輸出頻譜結(jié)果幫助讀者走通從采樣率配置到頻域分析的完整流程。資源包共 187 個(gè)文件約 4.85MB核心由 58 個(gè) H 頭文件與 47 個(gè) C 源碼文件組成另含工程構(gòu)建配置、鏈接腳本、編譯輸出和調(diào)試文件方便直接了解 ADC、DMA、定時(shí)器及 DSP 庫的協(xié)作關(guān)系。項(xiàng)目支持 512KHZ、256KHZ、128KHZ 三檔采樣頻率也可修改采樣點(diǎn)數(shù)與頻率適合開展采樣定理驗(yàn)證和頻譜分辨率觀察。已有 5782 人學(xué)習(xí)特別適合嵌入式與信號(hào)處理入門者、備賽學(xué)生以及需要搭建采集-FFT-串口輸出實(shí)驗(yàn)平臺(tái)的技術(shù)人員。拿到后可打開工程查看初始化、中斷配置與數(shù)據(jù)處理鏈路參考結(jié)構(gòu)與關(guān)鍵代碼快速移植到其他 STM32 項(xiàng)目中并擴(kuò)展顯示、存儲(chǔ)或上位機(jī)分析功能。 做STM32F407上的ADC采集DMA搬運(yùn)FFT頻譜分析這個(gè)組合聽起來像是教科書里平平無奇的一章但真正在項(xiàng)目里落地時(shí)我第一版代碼是拿定時(shí)器觸發(fā)ADC、在中斷里讀數(shù)據(jù)、攢夠1024個(gè)點(diǎn)再調(diào)用FFT結(jié)果頻譜圖出來差點(diǎn)讓我把板子扔了一個(gè)干干凈凈的1kHz正弦波譜線旁邊全是毛刺和底噪主峰旁邊的旁瓣拖得很長而且每次采出來的結(jié)果都在輕微漂移。折騰了好幾天才明白問題基本不在FFT算法本身而是在ADC采樣的“時(shí)間質(zhì)量”和“數(shù)據(jù)搬運(yùn)方式”上。這篇就專門聊透STM32F407的ADC、DMA、FFT三者如何正確配合以及每一個(gè)環(huán)節(jié)里那些不寫進(jìn)參考手冊(cè)的坑。1. 中斷采樣喂給FFT為什么頻譜會(huì)不忍直視1.1 采樣點(diǎn)的時(shí)間間隔被中斷延遲污染了先說一個(gè)容易被忽略的事實(shí)FFT的數(shù)學(xué)前提是輸入序列滿足“等時(shí)間間隔采樣”。這就像你拿節(jié)拍器給樂手打拍子節(jié)拍器必須每0.5秒敲一下不能有時(shí)候隔0.45秒、有時(shí)候隔0.62秒。ADC中斷采樣方式下的真實(shí)情況是每次轉(zhuǎn)換完成ADC硬件置EOC標(biāo)志觸發(fā)中斷CPU響應(yīng)中斷現(xiàn)場壓棧讀取DR寄存器存進(jìn)數(shù)組退出中斷。這一整套流程里只要有一個(gè)更高優(yōu)先級(jí)的中斷插進(jìn)來或者Flash讀指令時(shí)遇到總線等待下一次采樣的觸發(fā)就會(huì)被延后。在STM32F407跑到168MHz主頻時(shí)一次中斷響應(yīng)加數(shù)據(jù)搬遷本身只需要幾百納秒但系統(tǒng)里往往不止ADC一個(gè)中斷源。我當(dāng)時(shí)的工程里還掛著串口接收、定時(shí)器更新、按鍵掃描等中斷FFT采樣周期部分被這些中斷搶占后相鄰采樣點(diǎn)的時(shí)間間隔實(shí)際上是一個(gè)不均勻序列。FFT對(duì)這種非均勻采樣非常敏感直觀表現(xiàn)就是頻譜底噪整體抬高、主峰旁邊出現(xiàn)不規(guī)則的雜散分量而且看起來像是信號(hào)本身不干凈實(shí)際上是時(shí)間軸“臟”了。測量一下就能驗(yàn)證把輸入短接到GND理論上FFT結(jié)果應(yīng)該是全頻段接近0但中斷采樣方式下這段“0信號(hào)”的頻譜也會(huì)有一層噪聲底座原因就是采樣點(diǎn)在時(shí)間上的抖動(dòng)被轉(zhuǎn)化成了幅度上的隨機(jī)誤差。1.2 CPU被搬運(yùn)任務(wù)拖死算力全浪費(fèi)在搬數(shù)據(jù)上還有一個(gè)繞不開的問題中斷方式下ADC每采樣一個(gè)點(diǎn)CPU就要被中斷打斷一次。假設(shè)采樣率是100kHz那么1秒鐘就有10萬次中斷。每次中斷就算只占用幾百個(gè)周期累計(jì)起來也是一筆不小的開銷。最難受的是這些開銷全花在“把數(shù)據(jù)從外設(shè)寄存器搬到內(nèi)存數(shù)組”這種毫無技術(shù)含量的重復(fù)勞動(dòng)上。而FFT本身尤其是1024點(diǎn)、4096點(diǎn)這種規(guī)模CMSIS-DSP庫在F407上雖然只要幾百微秒但如果CPU一邊忙著搬數(shù)據(jù)、一邊響應(yīng)其他中斷FFT運(yùn)算可能被反復(fù)打斷不僅總耗時(shí)變長實(shí)時(shí)性也難以保證。更不要說產(chǎn)品里往往還要同時(shí)跑顯示刷新、通信協(xié)議、按鍵邏輯中斷采樣方案在這種場景下基本撐不住。用DMA替代中斷本質(zhì)上是把“搬運(yùn)工”的活交給專門的硬件通道去做。CPU只需要在緩沖區(qū)拿到一批數(shù)據(jù)后批量處理。這就是DMA方案在ADCFFT場景里成為默認(rèn)選擇的核心原因。1.3 數(shù)據(jù)形態(tài)也對(duì)不上ADC寄存器里是整數(shù)FFT要吃浮點(diǎn)序列這一點(diǎn)很少被新手注意到但排查起來非常隱蔽。ADC的DR寄存器是12位右對(duì)齊的uint16_t類型而CMSIS-DSP的浮點(diǎn)FFT函數(shù)輸入是float32_t數(shù)組。有人圖省事直接把uint16_t數(shù)組強(qiáng)轉(zhuǎn)成float指針傳給FFT函數(shù)結(jié)果數(shù)據(jù)解釋完全錯(cuò)亂。因?yàn)閒loat在內(nèi)存里的IEEE754格式和uint16_t的整數(shù)布局完全不是一回事。正確的做法是在DMA搬完數(shù)據(jù)后把每個(gè)uint16_t的ADC原始值先轉(zhuǎn)換為float并按需要線性映射到0.0~3.3V按電壓算或0.0~1.0按滿量程歸一化。這一步看起來簡單但它決定了后續(xù)所有幅值計(jì)算的物理含義我見過好幾個(gè)項(xiàng)目卡在“FFT結(jié)果數(shù)量級(jí)完全不對(duì)”上最后發(fā)現(xiàn)是這一步的類型映射寫錯(cuò)了。所以在整個(gè)數(shù)據(jù)通路里明確區(qū)分“ADC原始值”和“FFT輸入序列”是很重要的一件事。2. 把DMA加進(jìn)來數(shù)據(jù)流架構(gòu)與Buffer設(shè)計(jì)里最容易犯的錯(cuò)2.1 ADC、DMA、內(nèi)存三者的分工先梳理一下整體數(shù)據(jù)流HAL庫的HAL_ADC_Start_DMA函數(shù)啟動(dòng)ADC轉(zhuǎn)換同時(shí)啟動(dòng)DMA搬運(yùn)。ADC每完成一次轉(zhuǎn)換產(chǎn)生一個(gè)DMA請(qǐng)求DMA控制器將這個(gè)轉(zhuǎn)換結(jié)果從ADC的數(shù)據(jù)寄存器搬到內(nèi)存Buffer中全程不經(jīng)過CPU。當(dāng)搬運(yùn)次數(shù)達(dá)到設(shè)定的Buffer長度時(shí)DMA觸發(fā)傳輸完成中斷這時(shí)CPU才介入把Buffer里的數(shù)據(jù)取走做FFT處理。如果配置為單次模式DMA搬完設(shè)定長度后會(huì)自動(dòng)停止如果配置為循環(huán)模式DMA會(huì)“原地轉(zhuǎn)圈”Buffer被不斷刷新。FFT應(yīng)用里通常使用循環(huán)模式這樣采集是無限連續(xù)的CPU在任何時(shí)刻去讀Buffer拿到的都是最近一段時(shí)間的數(shù)據(jù)。這就像一條自動(dòng)傳送帶不斷把工件送到倉庫倉庫管理員只需要在傳送帶堆滿一批時(shí)過去取走不需要每送一個(gè)工件就跑一趟。值得注意的是ADC連續(xù)轉(zhuǎn)換模式加上DMA循環(huán)模式會(huì)讓Buffer里的數(shù)據(jù)一直在被DMA刷新。CPU讀取Buffer做FFT時(shí)如果正趕上DMA寫入同一塊緩沖區(qū)就會(huì)讀到一半是舊數(shù)據(jù)、一半是新數(shù)據(jù)FFT結(jié)果自然不對(duì)。這就引出了雙緩沖和乒乓操作的必要性稍后細(xì)說。2.2 Buffer長度到底是“點(diǎn)數(shù)”還是“點(diǎn)數(shù)×通道數(shù)”這里最容易想歪很多人初次配置時(shí)以為“我要做1024點(diǎn)FFT那DMA搬運(yùn)次數(shù)就設(shè)成1024”。這是在單通道場景下的正確思路。但如果開了ADC掃描模式比如同時(shí)采樣兩路信號(hào)DMA搬運(yùn)回來的數(shù)據(jù)是按通道順序交替排列的ch0的第一個(gè)點(diǎn)、ch1的第一個(gè)點(diǎn)、ch0的第二個(gè)點(diǎn)、ch1的第二個(gè)點(diǎn)……這種情況下DMA傳輸次數(shù)1024只代表了“兩個(gè)通道各采了512個(gè)點(diǎn)”拿這1024個(gè)數(shù)據(jù)硬塞給FFT兩個(gè)通道的波形混在一起頻譜完全是亂的。正確做法是先把需求想清楚如果你需要每個(gè)通道都做1024點(diǎn)FFT并且使用掃描模式DMA的傳輸長度應(yīng)該設(shè)置成FFT點(diǎn)數(shù)乘以通道數(shù)也就是1024×22048。然后在DMA中斷回調(diào)里按通道號(hào)進(jìn)行取模抽取偏移量為偶數(shù)的數(shù)據(jù)屬于ch0偏移量為奇數(shù)的數(shù)據(jù)屬于ch1。抽取出來的數(shù)據(jù)再分別裝入兩個(gè)float數(shù)組做FFT。我之前寫過一段提取邏輯大致套路如下#define FFT_POINTS 1024 #define ADC_CH_NUM 2 uint16_t adc_dma_buf[FFT_POINTS * ADC_CH_NUM]; float fft_input_ch0[FFT_POINTS]; float fft_input_ch1[FFT_POINTS]; for (uint16_t i 0; i FFT_POINTS; i) { fft_input_ch0[i] (float)adc_dma_buf[i * ADC_CH_NUM 0]; fft_input_ch1[i] (float)adc_dma_buf[i * ADC_CH_NUM 1]; }這段代碼雖然簡單但它把“物理通道”和“FFT譜線”正確對(duì)應(yīng)起來了比一頭扎進(jìn)FFT函數(shù)里調(diào)參數(shù)更重要。2.3 雙緩沖與半傳輸中斷采集和運(yùn)算的相位錯(cuò)開我剛才提到循環(huán)DMA模式會(huì)讓數(shù)據(jù)不斷刷新CPU在任意時(shí)刻讀取Buffer都可能讀到剛好被改寫一半的數(shù)據(jù)。解決這個(gè)問題最經(jīng)典的方案不是“加鎖”而是利用DMA的半傳輸中斷Half Transfer和傳輸完成中斷Transfer Complete實(shí)現(xiàn)乒乓操作。具體機(jī)制是DMA緩沖區(qū)被分成前后兩半前半段寫滿時(shí)觸發(fā)半傳輸中斷此時(shí)CPU處理前半段數(shù)據(jù)在這段時(shí)間里DMA繼續(xù)往后半段搬運(yùn)數(shù)據(jù)等后半段也寫滿了DMA從頭重新開始填充前半段同時(shí)觸發(fā)傳輸完成中斷CPU再去處理后半段數(shù)據(jù)。整個(gè)過程里CPU處理的永遠(yuǎn)是“上一次已經(jīng)填滿、且DMA暫時(shí)不會(huì)再寫”的那一半不會(huì)出現(xiàn)讀寫競爭。在STM32 HAL庫中對(duì)應(yīng)的回調(diào)函數(shù)是HAL_ADC_ConvHalfCpltCallback和HAL_ADC_ConvCpltCallback。我在實(shí)際項(xiàng)目里用標(biāo)志位通知主循環(huán)處理不在中斷里直接跑FFT避免長時(shí)間占用中斷上下文。volatile uint8_t fft_ready 0; volatile uint8_t half_flag 0; void HAL_ADC_ConvHalfCpltCallback(ADC_HandleTypeDef *hadc) { (void)hadc; half_flag 1; } void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef *hadc) { (void)hadc; half_flag 2; }主循環(huán)里判斷half_flag為1處理前半段為2處理后半段處理完清標(biāo)志。這個(gè)結(jié)構(gòu)保證了FFT輸入序列的完整性是長時(shí)間穩(wěn)定運(yùn)行的關(guān)鍵否則跑半小時(shí)后FFT結(jié)果就會(huì)莫名跳變。3. CubeMX配置里的三個(gè)“埋雷點(diǎn)”搬一次就停、數(shù)值亂跳、數(shù)組沒對(duì)齊3.1 DMA模式選Normal還是Circular決定數(shù)據(jù)是“搬一次”還是“持續(xù)搬”CubeMX里配置ADC的DMA時(shí)大多數(shù)人會(huì)隨手選一個(gè)Normal模式。Normal模式的含義是DMA傳輸完設(shè)定長度后傳輸自動(dòng)停止通道關(guān)閉。這時(shí)ADC雖然還在繼續(xù)轉(zhuǎn)換但轉(zhuǎn)換結(jié)果沒人搬了Buffer里的數(shù)據(jù)永遠(yuǎn)停留在第一次傳輸完的狀態(tài)?，F(xiàn)象就是串口打印第一次數(shù)據(jù)看起來正常之后打印出來的永遠(yuǎn)是同一批數(shù)據(jù)。解決方法是把DMA模式改成Circular。Circular模式下DMA傳輸計(jì)數(shù)遞減到0后自動(dòng)重裝初始值繼續(xù)下一次傳輸形成無限循環(huán)。這樣只要ADC不停轉(zhuǎn)換DMA就不停搬運(yùn)CPU隨時(shí)能拿到一份“剛剛更新過”的緩沖區(qū)數(shù)據(jù)。判斷是不是這個(gè)問題的快速方法在調(diào)試器里看DMA的NDTR寄存器。如果它停在一個(gè)固定值不動(dòng)而ADC狀態(tài)寄存器還在繼續(xù)轉(zhuǎn)換基本就是Normal模式?jīng)]跑了。改完模式后重啟工程N(yùn)DTR應(yīng)該來回滾動(dòng)數(shù)據(jù)也會(huì)持續(xù)刷新。3.2 采樣時(shí)間別用最小檔源阻抗會(huì)讓ADC讀數(shù)“飄”STM32F407的ADC是逐次逼近型SAR結(jié)構(gòu)內(nèi)部有一個(gè)采樣保持電容。在采樣階段ADC引腳通過內(nèi)部模擬開關(guān)給這個(gè)電容充電。如果開關(guān)閉合時(shí)間太短而信號(hào)源的輸出阻抗又比較高比如前級(jí)是幾十kΩ的電阻分壓網(wǎng)絡(luò)電容還沒來得及充滿開關(guān)就斷開了轉(zhuǎn)換結(jié)果自然偏低而且會(huì)隨著溫度、電源電壓輕微變化而隨機(jī)跳動(dòng)。F407的ADC采樣時(shí)間可以從3個(gè)周期一路調(diào)到480個(gè)周期。CubeMX里默認(rèn)配置經(jīng)常是3個(gè)周期這在信號(hào)源內(nèi)阻很低比如運(yùn)放直接驅(qū)動(dòng)時(shí)沒問題但如果前端是幾十kΩ的分壓電阻3個(gè)周期的采樣時(shí)間遠(yuǎn)遠(yuǎn)不夠。我個(gè)人的經(jīng)驗(yàn)是ADC時(shí)鐘21MHz時(shí)采樣時(shí)間至少取28個(gè)周期以上比較穩(wěn)妥信號(hào)源阻抗較高時(shí)甚至要拉到144個(gè)周期。當(dāng)然采樣時(shí)間變長會(huì)直接降低等效采樣率具體數(shù)值可以按下面這個(gè)關(guān)系估算采樣時(shí)間設(shè)置ADC時(shí)鐘21MHz下單個(gè)轉(zhuǎn)換周期理論上限采樣率約3周期15.5周期1.35Msps15周期27.5周期764ksps28周期40.5周期518ksps144周期156.5周期134ksps480周期492.5周期42.6ksps如果只是做音頻頻段的FFT分析幾十kHz到一兩百kHz的等效采樣率完全夠用與其冒險(xiǎn)用短采樣時(shí)間去追采樣率不如把穩(wěn)定性放在第一位。3.3 數(shù)據(jù)寬度和內(nèi)存對(duì)齊CMSIS-DSP函數(shù)有硬性要求ADC的DR寄存器有效位數(shù)是16位右對(duì)齊所以DMA搬運(yùn)的數(shù)據(jù)寬度必須配置為Half Word半字16位。如果配置成Word32位或者Byte8位搬回來的數(shù)據(jù)就會(huì)錯(cuò)位Word會(huì)把相鄰寄存器的內(nèi)容一起讀進(jìn)來Byte會(huì)丟失高8位。這還只是第一步。真正隱藏的坑在FFT輸入數(shù)組的定義上。CMSIS-DSP的arm_rfft_fast_f32等函數(shù)內(nèi)部會(huì)使用SIMD指令和雙字加載要求傳入的float數(shù)組必須是32位對(duì)齊的內(nèi)存地址。如果你在棧上定義一個(gè)普通局部數(shù)組編譯器默認(rèn)對(duì)齊可能只有4字節(jié)甚至更差調(diào)用FFT函數(shù)時(shí)輕則性能下降重則直接進(jìn)入HardFault。解決辦法是定義數(shù)組時(shí)加上對(duì)齊屬性比如在GCC/Keil環(huán)境下__ALIGNED(4) float fft_input[FFT_POINTS]; __ALIGNED(4) float fft_output[FFT_POINTS];CubeMX生成的代碼里如果使用DMA通常會(huì)自動(dòng)把DMA緩沖區(qū)設(shè)置為32位對(duì)齊但自己另外定義的FFT輸入輸出數(shù)組很容易漏掉這一點(diǎn)。還有另一個(gè)細(xì)節(jié)CMSIS庫的實(shí)數(shù)FFT函數(shù)要求輸入輸出數(shù)組必須互不重疊否則會(huì)破壞內(nèi)部狀態(tài)。4. 從ADC原始值到頻點(diǎn)幅值CMSIS-DSP的正確姿勢與幅度還原4.1 arm_rfft_fast_f32的輸入輸出布局別拿Matlab習(xí)慣來套ADC采樣得到的是一串實(shí)數(shù)序列。實(shí)數(shù)序列的FFT結(jié)果具有共軛對(duì)稱性所以CMSIS-DSP專門提供了arm_rfft_fast_f32來處理這種情況它比通用的復(fù)數(shù)FFT函數(shù)省了將近一半的內(nèi)存和計(jì)算量。但它的輸出布局非常容易讓人看錯(cuò)函數(shù)并不會(huì)直接輸出你想象的N個(gè)“頻率點(diǎn)幅值”而是輸出N個(gè)float在第0個(gè)和第N/2個(gè)位置放DC分量和奈奎斯特頻率分量的實(shí)數(shù)中間位置交替存放正頻率分量的實(shí)部和虛部。參考代碼結(jié)構(gòu)如下arm_rfft_fast_instance_f32 fft_inst; arm_rfft_fast_init_f32(fft_inst, FFT_POINTS); float fft_input[FFT_POINTS]; // 時(shí)域輸入由ADC原始值轉(zhuǎn)換而來 float fft_output[FFT_POINTS]; // 頻域輸出布局為實(shí)虛交替 arm_rfft_fast_f32(fft_inst, fft_input, fft_output, 0); // 手動(dòng)計(jì)算幅度0點(diǎn)直流除外 float mag_db; for (uint16_t k 0; k FFT_POINTS / 2; k) { float re, im; if (k 0 || k FFT_POINTS / 2) { re fft_output[k]; im 0.0f; } else { re fft_output[2 * k]; im fft_output[2 * k 1]; } mag_half[k] sqrtf(re * re im * im); }手寫遍歷取模其實(shí)也不復(fù)雜而且能順便把不同bin的實(shí)部/虛部對(duì)應(yīng)關(guān)系理清楚。等代碼跑通后再?zèng)Q定要不要換成arm_cmplx_mag_f32等批量函數(shù)。4.2 減均值、加窗與頻譜泄露FFT本質(zhì)上是對(duì)“N點(diǎn)序列在一個(gè)周期內(nèi)”做傅里葉分析它的隱含假設(shè)是這N點(diǎn)序列是周期信號(hào)的一個(gè)整周期。但實(shí)際采集的N個(gè)點(diǎn)往往首尾不連續(xù)相當(dāng)于在時(shí)域乘了一個(gè)矩形窗頻譜上就會(huì)有能量從真正的主瓣“漏”到旁瓣去表現(xiàn)為主峰旁邊拖著一串衰減的波紋。加窗的作用就是讓N點(diǎn)序列的首尾都平滑過渡到接近0減弱這種非線性截?cái)?。最常用的是Hann窗它能把旁瓣壓得很低代價(jià)是主瓣寬度增加一倍頻率分辨能力略微下降。在實(shí)際代碼中減均值這一步往往比加窗還要優(yōu)先如果不減去DC分量FFT結(jié)果的0Hz處會(huì)有一個(gè)很大的直流譜線而它在Hann窗作用下還會(huì)向鄰近頻點(diǎn)泄漏把低頻段的小信號(hào)都淹沒了。完整預(yù)處理邏輯一般是這樣的float mean 0.0f; for (uint16_t i 0; i FFT_POINTS; i) { mean fft_input[i]; } mean / FFT_POINTS; for (uint16_t i 0; i FFT_POINTS; i) { float w 0.5f - 0.5f * cosf(2.0f * PI * i / (FFT_POINTS - 1)); fft_input[i] (fft_input[i] - mean) * w; }注意Hann窗會(huì)讓信號(hào)的幅度乘上一個(gè)約等于0.5的系數(shù)相干增益所以在后面做幅值還原時(shí)要把這個(gè)系數(shù)補(bǔ)回來。4.3 幅值換算單邊譜、Hann窗恢復(fù)因子與頻率分辨率做完FFT取模之后得到的|X[k]|并不直接等于信號(hào)的幅度。因?yàn)镃MSIS的FFT變換沒有做歸一化N點(diǎn)FFT結(jié)果的量級(jí)大約是原始信號(hào)幅度的N/2倍正頻率部分。還原真實(shí)幅值的完整公式要分兩步第一步幅度歸一化|X[k]| / N這樣回退到原始信號(hào)幅度第二步單邊譜合并k0的DC分量不乘2其余k0的頻點(diǎn)要乘以2因?yàn)樨?fù)頻率部分的能量被折回正頻率第三步如果使用了Hann窗再乘上窗恢復(fù)系數(shù)2因?yàn)镠ann相干增益約為0.5。三部分合起來對(duì)k0的有效頻點(diǎn)恢復(fù)系數(shù)是4/N。代碼中可以寫float real_amp 2.0f * 2.0f / (float)FFT_POINTS * sqrtf(re * re im * im);頻率分辨率則由采樣率和FFT點(diǎn)數(shù)共同決定Δf fs / N。比如采樣率20kHz做1024點(diǎn)FFT每個(gè)bin對(duì)應(yīng)的頻率寬度大約是19.5Hz。這意味著兩個(gè)頻率差小于19.5Hz的信號(hào)會(huì)在同一個(gè)bin里疊加無法分辨。想做更精細(xì)的頻率分析要么降低采樣率要么增加點(diǎn)數(shù)二者要按具體場景權(quán)衡。舉個(gè)例子如果要分辨10Hz間隔的邊帶而信號(hào)最高頻率是5kHz那N至少要滿足fs/Δf 10000/10 1000向上取2的冪就選1024。如果把同樣的點(diǎn)數(shù)用在40kHz采樣率上分辨率就只有39Hz邊帶信息全糊掉了。這個(gè)取舍是FFT應(yīng)用中最需要對(duì)系統(tǒng)級(jí)需求有清晰認(rèn)知的地方。5. 實(shí)測翻車記錄五個(gè)真實(shí)故障的完整排查鏈路5.1 DMA只搬了一次就停之后數(shù)據(jù)全是重復(fù)的第一段現(xiàn)象串口打印ADC數(shù)據(jù)第一次打印看起來正常再打印發(fā)現(xiàn)數(shù)值完全沒變化像死機(jī)一樣。用調(diào)試器查看DMA的NDTR寄存器發(fā)現(xiàn)它固定在一個(gè)值上不動(dòng)。排查過程我首先檢查HAL_ADC_Start_DMA是否被反復(fù)調(diào)用結(jié)果發(fā)現(xiàn)只調(diào)用了一次這個(gè)沒問題。再查DMA初始化的模式參數(shù)發(fā)現(xiàn)CubeMX生成代碼里DMA_HandleTypeDef的Init.Mode字段是DMA_NORMAL問題就出在這里。原因Normal模式下DMA傳輸次數(shù)到達(dá)設(shè)定值后通道自動(dòng)關(guān)閉。ADC還在繼續(xù)轉(zhuǎn)換、繼續(xù)產(chǎn)生DMA請(qǐng)求但DMA通道已經(jīng)不再響應(yīng)數(shù)據(jù)自然不被搬運(yùn)。改為DMA_CIRCULAR后NDTR寄存器會(huì)開始循環(huán)變化數(shù)據(jù)恢復(fù)連續(xù)更新。這個(gè)案例看起來簡單但它很典型很多“程序跑著跑著數(shù)據(jù)不動(dòng)了”的問題不是程序邏輯死了而是外設(shè)配置模式用錯(cuò)了。排查時(shí)先檢查硬件外設(shè)狀態(tài)寄存器比反復(fù)看邏輯代碼高效得多。5.2 干凈信號(hào)源FFT卻一片尖峰底噪高出預(yù)期現(xiàn)象用信號(hào)發(fā)生器輸入一個(gè)非常干凈的1kHz正弦波FFT結(jié)果里除了1kHz主峰還在幾十Hz到幾百Hz區(qū)間出現(xiàn)一堆小尖峰底噪比預(yù)期高出一個(gè)數(shù)量級(jí)。排查過程我做了兩個(gè)實(shí)驗(yàn)先直接把ADC引腳短接到GND跑FFT底噪依然很高說明干擾在ADC通路內(nèi)部或參考電源上再把ADC引腳接到一個(gè)低噪基準(zhǔn)電壓源比如干凈的1.65V底噪降下來說明問題在輸入通路的前端或電源質(zhì)量上。最后在Vref引腳旁加強(qiáng)濾波電容并把該引腳的走線避開數(shù)字信號(hào)線頻譜干凈了不少。這里要特別說一下STM32F407的參考電壓。很多最小系統(tǒng)板上Vref直接接3.3V電源而這個(gè)3.3V同時(shí)給一堆數(shù)字芯片供電開關(guān)噪聲全都會(huì)通過Vref耦合進(jìn)ADC轉(zhuǎn)換結(jié)果。FFT對(duì)噪聲極其敏感這種電源噪聲在時(shí)域看起來不明顯在頻域卻會(huì)表現(xiàn)為寬頻底噪或尖峰。在追求頻譜精度的場合至少要給Vref加一個(gè)高性能LDO和低ESR電容模擬地和數(shù)字地單點(diǎn)連接這個(gè)投資非常值得。5.3 頻率整體偏移了幾個(gè)bin信號(hào)源頻率明明很準(zhǔn)現(xiàn)象輸入1kHz正弦波FFT峰值出現(xiàn)在約1050Hz的位置不是1000Hz而且偏差比例基本固定。排查過程先確認(rèn)信號(hào)發(fā)生器輸出頻率足夠準(zhǔn)然后懷疑采樣率的實(shí)際值和理論值不一致。F407的ADC時(shí)鐘來自APB2總線經(jīng)過預(yù)分頻后才能喂給ADC。如果APB2總線和ADC預(yù)分頻配置錯(cuò)了比如我以為ADC時(shí)鐘是21MHz實(shí)際卻是其他值那么真正的采樣率和代碼里計(jì)算用的采樣率就對(duì)不上頻率自然整體偏差。我用FFT主峰的已知頻率反推實(shí)際采樣率如果程序設(shè)定的fs是20kHzFFT顯示1kHz信號(hào)在1050Hz那么實(shí)際采樣率約等于20kHz×(1050/1000)21kHz。按這個(gè)反推值去核對(duì)時(shí)鐘樹發(fā)現(xiàn)確實(shí)是ADC預(yù)分頻配置問題。修正后主峰立刻回到1000Hz位置。這也引出一個(gè)很實(shí)用的標(biāo)定方法做頻率測量類產(chǎn)品時(shí)可以用一個(gè)已知精度的信號(hào)源輸入利用FFT峰值偏移反算實(shí)際采樣率再用軟件校準(zhǔn)系數(shù)修正比逐個(gè)寄存器排查快得多。5.4 單通道正常多通道數(shù)據(jù)串位現(xiàn)象兩個(gè)通道掃描模式采集ch0輸入的信號(hào)卻在ch1的頻譜里看到了或者兩個(gè)通道的頻譜混在一起。排查過程先確認(rèn)ADC掃描順序DMA搬運(yùn)的數(shù)據(jù)是按rank順序排列的ch0在ch1前。問題出在提取數(shù)據(jù)時(shí)直接用了前1024個(gè)點(diǎn)沒有按通道數(shù)取模。由于DMA緩沖區(qū)是兩路交替排列前1024個(gè)點(diǎn)里其實(shí)混合了兩路數(shù)據(jù)直接當(dāng)單序列FFT結(jié)果當(dāng)然不對(duì)。按前面2.2小節(jié)的代碼重新抽取后兩路頻譜恢復(fù)正常。這類問題在代碼里往往很難一眼看出來但只要在調(diào)試器里查看DMA緩沖區(qū)的排列規(guī)律很快就能定位。5.5 長時(shí)間運(yùn)行后FFT結(jié)果突然亂掉甚至HardFault現(xiàn)象系統(tǒng)連續(xù)運(yùn)行幾分鐘后FFT結(jié)果莫名其妙變成一片噪聲偶爾直接進(jìn)入HardFault。排查過程先檢查內(nèi)存越界。FFT的輸入輸出數(shù)組長度都是1024但如果手動(dòng)計(jì)算幅值的數(shù)組長度只分配了512訪問k1024的一半邊界時(shí)越界會(huì)覆蓋其他變量。排查之后發(fā)現(xiàn)真正的問題在于CPU處理Buffer時(shí)DMA正在往同一個(gè)Buffer里寫新數(shù)據(jù)導(dǎo)致FFT讀取的數(shù)據(jù)被撕裂。解決方案就是用前面說的雙緩沖乒乓結(jié)構(gòu)讓DMA和CPU錯(cuò)開時(shí)段操作不同的半?yún)^(qū)。改成乒乓結(jié)構(gòu)之后連續(xù)跑了一整夜頻譜一直穩(wěn)定。另外一個(gè)小經(jīng)驗(yàn)HAL_ADC_ConvHalfCpltCallback和HAL_ADC_ConvCpltCallback里不要直接做FFT運(yùn)算和打印它們會(huì)拖累中斷響應(yīng)可能影響ADC采樣時(shí)序。最好只做標(biāo)志位或半?yún)^(qū)索引記錄把重活丟給主循環(huán)或RTOS任務(wù)。6. 讓FFT結(jié)果更“能打”的后續(xù)手段如果只是做Demo前面這些已經(jīng)夠了。但要在實(shí)際產(chǎn)品里把頻譜當(dāng)數(shù)據(jù)源使用還有幾件事值得做。第一件事是抗混疊濾波。FFT的輸入信號(hào)里如果混入了超過采樣率一半的頻率成分這些高頻成分會(huì)發(fā)生混疊折疊回低頻區(qū)形成“幽靈譜線”。軟件層面沒法有效消除混疊必須在ADC輸入前端加RC低通濾波器或者用運(yùn)放搭一個(gè)有源低通把超出fs/2的成分提前濾掉。我做音頻采樣時(shí)習(xí)慣在ADC引腳前放一個(gè)一階RC截止頻率設(shè)置在fs/2附近對(duì)高頻噪聲的壓制立竿見影。第二件事是考慮FFT運(yùn)算的實(shí)時(shí)調(diào)度。如果采樣率是20kHz一個(gè)1024點(diǎn)FFT的窗口是51.2msCMSIS的arm_rfft_fast_f32在168MHz主頻下跑1024點(diǎn)大約只要幾十到一百多微秒算力綽綽有余。但如果還要同時(shí)做顯示刷新、SD卡存儲(chǔ)、通信協(xié)議建議把FFT任務(wù)放到低優(yōu)先級(jí)循環(huán)里采樣和搬運(yùn)則完全依賴DMA不要用阻塞方式等FFT結(jié)果。如果需要更高幀率可以適當(dāng)減少點(diǎn)數(shù)或者把兩個(gè)DMA緩沖區(qū)擴(kuò)展成4段環(huán)形。第三件事是校準(zhǔn)幅值。ADC的增益誤差和偏置誤差、前級(jí)放大電路的倍率、Vref的不精確都會(huì)讓FFT算出來的幅度與實(shí)際物理量對(duì)不上。在產(chǎn)品化時(shí)通常用標(biāo)準(zhǔn)信號(hào)源輸入幾個(gè)已知幅度和頻率的點(diǎn)擬合出一條幅值校準(zhǔn)曲線在軟件里做補(bǔ)償。這個(gè)過程屬于長期調(diào)校但一旦做完了你的FFT結(jié)果就不只是“看起來有譜”而是可以當(dāng)作測量數(shù)據(jù)用了。我現(xiàn)在的習(xí)慣是拿到任何STM32F407的FFT需求先把“時(shí)鐘樹到底怎么了”“DMA是循環(huán)還是單次”“FFT輸入是否對(duì)齊、是否加窗、是否減均值”這三件事查清楚再開始寫應(yīng)用邏輯。這三步看著基礎(chǔ)卻決定了后面所有結(jié)果的可靠性。如果你也正在被ADC采出來的數(shù)據(jù)喂給FFT后頻譜亂糟糟的問題折磨希望這篇記錄能讓你少走一點(diǎn)彎路。本文還有配套的精品資源點(diǎn)擊獲取