基于誤差可控的轉(zhuǎn)換變換心電信號去噪方法

摘" 要： 針對心電信號在去除噪聲時不可控誤差引起的幅值和形態(tài)變形問題，提出一種基于誤差可控的轉(zhuǎn)換變換心電信號去噪方法。通過轉(zhuǎn)換變換可以獲取心電信號的低頻分量和高頻分量：通過對特征波形和高頻噪聲處設(shè)置不同的誤差界可有效濾除高頻噪聲，并將特征波的重構(gòu)誤差控制在較小的范圍內(nèi)；對低頻分量進行cubic插值可以獲得更加光滑的基線漂移。從而可改善去噪引起的形態(tài)改變。在MIT?BIH數(shù)據(jù)庫的心電圖上進行實驗，結(jié)果表明：提出的方法可以有效濾除基線漂移和高頻噪聲，同時特征波形的幅值和形狀得到了很好的保持。與對比方法相比，去噪后心電信號的信噪比、均方根誤差、歸一化相關(guān)系數(shù)均得到改善。

關(guān)鍵詞： 心電信號去噪； 轉(zhuǎn)換變換； 基線漂移； 高頻噪聲； 信噪比； 均方根誤差； 歸一化相關(guān)系數(shù)

中圖分類號： TN911.7?34" " " " " " " " " " " " 文獻標(biāo)識碼： A" " " " " " " " " " " " "文章編號： 1004?373X（2024）17?0058?07

ECG signal denoising method based on shift transformation with controllable error

LI Xiaoyun1， 2， LI Tongliang3， ZHAO Huanyu1， 2， HUANG Shizhong1， 2

（1. Institute of Applied Mathematics， Hebei Academy of Sciences， Shijiazhuang 050081， China;

2. Information Security Authentication Technology Innovation Center of Hebei Province， Shijiazhuang 050081， China;

3. Institute of Biology， Hebei Academy of Sciences， Shijiazhuang 050081， China）

Abstract： In view of the amplitude and shape deformation caused by uncontrollable error in removing noise from electrocardiogram （ECG） signals， an ECG signal denoising method based on shift transformation with controllable error is presented. The low?frequency and high?frequency components of the ECG signals can be obtained by the shift transformation. By setting different error bounds on the characteristic waveform and high?frequency noise， the high?frequency noise can be filtered out effectively and the reconstruction error at the characteristic waveform can be kept within a small range. Smoother baseline drift can be obtained by implementing cubic interpolation for low?frequency components. On the basis of the above， the morphological changes caused by denoising can be improved. Experiments on the ECG of the database MIT?BIH show that the proposed method can remove baseline drift and high?frequency noise effectively， while maintaining the amplitude and shape of the characteristic waveform well. In comparison with the contrast methods， the signal?to?noise ratio （SNR）， the root mean square error （RMSE） and the normalized correlation coefficient （NCC） of the ECG signals after denoising reach satisfactory values.

Keywords： ECG signal denoising; shift transformation; baseline drift; high?frequency noise; SNR; RMSE; NCC

0" 引" 言

心電信號在采集過程中易受到各種噪聲干擾，這些噪聲可能會淹沒心電信號的有用特征，如特征波的幅值和形狀，這些特征通常與心律失常、心房病變等疾病有關(guān)。因此如何有效濾除各種噪聲，同時更多保留原始心電信號特征形態(tài)是醫(yī)學(xué)信息處理中需要解決的難題。

目前有各種去噪方法，如基于數(shù)字濾波[1?2]、中值濾波[3]、形態(tài)學(xué)處理[4]、神經(jīng)網(wǎng)絡(luò)[5?6]、小波變換[7?10]和經(jīng)驗?zāi)B(tài)分解[11?13]等方法。這些方法各有其特點，但也存在一些局限性。數(shù)字濾波器和中值濾波器在去除噪聲時，可能由于參數(shù)選取不當(dāng)導(dǎo)致信號嚴(yán)重失真；形態(tài)學(xué)濾波在處理基線漂移時可能會濾除有用信息；基于神經(jīng)網(wǎng)絡(luò)的去噪效果則依賴于訓(xùn)練數(shù)據(jù)集中噪聲的分布特征；小波變換雖然能夠去除高頻噪聲，但是在閾值處理時可能會導(dǎo)致特征波形誤差不可控，進而影響信號的幅值和形狀；與小波方法相比，經(jīng)驗?zāi)B(tài)分解方法可以自適應(yīng)地分解信號，但它仍然依賴于小波閾值去噪，無法避免不可控誤差導(dǎo)致的波形變形。

