本發(fā)明涉及生物醫(yī)學工程技術領域，特別是涉及一種全數(shù)字心電信號的噪聲檢測方法及系統(tǒng)。

背景技術：

心電信號在采集的過程中，不可避免地會包含許多噪聲，主要的有以下三種：基線干擾噪聲、工頻噪聲和肌電噪聲。這些噪聲對計算機自動診斷或是人工診斷均會帶來干擾，嚴重時甚至能導致錯誤的結論。因此在信號的采集中，往往加入噪聲去除算法，通常采用濾波的方法。但由于噪聲與心電信號的頻率本身有一定范圍的疊加，在濾除噪聲的同時不可避免地會對有用的心電信號產(chǎn)生影響。

同時，在心電信號的采集過程中，有的時候即便是加入了噪聲去除算法，由于噪聲干擾較大，采集的心電信號也無法作為實際有用信號進行分析，特別是對于經(jīng)驗不足的人員來說，由于對什么時候“噪聲小”的概念較為模糊，在信號的采集中可能會做多次重復的、無意義的工作。

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

針對噪聲消除會影響心電信號及信號采集時不能及時反映噪聲大小的問題，本發(fā)明提供一種全數(shù)字心電信號噪聲檢測方法，實現(xiàn)心電信號采集或是分析過程中的噪聲檢測。

本發(fā)明所采用的技術方案是：提供一種全數(shù)字心電信號的噪聲檢測方法，包括如下步驟：

s100、采集的模擬心電信號轉換后得到數(shù)字心電信號；

s200、檢測所述數(shù)字心電信號中的qrs波，將檢測結果進行qrs波起始點和終止點的計算，進行qrs波成分消除，并獲取當前基線值；

s300、對消除qrs波后的信號做高通濾波處理，并分別進行工頻噪聲檢測和肌電噪聲檢測。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測方法中，所述步驟s200還包括如下步驟：

s201、所述數(shù)字心電信號分別進行信號延時處理獲得延時信號，同時檢測所述數(shù)字心電信號中的qrs波；

s202、將檢測qrs波的結果進行qrs波起始點和終止點的計算，對所述延時信號進行qrs波成分消除；

s203、根據(jù)所述延時信號和qrs波起始點和終止點的計算結果，獲取當前基線值。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測方法中，所述步驟s202中的步驟“對所述延時信號進行qrs波成分消除”，具體包括：對qrs波的起始點到終止點部分采用插值方法消除，如公式(1)所示：

其中xoff和xon分別表示終止點和起始點的信號值，xoff和xon分別表示終止點和起始點的時間值，yn和xn分別表示當前需要插值的信號值和時間值。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測方法中，所述步驟s203中的所述當前基線值用于檢測基線噪聲的干擾情況。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測方法中，所述步驟s300中高通濾波處理的作用是濾除低頻成分，所述低頻成分包括基線、p波和t波；進行高通濾波處理后，分別進行如下處理：

s301、進行工頻噪聲干擾檢測后輸出工頻噪聲檢測結果；

s302、進行工頻濾波后將結果進行肌電噪聲檢測，并輸出肌電噪聲檢測結果。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測方法中，所述步驟s302中， 采用基于零極點對消的方式進行所述工頻濾波，其中的工頻頻率為60hz，采樣率為500hz，使用的轉移函數(shù)如公式(2)所示：

本發(fā)明還提供了一種全數(shù)字心電信號的噪聲檢測系統(tǒng)，包括：

檢測qrs波模塊，檢測數(shù)字心電信號中的qrs波；

計算模塊，根據(jù)所述檢測qrs波模塊的檢測結果進行qrs波起始點和終止點的計算，將結果分別送給消除qrs波模塊和基線值模塊；消除qrs波模塊，消除數(shù)字心電信號中的qrs波；基線值模塊，獲取當前基線值；

高通濾波模塊，進行高通濾波后將結果分別送給工頻噪聲檢測模塊和工頻濾波模塊；

工頻噪聲檢測模塊，進行工頻噪聲干擾檢測后輸出工頻噪聲檢測結果；

工頻濾波模塊，進行工頻濾波后將結果送給肌電噪聲檢測模塊輸出肌電噪聲檢測結果。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測系統(tǒng)中，還包括：

信號延時模塊，進行信號延時；消除qrs波模塊，接收由信號延時模塊和計算模塊的信息，消除數(shù)字心電信號中的qrs波；基線值模塊，接收由信號延時模塊和計算模塊的信息，獲取當前基線值。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測系統(tǒng)中，還包括：

基線干擾檢測模塊，根據(jù)所述基線值模塊輸出的所述當前基線值，進行基線干擾檢測，并輸出基線干擾檢測結果。

在本發(fā)明提供的全數(shù)字心電信號的噪聲檢測系統(tǒng)中，根據(jù)所述工頻噪聲檢測結果、肌電噪聲檢測結果、基線干擾檢測結果，控制濾波器的開關，和/或進行采集命令控制，實現(xiàn)心電信號的自動采集。

