簡(jiǎn)介面向 STM32F407 學(xué)習(xí)與開(kāi)發(fā)者以定時(shí)器觸發(fā) ADC 采樣用 DMA 搬運(yùn)數(shù)據(jù)后執(zhí)行 FFT并通過(guò)串口輸出頻譜結(jié)果幫助讀者走通從采樣率配置到頻域分析的完整流程。資源包共 187 個(gè)文件約 4.85MB核心由 58 個(gè) H 頭文件與 47 個(gè) C 源碼文件組成另含工程構(gòu)建配置、鏈接腳本、編譯輸出和調(diào)試文件方便直接了解 ADC、DMA、定時(shí)器及 DSP 庫(kù)的協(xié)作關(guān)系。項(xiàng)目支持 512KHZ、256KHZ、128KHZ 三檔采樣頻率也可修改采樣點(diǎn)數(shù)與頻率適合開(kāi)展采樣定理驗(yàn)證和頻譜分辨率觀察。已有 5782 人學(xué)習(xí)特別適合嵌入式與信號(hào)處理入門者、備賽學(xué)生以及需要搭建采集-FFT-串口輸出實(shí)驗(yàn)平臺(tái)的技術(shù)人員。拿到后可打開(kāi)工程查看初始化、中斷配置與數(shù)據(jù)處理鏈路參考結(jié)構(gòu)與關(guān)鍵代碼快速移植到其他 STM32 項(xiàng)目中并擴(kuò)展顯示、存儲(chǔ)或上位機(jī)分析功能。 做STM32F407上的ADC采集DMA搬運(yùn)FFT頻譜分析這個(gè)組合聽(tīng)起來(lái)像是教科書里平平無(wú)奇的一章但真正在項(xiàng)目里落地時(shí)我第一版代碼是拿定時(shí)器觸發(fā)ADC、在中斷里讀數(shù)據(jù)、攢夠1024個(gè)點(diǎn)再調(diào)用FFT結(jié)果頻譜圖出來(lái)差點(diǎn)讓我把板子扔了一個(gè)干干凈凈的1kHz正弦波譜線旁邊全是毛刺和底噪主峰旁邊的旁瓣拖得很長(zhǎng)而且每次采出來(lái)的結(jié)果都在輕微漂移。折騰了好幾天才明白問(wèn)題基本不在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í)間間隔被中斷延遲污染了先說(shuō)一個(gè)容易被忽略的事實(shí)FFT的數(shù)學(xué)前提是輸入序列滿足“等時(shí)間間隔采樣”。這就像你拿節(jié)拍器給樂(lè)手打拍子節(jié)拍器必須每0.5秒敲一下不能有時(shí)候隔0.45秒、有時(shí)候隔0.62秒。ADC中斷采樣方式下的真實(shí)情況是每次轉(zhuǎn)換完成ADC硬件置EOC標(biāo)志觸發(fā)中斷CPU響應(yīng)中斷現(xiàn)場(chǎng)壓棧讀取DR寄存器存進(jìn)數(shù)組退出中斷。這一整套流程里只要有一個(gè)更高優(yōu)先級(jí)的中斷插進(jìn)來(lái)或者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ī)則的雜散分量而且看起來(lái)像是信號(hào)本身不干凈實(shí)際上是時(shí)間軸“臟”了。測(cè)量一下就能驗(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è)繞不開(kāi)的問(wèn)題中斷方式下ADC每采樣一個(gè)點(diǎn)CPU就要被中斷打斷一次。假設(shè)采樣率是100kHz那么1秒鐘就有10萬(wàn)次中斷。每次中斷就算只占用幾百個(gè)周期累計(jì)起來(lái)也是一筆不小的開(kāi)銷。最難受的是這些開(kāi)銷全花在“把數(shù)據(jù)從外設(shè)寄存器搬到內(nèi)存數(shù)組”這種毫無(wú)技術(shù)含量的重復(fù)勞動(dòng)上。而FFT本身尤其是1024點(diǎn)、4096點(diǎn)這種規(guī)模CMSIS-DSP庫(kù)在F407上雖然只要幾百微秒但如果CPU一邊忙著搬數(shù)據(jù)、一邊響應(yīng)其他中斷FFT運(yùn)算可能被反復(fù)打斷不僅總耗時(shí)變長(zhǎng)實(shí)時(shí)性也難以保證。更不要說(shuō)產(chǎn)品里往往還要同時(shí)跑顯示刷新、通信協(xié)議、按鍵邏輯中斷采樣方案在這種場(chǎng)景下基本撐不住。用DMA替代中斷本質(zhì)上是把“搬運(yùn)工”的活交給專門的硬件通道去做。CPU只需要在緩沖區(qū)拿到一批數(shù)據(jù)后批量處理。這就是DMA方案在ADCFFT場(chǎng)景里成為默認(rèn)選擇的核心原因。1.3 數(shù)據(jù)形態(tài)也對(duì)不上ADC寄存器里是整數(shù)FFT要吃浮點(diǎn)序列這一點(diǎn)很少被新手注意到但排查起來(lái)非常隱蔽。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按滿量程歸一化。這一步看起來(lái)簡(jiǎn)單但它決定了后續(xù)所有幅值計(jì)算的物理含義我見(jiàn)過(guò)好幾個(gè)項(xiàng)目卡在“FFT結(jié)果數(shù)量級(jí)完全不對(duì)”上最后發(fā)現(xiàn)是這一步的類型映射寫錯(cuò)了。所以在整個(gè)數(shù)據(jù)通路里明確區(qū)分“ADC原始值”和“FFT輸入序列”是很重要的一件事。2. 把DMA加進(jìn)來(lái)數(shù)據(jù)流架構(gòu)與Buffer設(shè)計(jì)里最容易犯的錯(cuò)2.1 ADC、DMA、內(nèi)存三者的分工先梳理一下整體數(shù)據(jù)流HAL庫(kù)的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)過(guò)CPU。當(dāng)搬運(yùn)次數(shù)達(dá)到設(shè)定的Buffer長(zhǎng)度時(shí)DMA觸發(fā)傳輸完成中斷這時(shí)CPU才介入把Buffer里的數(shù)據(jù)取走做FFT處理。如果配置為單次模式DMA搬完設(shè)定長(zhǎng)度后會(huì)自動(dòng)停止如果配置為循環(huán)模式DMA會(huì)“原地轉(zhuǎn)圈”Buffer被不斷刷新。FFT應(yīng)用里通常使用循環(huán)模式這樣采集是無(wú)限連續(xù)的CPU在任何時(shí)刻去讀Buffer拿到的都是最近一段時(shí)間的數(shù)據(jù)。這就像一條自動(dòng)傳送帶不斷把工件送到倉(cāng)庫(kù)倉(cāng)庫(kù)管理員只需要在傳送帶堆滿一批時(shí)過(guò)去取走不需要每送一個(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ì)說(shuō)。2.2 Buffer長(zhǎng)度到底是“點(diǎn)數(shù)”還是“點(diǎn)數(shù)×通道數(shù)”這里最容易想歪很多人初次配置時(shí)以為“我要做1024點(diǎn)FFT那DMA搬運(yùn)次數(shù)就設(shè)成1024”。這是在單通道場(chǎng)景下的正確思路。但如果開(kāi)了ADC掃描模式比如同時(shí)采樣兩路信號(hào)DMA搬運(yùn)回來(lái)的數(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的傳輸長(zhǎng)度應(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。抽取出來(lái)的數(shù)據(jù)再分別裝入兩個(gè)float數(shù)組做FFT。我之前寫過(guò)一段提取邏輯大致套路如下#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]; }這段代碼雖然簡(jiǎn)單但它把“物理通道”和“FFT譜線”正確對(duì)應(yīng)起來(lái)了比一頭扎進(jìn)FFT函數(shù)里調(diào)參數(shù)更重要。2.3 雙緩沖與半傳輸中斷采集和運(yùn)算的相位錯(cuò)開(kāi)我剛才提到循環(huán)DMA模式會(huì)讓數(shù)據(jù)不斷刷新CPU在任意時(shí)刻讀取Buffer都可能讀到剛好被改寫一半的數(shù)據(jù)。解決這個(gè)問(wèn)題最經(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從頭重新開(kāi)始填充前半段同時(shí)觸發(fā)傳輸完成中斷CPU再去處理后半段數(shù)據(jù)。整個(gè)過(guò)程里CPU處理的永遠(yuǎn)是“上一次已經(jīng)填滿、且DMA暫時(shí)不會(huì)再寫”的那一半不會(huì)出現(xiàn)讀寫競(jìng)爭(zhēng)。在STM32 HAL庫(kù)中對(duì)應(yīng)的回調(diào)函數(shù)是HAL_ADC_ConvHalfCpltCallback和HAL_ADC_ConvCpltCallback。我在實(shí)際項(xiàng)目里用標(biāo)志位通知主循環(huán)處理不在中斷里直接跑FFT避免長(zhǎng)時(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輸入序列的完整性是長(zhǎng)時(shí)間穩(wěn)定運(yùn)行的關(guān)鍵否則跑半小時(shí)后FFT結(jié)果就會(huì)莫名跳變。3. CubeMX配置里的三個(gè)“埋雷點(diǎn)”搬一次就停、數(shù)值亂跳、數(shù)組沒(méi)對(duì)齊3.1 DMA模式選Normal還是Circular決定數(shù)據(jù)是“搬一次”還是“持續(xù)搬”CubeMX里配置ADC的DMA時(shí)大多數(shù)人會(huì)隨手選一個(gè)Normal模式。Normal模式的含義是DMA傳輸完設(shè)定長(zhǎng)度后傳輸自動(dòng)停止通道關(guān)閉。這時(shí)ADC雖然還在繼續(xù)轉(zhuǎn)換但轉(zhuǎn)換結(jié)果沒(méi)人搬了Buffer里的數(shù)據(jù)永遠(yuǎn)停留在第一次傳輸完的狀態(tài)?，F(xiàn)象就是串口打印第一次數(shù)據(jù)看起來(lái)正常之后打印出來(lái)的永遠(yuǎn)是同一批數(shù)據(jù)。解決方法是把DMA模式改成Circular。Circular模式下DMA傳輸計(jì)數(shù)遞減到0后自動(dòng)重裝初始值繼續(xù)下一次傳輸形成無(wú)限循環(huán)。這樣只要ADC不停轉(zhuǎn)換DMA就不停搬運(yùn)CPU隨時(shí)能拿到一份“剛剛更新過(guò)”的緩沖區(qū)數(shù)據(jù)。判斷是不是這個(gè)問(wèn)題的快速方法在調(diào)試器里看DMA的NDTR寄存器。如果它停在一個(gè)固定值不動(dòng)而ADC狀態(tài)寄存器還在繼續(xù)轉(zhuǎn)換基本就是Normal模式?jīng)]跑了。改完模式后重啟工程N(yùn)DTR應(yīng)該來(lái)回滾動(dòng)數(shù)據(jù)也會(huì)持續(xù)刷新。3.2 采樣時(shí)間別用最小檔源阻抗會(huì)讓ADC讀數(shù)“飄”STM32F407的ADC是逐次逼近型SAR結(jié)構(gòu)內(nèi)部有一個(gè)采樣保持電容。在采樣階段ADC引腳通過(guò)內(nèi)部模擬開(kāi)關(guān)給這個(gè)電容充電。如果開(kāi)關(guān)閉合時(shí)間太短而信號(hào)源的輸出阻抗又比較高比如前級(jí)是幾十kΩ的電阻分壓網(wǎng)絡(luò)電容還沒(méi)來(lái)得及充滿開(kāi)關(guān)就斷開(kāi)了轉(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í)沒(méi)問(wèn)題但如果前端是幾十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í)間變長(zhǎng)會(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位搬回來(lái)的數(shù)據(jù)就會(huì)錯(cuò)位Word會(huì)把相鄰寄存器的內(nèi)容一起讀進(jìn)來(lái)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庫(kù)的實(shí)數(shù)FFT函數(shù)要求輸入輸出數(shù)組必須互不重疊否則會(huì)破壞內(nèi)部狀態(tài)。4. 從ADC原始值到頻點(diǎn)幅值CMSIS-DSP的正確姿勢(shì)與幅度還原4.1 arm_rfft_fast_f32的輸入輸出布局別拿Matlab習(xí)慣來(lái)套ADC采樣得到的是一串實(shí)數(shù)序列。實(shí)數(shù)序列的FFT結(jié)果具有共軛對(duì)稱性所以CMSIS-DSP專門提供了arm_rfft_fast_f32來(lái)處理這種情況它比通用的復(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)換而來(lái) 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)序列的首尾都平滑過(guò)渡到接近0減弱這種非線性截?cái)?。最常用的是Hann窗它能把旁瓣壓得很低代價(jià)是主瓣寬度增加一倍頻率分辨能力略微下降。在實(shí)際代碼中減均值這一步往往比加窗還要優(yōu)先如果不減去DC分量FFT結(jié)果的0Hz處會(huì)有一個(gè)很大的直流譜線而它在Hann窗作用下還會(huì)向鄰近頻點(diǎn)泄漏把低頻段的小信號(hào)都淹沒(méi)了。完整預(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ǔ)回來(lái)。4.3 幅值換算單邊譜、Hann窗恢復(fù)因子與頻率分辨率做完FFT取模之后得到的|X[k]|并不直接等于信號(hào)的幅度。因?yàn)镃MSIS的FFT變換沒(méi)有做歸一化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。三部分合起來(lái)對(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里疊加無(wú)法分辨。想做更精細(xì)的頻率分析要么降低采樣率要么增加點(diǎn)數(shù)二者要按具體場(chǎng)景權(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í)測(cè)翻車記錄五個(gè)真實(shí)故障的完整排查鏈路5.1 DMA只搬了一次就停之后數(shù)據(jù)全是重復(fù)的第一段現(xiàn)象串口打印ADC數(shù)據(jù)第一次打印看起來(lái)正常再打印發(fā)現(xiàn)數(shù)值完全沒(méi)變化像死機(jī)一樣。用調(diào)試器查看DMA的NDTR寄存器發(fā)現(xiàn)它固定在一個(gè)值上不動(dòng)。排查過(guò)程我首先檢查HAL_ADC_Start_DMA是否被反復(fù)調(diào)用結(jié)果發(fā)現(xiàn)只調(diào)用了一次這個(gè)沒(méi)問(wèn)題。再查DMA初始化的模式參數(shù)發(fā)現(xiàn)CubeMX生成代碼里DMA_HandleTypeDef的Init.Mode字段是DMA_NORMAL問(wèn)題就出在這里。原因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ì)開(kāi)始循環(huán)變化數(shù)據(jù)恢復(fù)連續(xù)更新。這個(gè)案例看起來(lái)簡(jiǎn)單但它很典型很多“程序跑著跑著數(shù)據(jù)不動(dòng)了”的問(wèn)題不是程序邏輯死了而是外設(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í)。排查過(guò)程我做了兩個(gè)實(shí)驗(yàn)先直接把ADC引腳短接到GND跑FFT底噪依然很高說(shuō)明干擾在ADC通路內(nèi)部或參考電源上再把ADC引腳接到一個(gè)低噪基準(zhǔn)電壓源比如干凈的1.65V底噪降下來(lái)說(shuō)明問(wèn)題在輸入通路的前端或電源質(zhì)量上。最后在Vref引腳旁加強(qiáng)濾波電容并把該引腳的走線避開(kāi)數(shù)字信號(hào)線頻譜干凈了不少。這里要特別說(shuō)一下STM32F407的參考電壓。很多最小系統(tǒng)板上Vref直接接3.3V電源而這個(gè)3.3V同時(shí)給一堆數(shù)字芯片供電開(kāi)關(guān)噪聲全都會(huì)通過(guò)Vref耦合進(jìn)ADC轉(zhuǎn)換結(jié)果。FFT對(duì)噪聲極其敏感這種電源噪聲在時(shí)域看起來(lái)不明顯在頻域卻會(huì)表現(xiàn)為寬頻底噪或尖峰。在追求頻譜精度的場(chǎng)合至少要給Vref加一個(gè)高性能LDO和低ESR電容模擬地和數(shù)字地單點(diǎn)連接這個(gè)投資非常值得。5.3 頻率整體偏移了幾個(gè)bin信號(hào)源頻率明明很準(zhǔn)現(xiàn)象輸入1kHz正弦波FFT峰值出現(xiàn)在約1050Hz的位置不是1000Hz而且偏差比例基本固定。排查過(guò)程先確認(rèn)信號(hào)發(fā)生器輸出頻率足夠準(zhǔn)然后懷疑采樣率的實(shí)際值和理論值不一致。F407的ADC時(shí)鐘來(lái)自APB2總線經(jīng)過(guò)預(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í)鐘樹(shù)發(fā)現(xiàn)確實(shí)是ADC預(yù)分頻配置問(wèn)題。修正后主峰立刻回到1000Hz位置。這也引出一個(gè)很實(shí)用的標(biāo)定方法做頻率測(cè)量類產(chǎn)品時(shí)可以用一個(gè)已知精度的信號(hào)源輸入利用FFT峰值偏移反算實(shí)際采樣率再用軟件校準(zhǔn)系數(shù)修正比逐個(gè)寄存器排查快得多。5.4 單通道正常多通道數(shù)據(jù)串位現(xiàn)象兩個(gè)通道掃描模式采集ch0輸入的信號(hào)卻在ch1的頻譜里看到了或者兩個(gè)通道的頻譜混在一起。排查過(guò)程先確認(rèn)ADC掃描順序DMA搬運(yùn)的數(shù)據(jù)是按rank順序排列的ch0在ch1前。問(wèn)題出在提取數(shù)據(jù)時(shí)直接用了前1024個(gè)點(diǎn)沒(méi)有按通道數(shù)取模。由于DMA緩沖區(qū)是兩路交替排列前1024個(gè)點(diǎn)里其實(shí)混合了兩路數(shù)據(jù)直接當(dāng)單序列FFT結(jié)果當(dāng)然不對(duì)。按前面2.2小節(jié)的代碼重新抽取后兩路頻譜恢復(fù)正常。這類問(wèn)題在代碼里往往很難一眼看出來(lái)但只要在調(diào)試器里查看DMA緩沖區(qū)的排列規(guī)律很快就能定位。5.5 長(zhǎng)時(shí)間運(yùn)行后FFT結(jié)果突然亂掉甚至HardFault現(xiàn)象系統(tǒng)連續(xù)運(yùn)行幾分鐘后FFT結(jié)果莫名其妙變成一片噪聲偶爾直接進(jìn)入HardFault。排查過(guò)程先檢查內(nèi)存越界。FFT的輸入輸出數(shù)組長(zhǎng)度都是1024但如果手動(dòng)計(jì)算幅值的數(shù)組長(zhǎng)度只分配了512訪問(wèn)k1024的一半邊界時(shí)越界會(huì)覆蓋其他變量。排查之后發(fā)現(xiàn)真正的問(wèn)題在于CPU處理Buffer時(shí)DMA正在往同一個(gè)Buffer里寫新數(shù)據(jù)導(dǎo)致FFT讀取的數(shù)據(jù)被撕裂。解決方案就是用前面說(shuō)的雙緩沖乒乓結(jié)構(gòu)讓DMA和CPU錯(cuò)開(kāi)時(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)里如果混入了超過(guò)采樣率一半的頻率成分這些高頻成分會(huì)發(fā)生混疊折疊回低頻區(qū)形成“幽靈譜線”。軟件層面沒(méi)法有效消除混疊必須在ADC輸入前端加RC低通濾波器或者用運(yùn)放搭一個(gè)有源低通把超出fs/2的成分提前濾掉。我做音頻采樣時(shí)習(xí)慣在ADC引腳前放一個(gè)一階RC截止頻率設(shè)置在fs/2附近對(duì)高頻噪聲的壓制立竿見(jiàn)影。第二件事是考慮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算出來(lái)的幅度與實(shí)際物理量對(duì)不上。在產(chǎn)品化時(shí)通常用標(biāo)準(zhǔn)信號(hào)源輸入幾個(gè)已知幅度和頻率的點(diǎn)擬合出一條幅值校準(zhǔn)曲線在軟件里做補(bǔ)償。這個(gè)過(guò)程屬于長(zhǎng)期調(diào)校但一旦做完了你的FFT結(jié)果就不只是“看起來(lái)有譜”而是可以當(dāng)作測(cè)量數(shù)據(jù)用了。我現(xiàn)在的習(xí)慣是拿到任何STM32F407的FFT需求先把“時(shí)鐘樹(shù)到底怎么了”“DMA是循環(huán)還是單次”“FFT輸入是否對(duì)齊、是否加窗、是否減均值”這三件事查清楚再開(kāi)始寫應(yīng)用邏輯。這三步看著基礎(chǔ)卻決定了后面所有結(jié)果的可靠性。如果你也正在被ADC采出來(lái)的數(shù)據(jù)喂給FFT后頻譜亂糟糟的問(wèn)題折磨希望這篇記錄能讓你少走一點(diǎn)彎路。本文還有配套的精品資源點(diǎn)擊獲取