牛曉東 盧莉蓉 王鑒 韓星程 郭樹言 王黎明?

1) (中北大學(xué),信息探測(cè)與處理山西省重點(diǎn)實(shí)驗(yàn)室,太原 030051)

2) (長(zhǎng)治醫(yī)學(xué)院物理教研室,長(zhǎng)治 046000)

3) (長(zhǎng)治醫(yī)學(xué)院生物醫(yī)學(xué)工程系,長(zhǎng)治 046000)

心電圖(electrocardiogram,ECG)診斷心臟疾病的嚴(yán)格標(biāo)準(zhǔn),要求有效地消除噪聲并準(zhǔn)確地重建ECG信號(hào).經(jīng)驗(yàn)?zāi)Ｊ椒纸?empirical mode decomposition,EMD)方法重建ECG 信號(hào)中,模式混疊及重建采用模式分量的識(shí)別以經(jīng)驗(yàn)為基礎(chǔ),導(dǎo)致重建ECG 信號(hào)準(zhǔn)確度降低,且方法不具有自適應(yīng)和通用性.本文首先基于積分均值定理提出一種改進(jìn)的EMD 方法—積分均值模式分解(integral mean mode decomposition,IMMD)方法,經(jīng)5000 個(gè)高斯白噪聲樣本的蒙特卡羅法驗(yàn)證,IMMD 方法比EMD 具有更優(yōu)多分辨率分析能力,能夠有效地緩解模式混疊.其次,基于ECG 信號(hào)內(nèi)固有心動(dòng)物理特征量識(shí)別重建ECG 信號(hào)所采用的模式分量,具有現(xiàn)實(shí)物理意義,因此,方法具有自適應(yīng)和通用性.經(jīng)驗(yàn)證,提出方法重建47 例ECG 信號(hào)與原ECG 信號(hào)的相關(guān)系數(shù)中: 31 例優(yōu)于變分模式分解方法; 33 例優(yōu)于Haar 小波軟閾值法; 42 例優(yōu)于集總經(jīng)驗(yàn)?zāi)Ｊ椒纸夥椒? 45 例優(yōu)于EMD 方法.相關(guān)系數(shù)均值為0.8904,方差為0.0071,表現(xiàn)穩(wěn)定且最優(yōu).

1 引 言

心電圖(electrocardiogram,ECG)記錄的是心臟生理電活動(dòng)信號(hào),ECG 檢查一般通過ECG 中P,QRS,T 特征波波形、時(shí)長(zhǎng)及幅值識(shí)別診斷各種心臟疾病.例如P 波形狀高尖,且振幅 ≥ 0.25 mV判斷為右房肥大; P,R 間期延長(zhǎng) > 0.21 s 存在房室傳導(dǎo)阻滯等[1].因此,ECG 信號(hào)輕微失真,也可能導(dǎo)致心臟疾病誤診.然而,ECG 信號(hào)是一種非線性、非平穩(wěn)弱信號(hào),極易被噪聲污染.噪聲一般包括兩類: 50 Hz (在一些國(guó)家為60 Hz)電源線工頻及由呼吸、電極或人體運(yùn)動(dòng)引起的基線漂移(baseline wander,BW)等窄帶噪聲; 主要由肌肉收縮引起的肌電干擾(muscle artifact,MA)等寬帶噪聲[2].兩種噪聲常常疊加引起ECG 信號(hào)畸變,嚴(yán)重影響ECG 診斷評(píng)估心臟病準(zhǔn)確性,甚至導(dǎo)致診斷失效.嚴(yán)格的ECG 檢查標(biāo)準(zhǔn)量以及ECG 信號(hào)的噪聲特點(diǎn),要求有效消除噪聲后重建的ECG 信號(hào)具有高準(zhǔn)確度.