為了解決上述問題，本文提出一種基于誤差可控的轉(zhuǎn)換變換[14?15]心電信號去噪方法。通過在轉(zhuǎn)換變換中對信號的有用部分和噪聲部分設(shè)置不同的誤差，以保持心電信號的原始特征形態(tài)。此外，本文還使用cubic插值[16]算法來修正基線漂移，從而減小信號的變形。這種方法旨在去除噪聲的同時，更大限度地保留心電信號的生理特征，以便更準(zhǔn)確地進行后期的心電信號的分析和診斷。

1" 轉(zhuǎn)換變換

轉(zhuǎn)換變換是一種以最大誤差為準(zhǔn)則的數(shù)據(jù)壓縮算法，優(yōu)勢在于不僅可保證整體誤差可控，還能將每一個數(shù)據(jù)重構(gòu)誤差控制在給定的誤差界內(nèi)。數(shù)據(jù)壓縮和去噪的共同點是兩者都需要減少冗余信息并盡可能多地保留重要信息，而心電信號中的一些噪聲信號可視作冗余信息，因此本文提出了一種利用壓縮算法實現(xiàn)去噪的策略。

圖1是一個轉(zhuǎn)換變換的誤差樹[T]，原始數(shù)據(jù)集[D={d0，d1，d2，…，d7}={8，12，9，1，7，3，6，10}]，給定誤差界[Δ=3]。通過轉(zhuǎn)換變換可對[D]構(gòu)建一個小于原始數(shù)據(jù)量的小波概要[S={s0，s2，s3，s5}={7，2.5，]-[1.5，4}]，利用[S]進行重構(gòu)，可得重構(gòu)數(shù)據(jù)集[D={d0，d1，d2，…，d7}={9.5，9.5，8.5，0.5，5.5，5.5，8.5，8.5}]。由此可見，轉(zhuǎn)換變換對數(shù)據(jù)進行壓縮時可將每一數(shù)據(jù)點的重構(gòu)誤差控制在給定誤差界[Δ=3]內(nèi)。

轉(zhuǎn)換變換是一種多尺度的信號分析方法，可以將信號分解成不同尺度的高頻分量和低頻分量。轉(zhuǎn)換變換相當(dāng)于在Haar小波變換基礎(chǔ)上引入了一個誤差界，這使得每一尺度的低頻分量不再限定為一個數(shù)值，而是一個數(shù)據(jù)范圍，這也為下一尺度高頻分量轉(zhuǎn)變?yōu)?提供了更多可能，從而實現(xiàn)了冗余信息的去除。具體來說，轉(zhuǎn)換變換需要從誤差樹的葉子節(jié)點開始自下而上計算每一節(jié)點的轉(zhuǎn)換系數(shù)[sj]（[j=0，1，2，…，n-1]，[n]為原始數(shù)據(jù)量）。葉子節(jié)點[di]（[i=0，1，2，…，n-1]）需采用區(qū)間范圍[[di，di]]來替代原始數(shù)據(jù)。重構(gòu)后需保證每一點誤差[di-di≤Δ]，其中[di]是數(shù)據(jù)[di]的重構(gòu)數(shù)值，因此[[di，di]]可以表示為[[di-Δ，di+Δ]]，如圖1所示，[d0]可以表示為[5，11]。接下來，需要計算內(nèi)部節(jié)點處的轉(zhuǎn)換系數(shù)[sj]（[j=1，2，…，n-1]），即高頻分量。[dL]和[dR]分別代表內(nèi)部節(jié)點的左右葉子節(jié)點。[dL]和[dR]的區(qū)間范圍可以表示為[[dL，dL]]和[[dR，dR]]。內(nèi)部節(jié)點處的轉(zhuǎn)換系數(shù)[sj]可通過如下步驟得到[15]：

1） 當(dāng)[[dL，dL]?[dR，dR]=?]（[?]代表空集）時，表明兩個葉子節(jié)點的距離大于給定誤差界，此時[sj]的計算公式如下：

[sj=（dL+dL）-（dR+dR）4] （1）

然后更新對應(yīng)節(jié)點處的區(qū)間范圍[[dj，dj]]，即低頻分量的取值區(qū)間：

