一種多級導(dǎo)聯(lián)心電信號qrs波形識別方法
【專利摘要】本發(fā)明涉及一種心電信號QRS波形識別的方法。本發(fā)明首先對采集到的心電信號進行濾波，再對II級導(dǎo)聯(lián)心電信號的R波進行識別，依據(jù)II級導(dǎo)聯(lián)信號的R波的位置作為參考對其他多級導(dǎo)聯(lián)心電信號的R波進行識別，之后對多級導(dǎo)聯(lián)信號的Q波S波進行識別，以實現(xiàn)對多導(dǎo)聯(lián)心電信號QRS波形識別。本發(fā)明的有益效果為，在小波變換識別R波提出了新的方法，可以避免傳統(tǒng)極大值-極小值方法中極大值與極小值的匹配錯誤。本發(fā)明尤其適用于心電信號QRS波形識別。
【專利說明】一種多級導(dǎo)聯(lián)心電信號QRS波形識別方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種心電信號QRS波形識別的方法，尤其涉及一種基于小波變換的多級導(dǎo)聯(lián)心電信號QRS波形識別方法。

【背景技術(shù)】
[0002]由于我國老齡化的加劇以及現(xiàn)有的醫(yī)療資源分布不平衡，醫(yī)療健康領(lǐng)域在我國變得越來越重要。其中心血管疾病作為危害現(xiàn)代人身體健康的重要疾病，心血管疾病的預(yù)防與治療變得尤為重要。心電信號波形的識別作為心血管疾病治療的一個重要環(huán)節(jié)，其理論研究需要提供高正確率。心電信號波形識別的步驟往往是從QRS波形識別開始。
[0003]在傳統(tǒng)的心電信號檢測過程中，通常采用標(biāo)準(zhǔn)的十二級導(dǎo)聯(lián)形式、單級導(dǎo)聯(lián)形式以及其他改進的導(dǎo)聯(lián)形式。但在傳統(tǒng)的心電信號波形的識別方法中，往往僅識別II級導(dǎo)聯(lián)波形。然而有些心律失常信號II級導(dǎo)聯(lián)波形與正常信號II級導(dǎo)聯(lián)波形差異很小，卻在其他級導(dǎo)聯(lián)波形上具有明顯的差異。比如對于室性早搏信號，其在Vl導(dǎo)聯(lián)上的特征信息比II級導(dǎo)聯(lián)上的信號的特征信息更為明顯。因此多級導(dǎo)聯(lián)信號QRS波形識別對其后續(xù)的特征信號提取以及心律失常信號分類很有幫助。本發(fā)明所述多級導(dǎo)聯(lián)是指標(biāo)準(zhǔn)十二級導(dǎo)聯(lián)形式或其他改進的導(dǎo)聯(lián)形式中的包含II級導(dǎo)聯(lián)信號的多個導(dǎo)聯(lián)信號組合，如后述的MLII級導(dǎo)聯(lián)與Vl級導(dǎo)聯(lián)的組合、MLII級與V5級導(dǎo)聯(lián)的組合。所述MLII級導(dǎo)聯(lián)是指改進肢體II級導(dǎo)聯(lián)。
[0004]心電信號QRS波形識別主要有如下方法:基于濾波和閾值檢測的方法、基于小波變換的方法，以及基于神經(jīng)網(wǎng)絡(luò)的方法。其中濾波和閾值檢測方法是基于時域的分析，對于大量的心電信號難以選取一個適合所有心電信號的閾值。神經(jīng)網(wǎng)絡(luò)方法需要提前訓(xùn)練，而且其訓(xùn)練時間長，難以實際應(yīng)用。小波變換是時間和頻率的局部變換，能有效地從信號中提取時頻域信息，適合心電信號QRS波形識別。
[0005]由于用戶在進行心電信號檢測時，無法確定用戶的心電信號是正常心電信號，還是包含一種或多種心律失常心電信號。因此，心電信號QRS波形方法識別必須能夠識別出所有心律失常心電信號的QRS波形。各種心律失常心電信號在II級導(dǎo)聯(lián)信號中都是正向的，而在一些其他級導(dǎo)聯(lián)信號中并不一定是正向的。這給多級導(dǎo)聯(lián)心電信號QRS波形識別帶來了一些困難。
[0006]傳統(tǒng)的使用小波變換進行R波識別方法往往是在頻域上通過閾值來獲取極大值和極小值，再根據(jù)這個極大值和極小值獲取極大值-極小值對中的過零點。將這個過零點視為R波出現(xiàn)的位置。對有些波形，使用這種方法獲得的極大值個數(shù)與極小值個數(shù)不匹配，從而造成無法獲取的極大值-極小值對的過零點。而且如果其中的一個極大值與極小值匹配錯誤，可能造成接下來的其他的極大值與極小值匹配錯誤，從而形成大量的錯誤的極大值-極小值對的過零點。
[0007]心電信號Q波S波的識別存在一些不方便識別的情況:一些心電信號Q波或S波可能消失；室性早搏II級導(dǎo)聯(lián)信號的S波消失，其隨后緊跟著倒置的T波。在識別S波時可能把倒置的T波作為QRS復(fù)合波的一部分進行識別。
[0008]傳統(tǒng)的使用小波變換進行Q波S波識別的方法往往采用類似與R波識別的頻域極大值-極小值對過零點的方法。R波對應(yīng)的過零點的前面一個過零點是Q波，其后的過零點是S波。然而如上述的一些特殊情況這種傳統(tǒng)方法沒有考慮從而會造成錯誤的識別。
[0009]心電信號QRS波形識別需要在多級導(dǎo)聯(lián)信號中識別多種心律失常信號的QRS波形，且識別的方法要求準(zhǔn)確性高，從而對心電信號的識別提高了要求。


