張喜紅　王玉香

(亳州職業(yè)技術(shù)學(xué)院， 安徽 亳州 236800)

?

基于EMD算法的心電信號基線漂移去除方法研究

張喜紅王玉香

(亳州職業(yè)技術(shù)學(xué)院， 安徽 亳州 236800)

摘要：采用經(jīng)驗(yàn)?zāi)B(tài)分解的方法，解決心電信號基線漂移時(shí)的最優(yōu)分解層數(shù)問題，通過合理選擇EMD分解層數(shù)來改進(jìn)基線漂移的去除方法。借助Matlab的工具箱，對MIT-BIH心律失常數(shù)據(jù)庫中較為潔凈的心電數(shù)據(jù)段，人為疊加心電常見噪聲信號以產(chǎn)生實(shí)驗(yàn)仿真信號，將所構(gòu)建的仿真信號經(jīng)EMD(經(jīng)驗(yàn)?zāi)B(tài)分解)分解后，將分解所得的殘余項(xiàng)去除，從而實(shí)現(xiàn)心電信號的基線漂移去除。對比不同分解層數(shù)下基線漂移去除的效果，可看到EMD層數(shù)越大，殘余項(xiàng)所在層的基線漂移特征越加明顯。最后對數(shù)據(jù)庫中受到嚴(yán)重基線漂移干擾的109號心電數(shù)據(jù)段進(jìn)行基線漂移去除，驗(yàn)證本方法的可行性。

關(guān)鍵詞：經(jīng)驗(yàn)?zāi)B(tài)分解； 基線漂移； 心電信號

心電信號濾波方面的相關(guān)理論是醫(yī)用電子儀器專業(yè)學(xué)生的必修理論，為了更好地培養(yǎng)學(xué)生本領(lǐng)域的相關(guān)技能及知識，亳州職業(yè)技術(shù)學(xué)院在安徽省醫(yī)電省級特色與省級實(shí)習(xí)實(shí)訓(xùn)示范中心兩大項(xiàng)目建設(shè)過程中，自主研發(fā)了開放式心電示教實(shí)驗(yàn)實(shí)訓(xùn)平臺(tái)[1]。在心電信號采集過程中，難免會(huì)混入常見的多種噪聲信號，因此在分析心電信號時(shí)需對采集的心電信號進(jìn)行各種濾波處理。基線漂移是心電信號采集過程中最常見的噪聲信號之一，它是人體呼吸波動(dòng)和電極運(yùn)動(dòng)引起的低頻干擾，其頻率通常小于5 Hz。因心電信號自身包含較多有用的低頻信號，心電信號因基線漂移的存在將會(huì)淹沒其有用的信息，一定程度上影響了診斷的準(zhǔn)確性[2]。基于以上客觀背景，在Matlab平臺(tái)下研究采用經(jīng)驗(yàn)?zāi)B(tài)分解的方法，解決基線漂移時(shí)的相關(guān)問題。

1材料及方法

1.1材料

實(shí)驗(yàn)中，首先在MIT-BIH數(shù)據(jù)庫中截取一段高信噪比原始心電信號[3]。圖1所示為原始心電信號。圖2所示為仿真心電信號。在Matlab中，通過式(1)，構(gòu)建在采集環(huán)境相對惡劣、采集設(shè)備抗噪聲能力較差情況下的仿真心電信號：

y(n)=x(n)+0.2sin(100πt)+

0.4sin(πt)+g(t)

(1)

式中：x(n) —— 截取的原始心電信號；

g(t) —— 高斯白噪聲，模擬肌電干擾；

0.2 sin(100 πt) —— 模擬50 Hz 的工頻噪聲，其幅度為0.2；

0.4 sin(πt) —— 模擬基線漂移干擾，其幅度為0.4；

y(n) —— 仿真信號輸出。

圖1　原始心電信號

圖2　仿真心電信號

1.2方法

1.2.1EMD基本理論