[dj=max{dL-sj，dR+sj}] （2）

[dj=min{dL-sj，dR+sj}] （3）

2） 當(dāng)[[dL，dL]?[dR，dR]≠?]時，表明兩個葉子節(jié)點的距離在給定誤差界內(nèi)。此時轉(zhuǎn)換系數(shù)[sj]置0，從而實現(xiàn)壓縮，從去噪角度可看作對不重要信息的去除。此時[[dj，dj]]可表示為：

[[dj，dj]=[dL，dL]?[dR，dR]] （4）

以上步驟稱作一步轉(zhuǎn)換變換，迭代地將每一級節(jié)點的低頻分量進行一步轉(zhuǎn)換變換，直到最后的根節(jié)點。根節(jié)點處的轉(zhuǎn)換系數(shù)[s0]通常表示為[s1]對應(yīng)的內(nèi)部節(jié)點區(qū)間范圍的均值，這個過程就是轉(zhuǎn)換變換。數(shù)據(jù)重構(gòu)公式可表示為：

[d（S）i=sj∈path（di）δijsj] （5）

式中：[S]為[sj]的集合；[path（di）]代表誤差樹[T]中葉子節(jié)點[di]的所有祖先節(jié)點。當(dāng)[di∈sj]的左子樹時，[δij=1]；當(dāng)[di∈sj]的右子樹時，[δij=-1]。

為了與小波系數(shù)度量結(jié)果相統(tǒng)一，需要對轉(zhuǎn)換系數(shù)進行歸一化，本文將公式（1）更新為：

[sj=（dL+dL）-（dR+dR）4×2] （6）

將式（2）和式（3）更新為：

[dj=max{2dL-sj，2dR+sj}] （7）

[dj=min{2dL-sj，2dR+sj}] （8）

將重構(gòu)公式中參數(shù)[δij]更新為[±22]。

在轉(zhuǎn)換變換的過程中，心電信號被逐級分解。圖2為8級分解的示意圖，其中[Di]代表低頻分量，[Si]代表高頻分量，[i=1，2，…，8]。

2" 基于轉(zhuǎn)換變換的心電信號去噪算法

2.1" 算法流程圖

圖3是本文提出的基于轉(zhuǎn)換變換的心電信號去噪算法流程圖。該算法以轉(zhuǎn)換變換為基礎(chǔ)，根據(jù)心電信號采樣頻率與基線漂移信號頻率之間的關(guān)系，確定分解級數(shù)[L]；通過對特征波形和高頻噪聲處設(shè)定自適應(yīng)誤差界，可以去除高頻噪聲并保持特征波形不變；對第[L]級低頻分量進行cubic插值可得到基線漂移信號；數(shù)據(jù)重構(gòu)后剔除基線漂移信號，可以得到最終的去噪心電信號。

2.2" 誤差設(shè)置與高頻噪聲去除

圖4是對心電信號進行8級轉(zhuǎn)換分解的一個例子，此時[Δ=0]，表示不去除任何冗余信息。圖4中ECG是含有肌電噪聲的心電信號，[S1]～[S8]分別對應(yīng)每一級高頻分量，[D8]代表第8級低頻分量。從圖中可以看出，在高頻分量[S1]～[S3]中都存在一定程度類似白噪聲的噪聲分量。

通常小波閾值去噪會對含噪的分量設(shè)置閾值，將小于閾值的系數(shù)置0來去除高頻噪聲。而轉(zhuǎn)換變換是對原始數(shù)據(jù)設(shè)定誤差界（閾值），比較相鄰數(shù)據(jù)之差是否在誤差范圍內(nèi)，滿足條件的系數(shù)會直接置0。而肌電引起的高頻噪聲類似于白噪聲，幅值遠低于特征波，并且相鄰數(shù)值波動不大，因此可通過設(shè)置誤差界來去除高頻噪聲。

在轉(zhuǎn)換變換中通常設(shè)置固定的誤差界作用于每一個數(shù)據(jù)點，但是為了濾除心電信號中的高頻噪聲同時保持特征波形不變，需要對噪聲處和特征波形處設(shè)置不同的誤差界。本文誤差界的設(shè)置方法如下：首先將數(shù)據(jù)長度為[n]的心電信號分成不重疊的小段，每一小段的長度為step；再分別計算每一小段心電信號的標(biāo)準(zhǔn)差[stdj]，[j=1，2，…，g]，[g=nstep]，[]符號代表向下取整。[stdj]的計算公式如下：