【發(fā)明內(nèi)容】

[0010]本發(fā)明的目的，就是針對上述問題，提出一種提出了多級導(dǎo)聯(lián)心電信號QRS波形識別方法。
[0011]本發(fā)明的技術(shù)方案:一種多級導(dǎo)聯(lián)心電信號QRS波形識別方法，其特征在于，包括以下步驟:
[0012]a.濾波:對采集到的心電信號進行濾波；具體方法為:
[0013]al.使用Daubechies小波的5階小波db5對心電信號進行分解,得到第7尺度低頻信號對應(yīng)心電信號的基線漂移噪聲，將其濾除再進行重構(gòu)從而濾除基線漂移噪聲；
[0014]a2.使用symlet小波8階小波sym8對各層噪聲分別進行估計、調(diào)整進行硬閾值濾波，從而濾除工頻干擾噪聲和肌電干擾噪聲；
[0015]b.對II級導(dǎo)聯(lián)心電信號的R波進行識別:使用硬閾值進行極大值檢測，再進行RR間距的調(diào)整，之后進行時域調(diào)整識別R波；具體方法為:
[0016]bl.使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的高頻信號；
[0017]b2.設(shè)定固定的幅度閾值THR_amplitude，查找大于該幅度閾值的信號點并存儲為X_i;對X-1做差分，其差分結(jié)果為大于固定閾值的點之間的間距X-diff ；
[0018]b3.設(shè)置RR間距閾值THRfinteval，選擇大于該間距閾值的信號X_intCTval，X_inte?al是每段大于幅度閾值的信號的右側(cè)邊界點的位置，從而可以獲得右側(cè)邊界點X ok ；
[0019]b4.根據(jù)大于幅度閾值的信號的右側(cè)邊界位置X_intOTval來查找每段大于幅度閾值的信號的最大值，該最大值為高頻信號的極大值；
[0020]b5.根據(jù)極大值位置在時域中進行時域調(diào)整，具體為:在濾波后的心電信號中的極大值位置左右1/12秒的時間窗內(nèi)查找最大值，該最大值便是時域上的R波；
[0021]b6.在時域上再次進行RR間距閾值濾波，具體為:設(shè)定RR間距對步驟b5得到的R波進行濾波，濾波后的R波作為識別得到的R波；所述RR間距的范圍為0.6?I秒；
[0022]c.心電信號其他多級導(dǎo)聯(lián)R波識別:根據(jù)II級導(dǎo)聯(lián)中識別的R波的位置作為基礎(chǔ)再進行時域調(diào)整，獲得其他級導(dǎo)聯(lián)的R波識別；具體方法為:
[0023]Cl.檢測其他級導(dǎo)聯(lián)時間窗內(nèi)的最大值和最小值的絕對值，比較最大值與最小值的絕對值，如果最大值大于最小值的絕對值。若是，則進入步驟c2 ;若否，則進入步驟c3:
[0024]c2.判斷R波是正向的，找出相應(yīng)的最大值作為識別的R波；
[0025]c3.判斷R波是負向的，找出相應(yīng)的最小值作為識別的R波；
[0026]d.心電信號Q波S波識別:用小波分解獲取第4尺度的低頻信號，在低頻信號上查找R波波峰兩側(cè)的波谷，在波谷左右兩側(cè)查找極大值，之后進行時域調(diào)整識別Q波S波；具體方法為:
[0027]dl.使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的低頻信號；
[0028]d2.在低頻信號中給定的時間窗內(nèi)查找R波波峰左右兩側(cè)的最小值來識別R波波谷;
[0029]d3.在R波波谷左右兩側(cè)給定的時間窗內(nèi)取最大值，判斷該最大值是否在左右兩側(cè)邊沿上，若否，則進入步驟d3，若是，則在時間窗內(nèi)取極大值，所述極大值為從一側(cè)沿一個方向找到大于前后兩邊的幅度的位置；
[0030]d4.在濾波后的心電信號上進行時域調(diào)整，具體為在微小的時間窗上查找極大值，若查到到極大值，則通過極大值識別Q波S波；若不存在極大值，則判定Q波S波消失，并將消失的Q波S波定為R波波谷。
[0031]本發(fā)明的有益效果為，本發(fā)明對多級導(dǎo)聯(lián)的心電信號進行QRS波形識別而非僅僅對II級導(dǎo)聯(lián)波形進行識別，這樣可以為后續(xù)特征提取過程提供更全面的特征信息；本發(fā)明在小波變換識別R波提出了新的方法，可以避免傳統(tǒng)極大值-極小值方法中極大值與極小值的匹配錯誤；本發(fā)明的Q波S波識別可以有效識別Q波消失或S波的消失以及S波消失后緊跟T波等特殊情況，避免了上述情況對波形識別造成的誤判；本發(fā)明所提供的波形識別方法準(zhǔn)確率達到實際應(yīng)用水平，可以為后續(xù)的特征提取過程提供更準(zhǔn)確的特征信息。

