本發(fā)明涉及一種用于從心電圖波形中提取諸如心跳間隔(R-R間隔)的生物信息的心跳檢測方法和心跳檢測設(shè)備。

背景技術(shù)：

ECG(心電圖)波形是通過觀察和記錄心臟的電活動獲得的，并且是通過以一般方法將電極附著到身體表面來測量的。作為ECG波形導(dǎo)聯(lián)(lead)系統(tǒng)，即電極布置，存在使用肢體和胸部的各種類型。在心前導(dǎo)聯(lián)中的V3至V5導(dǎo)聯(lián)中，電極被布置在左胸部。在適于長時間監(jiān)測ECG波形的CC5導(dǎo)聯(lián)中，電極被布置在左胸部和右胸部的對稱位置處。這些導(dǎo)聯(lián)的優(yōu)點在于能夠獲得具有較大振幅的穩(wěn)定波形。

圖9示出了ECG波形的示例。在圖9中，縱坐標表示電位，橫坐標表示時間。ECG波形由連續(xù)的心跳波形形成，并且一個心跳波形由反映心房和心室的活動的分量形成，比如P波、Q波、R波、S波和T波。

已知從ECG波形獲得的諸如R-R間隔的生物信息是反映自主活動的指標。它對于評估自主功能以獲得日常生活中的ECG波形以及從檢測到的心跳中分析心跳波動的數(shù)據(jù)是有用的。此外，還能夠用于根據(jù)運動期間的心跳數(shù)據(jù)估計運動耐量，并進行優(yōu)化等。

作為常規(guī)心跳檢測方法，以下文獻是已知的。日本專利特開No.2002-78695公開了一種用于消除ECG波形的基線波動的布置。此外，日本專利特開No.2003-561公開了一種使用基于波峰和波谷之間的振幅的閾值來識別R波的布置。

在文獻“ECG Implementation on the TMS320C5515DSP Medical Development Kit(MDK)with the ADS1298ECG-FE”(德州儀器公司，<http：//www.ti.com/lit/an/sprabj1/sprabj1.pdf>，2011)中描述了一種基于通過計算ECG波形的一階導(dǎo)數(shù)而獲得的值的變化來獲得R-R間隔等的方法。在該心跳檢測方法中，更具體地，獲得第(n+1)個采樣值和第(n-1)個采樣值之差的絕對值，基于閾值檢測峰值，然后，將兩個峰之間的時間寬度設(shè)置為R-R間隔。

技術(shù)實現(xiàn)要素：

本發(fā)明解決的問題

然而，上述心跳檢測方法具有以下問題。也就是說，如果要記錄或分析日常生活中或運動期間的心跳數(shù)據(jù)，則與靜止狀態(tài)下的測量不同的是，由身體運動等引起的噪聲可能會被混入ECG波形中。

圖10和11示出了用于說明常規(guī)問題的曲線圖，其中橫坐標表示時間[ms]，縱坐標表示用數(shù)字值代替的電位[任意單位]。在圖10和圖11中，X表示ECG波形的采樣數(shù)據(jù)的示例，R表示要被檢測為心跳的R波。圖10和11中的ECG波形是相同的。

在圖10和圖11中，盡管在21,000ms之后觀測到正常的ECG波形圖案，但是在21,000ms之前噪聲分量疊加在原始ECG波形上。即使試圖從這樣的ECG波形來基于閾值檢測心跳，振幅波動也會很大，從而，難以設(shè)置適當(dāng)?shù)拈撝?。此外，即使根?jù)峰值來周期性地更新閾值，在包括噪聲的部分中閾值也增加。因此，不必要地錯過了沒有任何噪聲的部分中的心跳。