窄帶噪聲頻譜與ECG 信號(hào)頻譜不重疊,常見線性和非線性信號(hào)處理方法都可有效校正[3?5],而寬帶噪聲全頻域污染ECG 信號(hào),導(dǎo)致嚴(yán)重特征波失真,成為ECG 信號(hào)消噪重點(diǎn).迄今為止,已有多種非線性信號(hào)處理方法應(yīng)用于消除寬帶噪聲并重建ECG 信號(hào),但各種方法均存在缺點(diǎn),重建準(zhǔn)確度不令人滿意.其中,自適應(yīng)濾波需要額外地采集噪聲相關(guān)參考信號(hào)[6]; 奇異值分解(singular value decomposition,SVD)[7]和小波閾值法[8,9]中,閾值過小消噪不理想,過大將導(dǎo)致波形嚴(yán)重畸變,并且閾值不能自適應(yīng)噪聲水平的波動(dòng)[10].

經(jīng)驗(yàn)?zāi)Ｊ椒纸?empirical mode decomposition,EMD)是一種自適應(yīng)且數(shù)據(jù)驅(qū)動(dòng)的非線性信號(hào)處理方法[11],可以克服SVD 和小波閾值法缺乏自適應(yīng)的缺點(diǎn),成為近年來生理信號(hào)[12]、醫(yī)學(xué)圖像[13]、特別是ECG 信號(hào)重建的研究熱點(diǎn).然而,EMD 域內(nèi)重建ECG 信號(hào)存在兩個(gè)問題: 1)瞬時(shí)模態(tài)分量(intrinsic mode function,IMF)之間的模式混疊,導(dǎo)致ECG 信號(hào)重建不準(zhǔn)確.集總經(jīng)驗(yàn)?zāi)Ｊ椒纸?ensemble empirical mode decomposition,EEMD)[14]以及變分模式分解(variational mode decomposition,VMD)[15]能夠緩解模式混疊,并已應(yīng)用于ECG 信號(hào)消噪和重建[16,17].然而,嚴(yán)格意義上,對(duì)每一個(gè)ECG 信號(hào),EEMD 都需要預(yù)先多次測(cè)試并選擇合適的輔助噪聲幅度和數(shù)量;VMD 同樣需要預(yù)先多次測(cè)試并選擇合適的分解層數(shù)、平衡參數(shù)、噪聲松弛等參數(shù),因此,兩種方法都不屬于自適應(yīng)且不具有通用性.2) 重建ECG 信號(hào)采用的模式分量基于工作者經(jīng)驗(yàn)識(shí)別[18,19],不具有通用性.采用IMF 的香農(nóng)熵、平均值、方差等統(tǒng)計(jì)特征量的閾值識(shí)別IMF,閾值也基于經(jīng)驗(yàn)設(shè)置[20].

極值域模式分解(extremum field mean mode decomposition,EMMD)[21]是一種改進(jìn)的EMD 方法.EMMD 采 用 積 分 均 值 定 理(integral mean value theorem,IMVT)求信號(hào)均值,能夠增強(qiáng)分解信號(hào)能力,但EMMD 算法分解信號(hào)存在過分解問題.為了在EMD 域內(nèi)進(jìn)一步提高重建ECG 信號(hào)準(zhǔn)確度,同時(shí)保證重建方法具有自適應(yīng)和通用性,本文提出一種重建ECG 信號(hào)新方法.新方法首先改進(jìn)EMMD,保留IMVT 均值方法并解決算法導(dǎo)致的過分解問題.為了區(qū)別于EMMD,本文將改進(jìn)的EMMD 稱為積分均值模式分解(integral mean mode decomposition,IMMD).經(jīng)驗(yàn)證可以發(fā)現(xiàn),IMMD 比EMD 具有更優(yōu)的多分辨率分解能力,并且可以有效地緩解IMF 之間的模式混疊.其次,心電活動(dòng)或心肌細(xì)胞的除極、復(fù)極過程對(duì)應(yīng)一個(gè)完整的心動(dòng)周期,因此,ECG 信號(hào)呈現(xiàn)心動(dòng)的周期性波動(dòng),其IMF 分量同樣表現(xiàn)出心動(dòng)物理周期特征.本文采用心動(dòng)物理周期特征識(shí)別屬于ECG信號(hào)的IMF,識(shí)別方法符合ECG 信號(hào)的物理本質(zhì),并保證重建ECG 信號(hào)方法的自適應(yīng)及通用性.