【專利附圖】

【附圖說明】
[0032]圖1濾除噪聲前的心電信號；
[0033]圖2濾除噪聲后的心電信號；
[0034]圖3使用sym8小波單支重構(gòu)第四尺度高頻信號；
[0035]圖4大于幅度閾值的信號X i ；
[0036]圖5大于幅度閾值信號的右側(cè)邊界邊界點X ok ；
[0037]圖6高頻信號的極大值；
[0038]圖711級導(dǎo)聯(lián)和Vl級導(dǎo)聯(lián)的R波識別效果圖；
[0039]圖8100秒內(nèi)的R波識別的效果；
[0040]圖9正常信號的MLII級導(dǎo)聯(lián)和Vl級導(dǎo)聯(lián)的Q波和S波的識別；
[0041]圖10對正常信號的MLII級導(dǎo)聯(lián)和V5級導(dǎo)聯(lián)的Q波和S波的識別；
[0042]圖11多種心律異常信號的MLII級導(dǎo)聯(lián)和Vl即導(dǎo)聯(lián)識別。

【具體實施方式】
[0043]下面結(jié)合附圖對本發(fā)明進行詳細的描述
[0044]本發(fā)明的過程如下:首先對采集到的心電信號進行濾波，再對II級導(dǎo)聯(lián)心電信號的R波進行識別，依據(jù)II級導(dǎo)聯(lián)信號的R波的位置作為參考對其他多級導(dǎo)聯(lián)心電信號的R波進行識別，之后對多級導(dǎo)聯(lián)信號的Q波S波進行識別，以實現(xiàn)對多導(dǎo)聯(lián)心電信號QRS波形識別。
[0045]濾波過程描述如下:采集到的心電信號包含基線漂移噪聲、工頻干擾噪聲、肌電干擾噪聲等噪聲?；€漂移噪聲是低頻噪聲，工頻干擾噪聲和肌電干擾噪聲屬于高頻噪聲。本發(fā)明使用Daubechies小波的5階小波db5對心電信號進行分解，得到第7尺度低頻信號對應(yīng)心電信號的基線漂移噪聲，將其濾除再進行重構(gòu)得到濾除基線漂移噪聲的效果。濾除工頻干擾噪聲和肌電干擾噪聲是采用小波進行一維自動降噪，這里使用symlet小波8階小波sym8對各層噪聲分別進行估計、調(diào)整進行硬閾值濾波。采集到的原始心電信號如圖1所示，濾除噪聲后的心電信號如圖2所示。本文在濾除噪聲后的心電信號的基礎(chǔ)上進行QRS波形識別，以避免這些噪聲對QRS波形的識別造成影響。
[0046]心電信號II級導(dǎo)聯(lián)R波識別過程:本發(fā)明提出使用硬閾值進行極大值檢測，再進行RR間距的調(diào)整，之后進行時域調(diào)整。在時域調(diào)整中，主要是在給定的較小時間窗內(nèi)尋找最大值，以及根據(jù)RR間距進行再調(diào)整。
[0047]本發(fā)明沒有直接在濾除噪聲后的信號直接查找間距的最大值，是因為每條記錄的心電信號是對不同的人進行的檢測，對應(yīng)的波形的幅度值存在偏差。并且有些心律失常心電信號的T波的幅度大于R波的幅度。因此難以選擇一個大于所有記錄的T波又小于所有記錄的R波的閾值。
[0048]本發(fā)明使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的高頻信號。第四尺度高頻信號如圖3所示。II級導(dǎo)聯(lián)高頻信號僅在R波波峰處對應(yīng)的位置上幅度較大，在其他位置上幅度很小。因此可以設(shè)定一個固定的幅度閾值THR_amplitude，查找大于該幅度閾值的信號點并存儲為X-J如圖4所示)，對X_i做差分，差分結(jié)果就是大于固定閾值的點之間的間距X_diff。根據(jù)觀察，設(shè)置一個RR間距閾值THR_.1nteval，選擇大于該間距閾值的信號X-1ntCTval，1_-是每段大于幅度閾值的信號的右側(cè)邊界點的位置。從而可以獲得右側(cè)邊界點X λ(如圖5所示)。
[0049]具體的公式如下，其中ECG high_fraiumc;y為心電信號的第四尺度的高頻信號。
[0050]X j — find (ECG high-frequency^THR EmpIitutJe)
[0051]X diff = Xi+1 -X i, i = I...Length(Xi)-1
[0052]X interval — find (X diff>THR ER intervai)
[0053]X ok = Xi(X interval)
[0054]根據(jù)大于幅度閾值的信號的右側(cè)邊界位置X_inteval來查找每段大于幅度閾值的信號的最大值。該最大值便是高頻信號的極大值(如圖6所示)，由于極大值點與極大值-極小值過零點的間距很小，完全可以通過時域調(diào)整來彌補誤差，并且這樣可以避免獲取極大值-極小值過零點的過程中的繁瑣處理以及一些不當(dāng)?shù)奶幚韼淼娜毕荩虼吮景l(fā)明直接采用極大值而非極大值-極小值過零點。
[0055]根據(jù)圖6中的極大值位置在時域中進行時域調(diào)整:在濾除噪聲后的心電信號中的極大值位置左右1/12秒的時間窗內(nèi)查找最大值，該最大值便是時域上的R波。
[0056]由于在高頻信號中的間距閾值是為了選取大于幅度閾值的點，因此這里的RR間距閾值要比正常的RR間距小一些。因此可能存在將幅度大的T波識別為R波的誤判，因此這里在時域上再次進行RR間距閾值濾波。心率是表示心臟每分鐘跳動的次數(shù)，心臟每跳動一次都會帶來一個R波波峰。因此RR間距可以從心率中進行估算。人體的心率范圍是60拍/分到100拍/分，因此RR間距在0.6秒到I秒之間。因此這里設(shè)定0.6秒為RR間距。對上述估計的R波進行濾波。濾除后的R波作為識別得到的R波。
[0057]心電信號其他多級導(dǎo)聯(lián)R波識別過程如下:其他多級導(dǎo)聯(lián)信號可以根據(jù)相同的方法進行識別，但是這樣會增加計算量。而且對于一些特殊情況，可能出現(xiàn)其他級導(dǎo)聯(lián)識別的R波信號與II級導(dǎo)聯(lián)識別的R波數(shù)量不匹配。這會影響到后續(xù)特征信息的提取。因此，本發(fā)明采用的方法是對于其他導(dǎo)聯(lián)的心電信號采用II級導(dǎo)聯(lián)中識別的R波的位置作為基礎(chǔ)再進行時域調(diào)整，以獲得其他級導(dǎo)聯(lián)的R波識別。
[0058]心電信號其他級導(dǎo)聯(lián)R波識別的時域調(diào)調(diào)整與II級導(dǎo)聯(lián)中R波識別的時域調(diào)整的時域調(diào)整相似。只不過由于一些心律失常心電信號的其他級導(dǎo)聯(lián)的R波存在負向波。因此在進行時域調(diào)整時需要提前進行判斷:檢測其他級導(dǎo)聯(lián)時間窗內(nèi)的最大值，將該最大值與閾值進行比較，如果大于該閾值，則認(rèn)為R波是正向的，否則R波是負向的。再找出相應(yīng)的最大值或最小值作為識別的R波。
[0059]最終II級導(dǎo)聯(lián)和Vl級導(dǎo)聯(lián)的R波識別的結(jié)果如圖7所示。100秒內(nèi)的心電信號的R波識別如圖8所示。由此可見R波的識別率非常高，并且沒有出現(xiàn)傳統(tǒng)方法中由于幅度最大值與幅度最小值之間的偏差很大而造成的幅度小的心電信號的漏檢的情況。也沒有出現(xiàn)傳統(tǒng)的模極大值-極小值對的不匹配的情況，說明本發(fā)明提供的方法的程序的準(zhǔn)確率和魯棒性效果很好。
[0060]心電信號Q波S波識別過程:使用小波分解獲取第4尺度的低頻信號，在低頻信號上查找R波波峰兩側(cè)的波谷，在波谷左右兩側(cè)查找極大值，之后進行時域調(diào)整。
[0061]本發(fā)明使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的低頻信號。沒有直接在濾波后的信號上進行Q波S波識別是因為濾除噪聲后的心電信號仍然包含一些毛刺，這些毛刺在進行查找R波波谷和Q波S波時會造成誤判。
[0062]本發(fā)明在低頻信號中給定的時間窗內(nèi)查找R波波峰左右兩側(cè)的最小值來識別R波波谷。這里需要在一定的時間窗內(nèi)進行查找是因為有可能不存S波并跟著一個倒置的T波，因此可能將倒置的T波波峰誤識別為R波波谷。
[0063]在R波波谷左右兩側(cè)給定的時間窗內(nèi)取最大值，如果這個最大值在左右兩側(cè)邊沿上，則說明該最大值在P波/T波的邊沿上，已經(jīng)越過了 Q波/S波。對于上述情況在時間窗內(nèi)取極大值。取極大值即從一側(cè)沿一個方向找到大于前后兩邊的幅度的位置。
[0064]在低頻信號上找到Q波S波的大致位置后，在濾除噪聲后的心電信號上進行時域調(diào)整，即在在微小的時間窗上查找極大值。由于心電信號的波形可以通過R波來確定，因此其他級的心電信號的心電波形數(shù)目與II級導(dǎo)聯(lián)的心電信號波形數(shù)目相同。于是在其他級導(dǎo)聯(lián)信號上采用在II級導(dǎo)聯(lián)信號上相同的方法識別Q波S波，這里不再贅述。
[0065]為應(yīng)對前述中Q波S波消失，本發(fā)明取在極大值時，如果極大值不存在，則將消失的Q波S波定為R波波谷。為應(yīng)對T波倒置，本發(fā)明在上述的每步操作都給定一定的時間窗，使得S波的識別不會被誤識別到倒置的T波上。
[0066]有些心電信號的一些導(dǎo)聯(lián)的心電波形可能不存在完整的QRS波，如正常信號的Vl級導(dǎo)聯(lián)的波形。由于在獲得該信號前并不知道該信號時正常信號還是心律失常信號，因此無法依據(jù)不同的情況對信號做出不同的處理。這里本發(fā)明對這些沒有完整的QRS波的信號仍然采取相同的方法，雖然沒有完整的QRS波形，但是可以提供一些特征信息，并不影響后續(xù)的波形處理。
[0067]對正常信號的MLII級導(dǎo)聯(lián)和Vl級導(dǎo)聯(lián)Q波和S波的識別如圖9所示。對正常信號的MLII級導(dǎo)聯(lián)和V5級導(dǎo)聯(lián)的Q波S波的識別如圖10所示。對多種心律異常信號的MLII級導(dǎo)聯(lián)和V5即導(dǎo)聯(lián)識別如圖11所示。在這些圖像上最高點和最低點表示R波，R波兩側(cè)是R波波谷，在R波左側(cè)波谷的左側(cè)的點表示Q波，在R波右側(cè)波谷右側(cè)的點表示S波。如果在一個完整的心電圖波形中的點不足五個，說明相應(yīng)的Q波或S波消失，則將R波波谷作為Q波或S波。由圖11可見，標(biāo)記為5的室性早搏信號沒有將T波誤判為S波。這些圖形中的心電信號的QRS波形識別均正確。說明本發(fā)明能夠在QRS波形識別中取得很好的效果。
【權(quán)利要求】
1.一種多級導(dǎo)聯(lián)心電信號QRS波形識別方法，其特征在于，包括以下步驟: a.濾波:對采集到的心電信號進行濾波；具體方法為: al.使用Daubechies小波的5階小波db5對心電信號進行分解，得到第7尺度低頻信號對應(yīng)心電信號的基線漂移噪聲，將其濾除再進行重構(gòu)從而濾除基線漂移噪聲； a2.使用symlet小波8階小波sym8對各層噪聲分別進行估計、調(diào)整進行硬閾值濾波，從而濾除工頻干擾噪聲和肌電干擾噪聲； b.對II級導(dǎo)聯(lián)心電信號的R波進行識別:使用硬閾值進行極大值檢測，再進行RR間距的調(diào)整，之后進行時域調(diào)整識別R波；具體方法為: bl.使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的高頻信號；b2.設(shè)定固定的幅度閾值THR_amplitude，查找大于該幅度閾值的信號點并存儲為X_i，對X_i做差分，其差分結(jié)果為大于固定閾值的點之間的間距x_diff ； b3.設(shè)置RR間距閾值THR_KK_inte,val，選擇大于該間距閾值的信號X-1nterval，X-1nterval 是每段大于幅度閾值的信號的右側(cè)邊界點的位置，從而可以獲得右側(cè)邊界點X ok ； b4.根據(jù)大于幅度閾值的信號的右側(cè)邊界位置X_inteval來查找每段大于幅度閾值的信號的最大值，該最大值為高頻信號的極大值； b5.根據(jù)極大值位置在時域中進行時域調(diào)整，具體為:在濾波后的心電信號中的極大值位置左右1/12秒的時間窗內(nèi)查找最大值，該最大值便是時域上的R波； b6.在時域上再次進行RR間距閾值濾波，具體為:設(shè)定RR間距對步驟b5得到的R波進行濾波，濾波后的R波作為識別得到的R波；所述RR間距的范圍為0.6?I秒； c.心電信號其他多級導(dǎo)聯(lián)R波識別:根據(jù)II級導(dǎo)聯(lián)中識別的R波的位置作為基礎(chǔ)再進行時域調(diào)整，獲得其他級導(dǎo)聯(lián)的R波識別；具體方法為: Cl.檢測其他級導(dǎo)聯(lián)時間窗內(nèi)的最大值，判斷最大值是否大于閾值；若是，則進入步驟c2，若否，則進入步驟c3 ； c2.判斷R波是正向的，找出相應(yīng)的最大值作為識別的R波； c3.判斷R波是負向的，找出相應(yīng)的最小值作為識別的R波； d.心電信號Q波S波識別:用小波分解獲取第4尺度的低頻信號，在低頻信號上查找R波波峰兩側(cè)的波谷，在波谷左右兩側(cè)查找極大值，之后進行時域調(diào)整識別Q波S波；具體方法為: dl.使用sym8小波對一維小波系數(shù)進行單支重構(gòu)，獲取第四尺度的低頻信號；d2.在低頻信號中給定的時間窗內(nèi)查找R波波峰左右兩側(cè)的最小值來識別R波波谷；d3.在R波波谷左右兩側(cè)給定的時間窗內(nèi)取最大值，判斷該最大值是否在左右兩側(cè)邊沿上，若否，則進入步驟d3，若是，則在時間窗內(nèi)取極大值，所述極大值為從一側(cè)沿一個方向找到大于前后兩邊的幅度的位置； d4.在濾波后的心電信號上進行時域調(diào)整，具體為在微小的時間窗上查找極大值，若查到到極大值，則通過極大值識別Q波S波；若不存在極大值，則判定Q波S波消失，并將消失的Q波S波定為R波波谷。
【文檔編號】A61B5/0472GK104161510SQ201410426278
【公開日】2014年11月26日 申請日期:2014年8月27日 優(yōu)先權(quán)日:2014年8月27日
【發(fā)明者】竇衡, 孔飛, 張德波, 李玲 申請人:電子科技大學(xué)