基于信號質(zhì)量評估的可穿戴動(dòng)態(tài)心電監護

目前，心血管疾病已經(jīng)成為危害人類(lèi)健康的主要疾病之一。可穿戴心電技術(shù)是在人們日常穿戴的衣物中嵌入心電采集系統，使其在自然狀態(tài)下隨時(shí)隨地獲取心電數據，是人體心電實(shí)時(shí)監護的有效方法。但由于人體處于運動(dòng)狀態(tài)，心電信號受到的干擾大，難以正確處理和評估。目前國內外有很多學(xué)者對動(dòng)態(tài)心電進(jìn)行了研究，例如盛虎提出的新型動(dòng)態(tài)心電監護系統彌補了信號采集功能和分析功能脫離的不足，實(shí)現了心電數據實(shí)時(shí)分析，但是該系統成本較高[1]；國外有些學(xué)者將卡爾曼濾波應用于動(dòng)態(tài)心電信號檢測，并與自適應算法結合，以提高檢測的正確率，但由于算法比較復雜，不能進(jìn)行實(shí)時(shí)監控[2]；李橋等人將卡爾曼濾波應用于危重病人監護，并結合了信號質(zhì)量評估，但主要針對靜態(tài)生理信號，不適用于動(dòng)態(tài)心電監護[3]。

本文提出了一種新的基于信號質(zhì)量評估和卡爾曼濾波的可穿戴動(dòng)態(tài)心電監護的方法，該方法將信號質(zhì)量評估、卡爾曼濾波以及可穿戴技術(shù)相結合，實(shí)現了便攜式高可靠的長(cháng)時(shí)間人體心率評估。

1 可穿戴動(dòng)態(tài)心電特征分析

心電信號作為心臟電信號在人體體表的表現，信號微弱，頻率主要介于0.01 Hz~100 Hz之間，其中與心率評估密切相關(guān)的QRS波群主要集中在0~(58±19)Hz，90%的頻譜能量集中在0.25 Hz~35 Hz之間，而高頻心電信號的頻帶范圍則在100 Hz~1 000 Hz[4]。在動(dòng)態(tài)心電信號采集的過(guò)程中，心電信號干擾中的肌電干擾、基線(xiàn)漂移（小于5 Hz）和運動(dòng)偽跡（小于7 Hz）均比靜態(tài)情況下嚴重。目前醫學(xué)上運用最廣泛的電極是傳統的氯化銀電極，雖然這種電極采集到的信號穩定，但是對皮膚損傷較大，不適合長(cháng)期使用?？纱┐鲃?dòng)態(tài)心電采集宜采用對人體無(wú)損害、能夠長(cháng)時(shí)間使用的織物柔性電極，但織物電極信號的自身特點(diǎn)，需要在處理時(shí)采用不同方法。

本文通過(guò)集成在智能服裝上的織物電極獲取運動(dòng)狀態(tài)下的人體心電信號，并進(jìn)行特征分析，作為建立相應處理模型的依據。靜止狀態(tài)與運動(dòng)狀態(tài)下心電信號的功率譜如圖1所示。從圖中可以看出當頻率在0 Hz~7 Hz時(shí)，運動(dòng)狀態(tài)下的功率要遠大于靜止狀態(tài)下的功率，這是因為在人體運動(dòng)時(shí)電極與人體的相對運動(dòng)劇烈，使得運動(dòng)偽跡干擾的影響增大，另外由人體呼氣引起的基線(xiàn)漂移干擾也同時(shí)增大。在QRS波群集中的頻段，運動(dòng)狀態(tài)下的功率略大于靜止狀態(tài)下的功率，這是由于人體運動(dòng)以及肌肉緊張收縮而引起的肌電干擾所致。

經(jīng)實(shí)驗證明，在走路狀態(tài)和慢跑狀態(tài)下用織物電極采集的信號，雖然受到了一定的干擾，但是波形并沒(méi)有失真，完全能夠用來(lái)進(jìn)行后期的處理。

2 基于信號質(zhì)量評估和卡爾曼濾波的動(dòng)態(tài)心率估計模型