另一方面，還提供了一種基于ECG波形的導(dǎo)數(shù)來檢測心跳的方法。通常，在數(shù)據(jù)處理中將ECG波形作為離散數(shù)據(jù)序列來處理。因此，獲得ECG波形的導(dǎo)數(shù)等效于獲得ECG波形的時間差。在圖10中，通過在ECG波形的采樣數(shù)據(jù)X的每個時刻處繪制通過從5ms之后的值減去5ms之前的值所獲得的值(即，一階導(dǎo)數(shù))來獲得DF。在通常的ECG波形中，通過獲得時間差，可以使從R波到S波的突然改變變得顯著，從而便于心跳的檢測。然而，圖10中的一階導(dǎo)數(shù)值DF包括突然改變的噪聲分量。因此，噪聲分量可能被誤認為是心跳，且難以基于一階導(dǎo)數(shù)值來檢測心跳。在圖10中，N表示容易被誤認為心跳的噪聲分量。

在圖11中，通過在圖10中的一階導(dǎo)數(shù)值DF的每個時刻處繪制通過從5ms之后的值減去5ms之前的值獲得的值(即，二階導(dǎo)數(shù))來獲得DS。除了心跳的峰值之外，圖11中的二階導(dǎo)數(shù)值DS還包括容易被誤認為是心跳的噪聲分量的峰值。

本發(fā)明是為了解決上述問題而提出的，其目的在于提供一種心跳檢測方法和心跳檢測設(shè)備，其能夠從在ECG波形上疊加有由身體運動引起的噪聲的數(shù)據(jù)中正確地檢測心跳及其時間。

解決問題的方法

根據(jù)本發(fā)明，提供了一種心跳檢測方法。該心跳檢測方法包括：計算步驟，根據(jù)活體的心電圖波形的采樣數(shù)據(jù)序列，針對每個采樣時刻計算采樣數(shù)據(jù)的改變量和改變程度之一；乘法步驟，通過將時刻K處的采樣數(shù)據(jù)的改變量和改變程度之一乘以時刻K處的采樣數(shù)據(jù)以及比時刻K早預(yù)定時間t的時刻處的采樣數(shù)據(jù)之一，來針對每個采樣時刻計算乘積；峰值檢測步驟，檢測所述乘積的峰值；以及心跳時刻確定步驟，將所述乘積的峰值的時刻設(shè)置為心跳時刻。

根據(jù)本發(fā)明，還提供了一種心跳檢測設(shè)備。該心跳檢測設(shè)備包括：計算裝置，用于根據(jù)活體的心電圖波形的采樣數(shù)據(jù)序列，針對每個采樣時刻計算采樣數(shù)據(jù)的改變量和改變程度之一；乘法裝置，用于通過將時刻K處的采樣數(shù)據(jù)的改變量和改變程度之一乘以時刻K處的采樣數(shù)據(jù)以及比時刻K早預(yù)定時間t的時刻處的采樣數(shù)據(jù)之一，來針對每個采樣時刻計算乘積；峰值檢測裝置，用于檢測所述乘積的峰值；以及心跳時刻確定裝置，用于將所述乘積的峰值的時刻設(shè)置為心跳時刻。

本發(fā)明的效果

根據(jù)本發(fā)明，針對活體的心電圖波形的采樣數(shù)據(jù)序列中的每個采樣時刻，計算采樣數(shù)據(jù)的改變量或改變程度；通過將時刻K處的采樣數(shù)據(jù)的改變量或改變程度乘以時刻K處的采樣數(shù)據(jù)或比時刻K早預(yù)定時間t的時刻處的采樣數(shù)據(jù)，來針對每個采樣時刻計算乘積；檢測所述乘積的峰值；并將乘積的峰值的時刻設(shè)為心跳時刻。在本發(fā)明中，有可能強調(diào)從心跳導(dǎo)出的峰值分量，并甚至可以從在心電圖波形上疊加有由身體運動引起的噪聲的采樣數(shù)據(jù)序列中正確地檢測心跳。

附圖說明

圖1示出了用于解釋根據(jù)本發(fā)明的原理的圖。