2 方法理論

2.1 IMMD 方法

設(shè)xi(t)(ti≤t≤ti+1,i=1,2,··· ,n ?1) 為信號(hào)x(t) 相鄰兩個(gè)極值x(ti)和x(ti+1) 之間的局部,IMVT 求信號(hào)局部xi(t) 的均值為

mi(t)固定于局部中點(diǎn).基于(1)式求出所有局部均值點(diǎn)后,三次樣條插值所有均值點(diǎn)構(gòu)建均值曲線m(t).

原型模式函數(shù)(proto mode function,PMF)[22]即為

上述過程稱為一次均值篩選[11],PMF 記為PMF1.PMF1作為一個(gè)新信號(hào),迭代重復(fù)均值篩選過程k(k= 2,3,···)次,得到PMFk.當(dāng)PMFk滿足柯西篩選停止準(zhǔn)則[11]:

篩選停止,PMFk即為IMF1.將剩余分量r1(t)=x(t) – IMF1 作為新信號(hào),重復(fù)上述所有過程得到IMF2.同理,可得其余IMF 分量及剩余分量r(t).

除局部均值外,EMMD 方法利用相鄰兩個(gè)局部均值mi(t)和mi+1(t)加權(quán)求取二者之間極值點(diǎn)ti+1處均值,即:

其中,

然而,m(ti+1)攜帶數(shù)據(jù)信息是mi(t)和mi+1(t)所攜信息的重復(fù).并且均值點(diǎn)數(shù)增加,極易導(dǎo)致樣條插值構(gòu)造均值曲線的不平滑,最終造成信號(hào)過分解.因此,IMMD 方法僅采用局部均值經(jīng)三次樣條插值構(gòu)建均值曲線.

2.2 IMMD 方法的多分辨率分解信號(hào)能力

文獻(xiàn)[23]研究了EMD 分解5000 個(gè)長(zhǎng)度N=512,均值為0,方差為1 的高斯白噪聲數(shù)據(jù)樣本,得到如下結(jié)論: EMD 類似于小波變換Mallat 算法,等效于一個(gè)恒定品質(zhì)因子Q的二分(或二進(jìn))濾波器組,具有多分辨率分解信號(hào)能力.本節(jié)將采用同樣方法,分析IMMD 分解高斯白噪聲特性,以此說明方法分解信號(hào)的性能.

IMMD 分解每個(gè)高斯白噪聲樣本,得到至少10 個(gè)IMF.對(duì)所有樣本相應(yīng)階數(shù)IMF 求平均傅里葉功率譜,結(jié)果見圖1(a).IMF1 相當(dāng)于一個(gè)高通濾波器,其余IMF 等效于一組重疊帶通濾波器,且后一IMF 等效帶通濾波器中心頻率約是前一 IMF等效帶通濾波器中心頻率的2/3.IMMD 方法分解高斯白噪聲等效于濾波器組,且具有三分特性.

除IMF1 外,后一IMF 等效帶通濾波器頻帶寬度大致是前一IMF 的2/3,所以等效帶通濾波器組具有恒Q性質(zhì).因此,IMFs 的功率譜存在自相似性[23]:

圖1 IMMD 和EMD 方法分解高斯白噪聲的等效濾波器組特性 (a) IMMD (實(shí)線)和EMD (虛線)的IMF 分量的平 均 功 率 譜; (b) 基 于(5)式,IMMD (實(shí) 線)和EMD (虛線)的IMF 分量的平均功率譜坍縮重合Fig.1.Equivalent filter banks of IMMD and EMD decomposing Gauss white noises: (a) Averaged power spectra of IMFs from IMMD (solid curves) and EMD (dotted curves);(b) collapse and coincidence of the average power spectrum of IMFs from IMMD (solid curves) and EMD (dotted curves) based on Eq.(5).

