量儀:從PPG信號(hào)采集到算法實(shí)現(xiàn)全解析)
簡介本資源是一套完整的基于STM32F103的心率脈搏測(cè)量儀課程設(shè)計(jì)與畢業(yè)設(shè)計(jì)實(shí)現(xiàn)方案面向嵌入式初學(xué)者、電子類專業(yè)學(xué)生及單片機(jī)項(xiàng)目開發(fā)者解決便攜式生理信號(hào)采集、調(diào)理、處理與多端顯示的全流程工程實(shí)踐問題。壓縮包含1018個(gè)文件總大小24.97MB其中C源碼564個(gè)含信號(hào)濾波、ADC采樣、隊(duì)列緩存等核心算法、頭文件249個(gè)定義硬件接口與功能模塊、Keil工程文件uvprojx/uvoptx、Altium Designer原理圖與PCB源文件schdoc/pcbdoc、以及IAR/ARM數(shù)學(xué)庫靜態(tài)鏈接文件.a/.lib等覆蓋軟硬件協(xié)同開發(fā)全鏈路。已有93人學(xué)習(xí)下載資源提供可直接編譯運(yùn)行的完整工程包含OLED動(dòng)態(tài)顯示、藍(lán)牙無線傳輸、MATLAB波形聯(lián)調(diào)、心率異常告警及自動(dòng)休眠等實(shí)用功能代碼結(jié)構(gòu)清晰、注釋充分配套信號(hào)調(diào)理電路設(shè)計(jì)說明與時(shí)間間隔平均值算法實(shí)現(xiàn)細(xì)節(jié)具備強(qiáng)復(fù)現(xiàn)性與教學(xué)參考價(jià)值。1. 項(xiàng)目拆解心率脈搏測(cè)量儀到底要做什么先把這個(gè)項(xiàng)目說透?;赟TM32設(shè)計(jì)的心率脈搏測(cè)量儀本質(zhì)上是一個(gè)以單片機(jī)為核心、通過傳感器采集人體脈搏信號(hào)經(jīng)算法處理后實(shí)時(shí)顯示心率數(shù)值的嵌入式系統(tǒng)。這類項(xiàng)目在課程設(shè)計(jì)、電子設(shè)計(jì)競賽、個(gè)人DIY里都非常常見因?yàn)樗募夹g(shù)棧覆蓋了信號(hào)采集、模擬前端、ADC采樣、數(shù)字濾波、算法識(shí)別、顯示交互、低功耗設(shè)計(jì)等多個(gè)嵌入式核心知識(shí)點(diǎn)而且軟硬件都有做完之后能拿得出手的東西很多。拿到這套“系統(tǒng)資料”時(shí)大多數(shù)人第一反應(yīng)是去看原理圖和代碼但我建議你先反過來——先把需求拆清楚再回來看資料你才會(huì)知道每份文件為什么要這么設(shè)計(jì)。這個(gè)系統(tǒng)要解決的核心問題有三個(gè)準(zhǔn)確采集到脈搏信號(hào)人體脈搏信號(hào)是微伏到毫伏級(jí)別的微弱生理信號(hào)且混有大量噪聲不像讀按鍵電平那么簡單必須經(jīng)過放大和濾波。從信號(hào)中計(jì)算出心率要么檢測(cè)時(shí)域波形的峰值間隔要么做頻域分析最終得出每分鐘心跳次數(shù)單位是bpmbeats per minute。以直觀方式呈現(xiàn)結(jié)果OLED屏幕實(shí)時(shí)顯示脈搏波形和心率數(shù)值必要時(shí)還能通過按鍵切換功能、設(shè)置報(bào)警閾值甚至通過串口把數(shù)據(jù)發(fā)給上位機(jī)做進(jìn)一步分析。所以整份系統(tǒng)資料的架構(gòu)就很清晰了傳感器模塊采集信號(hào) → 模擬前端調(diào)理 → STM32內(nèi)置ADC采樣 → 數(shù)字濾波與算法識(shí)別 → OLED顯示/串口輸出。這是一條非常標(biāo)準(zhǔn)的生物醫(yī)學(xué)信號(hào)采集處理鏈路掌握了這套流程今后做血氧儀、心電儀、肌電儀本質(zhì)上都是同一套思路只是傳感器和算法不同。1.1 心率測(cè)量原理三種主流方案你怎么選在動(dòng)手之前必須先理解心率是怎么測(cè)出來的。目前主流的方式有三種方案原理傳感器類型成本算法復(fù)雜度適合場景光電容積脈搏波PPG用LED照射皮膚光電二極管接收透射/反射光強(qiáng)變化血液脈動(dòng)引起光強(qiáng)變化如MAX30102、PulseSensor低中DIY、可穿戴設(shè)備心電法ECG測(cè)量心臟電活動(dòng)引起體表電位差R波峰值間隔即心搏周期AD8232等模擬前端中中高醫(yī)療級(jí)、運(yùn)動(dòng)監(jiān)測(cè)壓電/應(yīng)變式脈搏傳感通過手腕處壓力變化感知脈動(dòng)壓電薄膜、應(yīng)變片中低中傳統(tǒng)中醫(yī)脈診儀等這套系統(tǒng)資料里如果是基于PulseSensor或者M(jìn)AX30102的走的都是PPG方案。我自己做這個(gè)項(xiàng)目時(shí)也優(yōu)先推薦PPG原因是它結(jié)構(gòu)簡單、傳感器模塊成熟、上手快。PulseSensor本質(zhì)上是一個(gè)集成了綠光LED和環(huán)境光消除電路的模擬器件直接輸出一個(gè)與脈搏相關(guān)的電壓信號(hào)接在ADC引腳上就能采。MAX30102則更進(jìn)一步把LED、光電檢測(cè)、AD轉(zhuǎn)換、濾波都集成到一個(gè)芯片里通過I2C接口和主控通信可以做到更多配置。你要是只想要一個(gè)簡單能跑通的系統(tǒng)PulseSensor是首選你要是想做得更接近現(xiàn)代可穿戴設(shè)備MAX30102更合適。兩者的核心都是PPG但數(shù)據(jù)讀取方式和處理邏輯差別很大下面會(huì)展開講。1.2 系統(tǒng)架構(gòu)與軟硬件職責(zé)劃分從頂層看系統(tǒng)的模塊劃分如下傳感器輸入模塊采集脈搏信號(hào)輸出模擬電壓或數(shù)字量。信號(hào)調(diào)理模塊對(duì)模擬信號(hào)做隔直、放大、濾波讓信號(hào)幅度落在ADC的最佳輸入范圍。主控模塊STM32負(fù)責(zé)ADC采樣、數(shù)據(jù)緩存、算法處理、顯示驅(qū)動(dòng)、按鍵掃描。人機(jī)交互模塊OLED顯示屏、按鍵、蜂鳴器可選。電源模塊為所有模塊提供穩(wěn)定工作電壓。這套架構(gòu)里最容易翻車的往往不是單片機(jī)代碼而是模擬部分的信號(hào)調(diào)理。很多人PCB回來了程序調(diào)通了傳感器也讀數(shù)了但心率就是測(cè)不準(zhǔn)——波形糊成一片或者數(shù)值亂跳問題多半出在模擬前端的設(shè)計(jì)而不是算法。2. 主控選型與核心資源分配為什么STM32是這塊板子的最優(yōu)解做這個(gè)項(xiàng)目入門者可能會(huì)用51、Arduino進(jìn)階一點(diǎn)的用ESP32但STM32始終是這個(gè)題目最常出現(xiàn)的主角。原因有三ADC性能足夠STM32的12位逐次逼近型ADC采樣率最高能做到幾Msps對(duì)心率信號(hào)帶寬只有0.5Hz~5Hz左右來說綽綽有余而且多通道可以同時(shí)采集速度、溫度等其他信號(hào)方便擴(kuò)展。DMA和定時(shí)器聯(lián)動(dòng)方便可以配置成定時(shí)器觸發(fā)ADC采樣、DMA自動(dòng)搬運(yùn)數(shù)據(jù)到內(nèi)存CPU全程不參與采樣過程中還能同時(shí)做按鍵掃描和顯示刷新這套機(jī)制對(duì)實(shí)時(shí)性要求較高的生理信號(hào)采集非常關(guān)鍵。資料和生態(tài)成熟標(biāo)準(zhǔn)庫、HAL庫、寄存器版本隨便選網(wǎng)上工程模板一大堆遇到問題很容易檢索到解決方法。很多人糾結(jié)該用哪個(gè)型號(hào)。如果你看的是系統(tǒng)資料里的BOM表大概率是STM32F103C8T6藍(lán)色Pill板最常見或STM32F407系列。以STM32F103C8T6為例它擁有64KB Flash、20KB SRAM跑一個(gè)心率和展示程序綽綽有余即使加上UCOS或者裸機(jī)狀態(tài)機(jī)也夠用。3個(gè)USART一個(gè)可以連藍(lán)牙模塊一個(gè)接PC上位機(jī)調(diào)試。2個(gè)I2C如果選用MAX30102I2C直接掛載。10個(gè)ADC通道單通道定時(shí)采樣完全夠用多通道還能擴(kuò)展測(cè)量體溫。2.1 選芯片前先估算資源占用我自己的習(xí)慣是寫代碼前先估算資源避免寫到一半發(fā)現(xiàn)Flash或SRAM不夠。以脈率檢測(cè)為例采樣率設(shè)為100Hz只測(cè)0.5~5Hz信號(hào)做20倍過采樣抗混疊效果比50Hz更穩(wěn)。每次采樣產(chǎn)生一個(gè)uint16_t如果做滑動(dòng)窗口緩存10秒數(shù)據(jù)所需內(nèi)存 100Hz × 10s × 2字節(jié) 2KB占20KB的10%沒問題。如果做FFT頻域分析取1024點(diǎn)FFT復(fù)數(shù)數(shù)組需要2×1024×4字節(jié) 8KB加上其他緩存也還是夠的。時(shí)域峰值檢測(cè)算法則省內(nèi)存得多只用幾個(gè)變量和一個(gè)小數(shù)組幾乎可以忽略不計(jì)。所以STM32F103C8T6的內(nèi)存完全夠用選擇它做設(shè)計(jì)不是硬件性能的妥協(xié)而是性價(jià)比的合理匹配。2.2 選型對(duì)比PulseSensor還是MAX30102繼續(xù)細(xì)化傳感器選型這里我做了一個(gè)對(duì)比表方便你參考對(duì)比項(xiàng)PulseSensorMAX30102輸出接口模擬電壓0~3.3VI2C數(shù)字輸出信號(hào)內(nèi)容原始PPG波形經(jīng)前端放大濾波內(nèi)部ADC轉(zhuǎn)換后的紅外/紅光數(shù)據(jù)供電電壓3.3V~5V典型3.3V1.8V~3.3V驅(qū)動(dòng)難度簡單ADC直接采集中需要配置寄存器、讀FIFO抗干擾中易受環(huán)境光和運(yùn)動(dòng)干擾較好有多級(jí)濾波和掩碼設(shè)置價(jià)格十幾元幾十元但有些模塊可能虛標(biāo)適合階段學(xué)原理、快速驗(yàn)證做產(chǎn)品原型、深入調(diào)優(yōu)這里要特別提醒一點(diǎn)MAX30102的模塊現(xiàn)在淘寶上很多但有相當(dāng)一部分是翻新芯片或者傳感器臟污的測(cè)量的紅外數(shù)據(jù)會(huì)出現(xiàn)明顯周期性噪聲而且個(gè)別模塊還貼錯(cuò)了電阻導(dǎo)致I2C地址不一致。購買時(shí)盡量選正規(guī)渠道買到之后先用邏輯分析儀或者示波器看一下I2C時(shí)序是否正確別一上來就套用網(wǎng)絡(luò)上的驅(qū)動(dòng)。為什么我用過之后還是建議初次上手先從PulseSensor開始因?yàn)樗哪M輸出讓你能直接看到PPG波形長什么樣這對(duì)你理解后續(xù)的濾波和算法是巨大幫助。示波器放在傳感器輸出端能看到一個(gè)隨心跳起伏的脈沖波形它幅值大、形態(tài)清晰你就能知道算法該往哪個(gè)方向調(diào)。而MAX30102的輸出是FIFO里的數(shù)據(jù)雖然也能調(diào)但你少了一層“看得見摸得著”的直覺。3. 硬件設(shè)計(jì)與信號(hào)調(diào)理的關(guān)鍵細(xì)節(jié)這一部分是目前絕大多數(shù)教程和資料里講得最水、但恰恰是最容易出問題的環(huán)節(jié)。如果你的PCB板從模塊飛線出來的信號(hào)很差下面這些細(xì)節(jié)一定要重新過一遍。3.1 模擬前端從傳感器到ADC之間經(jīng)歷了什么以PulseSensor模塊為例它內(nèi)部其實(shí)已經(jīng)做了一部分信號(hào)調(diào)理板上集成了一塊運(yùn)算放大器和濾波網(wǎng)絡(luò)輸出的電壓信號(hào)大致在1/2 VCC附近典型值1.65V疊加一個(gè)和脈搏同步的交流分量。但這個(gè)交流分量幅度很小大約只有幾十到幾百毫伏而且疊加了較大的直流偏置。如果你直接用STM32的ADC去采這個(gè)信號(hào)會(huì)面臨兩個(gè)問題直流偏置占用ADC動(dòng)態(tài)范圍STM32的ADC輸入范圍是0~3.3V如果信號(hào)偏置在1.65V交流幅值只有200mV那ADC的有效分辨范圍只用了大約12%意味著你的12位ADC實(shí)際有效位數(shù)下降到了約9~10位采樣精度損失不小。噪聲混疊工頻50Hz干擾、LED驅(qū)動(dòng)紋波、電磁干擾都會(huì)進(jìn)入測(cè)量鏈路如果不做濾波波形上會(huì)出現(xiàn)明顯的毛刺。所以推薦的做法是在STM32 ADC前端增加一級(jí)有源二階低通濾波器或者至少用RC低通截止頻率設(shè)置為5~10Hz左右同時(shí)用隔直電容去掉直流偏置再用運(yùn)放把交流信號(hào)放大到合適的幅度。具體參數(shù)我給出一種常見實(shí)現(xiàn)隔直電容10μF耦合傳感器信號(hào)。偏置電阻兩個(gè)10kΩ電阻分壓把運(yùn)放同相端偏置到1.65V。運(yùn)放LM358或MCP6002單電源供電工作在放大狀態(tài)。低通濾波用Sallen-Key拓?fù)浣刂诡l率 f_c 1 / (2πRC)如果R16kΩ、C1μF算出來約10Hz。當(dāng)然實(shí)際項(xiàng)目中很多人省略這級(jí)電路直接把PulseSensor輸出接到ADC。項(xiàng)目能跑但你得到的波形噪聲會(huì)大不少后續(xù)算法需要做大量濾波處理才能穩(wěn)定識(shí)別。我建議系統(tǒng)資料里如果有這部分說明一定要讀透如果資料里沒有你自己加一級(jí)也非常有必要。3.2 電源和PCB布局為什么心電信號(hào)老是跳心率測(cè)量系統(tǒng)對(duì)電源紋波和地線噪聲極其敏感。我這里不是玄學(xué)是實(shí)際踩坑總結(jié)出來的規(guī)律。首先ADC的參考電壓必須穩(wěn)定。STM32F103的VREF通常直接和VDD相連如果你的3.3V穩(wěn)壓輸出紋波大ADC的滿量程基準(zhǔn)就會(huì)跟著抖采樣數(shù)值就開始“漂”。建議在電源輸入端加鉭電容10μF以上和陶瓷電容0.1μF的組合形成低阻抗去耦網(wǎng)絡(luò)。其次是地線問題。傳感器信號(hào)屬于微弱模擬信號(hào)顯示器、蜂鳴器、按鍵這些數(shù)字外設(shè)的開關(guān)噪聲會(huì)通過地線串?dāng)_到模擬部分。推薦的PCB布線原則是單點(diǎn)接地模擬地和數(shù)字地在主電容處單點(diǎn)匯合不要大面積鋪銅直接連通。信號(hào)走線盡量遠(yuǎn)離高頻數(shù)字線尤其是OLED的I2C或SPI線不要和傳感器信號(hào)線平行走長距離。傳感器模塊和主控板之間用屏蔽線或短而粗的跳線如果使用杜邦線飛線盡量控制在10cm以內(nèi)。有朋友說我做的是雙面板就幾個(gè)元器件不做模擬地和數(shù)字地分區(qū)也沒事吧我一開始也是這么想的后來發(fā)現(xiàn)OLED刷新屏幕時(shí)波形上就出現(xiàn)一系列規(guī)律毛刺一查就是數(shù)字部分的瞬態(tài)電流通過公共地線耦合到模擬端了。后來把模擬前端獨(dú)立成一小塊用0Ω電阻做單點(diǎn)互聯(lián)問題立刻改善。4. 固件開發(fā)從ADC采樣到心率算法的完整鏈路固件開發(fā)是整個(gè)系統(tǒng)里工作量最大的部分也是系統(tǒng)資料中代碼文件最密集的地方。這里我不會(huì)照搬某一套代碼而是把核心邏輯和設(shè)計(jì)思路講清楚讓你拿到任何一份資料都能快速看懂。4.1 ADC采樣策略定時(shí)器觸發(fā)DMA搬運(yùn)讓CPU做更重要的活心率采樣對(duì)時(shí)序的要求是“每隔固定間隔采一個(gè)點(diǎn)”而不是“能采多少采多少”。如果用HAL庫里的HAL_ADC_GetValue()函數(shù)在循環(huán)里反復(fù)讀你會(huì)發(fā)現(xiàn)問題循環(huán)里還有OLED刷新和按鍵掃描執(zhí)行時(shí)間不固定導(dǎo)致采樣間隔抖動(dòng)波形時(shí)間軸變形。高頻采樣時(shí)反復(fù)軟觸發(fā)ADCCPU占用率高系統(tǒng)響應(yīng)變慢。所以我推薦的方案是定時(shí)器觸發(fā)ADC采樣 DMA自動(dòng)搬運(yùn)數(shù)據(jù)到內(nèi)存緩沖區(qū)。這樣配置好之后ADC完全由硬件定時(shí)觸發(fā)以穩(wěn)定的100Hz或200Hz采樣率持續(xù)工作DMA自動(dòng)把ADC轉(zhuǎn)換結(jié)果搬進(jìn)RAM數(shù)組CPU全程不參與只有在緩沖區(qū)滿或需要計(jì)算時(shí)才去取數(shù)據(jù)。具體思路如下以STM32F103 HAL庫為例配置一個(gè)定時(shí)器例如TIM2將其更新事件作為ADC觸發(fā)源設(shè)置周期為10ms對(duì)應(yīng)100Hz采樣率。配置ADC1的通道比如ADC1_IN1對(duì)應(yīng)PA1設(shè)置為定時(shí)器觸發(fā)模式分辨率12位單通道。配置DMA1通道1這里以ADC1的DMA輸出為例把ADC轉(zhuǎn)換值持續(xù)搬運(yùn)到uint16_t adc_buffer[512]數(shù)組。在main函數(shù)中使能ADC的DMA中斷或環(huán)形緩沖半滿中斷當(dāng)數(shù)據(jù)積累到一定長度后啟動(dòng)心率計(jì)算。這里插一句如果你用的是系統(tǒng)資料里的標(biāo)準(zhǔn)庫版本操作步驟是相同的只是API不同。不要糾結(jié)標(biāo)準(zhǔn)庫還是HAL庫關(guān)鍵是理解背后的硬件鏈路。用標(biāo)準(zhǔn)庫的話代碼量更少更直觀用HAL庫的話配置更規(guī)范便于移植到其他STM32型號(hào)。4.2 數(shù)字濾波先平滑再識(shí)別否則算法再牛也沒用經(jīng)過ADC采樣得到的原始PPG信號(hào)依然包含高頻噪聲、基線漂移和運(yùn)動(dòng)偽跡。如果直接對(duì)這個(gè)信號(hào)做峰值檢測(cè)結(jié)果會(huì)非常不穩(wěn)定心率可能一會(huì)兒60一會(huì)兒120。所以在算法之前一定要做數(shù)字濾波。常用的三種濾波手段滑動(dòng)平均濾波移動(dòng)平均對(duì)最近N個(gè)點(diǎn)求平均作為當(dāng)前輸出值。比如N5在100Hz采樣率下就是50ms的滑動(dòng)窗口能有效濾除高頻毛刺但對(duì)基線漂移無效。這是最簡單、最不容易出錯(cuò)的濾波方式。IIR低通濾波比如一階低通 y[n] a * x[n] (1 - a) * y[n-1]a取0.1~0.2相當(dāng)于低通截止頻率較低。計(jì)算簡單適合MCU實(shí)時(shí)處理。中值濾波對(duì)滑動(dòng)窗口排序取中值對(duì)脈沖噪聲比如突然的尖峰有很好的抑制作用但計(jì)算量稍大。我實(shí)際使用下來的組合是先做一階高通去除基線漂移再做滑動(dòng)平均去毛刺最后用一個(gè)30Hz左右的低通。當(dāng)然這個(gè)組合的具體系數(shù)需要根據(jù)你的采樣率調(diào)整。寫代碼時(shí)注意一點(diǎn)濾波器初始化時(shí)濾波器的狀態(tài)變量比如前一次的y值要先初始化為第一次采樣的值不要讓初始狀態(tài)從0開始否則前一兩秒數(shù)據(jù)會(huì)出現(xiàn)明顯跳變。4.3 心率算法時(shí)域峰值檢測(cè)的完整步驟心率測(cè)量的核心是從PPG波形中識(shí)別每一次心搏。業(yè)界通用的方法主要有兩類時(shí)域峰值檢測(cè)和頻域功率譜分析。我做這個(gè)項(xiàng)目時(shí)最穩(wěn)定、最容易調(diào)試的是改進(jìn)的時(shí)域峰值檢測(cè)法配合自適應(yīng)閾值。下面講一下它的具體流程這套流程在系統(tǒng)資料里通常也會(huì)出現(xiàn)只是可能寫得比較粗糙。第一步滑動(dòng)窗口劃分。設(shè)窗口長度10秒每1秒滑動(dòng)一次每次取最新10秒數(shù)據(jù)做分析。第二步計(jì)算閾值。對(duì)窗口內(nèi)數(shù)據(jù)求最大值maxVal和最小值minVal初始閾值threshold (maxVal minVal) / 2。但如果窗口內(nèi)有劇烈運(yùn)動(dòng)尖峰直接取中點(diǎn)可能誤判。所以用“自適應(yīng)系數(shù)”調(diào)節(jié)比如threshold minVal 0.5 * (maxVal - minVal)也可以引入前幾秒的閾值做加權(quán)平均讓閾值變化更平滑。第三步峰值判定。遍歷窗口內(nèi)每個(gè)采樣點(diǎn)如果當(dāng)前點(diǎn)的值大于閾值且大于前一個(gè)點(diǎn)且大于后一個(gè)點(diǎn)就認(rèn)為這是一個(gè)局部峰。但要考慮到PPG波形在心跳周期內(nèi)可能有重搏波dicrotic notch即第一峰值后可能還有一個(gè)較小的隆起。如果你把重搏波也當(dāng)成一個(gè)心搏心率就會(huì)翻倍。解決辦法是設(shè)置最小峰間距refractory period比如0.3秒也就是兩個(gè)峰之間至少間隔300ms對(duì)應(yīng)最大心率200bpm在局部峰判定時(shí)如果距離上一次有效峰不足300ms就丟棄。第四步計(jì)算心率。統(tǒng)計(jì)窗口內(nèi)有效峰個(gè)數(shù)peakCount那么HR peakCount × 60 / windowSeconds。例如10秒窗口內(nèi)檢測(cè)到12個(gè)峰HR 12 × 6 72bpm。第五步平滑輸出。為了防止顯示數(shù)值跳變太突兀我常用指數(shù)平滑displayHR alpha * currentHR (1 - alpha) * lastDisplayHRalpha取0.3~0.5。這套算法是我做這個(gè)項(xiàng)目時(shí)反復(fù)調(diào)出來的結(jié)果。需要注意的是不同人的脈搏波形形態(tài)差異很大手指按壓力度、皮膚顏色、環(huán)境光都會(huì)影響波形形態(tài)所以閾值不能是一成不變的固定值必須自適應(yīng)。自適應(yīng)閾值的實(shí)現(xiàn)細(xì)節(jié)往往決定了一套資料寫的程序質(zhì)量高不高。4.4 MAX30102和PulseSensor的算法差異如果你用的是MAX30102原理還是一樣的但有兩個(gè)關(guān)鍵差異MAX30102輸出的是紅外IR通道和紅光Red通道兩組數(shù)據(jù)每組又分多個(gè)數(shù)據(jù)位。通常用IR通道做心率檢測(cè)因?yàn)榧t外光受血氧飽和度變化影響小信號(hào)更穩(wěn)定。MAX30102內(nèi)部自帶環(huán)境光消除、信號(hào)放大和低通濾波但它的原始FIFO數(shù)據(jù)依然存在基線漂移和運(yùn)動(dòng)偽跡所以你仍然要做數(shù)字濾波和峰值檢測(cè)。另外MAX30102有個(gè)重要配置是采樣率和ADC范圍。采樣率建議設(shè)置為100Hz或200HzADC范圍設(shè)為409612位即可。如果兩個(gè)通道都開數(shù)據(jù)量會(huì)翻倍注意FIFO溢出問題。說到這就得提一個(gè)坑MAX30102的FIFO溢出標(biāo)志。如果你讀FIFO的速度跟不上寫入速度數(shù)據(jù)就會(huì)丟失心跳波形會(huì)出現(xiàn)斷層。因此要確保在主循環(huán)中及時(shí)讀取FIFO數(shù)據(jù)或者開啟中斷通過INT引腳通知主控。5. 顯示與交互OLED屏、按鍵和附加功能怎么設(shè)計(jì)一般系統(tǒng)資料的硬件設(shè)計(jì)都會(huì)配置一個(gè)I2C接口的OLED顯示屏常見0.96寸SSD1306驅(qū)動(dòng)外加兩個(gè)按鍵還可能有一個(gè)蜂鳴器或LED指示燈。這部分雖然不復(fù)雜但如果設(shè)計(jì)不合理會(huì)嚴(yán)重拖累系統(tǒng)體驗(yàn)。5.1 OLED顯示刷新策略不要阻塞在主循環(huán)0.96寸OLED分辨率128x64SSD1306驅(qū)動(dòng)通過I2C接口通信I2C速度一般設(shè)置在400kHz。很多人寫OLED驅(qū)動(dòng)時(shí)直接用HAL_I2C_Mem_Write一個(gè)像素一個(gè)像素地刷全屏刷新一次可能需要幾十毫秒甚至上百毫秒。如果主循環(huán)里還帶著心率算法20ms刷一次屏CPU基本就卡在I2C上了。我建議的顯示策略分區(qū)域刷新只刷新數(shù)值變化的地方心率數(shù)值、時(shí)間、狀態(tài)圖標(biāo)波形區(qū)域可以整體繪制但用局部更新代替全屏清屏。使用顯存緩沖在內(nèi)存中維護(hù)一個(gè)128x8字節(jié)1KB的顯存數(shù)組所有繪圖操作都在數(shù)組上做最后一次性刷新到屏上。1KB的SRAM對(duì)STM32F103來說完全能承受。主循環(huán)中設(shè)置刷新頻率OLED刷新限制在20~30Hz約30~50ms刷新一次不要和ADC采樣同步避免刷屏?xí)r影響采樣時(shí)序。實(shí)際操作時(shí)你會(huì)發(fā)現(xiàn)把OLED刷屏和ADC采樣放在同一個(gè)中斷優(yōu)先級(jí)里會(huì)出問題。我建議ADC的DMA中斷優(yōu)先級(jí)高于OLED刷新的任何操作OLED刷新放在主循環(huán)中通過標(biāo)志位控制。這樣采樣不會(huì)因?yàn)樗⑵帘谎舆t。5.2 按鍵與交互邏輯設(shè)計(jì)常見的交互功能有兩個(gè)一個(gè)是短按切換顯示模式波形心率、只有心率、歷史記錄另一個(gè)是長按開關(guān)蜂鳴器報(bào)警。按鍵掃描可以在ADC中斷里做也可以在主循環(huán)里用定時(shí)器掃描。要避免的是按鍵消抖和長按判定邏輯寫得太復(fù)雜導(dǎo)致主循環(huán)卡頓。一個(gè)實(shí)用技巧是每10ms掃描一次按鍵讀取引腳電平。維護(hù)一個(gè)按鍵狀態(tài)變量連續(xù)N次比如5次讀到同一電平才判定按鍵有效實(shí)現(xiàn)消抖。長按檢測(cè)則記錄按下開始時(shí)間超過1秒認(rèn)為長按。這套代碼邏輯非常簡單但非??煽?。我在實(shí)際測(cè)試中相比單獨(dú)用delay延時(shí)消抖這套非阻塞方式能保證按鍵響應(yīng)時(shí)不影響采樣和其他任務(wù)。5.3 擴(kuò)展功能蜂鳴器報(bào)警與上位機(jī)串口輸出不少系統(tǒng)資料還會(huì)加一個(gè)蜂鳴器報(bào)警功能用于心率過高或過低時(shí)提醒。這個(gè)功能在代碼層面很簡單只要在算法計(jì)算完心率后再做一個(gè)閾值判斷就行。但要注意心率閾值設(shè)置不能太窄否則稍一活動(dòng)就會(huì)被誤報(bào)。我做過測(cè)試靜息狀態(tài)下心率在70附近但手輕輕動(dòng)一下PPG信號(hào)產(chǎn)生的運(yùn)動(dòng)偽跡可能讓算法短暫誤判到150以上如果你閾值設(shè)在120就會(huì)誤報(bào)。解決方法是對(duì)心率結(jié)果再做一次合理性校驗(yàn)比如連續(xù)3次計(jì)算結(jié)果都在120以上才觸發(fā)報(bào)警且報(bào)警持續(xù)一定時(shí)間才發(fā)聲避免單個(gè)波動(dòng)觸發(fā)。串口輸出則更推薦加上方便調(diào)試。在開發(fā)時(shí)STM32通過串口把原始采樣波形和計(jì)算出的心率值發(fā)送給電腦串口助手立刻就能看到波形是否正常。這里有個(gè)小技巧發(fā)送數(shù)據(jù)的波特率用921600或460800避免波特率太低導(dǎo)致的數(shù)據(jù)傳輸瓶頸。當(dāng)然串口發(fā)送的字符盡量精簡比如用“123,72\n”這種格式方便上位機(jī)解析。6. 從資料包到實(shí)物我的整合經(jīng)驗(yàn)與避坑清單當(dāng)你拿到一套“系統(tǒng)資料”時(shí)里面通常有原理圖、PCB、源代碼、文檔和工程文件但有時(shí)候你會(huì)發(fā)現(xiàn)資料里的代碼跟你手里買的傳感器版本不匹配、原理圖版本號(hào)不統(tǒng)一甚至引腳定義有沖突。這時(shí)候最忌諱的就是直接燒錄代碼然后看效果一旦現(xiàn)象不對(duì)就會(huì)在幾個(gè)模塊之間來回猜非常浪費(fèi)時(shí)間。6.1 拿到資料后的第一件事核對(duì)引腳和版本我拿到這類資料后第一件事是先整理出一張“資源映射表”把STM32的每個(gè)引腳對(duì)應(yīng)哪個(gè)外設(shè)、哪個(gè)傳感器通道、哪個(gè)LED、哪個(gè)按鍵全部列出來。比如STM32引腳功能連接外設(shè)PA0ADC1_IN0PulseSensor輸出PB6I2C1_SCLOLEDPB7I2C1_SDAOLEDPA1GPIO_Input按鍵1PA2GPIO_Output蜂鳴器PA9/PA10USART1調(diào)試串口列出這張表后再對(duì)照原理圖和代碼里的引腳定義如果發(fā)現(xiàn)沖突立即用宏定義改掉而不是在程序里到處打補(bǔ)丁。這一步能節(jié)省你至少半天的排錯(cuò)時(shí)間。還有一點(diǎn)很多資料代碼是用標(biāo)準(zhǔn)庫寫的而你當(dāng)前的開發(fā)環(huán)境可能是HAL庫這時(shí)不要強(qiáng)行移植可以先跑通原有工程理解邏輯后再換自己熟悉的庫。6.2 調(diào)試過程的常見坑下面是我實(shí)際調(diào)試中遇到幾次的坑也是群里大家問得最多的問題。讀到的ADC值一直為0: 檢查ADC引腳是否被復(fù)用成了別的功能比如PB0/PC13可能有JTAG或RTC功能沖突還有檢查DMA配置和ADC掃描模式是否沖突。大部分人最終查出來是沒開ADC時(shí)鐘或者DMA通道配錯(cuò)。波形有規(guī)律的尖刺: 大概率是OLED刷新或串口發(fā)送引入了干擾注意電源去耦或者把顯示刷新頻率降低一點(diǎn)。心率突然翻倍或減半: 一定是重搏波誤檢或漏檢。回頭檢查最小峰間距設(shè)沒設(shè)閾值是不是沒有自適應(yīng)。手指按上去沒信號(hào): 一是傳感器位置不對(duì)PPG傳感器和指尖接觸要穩(wěn)定二是傳感器供電電壓設(shè)計(jì)錯(cuò)誤部分PulseSensor需要5V供電才能有好的信號(hào)幅度。程序進(jìn)了某個(gè)函數(shù)就卡死: 檢查是否有死循環(huán)等待某個(gè)外設(shè)標(biāo)志位特別是I2C設(shè)備沒應(yīng)答時(shí)會(huì)卡在HAL_I2C_Mem_Read里。給I2C加個(gè)超時(shí)是經(jīng)驗(yàn)之談。6.3 性能評(píng)估與優(yōu)化方向項(xiàng)目做出來如果能顯示心率、能報(bào)警、能顯示波形已經(jīng)達(dá)到了基本要求。但如果繼續(xù)優(yōu)化我會(huì)建議你按這個(gè)順序來優(yōu)化采樣率與濾波器系數(shù)通過實(shí)驗(yàn)找到適合自己傳感器的濾波參數(shù)把波形調(diào)得更干凈心率計(jì)算的穩(wěn)定性會(huì)明顯改善。算法升級(jí)將時(shí)域峰值檢測(cè)替換為基于小波變換或自適應(yīng)閾值結(jié)合模板匹配的方法對(duì)運(yùn)動(dòng)偽跡的魯棒性會(huì)更好。降低功耗如果后面裝電池需要用STM32的睡眠模式只在采樣時(shí)刻喚醒MAX30102可以降低采樣率或關(guān)閉LED從而降低功耗。增加數(shù)據(jù)記錄把心率數(shù)據(jù)存到Flash芯片中做成可回看的日志。6.4 個(gè)人體會(huì)與建議最后分享一點(diǎn)我個(gè)人的整體感受。這個(gè)項(xiàng)目好上手但真正做好、做穩(wěn)涉及的知識(shí)其實(shí)一點(diǎn)也不少。如果你只是照著資料燒錄代碼、看串口輸出那學(xué)到的東西有限。我建議你至少自己改三處把心率算法從頭寫一遍哪怕剛開始實(shí)現(xiàn)得比資料里的丑但寫一遍你才真正理解峰值檢測(cè)存在的那些坑。用示波器看一次傳感器輸出的原始信號(hào)再對(duì)比經(jīng)過你濾波程序之后的波形你會(huì)對(duì)數(shù)字濾波有極其直觀的印象。打斷一下自己的思維試著用MAX30102替換掉PulseSensor哪怕只是移植驅(qū)動(dòng)你也會(huì)把I2C、FIFO、寄存器配置這些知識(shí)點(diǎn)融會(huì)貫通。嵌入式系統(tǒng)里沒有“純軟件項(xiàng)目”和“純硬件項(xiàng)目”心率脈搏測(cè)量儀恰恰是把兩邊都串起來的那個(gè)典型。你把這套東西消化透了后面再做其他傳感器項(xiàng)目尤其是模擬信號(hào)處理類的會(huì)輕松很多。本文還有配套的精品資源點(diǎn)擊獲取