[stdj=1stepi=1stepdi-μ2] （9）

[μ=1stepi=1stepdi] （10）

式中[μ]為均值。在得到每一段的標(biāo)準(zhǔn)差后，需求得所有小段標(biāo)準(zhǔn)差的平均值（mstd）：

[mstd=1gj=1gstdj] （11）

據(jù)此可以得到一個閾值[thr=α×mstd]，其中[α]為調(diào)節(jié)因子。

根據(jù)閾值和每一段標(biāo)準(zhǔn)差之間的關(guān)系可判斷出該波段是屬于噪聲處還是特征波形處，將特征波形處的誤差界設(shè)定為0以保證重構(gòu)后數(shù)值不變，將噪聲處的誤差界設(shè)定為閾值以去除高頻噪聲。根據(jù)此閾值可設(shè)定每一段的自適應(yīng)誤差界[Δj]：

[Δj=0，stdjgt;thrthr，stdj≤thr] （12）

經(jīng)過大量實驗得出，當(dāng)step=8時通常能夠得到良好的去噪效果，調(diào)節(jié)因子[α]取值通常在1～3之間，當(dāng)[α]取值越大時噪聲去除越徹底，但是有用信號也被去除；當(dāng)[α]取值越小時噪聲去除能力越弱，因此需要根據(jù)待處理信號的噪聲水平相應(yīng)調(diào)整，以達到兼具去噪和保持特征的能力。

圖5是對圖4中含噪的心電信號設(shè)定自適應(yīng)誤差界后的轉(zhuǎn)換分解圖和重構(gòu)的心電圖，可以看到重構(gòu)的心電信號高頻噪聲被去除，相應(yīng)的高頻分量[S1～S3]中的噪聲也被去除。

2.3" 基線校正

通過轉(zhuǎn)換變換后信號被逐級分解，其低頻分量的頻率也逐級減半。通常基線漂移的頻帶處于0.05～1 Hz。因此，可根據(jù)心電信號的采樣頻率和基線漂移的頻率范圍來確定分解級數(shù)[L]。本文心電信號采樣頻率為360 Hz。當(dāng)對360 Hz的心電信號分解到第8級時，其低頻信號的頻帶范圍約為0～0.7 Hz，與基線漂移的頻帶范圍大致相同，因此可確定分解級數(shù)[L]=8。

經(jīng)過轉(zhuǎn)換變換后低頻分量數(shù)據(jù)量小于原始數(shù)據(jù)量，在此采用cubic插值[16]將該級低頻分量進行擴展，并且對插值后的數(shù)據(jù)規(guī)范化可得基線漂移信號，規(guī)范化計算公式如下：

[W=P（2）L] （13）

式中：[P]代表cubic插值后得到的信號；[W]代表對[P]進行規(guī)范化后的基線漂移信號。在獲得基線漂移后，利用重構(gòu)信號減去基線漂移信號即可得到去噪的心電信號。圖6是采用轉(zhuǎn)換變換去噪的過程圖，可以看出高頻噪聲和基線漂移被有效去除。

3" 仿真實驗

3.1" 實驗數(shù)據(jù)及評價指標(biāo)

本文實驗數(shù)據(jù)來源于麻省理工學(xué)院的MIT?BIH心率失常數(shù)據(jù)庫[17]，數(shù)據(jù)采樣頻率為360 Hz。本文對實驗數(shù)據(jù)進行了定性分析和定量分析。定性分析主要從主觀視覺方面進行判斷，定量分析采用SNR、RMSE、NCC來衡量[5，10，12]。SNR主要用于衡量降噪的程度，數(shù)值越大，說明信號中所含的噪聲越少，降噪效果就越好，反之，降噪效果越差。

[SNR=10lgi=1nX（i）2X（i）-Y（i）2] （14）

式中：[X]代表原始“干凈”的心電信號；[Y]代表去噪后的心電信號；[n]表示信號數(shù)據(jù)長度。

RMSE是從整體上來評價兩個信號之間的幅值差異，數(shù)值越小，表明兩個信號之間的幅值差異越小，反之越大。

[RMSE=1ni=1nX（i）-Y（i）2] （15）

NCC主要用于表征兩個特征曲線之間的相似程度，數(shù)值越接近1，表明兩個信號波形的特征越相似，即波形相似度越高，反之，越低。