圖2示出了用于解釋根據(jù)本發(fā)明的原理的圖。

圖3是用于解釋根據(jù)本發(fā)明的原理的圖。

圖4是用于解釋根據(jù)本發(fā)明的原理的圖。

圖5是示出了根據(jù)本發(fā)明的實施例的心跳檢測設(shè)備的布置的框圖。

圖6是用于解釋根據(jù)本發(fā)明的實施例的心跳檢測設(shè)備的差值計算單元和乘法單元的操作的流程圖。

圖7是用于解釋根據(jù)本發(fā)明的實施例的心跳檢測設(shè)備的峰值檢測單元和心跳時刻確定單元的操作的流程圖。

圖8是用于解釋根據(jù)本發(fā)明的實施例的心跳檢測設(shè)備的峰值檢測單元的另一操作的流程圖。

圖9是示出了心電圖波形的示例的圖。

圖10示出了用于解釋常規(guī)問題的圖。以及

圖11示出了用于解釋常規(guī)問題的圖。

具體實施方式

[發(fā)明原理]

圖1至圖4示出了用于解釋根據(jù)本發(fā)明的原理的圖。在圖1至圖4中，橫軸表示時間(ms)，縱軸表示用數(shù)字值代替的電位(任意單位)。

圖1示出了圖10中的部分的放大圖?？梢岳斫?，圖1中的R波的峰值R1在一階導(dǎo)數(shù)值DF的峰值D1之前10ms的時刻處。因此，為了強調(diào)峰值R1和D1，對通過將一階導(dǎo)數(shù)值DF乘以采樣值X之前10ms處的采樣值而獲得的數(shù)據(jù)序列進行追蹤，其中一階導(dǎo)數(shù)值DF是采樣值X之前和之后的值之間的差。

圖2示出了圖11中的部分的放大圖?？梢岳斫猓瑘D2中R波的峰R1處于與二階導(dǎo)數(shù)值DS的峰值D2相同的位置處。

圖3是通過繪制通過將圖10中每個時刻處的一階導(dǎo)數(shù)值DF乘以該時刻之前10毫秒處的采樣值而計算的值所獲得的圖。從圖3可以看出，在強調(diào)了與心跳對應(yīng)的峰值M1的同時，圖10中所示出的心跳之間的噪聲分量降低，并且可以基于閾值Th來檢測與心跳相對應(yīng)的峰值M1。

類似地，圖4是通過繪制通過將圖11中每個時刻處的二階導(dǎo)數(shù)值DS乘以圖11中的時刻處的采樣值X而計算的值所獲得的圖。從圖4中可以看出，噪聲分量降低，與圖3類似，可以基于閾值Th檢測與心跳相對應(yīng)的峰值M2。

如上所述，根據(jù)本發(fā)明，在ECG波形的采樣數(shù)據(jù)序列中，將給定時刻處的采樣值的改變量或改變程度乘以給定時刻處的采樣值或比所述給定時間早預(yù)定時間t的時刻處的采樣值，從而將所述乘積的峰值的時刻設(shè)置為心跳時刻。

ECG波形的采樣數(shù)據(jù)序列包括對應(yīng)于每次心跳的R波的峰值分量。此外，針對每次心跳，通過計算ECG波形的導(dǎo)數(shù)而獲得的數(shù)據(jù)序列包括與R波和S波之間的突然改變相對應(yīng)的峰值分量。也就是說，每個數(shù)據(jù)序列包括相同心跳節(jié)奏的峰值分量。如果通過將這些數(shù)據(jù)序列移位給定的時間寬度來使這些數(shù)據(jù)序列重疊，則峰值分量彼此同步。因此，通過在適當(dāng)?shù)臈l件下乘以數(shù)據(jù)序列，有可能強調(diào)從心跳中導(dǎo)出的峰值分量。另一方面，由身體運動等引起的噪聲的波動分量與心跳無關(guān)地出現(xiàn)，并且傾向于通過乘以數(shù)據(jù)序列而被平滑。因此，可以使得僅心跳顯著，從而便于檢測它們。