其中Sk′(f) 是第k′個(gè)IMF 的 功率 譜(Sk(f) 是 第k個(gè)IMF 的功率譜),k′ >k≥2 ;ρ是常數(shù),對(duì)于IMMD,ρ=3/2 ,對(duì)于EMD,ρ=2 .基于(5)式對(duì)平均功率譜標(biāo)準(zhǔn)化,所有帶通IMF 的功率譜坍縮并重合為一條曲線.對(duì)于EMD,帶通濾波器混疊較多低頻成分,關(guān)于中心頻率對(duì)稱度差,窄帶特性差,需要適當(dāng)調(diào)整ρ值為2.01[24].對(duì)于IMMD,等效帶通濾波器關(guān)于中心頻率對(duì)稱度及窄帶特性明顯優(yōu)于EMD,兩套等效帶通濾波器組的平均功率譜標(biāo)準(zhǔn)化結(jié)果如圖1(b)所示.

IMMD 分解高斯白噪聲等效濾波器組具有三分特性,等效濾波器個(gè)數(shù)多于EMD (圖1),濾波器組窄帶特性及恒Q性質(zhì)優(yōu)于EMD,因此,具有更好的多分辨率分析能力.除上述5000 個(gè)高斯白噪聲樣本數(shù)據(jù)的蒙特卡羅方法實(shí)驗(yàn)外,本課題組基于大量數(shù)據(jù)實(shí)驗(yàn),驗(yàn)證得到: 相比EMD,IMMD 具有更好的多分辨率分析能力,可以有效地緩解IMF分量之間模式混疊.

2.3 心動(dòng)物理特征識(shí)別IMF

ECG 信號(hào)是記錄人體心臟電活動(dòng)信號(hào),具有心動(dòng)物理周期特性,其IMF 分量將呈現(xiàn)心動(dòng)周期或心率(heart rate,HR)特征.

任一個(gè)IMF 分量可表達(dá)為[25]

其中φ(t) 為IMF 頻率調(diào)制,a(t) 為IMF 振幅調(diào)制.心動(dòng)周期特征通過兩種模式作用于IMF 分量:1) 屬于ECG 信號(hào)的低階IMF,心動(dòng)周期特征作用于調(diào)幅a(t) .此時(shí),IMF 分量的包絡(luò)具有心動(dòng)周期特性,包絡(luò)的頻譜中最大幅值對(duì)應(yīng)頻率等于HR.2) 屬于ECG 信號(hào)的高階IMF,心動(dòng)周期特征作用于調(diào)頻φ(t) .此時(shí),IMF 具有周期性心動(dòng)的諧波特性,頻譜中最大幅值對(duì)應(yīng)頻率為HR 的整數(shù)倍(通常為1—3 倍).噪聲屬性IMF 分量不能表現(xiàn)上述兩種心動(dòng)周期特性.

2.4 IMMD 域內(nèi)心動(dòng)物理特征識(shí)別模式分量重建ECG 方法

1) 含噪ECG 信號(hào)經(jīng)IMMD 分解為一套IMFs分量集合;

2) 求出每一個(gè)IMF 的上包絡(luò)功率譜,若譜中最大幅值對(duì)應(yīng)頻率等于HR,則IMF 被識(shí)別為ECG信號(hào)分量;

3) 對(duì)步驟2 中未被識(shí)別為ECG 信號(hào)分量的IMF,求出每一個(gè)IMF 功率譜,若譜中最大幅值對(duì)應(yīng)頻率等于HR 整數(shù)倍,則IMF 也被識(shí)別為ECG信號(hào)分量;

4) 采用上述所有被識(shí)別為ECG 信號(hào)分量的IMF重建ECG 信號(hào).

3 實(shí)驗(yàn)驗(yàn)證

3.1 105 ECG 信號(hào)[26]實(shí)驗(yàn)