準確的心率估計是可穿戴監護系統的最基本要求，R波檢測算法是獲得心率的最便利的方法，然而心率的檢測經(jīng)常受到干擾的影響而出現錯誤?？柭鼮V波算法簡(jiǎn)單，所需的數據量小，在沒(méi)有信號先驗知識的情況下仍能有效地估計信號趨勢的變化和噪聲干擾，是一種估計信號趨勢的有效方法。信號質(zhì)量指數 SQI(Signal Quality Index)可以實(shí)時(shí)地表征動(dòng)態(tài)心電信號的質(zhì)量，從而作為卡爾曼濾波器參數調節的依據，提高估計準確性。本文提出了一種新的基于信號質(zhì)量評估和卡爾曼濾波的動(dòng)態(tài)心率測量和評估的方法，首先對心電信號進(jìn)行R波提取并計算心率，接著(zhù)利用R波檢測和加速度計的結果來(lái)獲得運動(dòng)狀態(tài)下心電信號質(zhì)量指數SQI，然后根據SQI對卡爾曼濾波器參數進(jìn)行調節，最后對心率進(jìn)行重新估計，以獲得更為準確的結果。系統結構圖如圖2所示

2.1 動(dòng)態(tài)心電R波識別和心率測量

R波檢測算法必須具備準確性與快速性?xún)蓚€(gè)特點(diǎn)，而這兩者之間又存在著(zhù)矛盾。常用的算法有面積法、小波變換法、幅度法和斜率法等，這些算法都是基于心電信號中的R波具有幅度和斜率較大的特點(diǎn)實(shí)現的。其中面積法雖然準確率較高，但算法很復雜，而小波變換法速度較慢，不適合用于實(shí)時(shí)分析，因此本文采用的R波檢測算法是幅度法和斜率法。

幅度法算法簡(jiǎn)單、速度快，在干擾較小時(shí)有很高的準確率，但是這種算法比較容易受到心電信號中的大T波干擾，從而導致心率估計不準，如圖3所示。斜率法的抗干擾能力比幅度法強，準確率高，缺點(diǎn)是算法比較復雜，容易受高頻肌電干擾的影響。

識別出R波過(guò)后，本文以10 s為一個(gè)時(shí)間窗口，計算該窗口內R-R間隔的均值，并以此得到10 s內的心率測量值，公式為：

其中，N(k)表示以某一秒的起始時(shí)刻為中心左右各取5 s，在這10 s的窗口內兩種R波檢測算法檢測出的R波匹配的數目。NA表示幅度算法檢測出的R波的個(gè)數，NS表示斜率算法檢測出的R波的個(gè)數，S表示加速度計檢測出的人體運動(dòng)的劇烈程度，S的值介于0～1之間，當S等于0時(shí)表示人體處于靜止狀態(tài)，S的值越大表示人體運動(dòng)越劇烈。前面提到，幅度算法在干擾較小時(shí)對于R波的檢測有很高的準確性，而斜率算法的抗干擾能力要強于幅度算法。因此在體現人體運動(dòng)劇烈程度的S的值越小時(shí)，幅度算法檢測出的結果所占的權重越大；而當S 的值越大，即人體運動(dòng)程度越劇烈時(shí)，斜率算法檢測出的結果所占的權重越大。

由式（8）可知，心電信號質(zhì)量指數SQI的值介于0～1之間，接近或等于0表示心電信號質(zhì)量很差，而接近或等于1表示心電信號質(zhì)量很高。

2.4 基于SQI的卡爾曼濾波心率估計器參數調節

當心電信號由于人體運動(dòng)而受到干擾時(shí)，使用R波檢測算法得到的心率將出現一定誤差，本文提出在運動(dòng)狀態(tài)下對R波檢測算法得到的心率應用卡爾曼濾波進(jìn)行最佳估計，并研究了通過(guò)心電信號質(zhì)量指數SQI調節卡爾曼濾波器增益系數的方法。由于當心電信號在人體運動(dòng)時(shí)會(huì )受到干擾，這將使心電信號質(zhì)量指數發(fā)生變化，而這種干擾大小的變化表現為卡爾曼濾波器方程中測量噪聲協(xié)方差R的變化，因此根據SQI值調節測量噪聲協(xié)方差R：

由式（9）可知，當心電信號質(zhì)量相對較高(即SQI值較大)時(shí)，R的值相對較小，從而Kg(k)值較大，此時(shí)較多地由測量值Z(k)來(lái)調整心率值；反之，當心電信號質(zhì)量較低(即SQI值較小)時(shí)，R值相對較大，Kg(k)值減小，此時(shí)較多地由先驗估計來(lái)調整心率值，從而避免運動(dòng)干擾的影響。

2.5 卡爾曼濾波方程初始值的確定

在運用卡爾曼濾波進(jìn)行動(dòng)態(tài)心率估計之前，應確定濾波器方程中系數的初始值。假設人的心率不會(huì )發(fā)生突變，即k時(shí)刻的心率近似等于k-1時(shí)刻的心率，因此 A=1。另外，前面提到如果系統沒(méi)有控制量，U(k)可以為0。由于系統的測量值由R波檢測算法獲得，和心率直接對應，因此H=1。系統測量噪聲協(xié)方差 R=R0=1，該系數會(huì )隨著(zhù)心電信號的SQI值的變化而變化。系統過(guò)程協(xié)方差Q初值設為0.1，P(0|0)=1。

