[發明專利]心電信號R波的提取方法

有效

• 申請號: 201310633841.1
• 申請日: 2013-11-29
• 公開（公告）號: CN103584854A
• 公開（公告）日: 2014-02-19
• 發明（設計）人: 龐宇;李章勇;李捷;吳強;趙何婷;任斌斌;李澤穎
• 申請（專利權）人: 重慶海睿科技有限公司
• 主分類號: A61B5/0456
• 分類號: A61B5/0456
• 代理公司: 重慶為信知識產權代理事務所(普通合伙) 50216
• 代理人: 余錦曦

技術領域

本發明涉及到醫療器械技術領域，具體地說，是一種心電信號R波的提取方法。?

背景技術

QRS波是心電圖ECG（Eletrocardiogram）信號的主要特征，因此在Indian信號最重要的是檢測出QRS波。這是因為只有在確定QRS波群后下才能計算心率、檢測心率變異，并進一步檢測心電信號的其他細節。目前R波的檢測有以下幾種檢測方法：差分法、濾波器組法、小波變換法等。差分法具有直觀形象計算量小的優點，但是抗干擾能力弱；濾波器組法具有運算量小，抗干擾能力較強的優點，但是時常出現錯檢和漏檢情況；小波變換法對R波的提取效果很好，但是運算大，運算復雜。?

因此，需要一種能夠克服抗干擾能力弱、準確率低、運算量大、運算復雜等缺陷的提取方法。?

發明內容

針對現有技術的不足，本發明的目的之一是提供一種抗干擾能力強、準確率高、運算量小的心電信號R波的提取方法。?

為達到上述目的，本發明表述一種心電信號R波的提取方法，具體按照以下步驟進行：?

步驟1：獲取預定時間T內的心電信號，并進行低通濾波預處理；?

步驟2：對濾波后的心電信號提取第一個R波，并將第一個R波所在心電信號段作為初始模板；?

步驟3：對獲取的心電信號進行雙閾值比較和粗細匹配計算，提取預定時間T內其它R波。?

作為進一步描述，所述步驟1中以頻率f=400Hz進行采樣，將采樣獲取的心電信號按照截止頻率100Hz進行低通濾波；所述步驟2中獲取第一個R波的具體方法按照以下步驟進行：?

步驟2-1：從起始處取寬度W1=2s的心電信號段；?

步驟2-2：按照等式（a）對心電信號做一階前向差分運算，得到差分信號dif(k)，等式中f(k)為濾波后的心電信號，f(k+1)是f(k)的移位序列；?

dif(k)＝f(k)-f(k+1)???（a）?

步驟2-3：找到該信號段的前向差分的最大值dif(k)_max所在的位置，記為K0點；?

步驟2-4：以K0點為中心，取寬度W2=100ms的心電信號段，找到該信號段中的最大值點K1，即為第一個R波所在位置。?

作為進一步描述，所述步驟3中的雙閾值比較和粗細匹配具體按照以下步驟進行：?

步驟3-1：以上一個R波所在點的位置為中心，選取寬度W3=150ms的心電信號段作為模板信號；?

步驟3-2：設置t0的初值并滑動t0，取寬度為W3=150ms的心電?信號段作為匹配信號；?

步驟3-3：更新幅度閾值；?

步驟3-4：將匹配信號段中的最大值與幅度閾值比較，小于閾值的則進入步驟3-5作粗糙匹配，大于閾值的則進入步驟3-6作精細匹配；?

步驟3-5：按照等式（b）計算粗糙匹配度Q1，并將得到的粗糙匹配度Q1與閾值c1比較，若Q1大于c1的，則返回步驟3-3重新更新幅度閾值，若Q1小于c1的，則修改t0=t1，返回步驟3-2；?

Q1=(1-|Mk-ki|Mk)*(1-|Mc-ci|Mc)---(b)]]>