EMD即經(jīng)驗(yàn)?zāi)B(tài)分解，是由Huang等人在1998年提出的一種信號處理分解算法，主要適用于非平穩(wěn)信號的分解處理。人體生理信號多具有非平穩(wěn)信號的特征，生物醫(yī)學(xué)工程研究人員對此算法給予了較多關(guān)注[4]。EMD算法的主要過程就是從復(fù)雜的非平穩(wěn)信號中分解并提取本征模態(tài)函數(shù)(IMF)，所提取的本征模態(tài)函數(shù)(IMF)的定義是必須滿足2個(gè)條件：(1)在整個(gè)數(shù)據(jù)序列集中，含有相同數(shù)量的極值點(diǎn)與過零點(diǎn)，或者相差的個(gè)數(shù)應(yīng)為0個(gè)或1個(gè)；(2)在每個(gè)時(shí)間點(diǎn)，由信號序列局部極大值所形成的包絡(luò)和由信號序列局部極小值所形成的包絡(luò)均值必須為零。EMD分解是通過3次樣條插值，將信號變換為式(2)所示的形式[3]：

(2)

式中：x(t) —— 原始信號；

N—— 分解出的本征模態(tài)函數(shù)(IMF)個(gè)數(shù)；

rn(t)—— 分解得到的殘余項(xiàng)。

1.2.2實(shí)驗(yàn)方法及步驟

研究過程中，在Matlab 7.0平臺(tái)上，通過相應(yīng)的工具箱函數(shù),對所構(gòu)建的仿真信號進(jìn)行處理：

(1)將仿真信號y(n)通過MATLAB工具箱函數(shù)進(jìn)行EMD多個(gè)層次分解，對比多個(gè)層次分解基線漂移成分的分布情況，確定合理的EMD分解層數(shù)；

(2)對多個(gè)層進(jìn)行分解，分別去除最后一層殘余項(xiàng)分量；

(3)用上述步驟保留下來的IMF分量重構(gòu)ECG信號y(n)，對比多層分解、去除殘余項(xiàng)、重構(gòu)得到的新信號，比較基線漂移的去除效果。

2實(shí)驗(yàn)結(jié)果與分析

2.1EMD分解實(shí)驗(yàn)與分析

在Matlab 7.0平臺(tái)上，對以上所構(gòu)建的仿真心電信號，分別進(jìn)行5、6、7、8層EMD模態(tài)分解，得到圖3所示的結(jié)果。分析可知，1、2、3、4層主要包含了工頻噪聲、肌電噪聲和QRS波信號；伴隨著經(jīng)驗(yàn)?zāi)B(tài)分解層數(shù)的增加，肌電干擾噪聲能量趨于分散，分布于各分解層；基線漂移噪聲伴隨著EMD分解層的增加，能量在最后一層(即EMD分解所得殘余項(xiàng))的分布更加集中；P、T波信號及低頻信號主要分布在其他層。此外，在執(zhí)行EMD分解算法的過程中，分解層數(shù)越多，運(yùn)行時(shí)間越長。總之，從分析結(jié)果來看，分解層數(shù)越多，其他類型的噪聲及有用的PQRS心電信號能量分布越散亂，基線漂移越多集中于EMD分解殘余項(xiàng)。綜合考慮，進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)的7層分解后，將基線漂移所在殘余項(xiàng)去除。這樣既可有效去除基線漂移，又不影響后續(xù)濾波處理，且算法效率適中。

圖3　不同分解層次的EMD分解結(jié)果

2.2基線漂移去除實(shí)驗(yàn)與分析

基于以上實(shí)驗(yàn)結(jié)果，設(shè)計(jì)基線漂移去除方法：(1)對含有基線漂移的信號x(n)進(jìn)行7層EMD分解；(2)將7層EMD分解所得的基線漂移集中分布?xì)堄囗?xiàng)直接剔除掉；(3)用保留下來的IMF分量重構(gòu)ECG信號y(n)。按照此方法分別對構(gòu)建的仿真心電信號及MIT-BIH心律失常數(shù)據(jù)庫中受到嚴(yán)重基線漂移干擾的109號、203號心電數(shù)據(jù)段進(jìn)行基線漂移去除，所得結(jié)果分別如圖4、圖5、圖6所示?；贓MD的7層分解，直接去除殘余項(xiàng)分量，能夠?qū)λ鶚?gòu)建的仿真心電信號及MIT-BIH心律失常數(shù)據(jù)庫中受到嚴(yán)重基線漂移干擾的109號、203號心電數(shù)據(jù)段去除基線漂移；同時(shí)，重構(gòu)的心電信號并未失真，基本上保留了原來的概貌特征，有利于后續(xù)的心電信號處理,但肌電噪聲、工頻噪聲仍無法去除。