與現(xiàn)有技術相比，本發(fā)明的有益效果是：在采集過程中根據(jù)噪聲大小的不同，控制濾波器的開關，降低由于濾波器對信號的影響；同時根據(jù)噪聲大小進行采集命令控制，實現(xiàn)心電信號的自動采集，減少人工干預。

附圖說明

為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術中的技術方案，下面將對實施例或現(xiàn)有技術描述中所需要使用的附圖作簡單地介紹。顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實施例，對于本領域技術人員來講，在無需付出創(chuàng)造性勞動的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1是本發(fā)明全數(shù)字心電信號的噪聲檢測系統(tǒng)的總體框圖；

圖2是qrs波形檢測流程圖；

圖3是qrs波對消和高通濾波結果；

圖4是高通濾波的幅頻響應；

圖5是工頻噪聲檢測詳細設計圖；

圖6是低通濾波的幅頻響應；

圖7是正弦信號發(fā)生器原理圖；

圖8是工頻濾波幅頻響應；

圖9是肌電噪聲檢測詳細設計圖；

圖10是基線干擾檢測詳細設計圖。

具體實施方式

下面將結合實施例，對本發(fā)明實施例中的技術方案進行清楚、完整地描述。

本發(fā)明的主要創(chuàng)新點在于，提供一種全數(shù)字心電信號的噪聲檢測方法及系統(tǒng)，對心電信號綜合波定位并得到其起始點和終止點，然后消除qrs波，同時得到該心搏的基線值，該基線值用于檢測基線噪聲的干擾情況，同時對消除qrs波后的信號做濾波處理，進行工頻噪聲和肌電噪聲的檢測。

本發(fā)明的實施例和附圖中，frequency表示頻率，單位hz；magnitude表示幅度，單位db，即分貝。

下面結合附圖對本發(fā)明做進一步說明。圖1示出了本發(fā)明全數(shù)字心電信號的噪聲檢測系統(tǒng)的總體框圖，如圖1所示，采集的模擬心電信號經(jīng)“模/數(shù)”轉換后得到數(shù)字心電信號(即為圖1中的ecg信號9)，信號9分別進入信號 延時模塊10(信號延時)和檢測qrs波模塊12(檢測qrs波)，信號延時模塊10保證后面處理的信號同步，而檢測qrs波模塊12進行qrs波的檢測然后將結果送給計算模塊14(計算qrs波起始點和終止點)進行qrs波起始點和終止點的計算，將結果分別送給消除qrs波模塊20(消除qrs波)和基線值模塊22(獲取當前基線值)。基線值模塊22接收由信號延時模塊10和計算模塊14的信息獲取當前基線值，將結果送給基線干擾檢測模塊44(基線干擾檢測)輸出基線干擾檢測結果54。消除qrs波模塊20接收由信號延時模塊10和計算模塊14的信息，將原始ecg信號中的qrs波成分消除，然后送給高通濾波模塊30(高通濾波)，高通濾波模塊30進行高通濾波后將結果分別送給工頻噪聲檢測模塊40(工頻噪聲檢測)和工頻濾波模塊32(工頻濾波)，工頻噪聲檢測模塊40進行工頻噪聲干擾檢測后輸出工頻噪聲檢測結果50；工頻濾波模塊32進行工頻濾波后將結果送給肌電噪聲檢測模塊42(肌電噪聲檢測)輸出肌電噪聲檢測結果52。

圖2為檢測qrs波模塊12的流程圖，采用pan的方法([1]panj,tompkinswj.areal-timeqrsdetectionalgorithm[j].biomedicalengineering,ieeetransactionson,1985(3):230-236)，對心電信號進行濾波后，依次進行差分、平方和移動平均處理，最后尋找峰值，找到峰值后與閾值進行比較以判斷是否找到qrs波并更新閾值，同時將結果13送給計算模塊14計算qrs波的起始點和終止點，得到結果15。

圖3為消除qrs波模塊20和高通濾波模塊30為利用信息15對信號11處理的結果。圖中a為一個受噪聲干擾的心電信號波形，b為消除qrs波模塊20處理的結果，對qrs波的起始點到終止點部分采用插值方法消除，采用線性插值的方法避免對后面的噪聲檢測帶來影響，如公式(1)所示。

其中xoff和xon分別表示終止點和起始點的信號值，xoff和xon分別表示終止點和起始點的時間值，yn和xn分別表示當前需要插值的信號值和時間值。

圖4為高通濾波器的幅度響應，高通濾波器的主要作用是濾除基線、p波 和t波等低頻成分，采用切比雪夫ii型的iir進行濾波器設計；圖3中c為高通濾波的結果。