[NCC=i=1nX（i）Y（i）i=1nX（i）2i=1nY（i）2] （16）

3.2" 去噪實驗比較

為了驗證本文的去噪效果，圖7展示了不同方法下MIT?BIH數(shù)據(jù)庫108號心電圖的去噪效果，其中對比文獻[10]采用改進小波閾值進行去噪，文獻[12]利用集合經(jīng)驗?zāi)B(tài)分解將心電信號分解成一組本征模態(tài)（IMF）分量，然后分別對高頻分量和低頻分量所屬的IMF分量進行改進小波閾值濾波和中值濾波。從圖7虛線框所示可以看出：文獻[10]去噪后心電信號在突變處容易產(chǎn)生變形，且仍存在基線漂移；文獻[12]去除了基線漂移和高頻噪聲，但是基于小波閾值去噪的局限性，在特征波突變處能看到一定程度的幅值和形狀改變；本文方法由于對特征波形處設(shè)置了零誤差，其幅值和形狀得到了較好的保持。

由于MIT?BIH數(shù)據(jù)庫原始心電信號都存在一定噪聲，直接作為參考信號會影響評價指標(biāo)計算的準(zhǔn)確性。本文借鑒文獻[18]的方法，把去噪后的心電信號，如圖7中“本文去噪心電信號”作為參考信號，對其加噪得到圖8中的“加噪心電信號”，然后利用不同方法去噪進行比較，最后計算各評價指標(biāo)來衡量去噪效果。在此以20 dB高斯白噪聲模擬肌電噪聲，以[0.3×sin（2π×0.15t）]模擬基線漂移信號。

圖8展示了加噪心電信號在不同方法下的去噪效果，從虛線框所示依然能夠看出文獻[10]和文獻[12]在去噪后特征波突變處幅值和形狀均有改變，并且文獻[10]的基線漂移不能得到矯正。本文方法在去噪的基礎(chǔ)上能夠更多地保留原始特征，具有更高的保真度。此外，圖9對比了文獻[12]和本文獲取的基線漂移與模擬的基線漂移。從圖9可以看出，本文獲取的基線漂移突變更少，符合基線緩變的特點，也更接近模擬基線漂移，因此可進一步降低基線漂移矯正帶來的信號形態(tài)損失。通過定量計算可得文獻[10]、文獻[12]以及本文去噪后心電信號的SNR分別為-3.200 6 dB、10.117 5 dB、10.740 9 dB；RMSE分別為0.210 3 mV、0.045 4 mV、0.042 2 mV；NCC分別為0.542 7、0.953 1、0.961 2。可以看出本文算法去噪后心電信號SNR最大、RMSE最小、NCC最大，因此整體去噪能力最好。從基線漂移的定量計算結(jié)果可得，文獻[12]和本文獲得的SNR分別為13.853 9 dB和14.027 3 dB；RMSE分別為0.042 2 mV和0.041 4 mV；NCC分別為0.978 0和0.980 0，可以看出本文獲得的基線漂移各項評價指標(biāo)均優(yōu)于文獻[12]。因文獻[10]不能矯正基線，故不作比較。

為了進一步驗證本文方法的有效性，本文另選取了數(shù)據(jù)庫中的105號、106號心電信號進行實驗，實驗中分別疊加了[0.3×sin（2π×0.15t）]模擬的基線漂移和20 dB、15 dB、10 dB的高斯白噪聲，定量計算結(jié)果如表1和表2所示。

表1為不同噪聲下的去噪心電信號的各項指標(biāo)參數(shù)，表2為不同去噪方法獲得的基線漂移與模擬基線漂移之間的各項性能指標(biāo)，加粗部分為最優(yōu)結(jié)果。從兩個表中可以看出，本文的去噪方法各項指標(biāo)更優(yōu)，因此表明本文去噪方法有用信息損失更少，能夠更好地反映原始心電信號特征。

4" 結(jié)" 語