圖4　仿真信號和濾波后的結(jié)果

圖5　109號數(shù)據(jù)和濾波后的結(jié)果

圖6　203號數(shù)據(jù)和濾波后的結(jié)果

3結(jié)語

心電信號是生物醫(yī)學(xué)工程領(lǐng)域常見的典型非線性、非平穩(wěn)信號，在采集的過程中易受到各種噪聲干擾。因其自身頻帶分布范圍較寬，用傳統(tǒng)的傅里葉變換進(jìn)行噪聲消除效果欠佳，而EMD算法能有效處理非線性、非平穩(wěn)信號。本次研究針對基線漂移去除，以構(gòu)建的仿真心電信號為實(shí)驗(yàn)源數(shù)據(jù)，借助Matlab 工具平臺(tái)，分別進(jìn)行了EMD多層分解實(shí)驗(yàn)。經(jīng)分析發(fā)現(xiàn)基線漂移噪聲主要集中分布于最后一層的殘余分量上，且7層EMD分解的運(yùn)行速度和分布區(qū)分度都較好。分別對所構(gòu)建的仿真心電信號及MIT-BIH心律失常數(shù)據(jù)庫中受到嚴(yán)重基線漂移干擾的109號、203號心電數(shù)據(jù)段進(jìn)行基線漂移去除，結(jié)果有效可行。

參考文獻(xiàn)

[1] 張喜紅.開放式心電圖機(jī)教學(xué)平臺(tái)的設(shè)計(jì)[J].長江大學(xué)學(xué)報(bào)(自然科學(xué)版)，2013(34):43-46.

[2] 丁瑞,李國軍,王青.心電信號基線漂移去除方法研究[J].云南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，36(5):665-660.

[3] 高彩紅.心電信號臨床信息的識別研究[D].鎮(zhèn)江：江蘇大學(xué),2010：16-26.

[4] 張喜紅,楊清志,王玉香.基于HHT的心電信號濾波算法的研究[J].吉林化工學(xué)院學(xué)報(bào)，2014(2)：51-54.

The Method of ECG Baseline Wander Removal Based on EMD

ZHANGXihongWANGYuxiang

(Bozhou Vocational and Technical College, Bozhou Anhui 236800, China)

Abstract:Based on Empirical Mode Decomposition (EMD) algorithm, this paper focues on the research of new methods of removal ECG baseline wander, thereby improving the quality of the ECG signal. Using MATLAB toolbox for MIT-BIH arrhythmia database relatively clean ECG data segment, we firstly artificially superimposed common ECG signal noise signal generator simulation experiment; then the simulated signals were constructed by various levels of decomposition EMD by decomposing the residual term to remove the resulting realize ECG baseline wander removal. Comparing the different layers decompose under the effect of baseline wander removal, we obtained the results: EMD layer is large, the baseline drift which features more residual layer appeared. The MIT-BIH arrhythmia database is disturbed by the baseline drift of the 109 ECG data segment, and the baseline wander is removed. The feasibility of this method is also verified in this paper.

Key words:Empirical Mode Decomposition; baseline wander; ECG signal

收稿日期：2015-09-28

基金項(xiàng)目：2010年度安徽省教育廳省級質(zhì)量工程項(xiàng)目“醫(yī)用電子儀器與維護(hù)特色專業(yè)建設(shè)”(20101459)；2011年度安徽省教育廳質(zhì)量工程項(xiàng)目“醫(yī)用儀器示范實(shí)習(xí)實(shí)訓(xùn)中心建設(shè)”(2011131)

作者簡介：張喜紅(1983 — )，男，山西朔州人，碩士，亳州職業(yè)技術(shù)學(xué)院講師，技師，研究方向?yàn)樯镝t(yī)學(xué)工程。

中圖分類號：TP391

文獻(xiàn)標(biāo)識碼：A

文章編號：1673-1980(2016)02-0098-03