選取MIT-BIH 心律失常數(shù)據(jù)庫中105 ECG、壓力測(cè)試數(shù)據(jù)庫中BW 和MA 噪聲[27],疊加構(gòu)成含噪105 ECG 信號(hào),采樣頻率為360 Hz,采樣長(zhǎng)度為3.6×103,HR 等于1.4 Hz,如圖2 所示.可見,原105 ECG 信號(hào)的噪聲污染程度較低,便于本節(jié)下述實(shí)驗(yàn)中的定量分析.

圖3 含噪105 ECG 信號(hào)的IMF (藍(lán)色曲線)及其包絡(luò)(紅色曲線)Fig.3.IMF envelopes (red curves) and IMFs (blue curves) of noisy No.105 ECG signal.

圖2 實(shí)驗(yàn)采用信號(hào)波形 (a) 105 ECG 信號(hào); (b) BW 信號(hào); (c) MA 信號(hào); (d)含噪105 ECG 信號(hào)Fig.2.Waveform of signals used in experiment: (a) No.105 ECG signal; (b) BW signal; (c) MA signal; (d) the noisy No.105 ECG signal.

圖4 含噪105 ECG 信號(hào)IMFs 的功率譜(藍(lán)色曲線),及IMFs 包絡(luò)的功率譜(紅色曲線)Fig.4.Power spectra of envelopes of IMFs (red curves) and of IMFs (blue curves) of noisy No.105 ECG signal.

IMMD 分解含噪105 ECG 信號(hào)得到IMF 分量及其包絡(luò)如圖3 所示,對(duì)應(yīng)頻譜如圖4 所示.IMF1和IMF2 及其包絡(luò)的頻域中,幅值最高點(diǎn)處頻率不等于1.4 Hz,因此沒有表現(xiàn)出心動(dòng)周期特征作用IMF 的兩種模式(詳見2.3 節(jié)),被識(shí)別為MA 噪聲高頻分量; IMF3—8 包絡(luò)的頻域中,幅值最高點(diǎn)對(duì)應(yīng)頻率等于心率1.4 Hz,符合心動(dòng)周期特征作用IMF 的第一種模式; IMF9—12 的頻域中,幅值最高點(diǎn)對(duì)應(yīng)頻率分別等于心率的3,3,2,1 倍,識(shí)別為ECG 信號(hào)諧波分量,符合心動(dòng)周期特征作用IMF 的第二種模式; IMF13—r 沒有表現(xiàn)出心動(dòng)周期特征作用IMF 的兩種模式,統(tǒng)一歸類為低頻噪聲.產(chǎn)生的原因包括: MA 噪聲低頻分量,BW 噪聲分量,IMMD 方法端點(diǎn)效應(yīng)以及原105 ECG 信號(hào)的非零均值.經(jīng)上述心動(dòng)物理周期特征識(shí)別,IMF3—12 為ECG 信號(hào)分量,IMF1 和IMF2 為高頻噪聲,IMF13—r 為低頻噪聲.

為了顯示所提出方法的能力,將本文方法同近年來常用VMD、小波軟閾值法、EEMD 以及EMD重建105 ECG 信號(hào)進(jìn)行對(duì)比,利用這五種方法重建的105 ECG 信號(hào)如圖5 所示.由圖5 可見,五種方法都很好地消除了BW 噪聲.由于MA 噪聲寬頻特性,五種方法重建ECG 信號(hào)中仍然存在少量MA 噪聲分量,特別是Harr 小波軟閾值方法.EMD方法重建ECG 信號(hào)畸變最為嚴(yán)重,其次是EEMD.實(shí)際中,EEMD 對(duì)含噪105 ECG 信號(hào)的多次重建結(jié)果之間都有輕微不同,這是由輔助白噪聲的隨機(jī)性引起的[28].另外,VMD 方法中特征R 波的峰值失真比本文提出方法嚴(yán)重,例如,圖5(b)中第10個(gè)R 波波峰峰值損失14.5%,本文方法為3.2%.

