動(dòng)開發(fā)實(shí)戰(zhàn))
前陣子做健康監(jiān)測小項(xiàng)目我在網(wǎng)購模塊里挑了一片 MAX30100 血氧心跳傳感器以為接上 I2C、讀幾個(gè)寄存器就能把數(shù)據(jù)跑起來結(jié)果從 OpenHarmony 側(cè)驅(qū)動(dòng)注冊(cè)到 PPG 信號(hào)穩(wěn)定輸出前后折騰了不少時(shí)間。這個(gè)過程中踩過的坑、驗(yàn)證過的方法整理出來就是這篇文章。內(nèi)容覆蓋 MAX30100 的測量原理、HDF 驅(qū)動(dòng)框架下的驅(qū)動(dòng)實(shí)現(xiàn)、FIFO 數(shù)據(jù)讀取與緩沖設(shè)計(jì)、血氧和心率算法的邊界以及最終調(diào)試驗(yàn)證時(shí)遇到的一堆實(shí)際問題。適合有 C 語言基礎(chǔ)、想直接在 OpenHarmony 小設(shè)備上跑通這顆傳感器驅(qū)動(dòng)開發(fā)流程的朋友。1. 先搞清楚MAX30100在測什么血氧和心率背后的光學(xué)原理很多教程把 MAX30100 說得太簡單了好像就是讀寄存器、拿數(shù)據(jù)、算結(jié)果三步。但如果你不了解它到底在測量什么物理量后面遇到數(shù)據(jù)亂跳、數(shù)值不準(zhǔn)的時(shí)候根本不知道從哪查起。所以第一步先把這顆傳感器的物理邏輯講清楚。1.1 紅光與紅外光的吸收差異SpO2測出來的依據(jù)血氧飽和度的核心是氧合血紅蛋白HbO2和還原血紅蛋白Hb對(duì)光的吸收特性不同。簡單說氧合血紅蛋白對(duì) 660nm 左右的紅光吸收更弱、對(duì) 940nm 左右的紅外光吸收更強(qiáng)還原血紅蛋白正好相反。MAX30100 內(nèi)部集成了紅光 LED 和紅外 LED再加上一個(gè)光電接收管。它讓紅光和紅外交替發(fā)光接收管測量經(jīng)過組織反射回來的光強(qiáng)根據(jù)兩種波長的吸收比就能推斷血液中氧合血紅蛋白的比例。理解這一點(diǎn)對(duì)驅(qū)動(dòng)開發(fā)很重要。因?yàn)槟阕x出來的不是血氧百分比而是兩個(gè)原始 ADC 值紅光通道和紅外通道的反射光強(qiáng)。血氧值必須用這兩個(gè)通道的相對(duì)關(guān)系算出來只看任意一個(gè)通道是得不出結(jié)果的。這也就解釋了為什么調(diào)試時(shí)一定要同時(shí)觀察兩條波形而不是只看一條。1.2 動(dòng)脈搏動(dòng)產(chǎn)生AC分量心率信號(hào)藏在0.5%~5%的變化里心率測量的原理是光電容積脈搏波PPG。心臟每次搏動(dòng)動(dòng)脈血管里的血容量會(huì)周期性變化這會(huì)導(dǎo)致皮膚組織對(duì)光的吸收量發(fā)生微小波動(dòng)。于是反射光強(qiáng)出現(xiàn)一個(gè)很小的周期性起伏這個(gè)起伏的頻率就是心率。這里有個(gè)關(guān)鍵概念光強(qiáng)信號(hào)里包含著巨大的直流分量DC它來自皮膚、骨骼、靜脈血等固定組織對(duì)光的穩(wěn)定吸收數(shù)值可能是幾千而心跳引起的交流分量AC可能只占 DC 的 0.5% 到 5%也就是幾十到幾百 LSB。MAX30100 的 ADC 可以直接輸出這個(gè)疊加了 DC 和 AC 的原始值驅(qū)動(dòng)層要做的就是把這個(gè)信號(hào)完整、連續(xù)、不失真地采集回來后續(xù)用濾波把 DC 去掉把 AC 分離出來。注意AC 分量極其微弱這就是為什么硬件上很小的電源紋波都可能直接淹沒信號(hào)。我一開始把傳感器模塊直接掛在開發(fā)板的 3.3V 上波形里全是周期性的毛刺后來專門用一個(gè) LDO 單獨(dú)給傳感器供電并且在靠近芯片的 VIN 引腳并了 10μF 加 0.1μF 的去耦電容信號(hào)才干凈下來。1.3 硬件連接與上電時(shí)序別讓電源紋波毀掉PPGMAX30100 是標(biāo)準(zhǔn) I2C 接口7 位地址是 0x57對(duì)應(yīng) 8 位寫地址 0xAE、讀地址 0xAF。除了 VIN、GND、SDA、SCL還有一個(gè) INT 引腳可以在 FIFO 達(dá)到某個(gè)閾值時(shí)產(chǎn)生中斷。我的板子上最初沒有接 INT完全靠輪詢方式讀取100Hz 采樣率下完全夠用。上電時(shí)序值得單獨(dú)說芯片上電后寄存器不是立刻穩(wěn)定的我從實(shí)際調(diào)試中發(fā)現(xiàn)剛上電立刻讀設(shè)備 ID 很容易失敗。建議流程是先讓電源穩(wěn)定幾毫秒再讀一次 PART_ID 驗(yàn)證 I2C 通訊確認(rèn)通訊正常以后再寫配置寄存器。如果一上來就連續(xù)寫一串配置偶爾會(huì)出現(xiàn)前面的配置丟失、后面的配置生效的詭異情況排查起來很費(fèi)時(shí)間。還有一個(gè)容易忽略的點(diǎn)LED 脈沖瞬間電流不小如果 VIN 走線和其它大電流外設(shè)共用一條電源線PPG 信號(hào)里會(huì)混入明顯的工頻或開關(guān)噪聲。條件允許的話把它當(dāng)作一個(gè)模擬前端器件來對(duì)待電源單獨(dú)走、地線單點(diǎn)接。2. 在OpenHarmony這邊驅(qū)動(dòng)代碼掛在哪里HDF的落地路徑傳感器物理上通了接下來是 OpenHarmony 側(cè)的問題。OpenHarmony 的驅(qū)動(dòng)框架叫 HDFHarmonyOS Driver Framework它有一套統(tǒng)一的驅(qū)動(dòng)模型從設(shè)備描述、驅(qū)動(dòng)加載到對(duì)外服務(wù)接口都有規(guī)范。剛開始接觸會(huì)有點(diǎn)懵其實(shí)拆開看就三層驅(qū)動(dòng)入口、設(shè)備描述HCS 配置、平臺(tái)接口I2C 等。2.1 HDF驅(qū)動(dòng)生命周期Bind、Init、Release到底干了什么HDF 驅(qū)動(dòng)的生命周期很清晰。系統(tǒng)啟動(dòng)時(shí)HDF 框架會(huì)解析配置文件找到匹配的設(shè)備節(jié)點(diǎn)和驅(qū)動(dòng)模塊然后依次調(diào)用三個(gè)函數(shù)Bind、Init、Release。Bind 負(fù)責(zé)設(shè)備路徑與驅(qū)動(dòng)服務(wù)的綁定通常用來給上層應(yīng)用提供一個(gè)可訪問的接口路徑Init 做真正的硬件初始化比如解析配置參數(shù)、打開 I2C 控制器、配置傳感器寄存器、創(chuàng)建采集線程Release 則是卸載時(shí)釋放資源、關(guān)閉設(shè)備。對(duì)應(yīng)到 MAX30100 驅(qū)動(dòng)上大部分核心工作都放在 Init 里Bind 如果暫時(shí)不做服務(wù)接口直接返回成功即可。驅(qū)動(dòng)入口結(jié)構(gòu)體需要指定模塊名HDF 框架憑 moduleName 找到這個(gè)驅(qū)動(dòng)模塊。注冊(cè)入口的宏在舊版本里是HDF_INIT新版本有的分支改成了HDF_DRIVER_INIT。如果你編譯報(bào)這個(gè)宏找不到去系統(tǒng)已有的驅(qū)動(dòng)源碼里搜一下照著寫就行。這類版本差異不用背遇到時(shí)現(xiàn)場對(duì)照最省事。2.2 HCS設(shè)備描述讓驅(qū)動(dòng)參數(shù)脫離C代碼HDF 用的是 HCS 配置文件它是類樹狀結(jié)構(gòu)的描述語言。設(shè)備在哪里、匹配哪個(gè)驅(qū)動(dòng)屬性、使用哪條 I2C 總線、從機(jī)地址是多少這些參數(shù)都寫在 HCS 里驅(qū)動(dòng)代碼里不再硬編碼。我項(xiàng)目里的 device_info.hcs 大致長這樣root { sample_host { device_max30100 { deviceName max30100_device; deviceMatchAttr max30100_config; #include max30100_config.hcs } } }對(duì)應(yīng)的 max30100_config.hcs 里放驅(qū)動(dòng)自定義屬性root { max30100_config { i2cBusNum 1; i2cAddr 0x57; sampleRate 100; redLedCurrent 0x21; irLedCurrent 0x24; } }這樣做的價(jià)值在于換了一塊板子、I2C 控制器變了只要改 HCS 文件C 代碼一行不動(dòng)。我后面前后換了兩個(gè)硬件平臺(tái)這個(gè)設(shè)計(jì)幫我省了很多事。注意deviceMatchAttr的值要和驅(qū)動(dòng)代碼里讀取配置時(shí)使用的匹配屬性一致否則驅(qū)動(dòng)能找到模塊但拿不到自定義參數(shù)。2.3 驅(qū)動(dòng)源碼目錄與BUILD.gn的組織方式驅(qū)動(dòng)源碼我放在板級(jí)目錄下的 drivers 子目錄里結(jié)構(gòu)很簡單device/board/xxx/drivers/max30100/ ├── max30100.c ├── max30100.h ├── max30100_config.hcs └── BUILD.gnBUILD.gn 用 HDF 提供的驅(qū)動(dòng)模板編譯簡化版長這樣import(//build/ohos.gni) hdf_driver(max30100_driver) { sources [ max30100.c, ] include_dirs [ //drivers/framework/include, //drivers/framework/include/platform, ] deps [ //drivers/framework/stream/platform/i2c:libi2c, //drivers/framework/core/adapter:libhdf_adapter, ] }不同 OpenHarmony 版本的模板名稱、依賴路徑差別不小我每次換版本都會(huì)對(duì)照系統(tǒng)里現(xiàn)成的驅(qū)動(dòng) BUILD.gn 改。不用死記關(guān)鍵是理解依賴關(guān)系驅(qū)動(dòng)代碼需要包含 HDF 核心頭文件和 I2C 平臺(tái)接口頭文件。3. 驅(qū)動(dòng)開發(fā)實(shí)操從寄存器配置到FIFO數(shù)據(jù)流前面準(zhǔn)備做足了這一章進(jìn)入代碼。我會(huì)按實(shí)際開發(fā)順序來講先寫 I2C 讀寫封裝再做初始化配置最后設(shè)計(jì)連續(xù)采集的數(shù)據(jù)通路。寫驅(qū)動(dòng)和寫業(yè)務(wù)代碼不一樣每一步都要能獨(dú)立驗(yàn)證別指望一口氣寫完再統(tǒng)一調(diào)試那樣出了問題根本定位不到。3.1 驅(qū)動(dòng)入口與I2C讀寫封裝先定義設(shè)備結(jié)構(gòu)體和兩個(gè) I2C 傳輸函數(shù)。HDF 的 I2C 接口在舊版本里是I2cOpen、I2cTransfer、I2cClose如果你的系統(tǒng)版本用的是新的I2cCnRead、I2cCnWrite替換對(duì)應(yīng)調(diào)用即可封裝思路是一樣的#include hdf_device_desc.h #include hdf_log.h #include device_resource_if.h #include i2c_if.h #define MAX30100_I2C_ADDR 0x57 #define MAX30100_REG_PART_ID 0x0D typedef struct { DevHandle i2cHandle; uint8_t i2cAddr; uint32_t busNum; uint32_t sampleRate; uint8_t redLedCurrent; uint8_t irLedCurrent; } Max30100Device; static Max30100Device g_max30100; static int32_t Max30100WriteReg(uint8_t reg, uint8_t val) { struct I2cMsg msg; uint8_t buf[2]; buf[0] reg; buf[1] val; msg.addr g_max30100.i2cAddr; msg.flags 0; msg.len 2; msg.buf buf; return I2cTransfer(g_max30100.i2cHandle, msg, 1); } static int32_t Max30100ReadRegs(uint8_t reg, uint8_t *buf, uint16_t len) { struct I2cMsg msgs[2]; uint8_t regAddr reg; msgs[0].addr g_max30100.i2cAddr; msgs[0].flags 0; msgs[0].len 1; msgs[0].buf regAddr; msgs[1].addr g_max30100.i2cAddr; msgs[1].flags I2C_FLAG_READ; msgs[1].len len; msgs[1].buf buf; return I2cTransfer(g_max30100.i2cHandle, msgs, 2); }寫寄存器用先發(fā)寄存器地址、再發(fā)數(shù)據(jù)的一次寫事務(wù)讀用先寫寄存器地址、再發(fā)讀事務(wù)的兩段式傳輸。我在調(diào)試中發(fā)現(xiàn)千萬不要把多個(gè)寄存器寫操作一次性拼到一個(gè)事務(wù)里MAX30100 對(duì) I2C 時(shí)序挺敏感分開寫、讀回驗(yàn)證穩(wěn)定性明顯更好。然后是驅(qū)動(dòng)入口static int32_t Max30100Bind(struct HdfDeviceObject *device) { (void)device; return HDF_SUCCESS; } static int32_t Max30100Init(struct HdfDeviceObject *device) { // 解析HCS屬性、打開I2C、驗(yàn)證芯片、初始化寄存器、啟動(dòng)采集線程 } static void Max30100Release(struct HdfDeviceObject *device) { if (g_max30100.i2cHandle ! NULL) { I2cClose(g_max30100.i2cHandle); g_max30100.i2cHandle NULL; } } struct HdfDriverEntry g_max30100DriverEntry { .moduleVersion 1, .moduleName max30100_driver, .Bind Max30100Bind, .Init Max30100Init, .Release Max30100Release, }; HDF_INIT(g_max30100DriverEntry);Init 里第一件事是解析 HCS 屬性。不同版本接口名有差異我用的版本是GetUint32ElemValuestruct DeviceResourceIface *iface DeviceResourceIfaceGet(); if (iface NULL) { return HDF_FAILURE; } if (iface-GetUint32ElemValue(device-property, i2cBusNum, g_max30100.busNum, 0) ! HDF_SUCCESS) { HDF_LOGE(read i2cBusNum failed); return HDF_FAILURE; } if (iface-GetUint32ElemValue(device-property, i2cAddr, g_max30100.i2cAddr, 0) ! HDF_SUCCESS) { HDF_LOGE(read i2cAddr failed); return HDF_FAILURE; }拿到總線號(hào)和地址之后打開 I2C 控制器g_max30100.i2cHandle I2cOpen((uint16_t)g_max30100.busNum); if (g_max30100.i2cHandle NULL) { HDF_LOGE(I2cOpen failed); return HDF_FAILURE; }先讀一次 PART_ID確認(rèn)芯片真的在總線上再做后續(xù)初始化。這個(gè)先驗(yàn)證再配置的習(xí)慣幫我避免了很多無效調(diào)試。3.2 初始化寄存器設(shè)置分步驗(yàn)證比一次寫齊更重要初始化寄存器最忌諱網(wǎng)上抄一份配置一寫了之。每個(gè)項(xiàng)目的電源環(huán)境、LED 電流、采樣率需求都不一樣我一開始把幾組配置一起寫結(jié)果信號(hào)不是飽和就是太弱最后只能逐個(gè)寄存器回退著查。關(guān)鍵的寄存器有這幾個(gè)寄存器地址我的配置策略FIFO_CONFIG0x06開啟滾動(dòng)模式允許新數(shù)據(jù)覆蓋舊數(shù)據(jù)防止溢出停采樣本平均設(shè)到 4可以平滑信號(hào)毛刺MODE_CONFIG0x07配成 SpO2 模式讓紅光和紅外 LED 交替工作SPO2_CONFIG0x08采樣率配到 100Hz這個(gè)頻率對(duì)心率信號(hào)充足ADC 量程和脈寬按手上手冊(cè)的位定義來選LED1_PA0x09紅光 LED 電流初始給一個(gè)中間值后面按波形調(diào)LED2_PA0x0A紅外 LED 電流同樣先給中間值配置代碼大致這樣static int32_t Max30100Config(void) { int32_t ret; ret Max30100WriteReg(0x06, 0x50); // 滾動(dòng)使能 樣本平均4 if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x07, 0x03); // SpO2模式 if (ret ! HDF_SUCCESS) return ret; // 采樣率等參數(shù)按數(shù)據(jù)手冊(cè)的位定義設(shè)置這里不寫死 ret Max30100WriteReg(0x08, sampleRateConfig); if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x09, g_max30100.redLedCurrent); if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x0A, g_max30100.irLedCurrent); if (ret ! HDF_SUCCESS) return ret; return HDF_SUCCESS; }配置完成后我習(xí)慣立刻回讀一遍確認(rèn)寄存器真的寫進(jìn)去了。有些時(shí)候 I2C 傳輸返回成功但設(shè)備端因?yàn)闀r(shí)序問題沒有真正接收回讀是最直接的驗(yàn)證手段。LED 電流的選擇值得單獨(dú)說電流太小AC 分量被噪聲淹沒電流太大接收管飽和ADC 讀數(shù)頂?shù)綕M量程血氧比值完全失真。我調(diào)這個(gè)花了很長時(shí)間最終原則是讓兩個(gè)通道的靜態(tài) DC 值落在 ADC 量程的中間偏下位置比如 16 位 ADC 的話讓 DC 落在 2000 到 3000 左右然后在這個(gè)基礎(chǔ)上微調(diào) AC 幅度。3.3 FIFO讀取與環(huán)形緩沖采集線程怎么設(shè)計(jì)MAX30100 內(nèi)部有個(gè) FIFO傳感器采集的數(shù)據(jù)先存在 FIFO 里驅(qū)動(dòng)再去讀。在 SpO2 模式下每個(gè)樣本包含紅光和紅外兩個(gè) 16 位 ADC 值即 4 字節(jié)。讀取時(shí)先看寫指針和讀指針?biāo)愠霎?dāng)前可讀幾個(gè)樣本再連續(xù)讀取。#define MAX30100_REG_FIFO_WR_PTR 0x02 #define MAX30100_REG_FIFO_RD_PTR 0x04 #define MAX30100_REG_FIFO_DATA 0x05 static int32_t Max30100ReadPpg(uint16_t *redArray, uint16_t *irArray, uint8_t *count) { uint8_t wr; uint8_t rd; uint8_t raw[32]; uint8_t available; int32_t ret; int32_t i; ret Max30100ReadRegs(MAX30100_REG_FIFO_WR_PTR, wr, 1); if (ret ! HDF_SUCCESS) return ret; ret Max30100ReadRegs(MAX30100_REG_FIFO_RD_PTR, rd, 1); if (ret ! HDF_SUCCESS) return ret; available (uint8_t)(((wr 32 - rd) 0x1F) / 4); if (available 0) { *count 0; return HDF_SUCCESS; } if (available 8) { available 8; // 單次最多拿8個(gè)樣本避免長時(shí)間占用I2C } ret Max30100ReadRegs(MAX30100_REG_FIFO_DATA, raw, (uint16_t)(available * 4)); if (ret ! HDF_SUCCESS) return ret; for (i 0; i available; i) { uint8_t *p raw[i * 4]; redArray[i] (uint16_t)((p[0] 8) | p[1]); irArray[i] (uint16_t)((p[2] 8) | p[3]); } *count available; return HDF_SUCCESS; }這里有個(gè)小細(xì)節(jié)計(jì)算可用樣本數(shù)時(shí)((wr 32 - rd) 0x1F)是把寫指針和讀指針放在同一個(gè)環(huán)形空間里求差避免指針回繞時(shí)相減出負(fù)數(shù)或出大數(shù)。MAX30100 的 FIFO 深度是 32 個(gè)樣本單位SpO2 模式下每個(gè)樣本點(diǎn)占 4 個(gè)單元所以再除以 4 得到可讀樣本數(shù)。采集線程我用了最直白的輪詢方式每隔 10ms 調(diào)一次Max30100ReadPpg把數(shù)據(jù)寫進(jìn)一個(gè)環(huán)形緩沖區(qū)。環(huán)形緩沖大小我設(shè)置成 128 組樣本大約是 1 秒多鐘的數(shù)據(jù)量足夠應(yīng)用層做算法運(yùn)算。生產(chǎn)者和消費(fèi)者的邏輯很簡單typedef struct { uint16_t red[128]; uint16_t ir[128]; uint16_t head; uint16_t tail; } PpgRingBuffer;采集線程向 head 寫應(yīng)用層從 tail 讀。如果應(yīng)用層來不及消費(fèi)就覆蓋最舊的數(shù)據(jù)保證拿到的永遠(yuǎn)是最新的一段波。數(shù)據(jù)規(guī)模不大用數(shù)組實(shí)現(xiàn)比鏈表簡單可靠得多也不會(huì)因?yàn)閮?nèi)存分配產(chǎn)生碎片問題。3.4 把數(shù)據(jù)送給應(yīng)用層輪詢、中斷與服務(wù)接口的選擇數(shù)據(jù)最后怎么交給應(yīng)用層我前前后后試過三種方式。最簡單的是輪詢驅(qū)動(dòng)里的采集線程負(fù)責(zé)讀 FIFO應(yīng)用層按需調(diào)用查詢接口從環(huán)形緩沖里取數(shù)據(jù)。這種方式代碼少100Hz 采樣率下足夠穩(wěn)我就先用它跑通了整條鏈路。第二種是用 INT 引腳。MAX30100 的 INT 在 FIFO 數(shù)據(jù)達(dá)到一定量時(shí)拉低接到開發(fā)板的 GPIO 中斷上驅(qū)動(dòng)收到中斷以后再去讀。好處是功耗低、不需要固定周期喚醒適合做低功耗穿戴設(shè)備的場景。代價(jià)是要多申請(qǐng)一個(gè) GPIO 中斷而且中斷服務(wù)函數(shù)里不能做太重的事通常是置一個(gè)標(biāo)志喚醒采集線程去干活。第三種是做成 HDF 服務(wù)由上層應(yīng)用綁定服務(wù)路徑后調(diào)用。這是最規(guī)范的做法但牽扯到 HDF 服務(wù)接口定義、權(quán)限配置、應(yīng)用側(cè)跨進(jìn)程調(diào)用工程量會(huì)大不少。我的經(jīng)驗(yàn)是先把數(shù)據(jù)鏈路調(diào)通再考慮要不要上服務(wù)框架。4. 原始PPG到手之后算法與驅(qū)動(dòng)的分工邊界驅(qū)動(dòng)把連續(xù)的 PPG 原始數(shù)據(jù)送到環(huán)形緩沖之后接下來才是真正意義的數(shù)據(jù)處理。這一章重點(diǎn)講算法和驅(qū)動(dòng)的分工因?yàn)槲以谶@個(gè)環(huán)節(jié)走過彎路一開始想在驅(qū)動(dòng)里把血氧和心率全算完結(jié)果耦合太重、調(diào)參也麻煩后面全部推倒重來。4.1 原始PPG波形長什么樣先學(xué)會(huì)畫圖在寫任何算法之前先把數(shù)據(jù)畫出來看。最笨也最有效的方法是串口按 CSV 格式打印采樣序號(hào)、紅光 ADC、紅外 ADCPC 端用 Python 腳本讀取并畫波形。正常的 PPG 波形長這樣每個(gè)通道都是一個(gè)接近直流的基線基線值在幾千左右上面疊加小幅度的、周期性的起伏每個(gè)起伏對(duì)應(yīng)一次心跳。紅光通道和紅外通道的起伏應(yīng)該基本同步但幅度可能不同。如果你看到兩個(gè)通道的起伏不同步、形狀差異很大說明某個(gè)通道飽和了或者 LED 電流配比不對(duì)。這一步?jīng)]有捷徑。我見過很聰明的同學(xué)直接跳過畫圖開始寫血氧算法最后調(diào)了兩天也不知道是算法問題還是信號(hào)問題。信號(hào)都沒看明白算法就是盲人摸象。4.2 心率計(jì)算峰值檢測與帶通濾波心率算法我放在應(yīng)用層做驅(qū)動(dòng)層只保證數(shù)據(jù)連續(xù)、帶時(shí)間戳。算法思路分三步濾波、峰檢測、間隔換算。先對(duì)紅外通道做帶通濾波。心率范圍一般在 30 到 200 次每分鐘換算成頻率大約是 0.5Hz 到 3.3Hz用一個(gè)兩步的 IIR 濾波器把直流和低頻漂移濾掉同時(shí)把 5Hz 以上的高頻噪聲去掉。濾波之后基線被拉平剩下的是清晰的脈搏波。然后是峰值檢測。我最常用的是一階差分加滑動(dòng)窗口// 偽代碼尋找局部最大值并限制最小波峰間隔 if (current prev current next current threshold (now - lastBeatTime) minInterval) { heartRate 60000 / (now - lastBeatTime); lastBeatTime now; }這個(gè)思路簡單但在信號(hào)質(zhì)量不錯(cuò)的場景下很有效。minInterval用來防止把同一波形的兩個(gè)小毛刺誤判成兩次心跳我一般設(shè) 300ms對(duì)應(yīng)最大 200bpm。注意心率輸出不要每拍都更新正確做法是取過去 4 到 8 拍的間隔做滑動(dòng)平均否則讀數(shù)會(huì)跳動(dòng)得很厲害。4.3 SpO2計(jì)算R值公式和查表思路血氧計(jì)算比心率稍微繞一點(diǎn)但原理不復(fù)雜。對(duì)紅光和紅外兩個(gè)通道分別提取交流分量 AC 和直流分量 DC計(jì)算比值 RR (AC_RED / DC_RED) / (AC_IR / DC_IR)然后通過經(jīng)驗(yàn)公式或查表把 R 映射到 SpO2。標(biāo)準(zhǔn)模型下R 值在 0.4 附近對(duì)應(yīng)接近 100% 的血氧R 為 1.0 左右大約對(duì)應(yīng) 94%R 為 2.0 時(shí)大約降到 75% 附近。這個(gè)區(qū)間可以做個(gè)粗略對(duì)照R值參考SpO20.4約100%0.6約97%1.0約94%1.5約85%2.0約75%注意這張表只能作為演示示意不能直接用于醫(yī)療用途。不同膚色、不同 LED 波長偏差、不同佩戴松緊都會(huì)改變映射曲線做到產(chǎn)品級(jí)必須用自己的硬件和采集環(huán)境做多用戶標(biāo)定。實(shí)際計(jì)算中AC 不能直接取原始數(shù)據(jù)的振幅一般是用帶通濾波后的信號(hào)按一個(gè)滑動(dòng)窗口算均方根或平均峰峰值。窗口太長響應(yīng)慢窗口太短數(shù)值抖動(dòng)大我最終選了 4 秒窗口每 1 秒輸出一次平衡了更新速度和穩(wěn)定性。4.4 為什么我不建議在驅(qū)動(dòng)層做算法我知道有人會(huì)覺得既然驅(qū)動(dòng)里已經(jīng)有數(shù)據(jù)了順手就把心率算了多方便。我確實(shí)這么干過然后果斷放棄了。原因是算法開發(fā)是高頻迭代的濾波器的截止頻率、窗口長度、閾值參數(shù)幾乎每天都要調(diào)。如果這些邏輯塞在驅(qū)動(dòng)里每次調(diào)參都要重新編譯整個(gè)驅(qū)動(dòng)鏡像、燒錄、重啟設(shè)備效率極低。驅(qū)動(dòng)層的正確姿勢是穩(wěn)定采數(shù)、打時(shí)間戳、管理緩沖、提供干凈的接口。算法放應(yīng)用層可以單獨(dú)編譯、單獨(dú)調(diào)試甚至可以在 PC 上用同樣的算法處理串口傳出來的數(shù)據(jù)。這個(gè)架構(gòu)上的分界比優(yōu)化任何一個(gè)具體算法細(xì)節(jié)都更重要。5. 想拿到穩(wěn)定的數(shù)據(jù)我踩過并解決過的坑最后這一章我把整個(gè)開發(fā)周期里遇到的最典型的幾個(gè)問題列出來。這些問題在網(wǎng)上都能搜到零散答案但很少有文章把排查思路講清楚我盡量還原當(dāng)時(shí)的判斷過程。5.1 設(shè)備ID讀出來全FF先查時(shí)序再查I2C第一次上電讀 PART_ID 返回 0xFF也就是總線上根本沒應(yīng)答。我當(dāng)時(shí)第一反應(yīng)是模塊壞了換了一片還是同樣現(xiàn)象才開始系統(tǒng)排查。排查順序是這樣的先查 I2C 地址確認(rèn)發(fā)的是 0xAE 寫地址再查 HCS 配置里的總線號(hào)會(huì)不會(huì)配到了一條沒有設(shè)備的空總線上然后用邏輯分析儀看 SCL/SDA 波形發(fā)現(xiàn) SDA 根本沒有被拉低說明從設(shè)備沒有拉總線最后檢查模塊板上的上拉電阻發(fā)現(xiàn)模塊自帶上拉但因?yàn)榻恿艘桓容^長的杜邦線加上開發(fā)板內(nèi)部也有上拉組合起來電平反而異常。換用短線和獨(dú)立 4.7k 上拉后ID 正常讀回。還有一個(gè)容易忽略的原因剛上電立刻讀 ID芯片還沒準(zhǔn)備好。我后來在 Init 里先OsalMSleep(10)再讀成功率明顯提升。讀到 ID 值之后要和你手上的芯片批次資料核對(duì)不同批次的 ID 可能會(huì)不一樣不同型號(hào)的兼容芯片也可能返回不同數(shù)值。5.2 FIFO溢出或長期為空采樣率與讀出節(jié)奏的平衡調(diào)試中我遇到過兩個(gè)極端一個(gè)現(xiàn)象是 FIFO 溢出標(biāo)志一直置起另一個(gè)現(xiàn)象是讀出來的樣本永遠(yuǎn)是 0。這兩個(gè)問題其實(shí)是同一類問題采樣率和讀取節(jié)奏失配。溢出很好理解采樣率 800Hz但上層讀取線程每隔 50ms 才讀一次FIFO 早就滿了。如果沒有開滾動(dòng)模式芯片會(huì)停止采樣數(shù)據(jù)就斷了即使在滾動(dòng)模式數(shù)據(jù)也會(huì)大量覆蓋丟失。解決方法是把采樣率降到 100Hz讀取線程保證至少每 20ms 讀一次。我之前用消息隊(duì)列間接讀延遲漲到幾十毫秒后數(shù)據(jù)就開始出現(xiàn)空洞后來直接用一個(gè)高優(yōu)先級(jí)的采集線程輪詢問題消失。長期為空的原因則更隱蔽。起初我懷疑是寄存器配置問題反復(fù)檢查都沒發(fā)現(xiàn)錯(cuò)誤。后來發(fā)現(xiàn)采集線程創(chuàng)建成功但一直阻塞在等待一個(gè)應(yīng)用層才會(huì)發(fā)送的信號(hào)量上導(dǎo)致 FIFO 里的數(shù)據(jù)永遠(yuǎn)沒人讀。這是一個(gè)純粹的調(diào)度設(shè)計(jì)錯(cuò)誤和傳感器本身無關(guān)。這也提醒我驅(qū)動(dòng)里的采集線程和上層消費(fèi)邏輯盡量解耦采集線程只依賴自己的定時(shí)器不要依賴外部信號(hào)。5.3 SpO2數(shù)值亂跳信號(hào)質(zhì)量評(píng)估比算法更關(guān)鍵血氧值亂跳是幾乎所有健康監(jiān)測設(shè)備都會(huì)遇到的問題。剛開始我的 R 值算出來一會(huì)兒 0.3、一會(huì)兒 2.5血氧從 100% 直接掉到 75%但人的狀態(tài)完全正常。仔細(xì)觀察原始波形后才明白是手指按壓力度變化和微小晃動(dòng)導(dǎo)致的運(yùn)動(dòng)偽跡。解決這個(gè)問題算法上最高效的不是優(yōu)化公式而是加信號(hào)質(zhì)量評(píng)估。我加了三個(gè)判斷條件一是 R 值必須落在 0.3 到 2.5 的合理區(qū)間內(nèi)超出就不輸出二是紅光和紅外通道的波峰間隔要基本一致差異超過 20% 就視為信號(hào)不可靠三是 AC 分量占比不能太低AC 太小說明信號(hào)太弱通常是因?yàn)槭种笡]有貼緊或 LED 電流不足。經(jīng)過這三道閘以后穩(wěn)定坐姿下血氧讀數(shù)能連續(xù)幾分鐘保持誤差在一個(gè)很小的范圍內(nèi)。運(yùn)動(dòng)狀態(tài)下想做到完全穩(wěn)定那就要上加速度傳感器輔助做運(yùn)動(dòng)偽跡消除已經(jīng)超出這顆芯片本身的能力范圍了。5.4 調(diào)試方法串口日志加Python繪圖具體調(diào)試手段上我推薦一套投入小、見效快的組合驅(qū)動(dòng)把數(shù)據(jù)以文本幀格式從串口發(fā)出來PC 端用 Python 的 pyserial 接收再用 matplotlib 實(shí)時(shí)繪圖。串口輸出格式我做成這樣frame_no, red, ir 123, 2400, 2800 124, 2420, 2820驅(qū)動(dòng)里注意不要在 I2C 讀數(shù)據(jù)的過程中直接打印日志串口打印很慢會(huì)把采樣節(jié)奏拖垮。正確做法是先把本次讀到的數(shù)據(jù)存到環(huán)形緩沖到主循環(huán)外面再按批量輸出。Python 側(cè)我每次只畫最近 10 秒的波形這樣既能看細(xì)節(jié)又不會(huì)糊成一團(tuán)。畫圖的時(shí)候重點(diǎn)關(guān)注三個(gè)指標(biāo)DC 基線是否穩(wěn)定、AC 幅度是否在合理范圍、兩個(gè)通道的波峰是否同步。這三個(gè)指標(biāo)正常再去看心率血氧計(jì)算結(jié)果否則先回頭查硬件和配置。還有一個(gè)小技巧用串口輸出的樣本序號(hào)檢查丟幀。如果序號(hào)不連續(xù)說明驅(qū)動(dòng)讀取側(cè)有丟數(shù)據(jù)這時(shí)候不管算法多好最終結(jié)果都會(huì)受影響。先把序號(hào)補(bǔ)成全連續(xù)再接算法這條調(diào)試順序我后來一直遵守。如果問我在這個(gè)項(xiàng)目里最大的體會(huì)其實(shí)不是學(xué)會(huì)了怎么配 HDF、怎么寫 I2C 封裝而是意識(shí)到驅(qū)動(dòng)開發(fā)真正的難點(diǎn)是數(shù)據(jù)質(zhì)量而不是代碼框架。MAX30100 是一顆很便宜的傳感器但想拿到穩(wěn)定的 PPG 數(shù)據(jù)從電源、布局、LED 電流到讀取節(jié)奏每一個(gè)環(huán)節(jié)都在影響最終結(jié)果。我建議你動(dòng)手時(shí)先別急著寫算法把波形打印出來盯一兩個(gè)晚上信號(hào)看著順眼了心率血氧自然就出來了。后面有空我再寫一篇應(yīng)用層怎么通過 HDF 服務(wù)讀取驅(qū)動(dòng)數(shù)據(jù)以及把實(shí)時(shí)波形傳到手機(jī)端調(diào)試的做法。 OpenHarmony下MAX30100血氧心率傳感器驅(qū)動(dòng)開發(fā)實(shí)戰(zhàn)
前陣子做健康監(jiān)測小項(xiàng)目我在網(wǎng)購模塊里挑了一片 MAX30100 血氧心跳傳感器以為接上 I2C、讀幾個(gè)寄存器就能把數(shù)據(jù)跑起來結(jié)果從 OpenHarmony 側(cè)驅(qū)動(dòng)注冊(cè)到 PPG 信號(hào)穩(wěn)定輸出前后折騰了不少時(shí)間。這個(gè)過程中踩過的坑、驗(yàn)證過的方法整理出來就是這篇文章。內(nèi)容覆蓋 MAX30100 的測量原理、HDF 驅(qū)動(dòng)框架下的驅(qū)動(dòng)實(shí)現(xiàn)、FIFO 數(shù)據(jù)讀取與緩沖設(shè)計(jì)、血氧和心率算法的邊界以及最終調(diào)試驗(yàn)證時(shí)遇到的一堆實(shí)際問題。適合有 C 語言基礎(chǔ)、想直接在 OpenHarmony 小設(shè)備上跑通這顆傳感器驅(qū)動(dòng)開發(fā)流程的朋友。1. 先搞清楚MAX30100在測什么血氧和心率背后的光學(xué)原理很多教程把 MAX30100 說得太簡單了好像就是讀寄存器、拿數(shù)據(jù)、算結(jié)果三步。但如果你不了解它到底在測量什么物理量后面遇到數(shù)據(jù)亂跳、數(shù)值不準(zhǔn)的時(shí)候根本不知道從哪查起。所以第一步先把這顆傳感器的物理邏輯講清楚。1.1 紅光與紅外光的吸收差異SpO2測出來的依據(jù)血氧飽和度的核心是氧合血紅蛋白HbO2和還原血紅蛋白Hb對(duì)光的吸收特性不同。簡單說氧合血紅蛋白對(duì) 660nm 左右的紅光吸收更弱、對(duì) 940nm 左右的紅外光吸收更強(qiáng)還原血紅蛋白正好相反。MAX30100 內(nèi)部集成了紅光 LED 和紅外 LED再加上一個(gè)光電接收管。它讓紅光和紅外交替發(fā)光接收管測量經(jīng)過組織反射回來的光強(qiáng)根據(jù)兩種波長的吸收比就能推斷血液中氧合血紅蛋白的比例。理解這一點(diǎn)對(duì)驅(qū)動(dòng)開發(fā)很重要。因?yàn)槟阕x出來的不是血氧百分比而是兩個(gè)原始 ADC 值紅光通道和紅外通道的反射光強(qiáng)。血氧值必須用這兩個(gè)通道的相對(duì)關(guān)系算出來只看任意一個(gè)通道是得不出結(jié)果的。這也就解釋了為什么調(diào)試時(shí)一定要同時(shí)觀察兩條波形而不是只看一條。1.2 動(dòng)脈搏動(dòng)產(chǎn)生AC分量心率信號(hào)藏在0.5%~5%的變化里心率測量的原理是光電容積脈搏波PPG。心臟每次搏動(dòng)動(dòng)脈血管里的血容量會(huì)周期性變化這會(huì)導(dǎo)致皮膚組織對(duì)光的吸收量發(fā)生微小波動(dòng)。于是反射光強(qiáng)出現(xiàn)一個(gè)很小的周期性起伏這個(gè)起伏的頻率就是心率。這里有個(gè)關(guān)鍵概念光強(qiáng)信號(hào)里包含著巨大的直流分量DC它來自皮膚、骨骼、靜脈血等固定組織對(duì)光的穩(wěn)定吸收數(shù)值可能是幾千而心跳引起的交流分量AC可能只占 DC 的 0.5% 到 5%也就是幾十到幾百 LSB。MAX30100 的 ADC 可以直接輸出這個(gè)疊加了 DC 和 AC 的原始值驅(qū)動(dòng)層要做的就是把這個(gè)信號(hào)完整、連續(xù)、不失真地采集回來后續(xù)用濾波把 DC 去掉把 AC 分離出來。注意AC 分量極其微弱這就是為什么硬件上很小的電源紋波都可能直接淹沒信號(hào)。我一開始把傳感器模塊直接掛在開發(fā)板的 3.3V 上波形里全是周期性的毛刺后來專門用一個(gè) LDO 單獨(dú)給傳感器供電并且在靠近芯片的 VIN 引腳并了 10μF 加 0.1μF 的去耦電容信號(hào)才干凈下來。1.3 硬件連接與上電時(shí)序別讓電源紋波毀掉PPGMAX30100 是標(biāo)準(zhǔn) I2C 接口7 位地址是 0x57對(duì)應(yīng) 8 位寫地址 0xAE、讀地址 0xAF。除了 VIN、GND、SDA、SCL還有一個(gè) INT 引腳可以在 FIFO 達(dá)到某個(gè)閾值時(shí)產(chǎn)生中斷。我的板子上最初沒有接 INT完全靠輪詢方式讀取100Hz 采樣率下完全夠用。上電時(shí)序值得單獨(dú)說芯片上電后寄存器不是立刻穩(wěn)定的我從實(shí)際調(diào)試中發(fā)現(xiàn)剛上電立刻讀設(shè)備 ID 很容易失敗。建議流程是先讓電源穩(wěn)定幾毫秒再讀一次 PART_ID 驗(yàn)證 I2C 通訊確認(rèn)通訊正常以后再寫配置寄存器。如果一上來就連續(xù)寫一串配置偶爾會(huì)出現(xiàn)前面的配置丟失、后面的配置生效的詭異情況排查起來很費(fèi)時(shí)間。還有一個(gè)容易忽略的點(diǎn)LED 脈沖瞬間電流不小如果 VIN 走線和其它大電流外設(shè)共用一條電源線PPG 信號(hào)里會(huì)混入明顯的工頻或開關(guān)噪聲。條件允許的話把它當(dāng)作一個(gè)模擬前端器件來對(duì)待電源單獨(dú)走、地線單點(diǎn)接。2. 在OpenHarmony這邊驅(qū)動(dòng)代碼掛在哪里HDF的落地路徑傳感器物理上通了接下來是 OpenHarmony 側(cè)的問題。OpenHarmony 的驅(qū)動(dòng)框架叫 HDFHarmonyOS Driver Framework它有一套統(tǒng)一的驅(qū)動(dòng)模型從設(shè)備描述、驅(qū)動(dòng)加載到對(duì)外服務(wù)接口都有規(guī)范。剛開始接觸會(huì)有點(diǎn)懵其實(shí)拆開看就三層驅(qū)動(dòng)入口、設(shè)備描述HCS 配置、平臺(tái)接口I2C 等。2.1 HDF驅(qū)動(dòng)生命周期Bind、Init、Release到底干了什么HDF 驅(qū)動(dòng)的生命周期很清晰。系統(tǒng)啟動(dòng)時(shí)HDF 框架會(huì)解析配置文件找到匹配的設(shè)備節(jié)點(diǎn)和驅(qū)動(dòng)模塊然后依次調(diào)用三個(gè)函數(shù)Bind、Init、Release。Bind 負(fù)責(zé)設(shè)備路徑與驅(qū)動(dòng)服務(wù)的綁定通常用來給上層應(yīng)用提供一個(gè)可訪問的接口路徑Init 做真正的硬件初始化比如解析配置參數(shù)、打開 I2C 控制器、配置傳感器寄存器、創(chuàng)建采集線程Release 則是卸載時(shí)釋放資源、關(guān)閉設(shè)備。對(duì)應(yīng)到 MAX30100 驅(qū)動(dòng)上大部分核心工作都放在 Init 里Bind 如果暫時(shí)不做服務(wù)接口直接返回成功即可。驅(qū)動(dòng)入口結(jié)構(gòu)體需要指定模塊名HDF 框架憑 moduleName 找到這個(gè)驅(qū)動(dòng)模塊。注冊(cè)入口的宏在舊版本里是HDF_INIT新版本有的分支改成了HDF_DRIVER_INIT。如果你編譯報(bào)這個(gè)宏找不到去系統(tǒng)已有的驅(qū)動(dòng)源碼里搜一下照著寫就行。這類版本差異不用背遇到時(shí)現(xiàn)場對(duì)照最省事。2.2 HCS設(shè)備描述讓驅(qū)動(dòng)參數(shù)脫離C代碼HDF 用的是 HCS 配置文件它是類樹狀結(jié)構(gòu)的描述語言。設(shè)備在哪里、匹配哪個(gè)驅(qū)動(dòng)屬性、使用哪條 I2C 總線、從機(jī)地址是多少這些參數(shù)都寫在 HCS 里驅(qū)動(dòng)代碼里不再硬編碼。我項(xiàng)目里的 device_info.hcs 大致長這樣root { sample_host { device_max30100 { deviceName max30100_device; deviceMatchAttr max30100_config; #include max30100_config.hcs } } }對(duì)應(yīng)的 max30100_config.hcs 里放驅(qū)動(dòng)自定義屬性root { max30100_config { i2cBusNum 1; i2cAddr 0x57; sampleRate 100; redLedCurrent 0x21; irLedCurrent 0x24; } }這樣做的價(jià)值在于換了一塊板子、I2C 控制器變了只要改 HCS 文件C 代碼一行不動(dòng)。我后面前后換了兩個(gè)硬件平臺(tái)這個(gè)設(shè)計(jì)幫我省了很多事。注意deviceMatchAttr的值要和驅(qū)動(dòng)代碼里讀取配置時(shí)使用的匹配屬性一致否則驅(qū)動(dòng)能找到模塊但拿不到自定義參數(shù)。2.3 驅(qū)動(dòng)源碼目錄與BUILD.gn的組織方式驅(qū)動(dòng)源碼我放在板級(jí)目錄下的 drivers 子目錄里結(jié)構(gòu)很簡單device/board/xxx/drivers/max30100/ ├── max30100.c ├── max30100.h ├── max30100_config.hcs └── BUILD.gnBUILD.gn 用 HDF 提供的驅(qū)動(dòng)模板編譯簡化版長這樣import(//build/ohos.gni) hdf_driver(max30100_driver) { sources [ max30100.c, ] include_dirs [ //drivers/framework/include, //drivers/framework/include/platform, ] deps [ //drivers/framework/stream/platform/i2c:libi2c, //drivers/framework/core/adapter:libhdf_adapter, ] }不同 OpenHarmony 版本的模板名稱、依賴路徑差別不小我每次換版本都會(huì)對(duì)照系統(tǒng)里現(xiàn)成的驅(qū)動(dòng) BUILD.gn 改。不用死記關(guān)鍵是理解依賴關(guān)系驅(qū)動(dòng)代碼需要包含 HDF 核心頭文件和 I2C 平臺(tái)接口頭文件。3. 驅(qū)動(dòng)開發(fā)實(shí)操從寄存器配置到FIFO數(shù)據(jù)流前面準(zhǔn)備做足了這一章進(jìn)入代碼。我會(huì)按實(shí)際開發(fā)順序來講先寫 I2C 讀寫封裝再做初始化配置最后設(shè)計(jì)連續(xù)采集的數(shù)據(jù)通路。寫驅(qū)動(dòng)和寫業(yè)務(wù)代碼不一樣每一步都要能獨(dú)立驗(yàn)證別指望一口氣寫完再統(tǒng)一調(diào)試那樣出了問題根本定位不到。3.1 驅(qū)動(dòng)入口與I2C讀寫封裝先定義設(shè)備結(jié)構(gòu)體和兩個(gè) I2C 傳輸函數(shù)。HDF 的 I2C 接口在舊版本里是I2cOpen、I2cTransfer、I2cClose如果你的系統(tǒng)版本用的是新的I2cCnRead、I2cCnWrite替換對(duì)應(yīng)調(diào)用即可封裝思路是一樣的#include hdf_device_desc.h #include hdf_log.h #include device_resource_if.h #include i2c_if.h #define MAX30100_I2C_ADDR 0x57 #define MAX30100_REG_PART_ID 0x0D typedef struct { DevHandle i2cHandle; uint8_t i2cAddr; uint32_t busNum; uint32_t sampleRate; uint8_t redLedCurrent; uint8_t irLedCurrent; } Max30100Device; static Max30100Device g_max30100; static int32_t Max30100WriteReg(uint8_t reg, uint8_t val) { struct I2cMsg msg; uint8_t buf[2]; buf[0] reg; buf[1] val; msg.addr g_max30100.i2cAddr; msg.flags 0; msg.len 2; msg.buf buf; return I2cTransfer(g_max30100.i2cHandle, msg, 1); } static int32_t Max30100ReadRegs(uint8_t reg, uint8_t *buf, uint16_t len) { struct I2cMsg msgs[2]; uint8_t regAddr reg; msgs[0].addr g_max30100.i2cAddr; msgs[0].flags 0; msgs[0].len 1; msgs[0].buf regAddr; msgs[1].addr g_max30100.i2cAddr; msgs[1].flags I2C_FLAG_READ; msgs[1].len len; msgs[1].buf buf; return I2cTransfer(g_max30100.i2cHandle, msgs, 2); }寫寄存器用先發(fā)寄存器地址、再發(fā)數(shù)據(jù)的一次寫事務(wù)讀用先寫寄存器地址、再發(fā)讀事務(wù)的兩段式傳輸。我在調(diào)試中發(fā)現(xiàn)千萬不要把多個(gè)寄存器寫操作一次性拼到一個(gè)事務(wù)里MAX30100 對(duì) I2C 時(shí)序挺敏感分開寫、讀回驗(yàn)證穩(wěn)定性明顯更好。然后是驅(qū)動(dòng)入口static int32_t Max30100Bind(struct HdfDeviceObject *device) { (void)device; return HDF_SUCCESS; } static int32_t Max30100Init(struct HdfDeviceObject *device) { // 解析HCS屬性、打開I2C、驗(yàn)證芯片、初始化寄存器、啟動(dòng)采集線程 } static void Max30100Release(struct HdfDeviceObject *device) { if (g_max30100.i2cHandle ! NULL) { I2cClose(g_max30100.i2cHandle); g_max30100.i2cHandle NULL; } } struct HdfDriverEntry g_max30100DriverEntry { .moduleVersion 1, .moduleName max30100_driver, .Bind Max30100Bind, .Init Max30100Init, .Release Max30100Release, }; HDF_INIT(g_max30100DriverEntry);Init 里第一件事是解析 HCS 屬性。不同版本接口名有差異我用的版本是GetUint32ElemValuestruct DeviceResourceIface *iface DeviceResourceIfaceGet(); if (iface NULL) { return HDF_FAILURE; } if (iface-GetUint32ElemValue(device-property, i2cBusNum, g_max30100.busNum, 0) ! HDF_SUCCESS) { HDF_LOGE(read i2cBusNum failed); return HDF_FAILURE; } if (iface-GetUint32ElemValue(device-property, i2cAddr, g_max30100.i2cAddr, 0) ! HDF_SUCCESS) { HDF_LOGE(read i2cAddr failed); return HDF_FAILURE; }拿到總線號(hào)和地址之后打開 I2C 控制器g_max30100.i2cHandle I2cOpen((uint16_t)g_max30100.busNum); if (g_max30100.i2cHandle NULL) { HDF_LOGE(I2cOpen failed); return HDF_FAILURE; }先讀一次 PART_ID確認(rèn)芯片真的在總線上再做后續(xù)初始化。這個(gè)先驗(yàn)證再配置的習(xí)慣幫我避免了很多無效調(diào)試。3.2 初始化寄存器設(shè)置分步驗(yàn)證比一次寫齊更重要初始化寄存器最忌諱網(wǎng)上抄一份配置一寫了之。每個(gè)項(xiàng)目的電源環(huán)境、LED 電流、采樣率需求都不一樣我一開始把幾組配置一起寫結(jié)果信號(hào)不是飽和就是太弱最后只能逐個(gè)寄存器回退著查。關(guān)鍵的寄存器有這幾個(gè)寄存器地址我的配置策略FIFO_CONFIG0x06開啟滾動(dòng)模式允許新數(shù)據(jù)覆蓋舊數(shù)據(jù)防止溢出停采樣本平均設(shè)到 4可以平滑信號(hào)毛刺MODE_CONFIG0x07配成 SpO2 模式讓紅光和紅外 LED 交替工作SPO2_CONFIG0x08采樣率配到 100Hz這個(gè)頻率對(duì)心率信號(hào)充足ADC 量程和脈寬按手上手冊(cè)的位定義來選LED1_PA0x09紅光 LED 電流初始給一個(gè)中間值后面按波形調(diào)LED2_PA0x0A紅外 LED 電流同樣先給中間值配置代碼大致這樣static int32_t Max30100Config(void) { int32_t ret; ret Max30100WriteReg(0x06, 0x50); // 滾動(dòng)使能 樣本平均4 if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x07, 0x03); // SpO2模式 if (ret ! HDF_SUCCESS) return ret; // 采樣率等參數(shù)按數(shù)據(jù)手冊(cè)的位定義設(shè)置這里不寫死 ret Max30100WriteReg(0x08, sampleRateConfig); if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x09, g_max30100.redLedCurrent); if (ret ! HDF_SUCCESS) return ret; ret Max30100WriteReg(0x0A, g_max30100.irLedCurrent); if (ret ! HDF_SUCCESS) return ret; return HDF_SUCCESS; }配置完成后我習(xí)慣立刻回讀一遍確認(rèn)寄存器真的寫進(jìn)去了。有些時(shí)候 I2C 傳輸返回成功但設(shè)備端因?yàn)闀r(shí)序問題沒有真正接收回讀是最直接的驗(yàn)證手段。LED 電流的選擇值得單獨(dú)說電流太小AC 分量被噪聲淹沒電流太大接收管飽和ADC 讀數(shù)頂?shù)綕M量程血氧比值完全失真。我調(diào)這個(gè)花了很長時(shí)間最終原則是讓兩個(gè)通道的靜態(tài) DC 值落在 ADC 量程的中間偏下位置比如 16 位 ADC 的話讓 DC 落在 2000 到 3000 左右然后在這個(gè)基礎(chǔ)上微調(diào) AC 幅度。3.3 FIFO讀取與環(huán)形緩沖采集線程怎么設(shè)計(jì)MAX30100 內(nèi)部有個(gè) FIFO傳感器采集的數(shù)據(jù)先存在 FIFO 里驅(qū)動(dòng)再去讀。在 SpO2 模式下每個(gè)樣本包含紅光和紅外兩個(gè) 16 位 ADC 值即 4 字節(jié)。讀取時(shí)先看寫指針和讀指針?biāo)愠霎?dāng)前可讀幾個(gè)樣本再連續(xù)讀取。#define MAX30100_REG_FIFO_WR_PTR 0x02 #define MAX30100_REG_FIFO_RD_PTR 0x04 #define MAX30100_REG_FIFO_DATA 0x05 static int32_t Max30100ReadPpg(uint16_t *redArray, uint16_t *irArray, uint8_t *count) { uint8_t wr; uint8_t rd; uint8_t raw[32]; uint8_t available; int32_t ret; int32_t i; ret Max30100ReadRegs(MAX30100_REG_FIFO_WR_PTR, wr, 1); if (ret ! HDF_SUCCESS) return ret; ret Max30100ReadRegs(MAX30100_REG_FIFO_RD_PTR, rd, 1); if (ret ! HDF_SUCCESS) return ret; available (uint8_t)(((wr 32 - rd) 0x1F) / 4); if (available 0) { *count 0; return HDF_SUCCESS; } if (available 8) { available 8; // 單次最多拿8個(gè)樣本避免長時(shí)間占用I2C } ret Max30100ReadRegs(MAX30100_REG_FIFO_DATA, raw, (uint16_t)(available * 4)); if (ret ! HDF_SUCCESS) return ret; for (i 0; i available; i) { uint8_t *p raw[i * 4]; redArray[i] (uint16_t)((p[0] 8) | p[1]); irArray[i] (uint16_t)((p[2] 8) | p[3]); } *count available; return HDF_SUCCESS; }這里有個(gè)小細(xì)節(jié)計(jì)算可用樣本數(shù)時(shí)((wr 32 - rd) 0x1F)是把寫指針和讀指針放在同一個(gè)環(huán)形空間里求差避免指針回繞時(shí)相減出負(fù)數(shù)或出大數(shù)。MAX30100 的 FIFO 深度是 32 個(gè)樣本單位SpO2 模式下每個(gè)樣本點(diǎn)占 4 個(gè)單元所以再除以 4 得到可讀樣本數(shù)。采集線程我用了最直白的輪詢方式每隔 10ms 調(diào)一次Max30100ReadPpg把數(shù)據(jù)寫進(jìn)一個(gè)環(huán)形緩沖區(qū)。環(huán)形緩沖大小我設(shè)置成 128 組樣本大約是 1 秒多鐘的數(shù)據(jù)量足夠應(yīng)用層做算法運(yùn)算。生產(chǎn)者和消費(fèi)者的邏輯很簡單typedef struct { uint16_t red[128]; uint16_t ir[128]; uint16_t head; uint16_t tail; } PpgRingBuffer;采集線程向 head 寫應(yīng)用層從 tail 讀。如果應(yīng)用層來不及消費(fèi)就覆蓋最舊的數(shù)據(jù)保證拿到的永遠(yuǎn)是最新的一段波。數(shù)據(jù)規(guī)模不大用數(shù)組實(shí)現(xiàn)比鏈表簡單可靠得多也不會(huì)因?yàn)閮?nèi)存分配產(chǎn)生碎片問題。3.4 把數(shù)據(jù)送給應(yīng)用層輪詢、中斷與服務(wù)接口的選擇數(shù)據(jù)最后怎么交給應(yīng)用層我前前后后試過三種方式。最簡單的是輪詢驅(qū)動(dòng)里的采集線程負(fù)責(zé)讀 FIFO應(yīng)用層按需調(diào)用查詢接口從環(huán)形緩沖里取數(shù)據(jù)。這種方式代碼少100Hz 采樣率下足夠穩(wěn)我就先用它跑通了整條鏈路。第二種是用 INT 引腳。MAX30100 的 INT 在 FIFO 數(shù)據(jù)達(dá)到一定量時(shí)拉低接到開發(fā)板的 GPIO 中斷上驅(qū)動(dòng)收到中斷以后再去讀。好處是功耗低、不需要固定周期喚醒適合做低功耗穿戴設(shè)備的場景。代價(jià)是要多申請(qǐng)一個(gè) GPIO 中斷而且中斷服務(wù)函數(shù)里不能做太重的事通常是置一個(gè)標(biāo)志喚醒采集線程去干活。第三種是做成 HDF 服務(wù)由上層應(yīng)用綁定服務(wù)路徑后調(diào)用。這是最規(guī)范的做法但牽扯到 HDF 服務(wù)接口定義、權(quán)限配置、應(yīng)用側(cè)跨進(jìn)程調(diào)用工程量會(huì)大不少。我的經(jīng)驗(yàn)是先把數(shù)據(jù)鏈路調(diào)通再考慮要不要上服務(wù)框架。4. 原始PPG到手之后算法與驅(qū)動(dòng)的分工邊界驅(qū)動(dòng)把連續(xù)的 PPG 原始數(shù)據(jù)送到環(huán)形緩沖之后接下來才是真正意義的數(shù)據(jù)處理。這一章重點(diǎn)講算法和驅(qū)動(dòng)的分工因?yàn)槲以谶@個(gè)環(huán)節(jié)走過彎路一開始想在驅(qū)動(dòng)里把血氧和心率全算完結(jié)果耦合太重、調(diào)參也麻煩后面全部推倒重來。4.1 原始PPG波形長什么樣先學(xué)會(huì)畫圖在寫任何算法之前先把數(shù)據(jù)畫出來看。最笨也最有效的方法是串口按 CSV 格式打印采樣序號(hào)、紅光 ADC、紅外 ADCPC 端用 Python 腳本讀取并畫波形。正常的 PPG 波形長這樣每個(gè)通道都是一個(gè)接近直流的基線基線值在幾千左右上面疊加小幅度的、周期性的起伏每個(gè)起伏對(duì)應(yīng)一次心跳。紅光通道和紅外通道的起伏應(yīng)該基本同步但幅度可能不同。如果你看到兩個(gè)通道的起伏不同步、形狀差異很大說明某個(gè)通道飽和了或者 LED 電流配比不對(duì)。這一步?jīng)]有捷徑。我見過很聰明的同學(xué)直接跳過畫圖開始寫血氧算法最后調(diào)了兩天也不知道是算法問題還是信號(hào)問題。信號(hào)都沒看明白算法就是盲人摸象。4.2 心率計(jì)算峰值檢測與帶通濾波心率算法我放在應(yīng)用層做驅(qū)動(dòng)層只保證數(shù)據(jù)連續(xù)、帶時(shí)間戳。算法思路分三步濾波、峰檢測、間隔換算。先對(duì)紅外通道做帶通濾波。心率范圍一般在 30 到 200 次每分鐘換算成頻率大約是 0.5Hz 到 3.3Hz用一個(gè)兩步的 IIR 濾波器把直流和低頻漂移濾掉同時(shí)把 5Hz 以上的高頻噪聲去掉。濾波之后基線被拉平剩下的是清晰的脈搏波。然后是峰值檢測。我最常用的是一階差分加滑動(dòng)窗口// 偽代碼尋找局部最大值并限制最小波峰間隔 if (current prev current next current threshold (now - lastBeatTime) minInterval) { heartRate 60000 / (now - lastBeatTime); lastBeatTime now; }這個(gè)思路簡單但在信號(hào)質(zhì)量不錯(cuò)的場景下很有效。minInterval用來防止把同一波形的兩個(gè)小毛刺誤判成兩次心跳我一般設(shè) 300ms對(duì)應(yīng)最大 200bpm。注意心率輸出不要每拍都更新正確做法是取過去 4 到 8 拍的間隔做滑動(dòng)平均否則讀數(shù)會(huì)跳動(dòng)得很厲害。4.3 SpO2計(jì)算R值公式和查表思路血氧計(jì)算比心率稍微繞一點(diǎn)但原理不復(fù)雜。對(duì)紅光和紅外兩個(gè)通道分別提取交流分量 AC 和直流分量 DC計(jì)算比值 RR (AC_RED / DC_RED) / (AC_IR / DC_IR)然后通過經(jīng)驗(yàn)公式或查表把 R 映射到 SpO2。標(biāo)準(zhǔn)模型下R 值在 0.4 附近對(duì)應(yīng)接近 100% 的血氧R 為 1.0 左右大約對(duì)應(yīng) 94%R 為 2.0 時(shí)大約降到 75% 附近。這個(gè)區(qū)間可以做個(gè)粗略對(duì)照R值參考SpO20.4約100%0.6約97%1.0約94%1.5約85%2.0約75%注意這張表只能作為演示示意不能直接用于醫(yī)療用途。不同膚色、不同 LED 波長偏差、不同佩戴松緊都會(huì)改變映射曲線做到產(chǎn)品級(jí)必須用自己的硬件和采集環(huán)境做多用戶標(biāo)定。實(shí)際計(jì)算中AC 不能直接取原始數(shù)據(jù)的振幅一般是用帶通濾波后的信號(hào)按一個(gè)滑動(dòng)窗口算均方根或平均峰峰值。窗口太長響應(yīng)慢窗口太短數(shù)值抖動(dòng)大我最終選了 4 秒窗口每 1 秒輸出一次平衡了更新速度和穩(wěn)定性。4.4 為什么我不建議在驅(qū)動(dòng)層做算法我知道有人會(huì)覺得既然驅(qū)動(dòng)里已經(jīng)有數(shù)據(jù)了順手就把心率算了多方便。我確實(shí)這么干過然后果斷放棄了。原因是算法開發(fā)是高頻迭代的濾波器的截止頻率、窗口長度、閾值參數(shù)幾乎每天都要調(diào)。如果這些邏輯塞在驅(qū)動(dòng)里每次調(diào)參都要重新編譯整個(gè)驅(qū)動(dòng)鏡像、燒錄、重啟設(shè)備效率極低。驅(qū)動(dòng)層的正確姿勢是穩(wěn)定采數(shù)、打時(shí)間戳、管理緩沖、提供干凈的接口。算法放應(yīng)用層可以單獨(dú)編譯、單獨(dú)調(diào)試甚至可以在 PC 上用同樣的算法處理串口傳出來的數(shù)據(jù)。這個(gè)架構(gòu)上的分界比優(yōu)化任何一個(gè)具體算法細(xì)節(jié)都更重要。5. 想拿到穩(wěn)定的數(shù)據(jù)我踩過并解決過的坑最后這一章我把整個(gè)開發(fā)周期里遇到的最典型的幾個(gè)問題列出來。這些問題在網(wǎng)上都能搜到零散答案但很少有文章把排查思路講清楚我盡量還原當(dāng)時(shí)的判斷過程。5.1 設(shè)備ID讀出來全FF先查時(shí)序再查I2C第一次上電讀 PART_ID 返回 0xFF也就是總線上根本沒應(yīng)答。我當(dāng)時(shí)第一反應(yīng)是模塊壞了換了一片還是同樣現(xiàn)象才開始系統(tǒng)排查。排查順序是這樣的先查 I2C 地址確認(rèn)發(fā)的是 0xAE 寫地址再查 HCS 配置里的總線號(hào)會(huì)不會(huì)配到了一條沒有設(shè)備的空總線上然后用邏輯分析儀看 SCL/SDA 波形發(fā)現(xiàn) SDA 根本沒有被拉低說明從設(shè)備沒有拉總線最后檢查模塊板上的上拉電阻發(fā)現(xiàn)模塊自帶上拉但因?yàn)榻恿艘桓容^長的杜邦線加上開發(fā)板內(nèi)部也有上拉組合起來電平反而異常。換用短線和獨(dú)立 4.7k 上拉后ID 正常讀回。還有一個(gè)容易忽略的原因剛上電立刻讀 ID芯片還沒準(zhǔn)備好。我后來在 Init 里先OsalMSleep(10)再讀成功率明顯提升。讀到 ID 值之后要和你手上的芯片批次資料核對(duì)不同批次的 ID 可能會(huì)不一樣不同型號(hào)的兼容芯片也可能返回不同數(shù)值。5.2 FIFO溢出或長期為空采樣率與讀出節(jié)奏的平衡調(diào)試中我遇到過兩個(gè)極端一個(gè)現(xiàn)象是 FIFO 溢出標(biāo)志一直置起另一個(gè)現(xiàn)象是讀出來的樣本永遠(yuǎn)是 0。這兩個(gè)問題其實(shí)是同一類問題采樣率和讀取節(jié)奏失配。溢出很好理解采樣率 800Hz但上層讀取線程每隔 50ms 才讀一次FIFO 早就滿了。如果沒有開滾動(dòng)模式芯片會(huì)停止采樣數(shù)據(jù)就斷了即使在滾動(dòng)模式數(shù)據(jù)也會(huì)大量覆蓋丟失。解決方法是把采樣率降到 100Hz讀取線程保證至少每 20ms 讀一次。我之前用消息隊(duì)列間接讀延遲漲到幾十毫秒后數(shù)據(jù)就開始出現(xiàn)空洞后來直接用一個(gè)高優(yōu)先級(jí)的采集線程輪詢問題消失。長期為空的原因則更隱蔽。起初我懷疑是寄存器配置問題反復(fù)檢查都沒發(fā)現(xiàn)錯(cuò)誤。后來發(fā)現(xiàn)采集線程創(chuàng)建成功但一直阻塞在等待一個(gè)應(yīng)用層才會(huì)發(fā)送的信號(hào)量上導(dǎo)致 FIFO 里的數(shù)據(jù)永遠(yuǎn)沒人讀。這是一個(gè)純粹的調(diào)度設(shè)計(jì)錯(cuò)誤和傳感器本身無關(guān)。這也提醒我驅(qū)動(dòng)里的采集線程和上層消費(fèi)邏輯盡量解耦采集線程只依賴自己的定時(shí)器不要依賴外部信號(hào)。5.3 SpO2數(shù)值亂跳信號(hào)質(zhì)量評(píng)估比算法更關(guān)鍵血氧值亂跳是幾乎所有健康監(jiān)測設(shè)備都會(huì)遇到的問題。剛開始我的 R 值算出來一會(huì)兒 0.3、一會(huì)兒 2.5血氧從 100% 直接掉到 75%但人的狀態(tài)完全正常。仔細(xì)觀察原始波形后才明白是手指按壓力度變化和微小晃動(dòng)導(dǎo)致的運(yùn)動(dòng)偽跡。解決這個(gè)問題算法上最高效的不是優(yōu)化公式而是加信號(hào)質(zhì)量評(píng)估。我加了三個(gè)判斷條件一是 R 值必須落在 0.3 到 2.5 的合理區(qū)間內(nèi)超出就不輸出二是紅光和紅外通道的波峰間隔要基本一致差異超過 20% 就視為信號(hào)不可靠三是 AC 分量占比不能太低AC 太小說明信號(hào)太弱通常是因?yàn)槭种笡]有貼緊或 LED 電流不足。經(jīng)過這三道閘以后穩(wěn)定坐姿下血氧讀數(shù)能連續(xù)幾分鐘保持誤差在一個(gè)很小的范圍內(nèi)。運(yùn)動(dòng)狀態(tài)下想做到完全穩(wěn)定那就要上加速度傳感器輔助做運(yùn)動(dòng)偽跡消除已經(jīng)超出這顆芯片本身的能力范圍了。5.4 調(diào)試方法串口日志加Python繪圖具體調(diào)試手段上我推薦一套投入小、見效快的組合驅(qū)動(dòng)把數(shù)據(jù)以文本幀格式從串口發(fā)出來PC 端用 Python 的 pyserial 接收再用 matplotlib 實(shí)時(shí)繪圖。串口輸出格式我做成這樣frame_no, red, ir 123, 2400, 2800 124, 2420, 2820驅(qū)動(dòng)里注意不要在 I2C 讀數(shù)據(jù)的過程中直接打印日志串口打印很慢會(huì)把采樣節(jié)奏拖垮。正確做法是先把本次讀到的數(shù)據(jù)存到環(huán)形緩沖到主循環(huán)外面再按批量輸出。Python 側(cè)我每次只畫最近 10 秒的波形這樣既能看細(xì)節(jié)又不會(huì)糊成一團(tuán)。畫圖的時(shí)候重點(diǎn)關(guān)注三個(gè)指標(biāo)DC 基線是否穩(wěn)定、AC 幅度是否在合理范圍、兩個(gè)通道的波峰是否同步。這三個(gè)指標(biāo)正常再去看心率血氧計(jì)算結(jié)果否則先回頭查硬件和配置。還有一個(gè)小技巧用串口輸出的樣本序號(hào)檢查丟幀。如果序號(hào)不連續(xù)說明驅(qū)動(dòng)讀取側(cè)有丟數(shù)據(jù)這時(shí)候不管算法多好最終結(jié)果都會(huì)受影響。先把序號(hào)補(bǔ)成全連續(xù)再接算法這條調(diào)試順序我后來一直遵守。如果問我在這個(gè)項(xiàng)目里最大的體會(huì)其實(shí)不是學(xué)會(huì)了怎么配 HDF、怎么寫 I2C 封裝而是意識(shí)到驅(qū)動(dòng)開發(fā)真正的難點(diǎn)是數(shù)據(jù)質(zhì)量而不是代碼框架。MAX30100 是一顆很便宜的傳感器但想拿到穩(wěn)定的 PPG 數(shù)據(jù)從電源、布局、LED 電流到讀取節(jié)奏每一個(gè)環(huán)節(jié)都在影響最終結(jié)果。我建議你動(dòng)手時(shí)先別急著寫算法把波形打印出來盯一兩個(gè)晚上信號(hào)看著順眼了心率血氧自然就出來了。后面有空我再寫一篇應(yīng)用層怎么通過 HDF 服務(wù)讀取驅(qū)動(dòng)數(shù)據(jù)以及把實(shí)時(shí)波形傳到手機(jī)端調(diào)試的做法。