《一种基于心电信号的心率计算方法.pdf》由会员分享，可在线阅读，更多相关《一种基于心电信号的心率计算方法.pdf（6页珍藏版）》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201710495246.4 (22)申请日 2017.06.26 (71)申请人 四川长虹电器股份有限公司 地址 621000 四川省绵阳市高新区绵兴东 路35号 (72)发明人 周杨 张雅勤 刘传银 张长安 朱文涛 (74)专利代理机构 四川省成都市天策商标专利 事务所 51213 代理人 蒋金梅 李洁 (51)Int.Cl. A61B 5/0245(2006.01) A61B 5/0456(2006.01) (54)发明名称 一种基于心电信号的心率计算方法 (57)摘要 本。
2、发明公开了一种基于心电信号的心率计 算方法， 涉及医学电子信息领域。 本发明提供的 方法包括以下步骤： 第一步： 获取心电信号的差 分绝对值信号； 第二步： 搜索峰值点； 第三步： 定 位R波； 第四步： 计算心率。 该方法从采集的一段 连续心电信号序列中提取特征点， R波顶点， 并计 算出相邻特征点间的时间间隔， 再计算这些时间 间隔的均值， 利用该均值即可换算得到一定时间 长度内的平均心率。 该方法可以快速准确地从心 电信号中自动计算得到心率值。 权利要求书1页 说明书3页 附图1页 CN 107095665 A 2017.08.29 CN 107095665 A 1.一种基于心电信号的心。
3、率计算方法， 其特征在于， 包括以下步骤： 第一步： 获取心电信号的差分绝对值信号 给定一段时长L秒的肢体导联I导联心电信号x1(n)， 其中n1， 1、 2、 、 f*L， f为信号采 样频率； 计算其差分信号d_x1(n)x1(n+k)x1(n)， 其中n1、 2、 、 f*Lk， k为小于n的一 个正整数； 对差分信号d_x1(n)中各个元素取绝对值， 得到差分绝对值信号abs_d_x1(n)； 第二步： 搜索峰值点 对差分绝对值信号abs_d_x1从起始位置点开始搜索， 当遇到幅值大于thres1的点p1， 则在区间p1， p2上搜索幅值大于thres1的点的数量cnt1， 其中p2、。
4、 p1为点的索引值， 且p2 与p1的时长间距为0.1秒， 如果cnt1/(p2p1)大于阈值thres2， thres20.5， 则将p1,p2 区间上幅值最大点定义为一个峰值点； 然后， 从p2+1点继续采用上述方法向前搜索， 如此循 环直到差分绝对值信号abs_d_x1结束位置时停止搜索； 此过程将获得若干个峰值点的集合 Aa1,a2,， a1、 a2分别表示各峰值点； 第三步： 定位R波 在心电信号x1(n)上以集合A中各峰值点为中心的邻域， 即区间aid， ai+d上搜索最 大幅值点， 这些点即为R波顶点； 其中ai表示集合Aa1,a2,中任意一峰值点的相应索 引值， d为正整数； 。
5、第四步： 计算心率 计算相邻两R波顶点的时间间隔T， 由T可以换算得到一个瞬时心率； 对所有瞬时心率求 平均， 即可得到该时长L秒内的平均心率。 2.根据权利要求1所述的一种基于心电信号的心率计算方法， 其特征在于： 第一步中， L 2。 3.根据权利要求1所述的一种基于心电信号的心率计算方法， 其特征在于： 第二步中， 令阈值thres14*m， 其中m为差分绝对值信号abs_d_x1的均值。 4.根据权利要求1所述的一种基于心电信号的心率计算方法， 其特征在于： 第三步中， d 为0.1*f向下取整的结果， f为信号采样频率。 5.根据权利要求1所述的一种基于心电信号的心率计算方法， 其特。
6、征在于： 第四步中， 瞬时心率为60/T。 6.根据权利要求1所述的一种基于心电信号的心率计算方法， 其特征在于： 第一步中， k 为0.01*f向下取整的结果， f为信号采样频率。 权 利 要 求 书 1/1 页 2 CN 107095665 A 2 一种基于心电信号的心率计算方法 技术领域 0001 本发明涉及医学电子信息领域， 具体涉及一种基于心电信号的心率计算方法。 背景技术 0002 心电图是临床最常用的检查之一， 应用广泛， 包括帮助诊断心律失常、 心肌缺血、 心肌梗死等。 心电图记录的是随心动周期变化的体表特定位置的电位差。 心电信号是典型 的人体电信号， 人体电信号本质是两点的。
7、电位差信号。 直接加电压于身体并通过一定的导 联方式就可以观察到心电信号。 导联方式即输入导线与电极放置在机体特定的测试部位 (正输入端)、 参比部位(负输入端)和接地部位的连接方式。 在心电图学中， 双极肢体导联又 称标准I、 II、 III导联， 它是以两肢体间的电位差作为所获取的体表心电。 由于双极肢体导 联系统获取的心电图中有较大的R波， 因此利用双极肢体导联就能够可靠地识别每次心跳 并进行节律分析。 0003 如果能够从某一导联的心电信号中准确识别出所有R波， 那么根据相邻R波的时间 间隔就可以很容易得到心率值。 显然利用双极肢体导联(如I导联)获取的心电信号R波幅值 较大的特点， 。
8、可以方便快捷识别其中的R波。 0004 由于心电信号容易因呼吸、 肢体运动等产生基线漂移， 以及易受肌电信号干扰， 因 此现有技术直接在原始心电信号上自动准确寻找R波有较大难度。 发明内容 0005 本发明的目的是针对上述现有技术中不足， 提出基于心电信号的差分信号搜寻R 波的心率计算方法， 用于解决上述问题。 0006 为了达到上述的技术效果， 本发明采取以下技术方案： 一种基于心电信号的心率 计算方法， 包括以下步骤： 0007 第一步： 获取心电信号的差分绝对值信号 0008 给定一段时长L秒的肢体导联I导联心电信号x1(n)， 其中n1， 1、 2、 、 f*L， f为信 号采样频率；。
9、 计算其差分信号d_x1(n)x1(n+k)x1(n)， 其中n1、 2、 、 f*Lk， k为小于n 的一个正整数； 对差分信号d_x1(n)中各个元素取绝对值， 得到差分绝对值信号abs_d_x1 (n)； 0009 第二步： 搜索峰值点 0010 对差分绝对值信号abs_d_x1从起始位置点开始搜索， 当遇到幅值大于thres1的点 p1， 则在区间p1， p2上搜索幅值大于thres1的点的数量cnt1， 其中p2、 p1为点的索引值， 且 p2与p1的时长间距为0.1秒， 如果cnt1/(p2p1)大于阈值thres2， thres20.5， 则将p1, p2区间上幅值最大点定义为一。
10、个峰值点； 然后， 从p2+1点继续采用上述方法向前搜索， 如 此循环直到差分绝对值信号abs_d_x1结束位置时停止搜索； 此过程将获得若干个峰值点的 集合Aa1,a2,， a1、 a2分别表示各峰值点； 0011 第三步： 定位R波 说 明 书 1/3 页 3 CN 107095665 A 3 0012 在心电信号x1(n)上以集合A中各峰值点为中心的邻域， 即区间aid， ai+d上搜 索最大幅值点， 这些点即为R波顶点； 其中ai表示集合Aa1,a2,中任意一峰值点的相 应索引值， d为正整数； 0013 第四步： 计算心率 0014 计算相邻两R波顶点的时间间隔T， 由T可以换算得到。
11、一个瞬时心率； 对所有瞬时心 率求平均， 即可得到该时长L秒内的平均心率； 其中T的单位为秒。 0015 进一步的技术方案是： 在第一步中， L2。 0016 进一步的技术方案是： 第二步中， 令阈值thres14*m， 其中m为差分绝对值信号 abs_d_x1的均值。 0017 进一步的技术方案是： 第三步中， d为0.1*f向下取整的结果， f为信号采样频率。 0018 进一步的技术方案是： 第四步中， 瞬时心率为60/T。 0019 进一步的技术方案是： 第一步中， k为0.01*f向下取整的结果， f为信号采样频率。 0020 本发明与现有技术相比， 具有以下的有益效果： 本发明提供的。
12、方法从采集的一段 连续心电信号序列中提取特征点(如R波顶点)， 并计算出相邻特征点间的时间间隔， 再计算 这些时间间隔的均值， 利用该均值即可换算得到一定时间长度内的平均心率。 该方法可以 快速准确地从心电信号中自动计算得到心率值。 附图说明 0021 图1是本发明实施例心电信号x1(n)的示意图； 0022 图2是本发明实施例差分绝对值信号abs_d_x1(n)的示意图。 具体实施方式 0023 下面结合本发明的实施例对本发明作进一步的阐述和说明。 0024 实施例： 0025 一种基于心电信号的心率计算方法， 第一步是获取心电信号的差分绝对值信号。 如图1所示， 已知一段长度L3秒的人体肢。
13、体导联I导联心电信号x1(n)， n1、 2、 、 f*L， 其 中信号采样频率f512。 其差分信号d_x1(n)x(n+k)x(n)， n1、 2、 、 f*Lk， k为0.01*f 向下取整的结果即5。 对d_x1(n)中各个元素取绝对， 得到差分绝对值信号abs_d_x1(n)， 如 图2所示。 0026 第二步， 搜索峰值点。 首先计算abs_d_x1(n)的均值m为161， 取阈值thres14*m， 即644。 阈值thres2取0.5。 其次， 如图2所示， 对信号abs_d_x1(n)从起始位置点abs_d_x1(1) 开始搜索， 当遇到幅值大于thres1的点p1时， 则在。
14、区间p1,p2上搜索幅值大于thres1的点 的数量cnt1， 其中p2、 p1为点的索引， 且p2与p1的时长间距为0.1秒， 如果cnt1/(p2p1)大于 阈值thres2， 则将p1,p2区间上幅值最大点定义为一个峰值点。 然后， 从p2+1点继续采用 上述方法向前搜索， 如此循环直到信号abs_d_x1(n)结束位置时停止搜索。 此过程将获得若 干个峰值点的集合Aa1,a2,， a1、 a2等表示各个峰值点。 0027 第三步， 定位R波。 在心电信号x1(n)上的以集合A中各峰值点为中心的邻域即区间 aid,ai+d上搜索最大幅值点， 这些点即为R波顶点， 如图1所示。 其中， d。
15、为0.1*f向下取整 的结果即51， ai表示集合Aa1,a2,中任意一峰值点的相应索引值。 说 明 书 2/3 页 4 CN 107095665 A 4 0028 第四步， 计算心率。 计算相邻两R波顶点的时间间隔T， 即0.84和0.86， 依据瞬时心 率r60/T， 可得对应的瞬时心率71.4和70.6； 对所有瞬时心率求平均， 即可得到该时长L内 的平均心率即70。 0029 可以理解的是， 以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施 方式， 然而本发明并不局限于此。 对于本领域内的普通技术人员而言， 在不脱离本发明的精 神和实质的情况下， 可以做出各种变型和改进， 这些变型和改进也视为本发明的保护范围。 说 明 书 3/3 页 5 CN 107095665 A 5 图1 图2 说 明 书 附 图 1/1 页 6 CN 107095665 A 6 。