采用重建ECG 與原ECG 信號(hào)的相關(guān)系數(shù)R定量描述重建準(zhǔn)確度,信噪比(signal-to-noise ratio,SNR)及均方誤差(mean square error,MSE)定量描述消噪能力,五種方法重建105 ECG 信號(hào)R,SNR,MSE 的值如表1 所列.可 見,提出方法的R 和SNR 值最大,MSE 值最小,重建105 ECG 信號(hào)和消噪能力最優(yōu).

3.2 47 例含噪ECG 信號(hào)實(shí)驗(yàn)

為了進(jìn)一步驗(yàn)證提出方法的有效性以及普遍適用性,選取MIT-BIH 心律失常數(shù)據(jù)庫中其余ECG 信號(hào)(心律失常ECG 信號(hào)對(duì)應(yīng)類型見表2),分別疊加BW 和MA 噪聲,構(gòu)成47 例含噪ECG信號(hào)(除去周期性極差的232 號(hào)),并采用上述五種方法重建47 例原ECG 信號(hào).

圖5 原105 ECG 信號(hào)(藍(lán)色點(diǎn)虛線)與由五種方法重建的105 ECG 信號(hào)(紅色實(shí)線) (a) 本文方法; (b) VMD; (c) Haar 小波軟閾值; (d) EMMD; (e) EMDFig.5.Original No.105 ECG signal (blue dotted curves) and the No.105 ECG signals (red solid curves) reconstructed by 5 methods: (a) The proposed method; (b) VMD; (c) Haar wavelet with soft threshold; (d) EEMD; (e) EMD.

表1 五種方法重建105 ECG 信號(hào)的特征量值Table 1.Characteristic values of No.105 ECG signals reconstructed by 5 methods.

表2 實(shí)驗(yàn)采用的部分ECG 信號(hào)對(duì)應(yīng)心律失常類型Table 2.Type of arrhythmia corresponding to some ECG signals used in the experiment.

由于大多數(shù)原ECG 信號(hào)本身含有一定量的低頻和高頻噪聲,嚴(yán)重影響SNR 和MSE 值的準(zhǔn)確度,因此,僅采用相關(guān)系數(shù)R對(duì)五種方法進(jìn)行對(duì)比評(píng)估.利用五種方法重建5 組47 例ECG 信號(hào)與原ECG 信號(hào)的相關(guān)系數(shù)柱狀圖如圖6 所示.其中,本文提出方法有31 例R值優(yōu)于VMD,33 例R值優(yōu)于Haar 小波,42 例R值優(yōu)于EEMD,45 例R值優(yōu)于EMD.

采用五種方法重建47 個(gè)ECG 信號(hào)的相關(guān)系數(shù)R值的統(tǒng)計(jì)盒形圖如圖7 所示,平均值及方差見表3.表3 中,本文提出方法對(duì)應(yīng)相關(guān)系數(shù)均值最高,方差僅大于小波閾值.由圖7 可見,本文提出方法明顯優(yōu)于其余四種方法.本文提出方法重建ECG 信號(hào)最為穩(wěn)定、準(zhǔn)確,其次為VMD,Haar小波閾值法,EEMD,EMD.

表3 五種方法重建47 個(gè)ECG 信號(hào)的R 值的均值與方差Table 3.Means and variances of R values of 47 ECG signals reconstructed by 5 methods.

4 討論與結(jié)論

相比經(jīng)驗(yàn)性包絡(luò)相減法求均值(EMD,EEMD方法內(nèi)),本文提出的IMMD 方法采用IMVT 求均值,具有嚴(yán)格的均值數(shù)學(xué)基礎(chǔ),因此,IMF 能夠從信號(hào)中被準(zhǔn)確地篩選剝離.經(jīng)5000 個(gè)高斯白噪聲樣本的蒙特卡羅法統(tǒng)計(jì)驗(yàn)證,IMMD 方法具有比EMD,EEMD (同EMD 方 法 一 樣,EEMD 具 有恒Q二分濾波器組特性[29])更優(yōu)的多分辨率分解信號(hào)能力,可以有效緩解IMF 之間的模式混疊.