針對心電信號去噪過程中誤差不可控、形態(tài)易變形的問題，本文提出一種基于誤差可控的轉(zhuǎn)換變換去噪方法。通過分析心電信號采樣頻率與基線漂移信號頻率之間的關(guān)系，給出了合適的分解級數(shù)[L]。在分解過程中通過對噪聲處設(shè)置較大誤差界來去除高頻噪聲，對特征波設(shè)置零誤差界來保持信號的原始形態(tài)。通過對第[L]級低頻分量進行cubic插值并規(guī)范化可獲得平滑的基線漂移信號，降低基線校正造成的信號失真。實驗結(jié)果表明，本文方法可以有效去除基線漂移以及高頻噪聲，并保持較高的波形保真度。但是由于本文的轉(zhuǎn)換變換方法本質(zhì)上是一種擴展的Haar小波變換，因此數(shù)據(jù)重構(gòu)后曲線不夠光滑，這也是后期需要進一步改進的地方。

注：本文通訊作者為黎彤亮。

參考文獻

[1] UMA A， KALPANA P. ECG noise removal using modified distributed arithmetic based finite impulse response filter [J]. Journal of medical imaging and health informatics， 2021， 11（5）： 1444?1452.

[2] HESAR H D， MOHEBBI M. An adaptive Kalman filter bank for ECG denoising [J]. IEEE journal of biomedical and health informatics， 2021， 25（1）： 13?21.

[3] PAUL A， PANJA S， DAS N， et al. Median filter based noise reduction and QRS detection in ECG signal [C]// Proceedings of International Conference on Industrial Instrumentation and Control. Heidelberg， Germany： Springer， 2022： 67?76.

[4] TAOULI S A， BEREKSI?REGUIG F. Noise and baseline wandering suppression of ECG signals by morphological filter [J]. Journal of medical engineering amp; technology， 2010， 34（2）： 87?96.

[5] ROMERO F P， PINOL D C， VAZQUEZ?SEISDEDOS C R. DeepFilter： An ECG baseline wander removal filter using deep learning techniques [J]. Biomedical signal processing and control， 2021， 70： 102992.

[6] SINGH P， PRADHAN G. A new ECG denoising framework using generative adversarial network [J]. IEEE/ACM transactions on computational biology and bioinformatics， 2021， 18（2）： 759?764.

[7] YANG H， WEI Z Q. An effective morphological?stabled denoising method for ECG signals using wavelet?based techniques [J]. International journal of biomedical engineering and technology， 2022， 39（3）： 263?282.

[8] 顧旋，張偉，呂珊珊，等.基于改進小波閾值算法的ECG肌電信號去噪研究[J].中國醫(yī)學(xué)物理學(xué)雜志，2023，40（2）：212?219.

[9] HOUAMED I， SAIDI L， SRAIRI F. ECG signal denoising by fractional wavelet transform thresholding [J]. Research on biomedical engineering， 2020， 36（3）： 349?360.

[10] 吳葉麗，行鴻彥，李瑾，等.改進閾值函數(shù)的小波去噪算法[J].電子測量與儀器學(xué)報，2022，36（4）：9?16.

[11] ZHANG D， WANG S， LI F， et al. An efficient ECG denoising method based on empirical mode decomposition， sample entropy， and improved threshold function [EB/OL]. [2020?12?22]. https：//onlinelibrary.wiley.com/doi/10.1155/2020/8811962.

[12] 尹麗，陳富民，張琦，等.采用集合經(jīng)驗?zāi)B(tài)分解和改進閾值函數(shù)的心電自適應(yīng)去噪方法[J].西安交通大學(xué)學(xué)報，2020，54（1）：101?107.

[13] 蔡幫貴，朱雨男，王彪，等.EEMD與改進小波閾值結(jié)合應(yīng)用于心電信號去噪研究[J].現(xiàn)代電子技術(shù)，2023，46（12）：137?140.

[14] PANG C Y， ZHANG Q， ZHOU X F， et al. Computing unrestricted synopses under maximum error bound [J]. Algorithmica， 2013， 65（1）： 1?42.

[15] GUO X Y， LI T L， LI X Y， et al. An efficient multidimensional [L∞] wavelet method and its application to approximate query processing [J]. World wide web， 2021， 24（1）： 105?133.

[16] DYER S A， DYER J S. Cubic?spline interpolation.1 [J]. IEEE instrumentation amp; measurement magazine， 2001， 4（1）： 44?46.

[17] MOODY G B， MARK R G. The impact of the MIT?BIH arrhythmia database [J]. IEEE engineering in medicine biology magazine， 2001， 20（3）： 45?50.

[18] 耿讀艷，趙寧，龐耕，等.改進EWT算法的心電信號去噪研究[J].電子測量技術(shù)，2021，44（23）：30?35.