心電信號噪聲的數(shù)字濾波研究

1、技術(shù)創(chuàng)新中文核心期刊微計算機信息(管控一體化2008年第24卷第6-3期360元/年郵局訂閱號:82-946現(xiàn)場總線技術(shù)應(yīng)用200例軟件時空心電信號噪聲的數(shù)字濾波研究Research of Digital Filtering Approach to ECG Signal Noise(北京工商大學(xué)陳天華韓力群鄭彧CHEN Tian-huaHAN Li-qun ZHENG Yu摘要:心血管疾病的正確診斷依賴于心電信號的準確獲取,提出了一種消除工頻干擾和基線漂移的數(shù)字濾波方法,該方法運算量小、設(shè)計簡單、易于實現(xiàn)。實驗表明該方法能有效濾除工頻干擾和抑制基線漂移,獲得真實的心電信號。關(guān)鍵詞:心電信號;數(shù)

2、字濾波;基線漂移;噪聲消除;數(shù)字信號處理中圖分類號:TN911文獻標識碼:A Abstract:A correct diagnosis of cardiac vascular disease rely on the precise ECG signal,this paper presented a new kind of design method of digital filter in order to eliminate the 50Hz noise and baseline drifts.The method have the characteristics of small calc

3、ula-tion amount,simple and realization easily.the experiment shows that it is effective in eliminating 50Hz noise and baseline drifts,es-pecially in strong noise,the real ECG signal can be obtained.Key words:ECG;Digital filter;baseline drifts;De-noise;Digital signal processing文章編號:1008-0570(200806-3

4、-0252-03現(xiàn)代醫(yī)學(xué)表明,心電信號(ECG含有臨床診斷心血管疾病的大量信息,ECG 的檢測與分析臨床診斷中具有重要價值,是了解心臟的功能與狀況、輔助診斷心血管疾病、評估各種治療方法有效性的重要手段。但由于實際檢測工況的非理想,ECG 信號中往往含有工頻噪聲及電極極化等引起的各種隨機噪聲。噪聲的存在降低了診斷的準確性,其中影響最大的是50Hz 工頻干擾和基線漂移噪聲。因此,在ECG 信號檢測過程中,如何抑制工頻干擾和基線漂移是必須解決的問題。1數(shù)字濾波系統(tǒng)心電信號是微弱低頻人體生理電信號,通常頻率在0.05100Hz,幅值不超過4mv,通過安裝在皮膚表面的電極來獲取。由于人體是一個復(fù)雜的生命

5、系統(tǒng),存在50H 工頻干擾及基線漂移等其他生理電信號的干擾。噪聲可能會影響到醫(yī)生的臨床診斷,因此,需對心電信號進行濾波,即必須做好前端數(shù)據(jù)采集的軟硬件設(shè)計以保證心電數(shù)據(jù)的可靠和準確。傳統(tǒng)醫(yī)療設(shè)備分別采用50Hz 帶阻濾波器和RC 高通濾波器濾除工頻干擾和基線漂移。但帶阻濾波器電路復(fù)雜,其特性對元器件的精度敏感,而基線漂移本質(zhì)上是一種緩慢變化的低頻信號,采用RC 濾波器很難將高通濾波器的過渡帶做得十分陡峭,基線漂移補償效果不理想。因此,模擬方法往往不太容易獲得很好的特性。數(shù)字濾波方法有具有許多優(yōu)點,ECG 數(shù)字濾波系統(tǒng)組成如圖1所示,來自各電極的多路心電微弱信號經(jīng)多路輸入緩沖器緩沖放大,經(jīng)導(dǎo)聯(lián)選

6、擇電路進行選擇后,由前置放大電路放大,進行電壓放大以適應(yīng)A/D 轉(zhuǎn)換的幅度要求,然后進行數(shù)字濾波并輸出心電信號。數(shù)字濾波用于心電信號消噪,不僅能提高儀器設(shè)備的性能,而且對于不同的使用環(huán)境(例如對于不同國家的50Hz 或60Hz 供電條件,只需重新設(shè)置軟件參數(shù)即可,大大降低了醫(yī)療設(shè)備硬件復(fù)雜程度,降低了產(chǎn)品成本、提高了通用性。圖1心電信號數(shù)據(jù)采集與處理框圖2工頻干擾濾波對于心電信號中的工頻干擾,簡單而有效的方法是采用梳狀濾波器濾波,這種FIR 濾波器簡單、容易實現(xiàn)、濾波效果好,節(jié)數(shù)為N 的梳狀濾波器的系統(tǒng)函數(shù)為:(1其中,是由N 節(jié)延時單元組成的梳狀濾波器,H 1(z的幅值響應(yīng)由許多頻率間隔相同

7、的通帶和阻帶組成,它只許一些特定頻率范圍的信號通過而阻止另一些特定頻率信號通過。H(z的頻率響應(yīng)為:(2根據(jù)式(2可知,梳狀濾波器是具有線性相位的FIR 低通濾波器,相移=(N-1/2,為使濾波器的直流增益為1,可以考慮在H(z中增加增益因子1/N,歸一化后的幅頻響應(yīng)為:(3N=10時,幅頻響應(yīng)如圖2所示,在頻率點k =2k/N,k=0,1,2,N-1處的幅值為0,因而在k =2k/N 附近形成了多個阻帶,對這些頻率的信號具有很好的抑制作用。由于梳妝濾波器的系數(shù)相等且都為1,因而容易實現(xiàn),常用于對醫(yī)療儀器中的特定信號進行濾波處理。陳天華:碩士副教授基金項目:北京市教委科技發(fā)展基金(KM2007

8、10011010252- 郵局訂閱號:82-946360元/年技術(shù)創(chuàng)新軟件時空PLC 技術(shù)應(yīng)用200例您的論文得到兩院院士關(guān)注圖2線性相位FIR 濾波器的幅頻響應(yīng)特性根據(jù)圖2可得第一旁瓣峰值衰減約為20dB,若不能滿足要求,可對式(1進行修正:(4其頻率響應(yīng)為: (5衰減速度為修正前的k 倍,k=2時修正濾波器的幅頻響應(yīng)如圖3所示,其衰減性能更好。圖3經(jīng)修正以后的梳狀濾波器幅頻響應(yīng)3基線漂移的抑制抑制基線漂移方法很多,各有特點?；趻佄锞€的擬合基線算法結(jié)構(gòu)復(fù)雜,采用單片機的嵌入式系統(tǒng)運算速度難以保證實時性要求,基于小波的自適應(yīng)濾波抑制ECG 基線漂移,運算較復(fù)雜,不適合用于小系統(tǒng)的監(jiān)護。本文采

9、用一種自適應(yīng)快速線性擬合來抑制基線漂移的方法,通過多點采樣,用最小二乘法擬合出基線的變化趨勢,經(jīng)修正得到ECG 信號。3.1QRS 波群的探測理論上,兩個連續(xù)QRS 波群之間的T-P 段代表了實際的ECG 基線水平。根據(jù)每一個心動周期的P-R 段特征點,在兩個特征點之間進行插值,然后再進行曲線擬合,即將插值和曲線擬合結(jié)合起來,就可以得到基線漂移曲線,原始信號減去基線漂移曲線即為ECG 信號。因此,采用分析斜率、幅度和寬度的方法來識別QRS 波群。根據(jù)ECG 的功率譜,QRS 波群的能量主要集中在15Hz 頻率左右,為消除ECG 中其它頻率成分對QRS 波群檢測的影響,需設(shè)計一個中心頻率為15H

10、Z 的帶通濾波器,ECG 信號通過該濾波器時對R 波以外的頻率成份進行了衰減。當(dāng)采樣頻率為500Hz 時,該濾波器的傳遞函數(shù)如下:(6由于R 波含有豐富的高頻諧波成分,為突出高頻諧波分量,采用如下微分器:經(jīng)處理后的信號通過帶通濾波器和微分器后,P 波和T 波都有顯著的衰減,相應(yīng)的QRS 波群的峰值進一步加強,為消除微分器處理后散粒隨機噪聲產(chǎn)生的干擾,可以采用式(7所表示的積分器進行修正:(7然后,采用斜率跟蹤法進行R 波識別,找出與QRS 波群檢測特征點前后60ms 相當(dāng)?shù)拈g期來識別R 波。當(dāng)采集波形的斜率變化超過設(shè)定的閾值時,斜率變化最大處即為R 波所在的位置。實驗的采樣時間為4.8秒,采樣

11、點數(shù)為2048個點。許多異常的QRS 波群有高幅度和上升不迅速的特點,為提高以R 波斜率檢測QRS 波群的可靠性,可采用移動窗口積分法驗證以保證可靠檢測。為便于進行積分運算,先將采樣點進行逐點平方運算,并進行移動窗口積分,其差分方程如下:N 為移動窗口寬度的采樣點數(shù)窗口寬度應(yīng)當(dāng)盡可能與QRS 波群寬度相近。若窗口太大,積分波形將把QRS 波群和T 波融合在一起,若窗口尺度太小,一個QRS 波群在輸出中可能會產(chǎn)生幾個峰值。3.2自適應(yīng)線性擬合在抑制基線漂移的算法中,須首先確定基線的變化趨勢,為此應(yīng)在兩個心動周期中找到一段反映ECG 基線變化趨勢的曲線。QRS 波群是心電信號的高頻部分,其中R 波

12、具有波形陡峭、幅度大、寬度窄、變化趨勢明顯等特點,因而以R 波為基準點尋找代表基線的時段。依次取每個檢測的R 波峰值點作為參考基準點,利用回溯法求出每個TP 間期。確定了曲線區(qū)間就可采用自適應(yīng)快速線性擬合法,擬合出基線的變化趨勢。采用最小二乘法計算各區(qū)間線性擬合的斜率和截距,擬合出基線的變化趨勢。對于觀測點(x i ,y i i=1,2n,若y 與x 具有線性模型,則任意一對觀測值(x i ,y i 存在如下關(guān)系:式中,i 為理想狀態(tài)和觀測值y i 的誤差,令:顯然,f(k,b越小,結(jié)果就越精確,因此根據(jù)最小二乘法可推導(dǎo)出(8若斜率k>0表示往上漂,k<0往下漂,k=0則沒有漂移。

13、根據(jù)分段修正所求出的k 值,判斷基線的變化趨勢,從而修正每段基線漂移值,得到真實穩(wěn)定的ECG 信號。4濾波效果圖4表明,未經(jīng)濾波的心電信號存在較大的工頻噪聲,若取采樣頻率500Hz,梳狀濾波器節(jié)數(shù)N=10,經(jīng)濾波后心電信號的干擾被基本濾除。圖5為ECG 信號頻譜,濾波前ECG 信號存在較強的工頻干擾,濾波后頻譜中已經(jīng)沒有50Hz 成分。數(shù)字濾波對253-技術(shù)創(chuàng)新中文核心期刊微計算機信息(管控一體化2008年第24卷第6-3期360元/年郵局訂閱號:82-946現(xiàn)場總線技術(shù)應(yīng)用200例軟件時空消除50Hz 工頻及倍頻干擾具有很好的效果。圖4-1含有工頻噪聲的心電信號圖4-2FIR 梳狀濾波處理后

14、的心電信號圖5-1濾波前ECG 頻譜圖5-2濾波后ECG 頻譜圖6-1是存在基線漂移(上漂的ECG 信號,采用自適應(yīng)擬合算法抑制效果如圖6-2所示;圖7-1是存在下漂的ECG 信號,抑效果如圖7-2所示。實驗表明,對于不同形式的基線漂移,該方法能很好的抑制ECG 信號基線漂移。與其他方法相比,自適應(yīng)擬合算法簡單、快速、易于實現(xiàn)、抑制效果好、信號失真度小。在嵌入式系統(tǒng)和對實時性要求較高的監(jiān)護儀器中具有實際應(yīng)用價值,也可用于其它信號檢測儀器中。圖6-1基線上漂的心電信號5結(jié)論采用數(shù)字濾波方法消除心電信號噪聲已經(jīng)成為心電信號處理的重要內(nèi)容。本文的創(chuàng)新點是用梳狀濾波器和自適應(yīng)擬合算法有效消除了心電信號

15、中的工頻噪聲和基線漂移,獲得真實ECG 信號。本方法設(shè)計簡單、運算量小,在抑制工頻干擾和基線漂移的同時,能很好的保持原有信號,降低了成本、提高了產(chǎn)品通用性。圖6-2抑制基線漂移后的心電信號圖7-1基線下漂的心電信號圖7-2抑制基線漂移后的心電信號參考文獻:1張立新,丁北生等.用IIR 算法消除ECG 基線漂移J.天津大學(xué)學(xué)報.1998.1;31(1:120124.2王建軍,陳日新等.數(shù)字信號處理在醫(yī)學(xué)上的應(yīng)用J.計算技術(shù)與自動化.2000.6.19(2:4547.3付宜利,王國光,靳保.基于FPGA 的高速實時FFT 處理器設(shè)計J.微計算機信息.2007.2,23(2:194195.4Willis J.Tompkins.Biomedical Digital Signal Processing M,New Jersey:Prentice Hall,2712725高清維,李海鷹等.基于平穩(wěn)小波變換的心電信號噪聲消除方法J.電子學(xué)報.2003.2;31(2:238240.作者簡介:陳天華,男,湖南長沙人,碩士,南京航空航天大學(xué)(本科及研究生畢業(yè),副教授,主要從事智能信息處理及計算機控制等方面的研究,近年來先后在國內(nèi)外重要刊物及EI 核心源刊上發(fā)表論文30多篇。Biography:Chen Tianhua,male,Changsha,