圖6 采用五種方法重建47 個(gè)ECG 信號(hào)的R 值柱狀圖Fig.6.Histogram of R values of 47 ECG signals reconstructed by 5 methods.

圖7 采用五種方法重建47個(gè)ECG信號(hào)的相關(guān)系數(shù)R 值的統(tǒng)計(jì)盒子圖Fig.7.Five box-plots of R values of 47 ECG signals reconstructed by 5 methods.

另外,ECG 信號(hào)具有心動(dòng)周期、心率等波動(dòng)物理特征,ECG 信號(hào)的IMF 分量同樣具有心動(dòng)周期或HR 的特性.因此,本文方法中使用心動(dòng)周期或HR 識(shí)別ECG 信號(hào)的IMF 分量,符合ECG 的物理本質(zhì)特性.EMD,EEMD 和VMD 域內(nèi)ECG 重建一般采用QRS 特征波經(jīng)驗(yàn)性識(shí)別IMF,適用于時(shí)域中具有顯著QRS 波動(dòng)的IMF 識(shí)別(例如圖3中IMF4).如果存在: 1) 時(shí)域中QRS 波動(dòng)不顯著但屬于ECG 信號(hào)的低階IMF (例如圖3 中IMF3);2) 時(shí)域完全沒有QRS 特征但屬于ECG 信號(hào)的高階IMF (例如圖3 中IMF12),該方法識(shí)別錯(cuò)誤.小波閾值法中固定閾值不能自適應(yīng)小波系數(shù)不同局部?jī)?nèi)噪聲水平波動(dòng),因此,重建ECG 信號(hào)可能存在局部消噪不理想(例如圖5(c)).所以,本文提出方法,能夠比上述四種方法進(jìn)一步提高ECG信號(hào)重建準(zhǔn)確度,且方法具有自適應(yīng)和通用性.經(jīng)47 例ECG 信號(hào)實(shí)驗(yàn)驗(yàn)證,提出的方法重建ECG與原ECG 信號(hào)相似度平均值為0.8904,方差為0.0071,且有31 例相似度值優(yōu)于VMD,33 例相似度值優(yōu)于Haar 小波,42 例相似度值優(yōu)于EEMD,45 例相似度值優(yōu)于EMD.本文提出方法重建ECG 信號(hào)能力體現(xiàn)了物理特征或現(xiàn)象在生物電信號(hào)處理中具有重要作用.

本方法中,建議通過含噪ECG 信號(hào)IMF 分量包絡(luò)頻譜圖,經(jīng)驗(yàn)確定HR.HR 可以通過嚴(yán)格的RR 波間隔得到,但需要大量算法實(shí)現(xiàn).實(shí)際上,人心率一般為1—1.7 Hz (60—100 次/min).對(duì)于采集良好的ECG 信號(hào),其IMF 分量包絡(luò)的頻譜圖中,1—2 Hz 之間的幅度最大值對(duì)應(yīng)頻率即為HR,并且它常常出現(xiàn)在多個(gè)IMF 分量包絡(luò)的頻譜中,容易辨別(如圖4 所示).本文選取的47 個(gè)ECG信號(hào),都是基于上述經(jīng)驗(yàn)方法確定HR.

提出的方法具有一定局限.對(duì)于極其特殊的232 ECG 信號(hào)(竇性心動(dòng)過緩、一級(jí)房室傳導(dǎo)阻滯和頻繁異位心房運(yùn)動(dòng),ECG 信號(hào)停頓持續(xù)長(zhǎng)達(dá)6 s)重建,本文提出方法失效.關(guān)于心動(dòng)周期特性極差的ECG 信號(hào)重建,以及從物理現(xiàn)象本質(zhì)探索ECG信號(hào)處理,是本課題組下一步的工作.