適當(dāng)?shù)貙⒈硎締畏宀ㄐ蔚囊浑A導(dǎo)數(shù)或二階導(dǎo)數(shù)用作用于觀察ECG波形的變化以及R波和S波之間的突然改變的值。在ECG波形的采樣數(shù)據(jù)序列中，時刻K處的采樣值的一階導(dǎo)數(shù)值是通過從時刻(k+W)處的采樣值減去時刻(K-W)處的采樣值而獲得的(W是例如5ms)。如果使用時刻K處的采樣值的改變量，即，如果將一階導(dǎo)數(shù)值作為乘以采樣值的值，則將比時刻K早預(yù)定時間t的時刻(K-t)處的采樣值乘以一階導(dǎo)數(shù)值。一階導(dǎo)數(shù)值的峰值出現(xiàn)在ECG波形的R波的峰值之后約10至12ms處。因此，如果使用一階導(dǎo)數(shù)值，則將預(yù)定時間t設(shè)定為滿足10ms≤t≤12ms。

通過從時刻(K+W)處的采樣值的一階導(dǎo)數(shù)值減去時刻(k-W)處的采樣值的一階導(dǎo)數(shù)值，獲得時刻K處的采樣值的二階導(dǎo)數(shù)值。如果使用時刻K處的采樣值的改變程度，即，如果將二階導(dǎo)數(shù)值作為乘以采樣值的值，則將時刻K處的采樣值或比時刻K早預(yù)定時間t的時刻(K-t)處的采樣值乘以二階導(dǎo)數(shù)值。二階導(dǎo)數(shù)值的峰值出現(xiàn)在ECG波形的R波的峰值之后大約0到1ms處。因此，如果使用二階導(dǎo)數(shù)值，則將預(yù)定時間t設(shè)定為滿足0ms＜t≤1ms。

當(dāng)心室肌細胞的去極化從內(nèi)膜進展到外膜時，通過捕獲電流以及電位的變化來獲得ECG波形的R波和S波，且電流的移動和速度幾乎不受個體差異(例如，身體類型差別)的影響。因此，峰值之間的時間間隔幾乎是恒定的。

[實施例]

下文將參考附圖描述本發(fā)明的實施例。圖5是示出了根據(jù)本發(fā)明的實施例的心跳檢測設(shè)備的布置的框圖。心跳檢測設(shè)備包括心電圖儀1、存儲單元2、差值計算單元3(計算裝置)、乘法單元4(乘法裝置)、峰值檢測單元5(峰值檢測裝置)和心跳時刻確定單元6(心跳時刻確定裝置)。

下面將描述根據(jù)實施例的心跳檢測方法。在本說明書中，將對檢測心跳和計算心跳的心跳時刻的過程進行說明。通過針對ECG波形數(shù)據(jù)的周期重復(fù)計算心跳時刻，連續(xù)獲得心跳時刻的順序數(shù)據(jù)，且可以根據(jù)該順序數(shù)據(jù)計算心跳波動的指標。

在本實施例中，X(i)表示通過對ECG波形進行采樣而獲得的數(shù)據(jù)序列，其中，i(i＝1、2、…)表示指派給一個采樣數(shù)據(jù)的數(shù)。當(dāng)然，隨著數(shù)i變大，采樣時刻越晚。此外，a表示將當(dāng)獲得采樣數(shù)據(jù)X(i)的一階導(dǎo)數(shù)值時的時間間隔的一半(上述W)除以采樣間隔所獲得的整數(shù)，b表示將當(dāng)采樣數(shù)據(jù)X(i)乘以其一階導(dǎo)數(shù)值時提供的恒定時間差t除以采樣間隔所獲得的整數(shù)，且Th表示用于獲得采樣數(shù)據(jù)X(i)與其一階導(dǎo)數(shù)值的乘積的峰值的閾值。