高通濾波后，信號31進入工頻噪聲檢測模塊40進行工頻噪聲的檢測。圖5為工頻噪聲檢測模塊40的詳細設計圖，工頻噪聲檢測模塊402正弦信號發(fā)生器，分別發(fā)出50/60hz的正弦和余弦信號(正弦信號發(fā)生器發(fā)出的信號頻率與ecg信號9的工頻信號一致)，其中信號400為余弦信號，信號401為正弦信號。信號31分別與信號400和信號401做乘積后分別通過低通濾波器403和低通濾波器404，兩個濾波器采用相同的設計方法，均為截止頻率為0.2的iir(infiniteimpulseresponse，無限脈沖響應)低通濾波器，其幅頻響應如圖6所示。低通濾波后，信號分別進入工頻噪聲檢測模塊405和工頻噪聲檢測模塊406進行平方后進行求和，將結果送往工頻噪聲檢測模塊407與閾值進行比較，將比較結果送往工頻噪聲檢測模塊408(工頻噪聲狀態(tài))，若求和值大于閾值，則工頻噪聲檢測模塊408的輸出為工頻噪聲存在，否則輸出為工頻噪聲不存在。

圖7為工頻噪聲檢測模塊402正弦信號發(fā)生器的實現(xiàn)形式，采用一個類似iir的方法，模塊450、模塊456和模塊458為常量系數(shù)，表示輸入信號經(jīng)過模塊后會與模塊系數(shù)相乘后輸出；模塊452和模塊454表示信號延時一個單位。

圖8為工頻濾波模塊32的幅頻響應圖，采用基于零極點對消的方式設計工頻濾波器，其中的工頻頻率為60hz，采樣率為500hz，濾波器的轉移函數(shù)如公式(2)所示。

信號31經(jīng)工頻濾波模塊32工頻濾波后，將結果信號33送往肌電噪聲檢測模塊42進行肌電噪聲檢測。圖9為肌電噪聲檢測模塊42的詳細設計圖。輸入信號33首先經(jīng)過肌電噪聲檢測模塊420進行絕對值運算，然后將結果460送往肌電噪聲檢測模塊421進行樣本延時操作，肌電噪聲檢測模塊421的樣本延時長度為50ms并輸出結果461，信號461與460進入肌電噪聲檢測模塊 422進行求和操作輸出結果462，信號462與463進入肌電噪聲檢測模塊423進行求和操作后輸出結果464，其中信號463為信號464進過延時肌電噪聲檢測模塊424進行一個采樣點延時后的結果。信號464進入肌電噪聲檢測模塊425與設定閾值進行比較，將結果送給肌電噪聲檢測模塊426，肌電噪聲檢測模塊426根據(jù)肌電噪聲檢測模塊425的輸出最后的肌電噪聲狀態(tài)52，若信號464大于閾值，則狀態(tài)52為肌電噪聲存在，否則為肌電噪聲不存在。

圖10為基線干擾檢測模塊44的詳細設計圖。在基線值模塊22獲取當前基線值信號35后，將結果送往基線干擾檢測模塊440，基線干擾檢測模塊440判斷零電位是否初始化完成，若未進行初始化，則進入基線干擾檢測模塊442進行零電位的初始化，然后將結果送往基線干擾檢測模塊445，基線干擾檢測模塊442在進行零電位初始化后，將零電位初始化狀態(tài)置為真，使得下次進入基線干擾檢測模塊440的結果判斷為真；若已經(jīng)初始化完成，則進入基線干擾檢測模塊441，將當前基線值與零電位做差后取絕對值，然后將結果與閾值比較，若小于閾值，則進入基線干擾檢測模塊444進行零電位更新，然后將結果送往基線干擾檢測模塊445，若基線干擾檢測模塊443比較結果大于閾值，直接將結果送往基線干擾檢測模塊445?；€干擾檢測模塊445根據(jù)輸入的信號進行基線干擾狀態(tài)的判斷并輸出結果54，若由基線干擾檢測模塊442或基線干擾檢測模塊444輸入的結果，則狀態(tài)54為基線干擾不存在；若由基線干擾檢測模塊443經(jīng)比較后直接輸入基線干擾檢測模塊445，則狀態(tài)54為基線干擾存在。

本發(fā)明的全數(shù)字心電信號的噪聲檢測系統(tǒng)，根據(jù)工頻噪聲檢測結果(50)、肌電噪聲檢測結果(52)、基線干擾檢測結果(54)，即根據(jù)采集過程中根據(jù)噪聲大小的不同，控制濾波器的開關，降低由于濾波器對信號的影響；同時根據(jù)噪聲大小進行采集命令控制，實現(xiàn)心電信號的自動采集，減少人工干預。

需要特殊說明的是，以上技術方案僅用于說明本發(fā)明而不用于限制本發(fā)明的范圍。此外，在閱讀了本發(fā)明的內容之后，本領域技術人員可以對本發(fā)明做出改動或修改，這些等價形式同樣在本申請所附權利要求書所限定的范圍之內。