3 結果和討論

3.1 評價(jià)數據庫

本文所采用的評價(jià)數據庫均由埃德ML870 PowerLab 8/30數據采集系統獲得，該系統通過(guò)屏蔽線(xiàn)與傳統的氯化銀電極、集成在智能服裝上的織物電極相連，采集人體在靜止、走路、慢跑等狀態(tài)下的生理信號，采樣頻率為1 000 Hz，各10 min時(shí)間，然后經(jīng)過(guò)放大、濾波和整形等處理，最終在終端上進(jìn)行實(shí)時(shí)顯示并保存心電數據。

3.2 結果

為了驗證本文提出的模型，首先選取了2 min標準心電信號，并對其進(jìn)行心率計算，如圖4中（a）所示，以此結果作為之后的評價(jià)標準。隨后在該心電信號的40 s~60 s數據段添加干擾，添加了干擾后的心電信號如圖4(b)所示，并用兩種R波檢測算法和基于卡爾曼濾波的估計算法對其進(jìn)行心率估計。

由標準心電信號計算所得的心率如圖5(a)所示，從圖中可以看出，因為心電信號的質(zhì)量較高，兩種R波檢測算法的準確率很高，所以在這2 min內得到的心率比較穩定，沒(méi)有心率突變的情況發(fā)生。

在添加了干擾過(guò)后，由R波檢測算法得到的心率（測量值）和由基于卡爾曼濾波的心率估計算法得到的心率（估計值）如圖5（b）所示。從圖中可以看出，由于 40 s~60 s的心電信號數據段受到了干擾，因此使得R波檢測算法出現了較大的誤差，出現了心率突變的情況。而由基于卡爾曼濾波的心率估計算法得到的心率值并沒(méi)有因為干擾受到太大的影響，心率仍然保持平穩。因此，基于卡爾曼濾波的心率估計算法具有很強的抗干擾性，是一種較為理想的動(dòng)態(tài)心率估計算法。

為了驗證該模型能夠用于動(dòng)態(tài)心電監護，同樣選取了2 min慢跑時(shí)的心電信號來(lái)進(jìn)行測試。

測試結果如圖6所示，由圖中可知，直接由R波檢測算法得到的心率在某些時(shí)間段會(huì )因人體運動(dòng)產(chǎn)生的干擾而出現較大的誤差，而由卡爾曼濾波的心率估計算法得到的結果卻有明顯的改善，心率能夠保持平穩，不會(huì )因為運動(dòng)干擾而出現大的誤差。

利用幅值法和斜率法雖然能夠計算心率，但心電信號中的各種干擾會(huì )影響計算結果，信號的信噪比越低，心率計算的誤差會(huì )越大。本文提出的新的基于信號質(zhì)量評估的卡爾曼濾波估計算法，可以顯著(zhù)改善心率估計的誤差。

針對可穿戴動(dòng)態(tài)心電監護系統需要嵌入服裝、柔性化和低功耗的要求，本文進(jìn)行了基于信號質(zhì)量評估和卡爾曼濾波的可穿戴動(dòng)態(tài)心電監護系統的設計，針對可穿戴心電信號處理系統抗干擾的要求，研究了人體生理信號質(zhì)量評估方法，提出了基于信號質(zhì)量指數的卡爾曼濾波器的設計，提高了心率估計的準確性，實(shí)現了可穿戴的動(dòng)態(tài)心電監護，并最終通過(guò)實(shí)際的測試表明了該設計的有效性和可靠性。

參考文獻

[1] 盛虎．新型動(dòng)態(tài)心電監護系統的研究[D]．大連：大連交通大學(xué)，2006．

[2] ARNOLD M, MILNER X H R, WITTE H，et al.Adaptive AR modeling of nonstationary time series by means of Kalman filtering[J]．IEEE Transactions on Biomedical Engineering, 1998,45(5)：553-562．

[3] 李橋．危重病人生命體征信號質(zhì)量評估與分析[D]．濟南：山東大學(xué)，2008．

[4] 石金蘭．面型智能服裝的亞健康評估系統的研究[D]．上海：東華大學(xué)，2008．

[5] 敬喜．卡爾曼濾波器及其應用基礎[M]．北京：國防工業(yè)出版社，1973．

更多醫療電子信息請關(guān)注：21ic醫療電子頻道