心電圖儀1測量活體(人體)(未示出)的ECG波形，并輸出ECG波形的采樣數(shù)據(jù)序列X(i)。此時，心電圖儀1通過將采樣時刻信息添加到每個采樣數(shù)據(jù)來輸出數(shù)據(jù)序列。注意到，測量ECG波形的實際方法是公知技術(shù)，并且將省略其詳細描述。

存儲單元2存儲從心電圖儀1輸出的ECG波形的采樣數(shù)據(jù)列X(i)和所述采樣時刻信息。

圖6是用于說明差值計算單元3和乘法單元4的操作的流程圖。差值計算單元3計算每個采樣時刻的采樣數(shù)據(jù)X(i)的一階導(dǎo)數(shù)值(X(i+a)-X(i-a))(圖6的步驟S100)。

乘法單元4針對每個采樣時刻計算采樣數(shù)據(jù)X(i)的一階導(dǎo)數(shù)值(X(i+a)-X(i-a))與比采樣數(shù)據(jù)X(i)早預(yù)定時間t的時刻處的采樣數(shù)據(jù)X(i-b)的乘積(X(i+a)-X(i-a))×X(i-b)(圖6的步驟S101)。

圖7是用于說明峰值檢測單元5和心跳時刻確定單元6的操作的流程圖。峰值檢測單元5搜索乘積(X(i+a)-X(i-a))×X(i-b)超過閾值Th的點。

首先，峰值檢測單元5將用于連續(xù)讀出采樣數(shù)據(jù)序列X(i)的數(shù)(計數(shù)器變量)i設(shè)置為初始值(在本示例中為n)(圖7的步驟S1)。接著，峰值檢測單元5將由乘法單元4計算的乘積(X(i+a)-X(i-a))×X(i-b)與閾值Th進行比較(圖7的步驟S2)。

如果乘積(X(i+a)-X(i-a))×X(i-b)等于或大于閾值Th(步驟S2中的“否”)，則峰值檢測單元5確定在由i指示的時刻附近不存在乘積的峰值，設(shè)置i＝i+1(圖7的步驟S3)，并返回到步驟S2。重復(fù)步驟S2和S3中的處理，直到乘積(X(i+a)-X(i-a))×X(i-b)變得小于閾值Th為止。

如果在步驟S2中確定乘積(X(i+a)-X(i-a))×X(i-b)小于閾值Th，則峰值檢測單元5確定在由i指示的時刻附近存在乘積的峰值，并前進到步驟S4和后續(xù)步驟中的指明峰值位置的過程。

峰值檢測單元5臨時將峰值P設(shè)置為(X(i+a)-X(i-a))×X(i-b)并對其進行存儲(圖7的步驟S4)。峰值檢測單元5將用于檢測乘積的峰值的計數(shù)器變量j和指示峰值位置的變量k設(shè)置為1(圖7的步驟S5)。

峰值檢測單元5將由乘法單元4計算的乘積(X(i+a+j)-X(i-a+j))×X(i-b+j)與當(dāng)前峰值P進行比較(圖7的S6)。如果乘積(X(i+a+j)-X(i-a+j))×X(i-b+j)等于或大于峰值P(步驟S6中為否)，則該處理前進到步驟S9，而不改變峰值P和變量k的值。

如果乘積(X(i+a+j)-X(i-a+j))×X(i-b+j)小于當(dāng)前峰值P，則峰值檢測單元5將峰值P更新為(X(i+a+j)-X(i-a+j))×X(i-b+j)，并對其進行存儲(圖7的步驟S7)。峰值檢測單元5將變量k替換為j(圖7的步驟S8)，并前進到步驟S9。在步驟S9中，峰值檢測單元5確定計數(shù)器變量j是否不超過預(yù)定值jmax，其中所述預(yù)定值jmax指定獲得峰值的范圍。

如果計數(shù)器變量j不超過預(yù)定的jmax，則峰值檢測單元5設(shè)置j＝j(luò)+1(圖7的步驟S10)，并且返回到步驟S6。重復(fù)步驟S6至S10中的處理，直到計數(shù)器變量j超過預(yù)定值jmax為止。如果計數(shù)器變量j超過預(yù)定值jmax，則峰值檢測單元5基于計數(shù)器變量j結(jié)束對峰值P的搜索。此時，最新的峰值P的時刻(即，由(i+k)指示的時刻)是心跳時刻的候選。

接著，心跳時刻確定單元6確定檢測出的心跳時刻是否正確，并選擇性地固定心跳時刻。

首先，心跳時刻確定單元6判斷由(i+k)指示的時刻T與之前檢測的心跳時刻T_(-1)之差是否等于或大于預(yù)定時間(圖7的步驟S11)。如果時刻T與前一心跳時刻T_(-1)之差小于預(yù)定時間，則心跳時刻確定單元6丟棄由(i+k)指示的時刻T，不采用它作為心跳時刻，并返回到步驟S3。

心跳間隔存在一般正常值范圍。如果檢測到與該范圍相比非常短的心跳間隔，則很有可能將由身體運動等而疊加在心電圖波形上的噪聲錯誤地識別為心跳。通過設(shè)置令乘積的檢測峰值的時刻T與前一心跳時刻T_(-1)之差等于或大于預(yù)定時間這一條件，可以防止由噪聲等引起的錯誤檢測。

此外，心跳時刻確定單元6確定當(dāng)由(i+k)所指示的時刻T被認為是心跳時刻時的心跳間隔(T-T_(-1))是否沒有以預(yù)定速率或更大速率從前一心跳間隔(T_(-1)-T_(-2))增加(圖7的步驟S12)。如果心跳間隔的增加速率(T-T_(-1))/(T_(-1)-T_(-2))等于或高于預(yù)定值，則心跳時刻確定單元6確定心跳間隔以預(yù)定速率或更大速率增加，并丟棄由(i+k)指示的時刻T，而不采用它作為心跳時刻，從而返回到步驟S3。

如果對給定心跳的檢測失敗，則獲得作為給定心跳前后的心跳之間的心跳間隔的數(shù)據(jù)指示大約是實際值的兩倍的值，且不適于用于評價自主功能等。通過設(shè)置令檢測到的心跳間隔沒有以預(yù)定速率或更大速率增加的這一條件，有可能從生物信息的分析目標中排除心跳檢測失敗的錯誤數(shù)據(jù)。

如果由(i+k)指示的時刻T與前一心跳時刻T_(-1)之差等于或大于預(yù)定時間，且心跳間隔的增加率(T-T_(-1))/T_(-1)-T_(-2))小于預(yù)定值，則心跳時刻確定單元6將由(i+k)所指示的時刻T作為心跳時刻(圖7的步驟S13)。

在步驟S13結(jié)束之后，設(shè)置i＝i+1，且該處理返回到步驟S2。這開始了對下一個心跳的檢測。備選地，可以跳過數(shù)值與預(yù)定時間相對應(yīng)的i的值，且該處理可以返回到步驟S2，其中所述預(yù)定時間比將要檢測的心跳間隔的最小值更短。通過重復(fù)步驟S2至S13中的處理，可以獲得心跳時刻的順序數(shù)據(jù)，并從順序數(shù)據(jù)獲得心跳波動的指標。

注意，可以基于到目前為止獲得的峰值P的平均值等來周期性地更新閾值Th。圖8示出了這種情況中的流程圖。如果在步驟S9中確定為“是”，則峰值檢測單元5將通過將在測量開始之后檢測到的峰值P的均值乘以預(yù)定系數(shù)r而獲得的值更新為最新閾值Th(圖8的步驟S14)。當(dāng)要檢測第m次心跳時，已經(jīng)檢測了(m-1)個峰值P，峰值P是最新檢測到的。因此，總共獲得m個峰值P。因此，將m個峰值P的平均值乘以系數(shù)r。注意，當(dāng)在步驟S11或S12中確定為“否”時的峰值不被用于計算閾值Th。

在ECG波形中，每個分量的幅度或改變程度根據(jù)個體差異或電極附著的方式而變化。通過基于到目前為止獲得的峰值P的平均值來設(shè)置用于檢測乘積的峰值的閾值Th，能夠降低由個體差異等引起的波形變化的影響。圖8中除了步驟S14的處理與參考圖7描述的處理相同。

通過將根據(jù)本實施例的心跳檢測方法應(yīng)用于獲得較大R波和較深S波的ECG導(dǎo)聯(lián)的ECG波形(例如，V3到V5導(dǎo)聯(lián)之一的ECG波形)，可以獲得顯著效果。特別優(yōu)選的是，將該方法應(yīng)用于CC5導(dǎo)聯(lián)或通常用于在日常生活中獲得ECG波形的類似導(dǎo)聯(lián)的ECG波形。

在本實施方式的心跳檢測方法中，能夠獲得正確心跳時刻的數(shù)據(jù)序列，并能夠基于該數(shù)據(jù)序列獲得心跳波動的高可靠性指標。

注意，本實施例已經(jīng)說明了使用采樣數(shù)據(jù)的一階導(dǎo)數(shù)值來檢測心跳的情況。然而，如上所述，可以使用二階導(dǎo)數(shù)值。在這種情況下，差值計算單元3針對每個采樣時刻計算采樣數(shù)據(jù)X(i)的二階導(dǎo)數(shù)值(圖6的步驟S100)。當(dāng)用DF(i)＝(X(i+a)-X(i-a))表示采樣數(shù)據(jù)X(i)的一階導(dǎo)數(shù)值時，采樣數(shù)據(jù)X(i)的二階導(dǎo)數(shù)值表示為(DF(i+a)-DF(i-a))。

因此，乘法單元4針對每個采樣時刻計算采樣數(shù)據(jù)X(i)的二階導(dǎo)數(shù)值(DF(i+a)-DF(i-a))與采樣數(shù)據(jù)X(i)的乘積(DF(i+a)-DF(i-a))×X(i)，或二階導(dǎo)數(shù)值(DF(i+a)-DF(i-a))與比采樣數(shù)據(jù)X(i)早預(yù)定時間t的時刻處的采樣數(shù)據(jù)X(i-b)的乘積(DF(i+a)-DF(i-a))×X(i-b)(圖6的步驟S101)。

在圖7和圖8的描述中，由(DF(i+a)-DF(i-a))×X(i)或(DF(i+a)-DF(i-a))×X(i-b)來代替乘積(X(i+a)-X(i-a))×X(i-b)，并由DF(i+a+j)-DF(i-a+j))×X(i+j)或DF(i+a+j)-DF(i-a+j))×X(i-b+j)來代替乘積(X(i+a+j)-X(i-a+j))×X(i-b+j)。因此，可以使用采樣數(shù)據(jù)的二階導(dǎo)數(shù)值來檢測心跳。

本實施例中所述的存儲單元2、差值計算單元3、乘法單元4、峰值檢測單元5和心跳時刻確定單元6可以由包括CPU(中央處理器)、存儲設(shè)備、接口以及控制這些硬件資源的程序的計算機來實現(xiàn)。CPU根據(jù)存儲在存儲設(shè)備中的程序執(zhí)行本實施例中描述的處理。

[工業(yè)實用性]

本發(fā)明適用于檢測活體的心跳的技術(shù)。

附圖標記的解釋

1心電圖儀；2存儲單元；3差值計算單元；4乘法單元；5峰值檢測單元；6心跳時刻確定單元。