《基于智能诊断的远程无线心电监护系统及特征提取方法.pdf》由会员分享，可在线阅读，更多相关《基于智能诊断的远程无线心电监护系统及特征提取方法.pdf（17页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103371814 A (43)申请公布日 2013.10.30 CN 103371814 A *CN103371814A* (21)申请号 201210108824.1 (22)申请日 2012.04.14 A61B 5/0402(2006.01) A61B 5/0205(2006.01) A61B 5/0456(2006.01) A61B 5/0452(2006.01) (71)申请人 兰州大学 地址 730000 甘肃省兰州市天水南路 222 号 (72)发明人 马义德 张少华 田毅 刘莉 (74)专利代理机构 兰州振华专利代理有限责任 公司 62102 代理人。
2、 董斌 (54) 发明名称 基于智能诊断的远程无线心电监护系统及特 征提取方法 (57) 摘要 基于智能诊断的远程无线心电监护系统及特 征提取方法， 由心电采集子系统、 血压采集子系 统、 智能诊断单元、 GPS 定位子系统、 人机交互子 系统、 无线传输子系统、 存储器单元与系统处理器 单元组成 ； 用户佩戴本系统后， 可以对其心电信 号进行实时监测并对心电波形实时显示， 当检测 到心脏工作异常时， 系统处理器将控制血压采集 子系统和 GPS 定位子系统获取用户当前血压值和 位置坐标， 并利用无线传输子系统将以上信息发 送给与用户相关联的各监护终端上 ； 特征提取方 法是采用基于阈值和奇点检。
3、测相结合的方法来实 现对心电信号中 R 波和 S-T 段的定位以及其他心 电特征的提取。 (51)Int.Cl. 权利要求书 2 页 说明书 7 页 附图 7 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书7页 附图7页 (10)申请公布号 CN 103371814 A CN 103371814 A *CN103371814A* 1/2 页 2 1. 基于智能诊断的远程无线心电监护系统， 包括 ： 心电采集子系统、 血压采集子系统、 智能诊断单元、 GPS 定位子系统、 人机交互子系统、 无线传输子系统、 存储器单元以及系统处 理器单元， 其特征是将智能。
4、诊断功能和无线传输功能相结合应用于心电监护设备中， 其中 ： 电极片采集到的心电信号经 A/D 转换模块转换为数字信号后经过预处理和特征提取后交 给智能诊断单元， 同时血压采集子系统将采集到的血压信号送给处理器作为心电特征诊断 的辅助信息， 智能诊断单元得出的最终诊断结果交给处理器进行存储， 并在液晶显示屏上 实时显示心电波形和诊断结果， GPS 定位子系统定位患者当前的精确位置， 处理器将诊断结 果和位置信息通过无线传输子系统发送给监测端。 2. 根据权利要求 1 所述的基于智能诊断的远程无线心电监护系统， 其特征在于心电采 集子系统由心电信号调理电路、 小波滤波器和心电特征提取模块构成， 。
5、其中小波滤波器采 用 coif5 小波滤波。 3. 根据权利要求 1 所述的基于智能诊断的远程无线心电监护系统， 其特征是 GPS 定位 子系统由串口接口芯片、 GPS 模块、 GPS 外置天线组成。 4. 根据权利要求 1 所述的基于智能诊断的远程无线心电监护系统， 其特征是智能诊断 单元中的特征提取功能单元， 用于检测 R 波、 平均心率 Rate 与 HRV 间期。 5. 根据权利要求 1 所述的基于智能诊断的远程无线心电监护系统， 其特征在于无线传 输子系统包括串口转换芯片、 GSM/CDMA 调制解调器、 外置天线， 利用无线传输子系统可以将 本设备以移动终端的形式挂载到公共网络上。。
6、 6. 基于智能诊断的远程无线心电监护系统的特征提取方法， 其特征是使用 R 波检测算 法， 采用基于阈值和奇点检测相结合的方法来实现对心电信号中 R 波和 S-T 段的定位以及 其他心电特征的提取， 其步骤为 ： (1) 设定一个时间窗， 在该时间段内采用二次差分方法， 寻找波形奇异点 ； 对长度为 N( 取 N 4096) 的数据进行奇点检测 ； 奇点检测采用的算法为对一个采样周期的心电信号数据， 按照公式 (1) 至 (3) 进行 diff(sign(diff(N) 运算 ： 其中diff为信号差分， sign为符号函数， N为一个采样周期的心电信号， 运算结果为-2 的点即为心电信号的。
7、极大值点 ； 保留得到的所有极大值在一个数组 A 内， 并记录其对应波形位置 ； 取阈值Rth， 与上述数组A的值比较， 保留大于Rth的极大值在数组B内， 并记录其对应 权 利 要 求 书 CN 103371814 A 2 2/2 页 3 波形位置 ； 认为 B 中数据为检测到的 R 波幅值， 其对应位置为 R 波位置 ； (2) 自适应阈值 Rth 设置 ： 选取的N点ECG滤波之后的数据， 做奇异点检测之后， 寻找奇异点中的极大值信号幅值 范围， 按公式 (4)， 均分为 15 份， 并将每份幅值记录在 Th 数组中 ； 设积分投影函数为 ： 其中 ： I i 15 且 i Z (6) 。
8、即做极大值点在划分的各幅值段 Th(i) 上的积分投影 ； 选取零分布与非零分布的交接点处 i 值， 计算 Th(i)， 确定 Rth 为 Th(i) ； 阈值确定出 R 波位置和幅值之后， 计算平均心率 Rate， 单位 ： 次 / 分， Rate 60*Nr/(nr(end)-nr(1)/fs) (8) 其中， Nr 是固定时间窗中测得的 R 波个数， nr(end) 是固定时间窗内最后一个 R 波的位 置， nr(1) 是固定时间窗内第一个 R 波的位置， fs是数据采样率 ； 根据平均心率加入相关误检、 漏检更正， 最终修正出正确 R 波位置和平均心率 Rate ； 然后用公式 (9)。
9、， 计算 HRV 间期。 HRV RR(i+1)-RR(i) (9) 其中 RR 为相邻两个 R 波间期。 综上算法， 可得到后续诊断所需特征参数 ： Rate 和 H RV ； (3) 根据此指标通过查表的方式将提取到的心电特征与 MIT-BI H 数据库分析总结得 出诊断结果。 权 利 要 求 书 CN 103371814 A 3 1/7 页 4 基于智能诊断的远程无线心电监护系统及特征提取方法 技术领域 0001 本发明涉及基于智能诊断的远程无线心电监护系统。 背景技术 0002 目前心脏病诊断主要以心电图和心脏彩超为主， 心电图主要应用于心律失常、 心 肌缺血等功能性心脏病诊断 ； 心。
10、脏彩超则对风湿性心脏病、 先天性心脏病等结构异常的心 脏病才有用。 由于功能性疾病往往具有病情隐蔽、 发病突然、 危险性高、 持续时间短等特征， 因此很难把握最佳的诊断时间， 发病时病人往往也很难进行有效的自救措施。 因此， 病人的 心电特征对于突发性心脏病的诊断和治疗有非常重大的指导作用。 0003 目前医疗机构较常用的心电采集设备主要包括三种 ： 1、 心电监护仪， 这类设备可 以检测患者的心搏频率、 呼吸、 血压、 脉搏、 血氧饱和度等生命特征， 在一些重症病房非常常 见。设备采集到病人的以上生理指标， 在屏幕上实时显示， 对术后或者危重病人的看护有 重要意义。但这类设备采集的心电信号多。
11、数是单导联信号， 往往只能反应病人心脏是否在 正常跳动， 对于心脏病的诊断， 没有太多的参考意义。2、 心电图机， 目前较常用的导联数是 十二导联， 这类心电采集设备可以采集较完整的心电信号， 能够比较全面的反应多种心脏 功能性异常。 但这类设备一般结构庞大， 每次只能记录很短时间的心电信号， 需要专业人员 操作， 因此只有在医院等大型的医疗机构才能采用这种设备进行检测。很多心脏病具有发 病突然、 持续时间短等特点， 往往会出现病人感觉到异常时前往医院进行诊断， 待到达医院 时病情得到缓解， 在心电图机上不能发现明显的病变特征的情况， 长期住院观察又将花费 大量的人力财力， 很多患者难以承受。。
12、 因此， 这类心电采集设备在实际应用中还存在很多弊 端。3、 动态心电采集设备 (Holter)， 这类设备主要是为了解决普通心电图机难以完成全天 候监测的问题设计的， 这类设备体积小， 便于随身佩戴， 主要包括心电数据的采集和存储两 部分， 大容量的存储设备将 24 小时的心电信号进行存储， 然后将这些信号传送给计算机进 行初步筛选， 将可能存在异常的部分送给医生进行经验判断。这类设备在一定程度上解决 了普通心电图机存在的问题， 但仍存在很多亟待解决的问题。 首先是这类设备非常昂贵， 需 要长时间佩戴的患者很多难以承受其巨大的花费。其次这类设备智能化水平较低， 需要经 过多道工序才能筛选出最。
13、终用来诊断需要的病变信号， 耗费大量的人力、 物力。 而且这种设 备对于一些急性、 突发使人失去自救能力的致命发病情况也束手无策， 不能有效的缩短救 援所需时间。因此还不能完全满足心脏病患者监护的全部需求。 发明内容 0004 本发明的目的是提供一种基于智能诊断的远程无线心电监护系统及特征提取方 法。 0005 本发明是基于智能诊断的远程无线心电监护系统及特征提取方法， 其远程无线心 电监护系统， 包括心电采集子系统、 血压采集子系统、 智能诊断单元、 GPS 定位子系统、 人机 交互子系统、 无线传输子系统、 存储器单元以及系统处理器单元， 其特征是将智能诊断功能 说 明 书 CN 1033。
14、71814 A 4 2/7 页 5 和无线传输功能相结合应用于心电监护设备中， 其中 ： 电极片采集到的心电信号经 A/D 转 换模块转换为数字信号后经过预处理和特征提取后交给智能诊断单元， 同时血压采集子系 统将采集到的血压信号送给处理器作为心电特征诊断的辅助信息， 智能诊断单元得出的最 终诊断结果交给处理器进行存储， 并在液晶显示屏上实时显示心电波形和诊断结果， GPS 定 位子系统定位患者当前的精确位置， 处理器将诊断结果和位置信息通过无线传输子系统发 送给监测端。 0006 基于智能诊断的远程无线心电监护系统的特征提取方法， 其特征是使用 R 波检测 算法， 采用基于阈值和奇点检测相结。
15、合的方法来实现对心电信号中 R 波和 S-T 段的定位以 及其他心电特征的提取， 其步骤为 ： 0007 (1) 设定一个时间窗， 在该时间段内采用二次差分方法， 寻找波形奇异点 ； 对长度 为 N( 取 N 4096) 的数据进行奇点检测 ； 0008 奇点检测采用的算法为对一个采样周期的心电信号数据， 按照公式 (1) 至 (3) 进 行 diff(sign(diff(N) 运算 ： 0009 0010 0011 0012 其中diff为信号差分， sign为符号函数， N为一个采样周期的心电信号， 运算结果 为 -2 的点即为心电信号的极大值点 ； 0013 保留得到的所有极大值在一个数。
16、组 A 内， 并记录其对应波形位置 ； 0014 取阈值Rth， 与上述数组A的值比较， 保留大于Rth的极大值在数组B内， 并记录其 对应波形位置 ； 0015 认为 B 中数据为检测到的 R 波幅值， 其对应位置为 R 波位置 ； 0016 (2) 自适应阈值 Rth 设置 ： 0017 选取的N点ECG滤波之后的数据， 做奇异点检测之后， 寻找奇异点中的极大值信号 幅值范围， 按公式 (4)， 均分为 15 份， 并将每份幅值记录在 Th 数组中 ； 0018 0019 设积分投影函数为 ： 0020 说 明 书 CN 103371814 A 5 3/7 页 6 0021 其中 ： 00。
17、22 1 i 15 且 i Z (6) 0023 0024 即做极大值点在划分的各幅值段 Th(i) 上的积分投影 ； 0025 选取零分布与非零分布的交接点处 i 值， 计算 Th(i)， 确定 Rth 为 Th(i) ； 0026 阈值确定出 R 波位置和幅值之后， 计算平均心率 Rate， 单位 ： 次 / 分， 0027 Rate 60*Nr/(nr(end)-nr(1)/fs) (8) 0028 其中， Nr 是固定时间窗中测得的 R 波个数， nr(end) 是固定时间窗内最后一个 R 波 的位置， nr(1) 是固定时间窗内第一个 R 波的位置， fs是数据采样率 ； 0029 。
18、根据平均心率加入相关误检、 漏检更正， 最终修正出正确 R 波位置和平均心率 Rate ； 0030 然后用公式 (9)， 计算 HRV 间期。 0031 HRV RR(i+1)-RR(i) (9) 0032 其中 RR 为相邻两个 R 波间期。 0033 综上算法， 可得到后续诊断所需特征参数 ： Rate 和 H RV ； 0034 (3) 根据此指标通过查表的方式将提取到的心电特征与 MIT-BI H 数据库分析总 结得出诊断结果。 0035 本发明针对现有的心电监护设备在对心脏病患者的诊断和监护过程中存在的问 题提供了一种将智能诊断和远程监护合为一体的高效的、 稳定的监护系统。该系统可。
19、以实 时监测病人的心电、 血压变化情况， 并进行综合的智能诊断， 当有病变特征出现时及时的通 过无线网络向病人家属或医生发送诊断报告， 并可通过 GPS 定位系统实时定位患者所处位 置， 缩短救援所需时间， 减少由于救援不力造成患者死亡的悲剧。 附图说明 0036 图 1 本发明整体系统框图， 图 2(a) 本发明心电采集子系统架构图， 图 2(b) 本发明 心电采集子系统智能诊断部分流程图， 图3(a)本发明血压采集子系统架构图， 图3(b)本发 明血压采集子系统流程图， 图 4(a) 本发明 GPS 定位子系统架构图， 图 4(b) 本发明 GPS 定位 子系统流程图， 图5(a)本发明人。
20、机交互子系统架构图， 图5(b)本发明人机交互子系统流程 图， 图 6(a) 本发明无线传输子系统架构图， 图 6(b) 本发明无线传输子系统流程图， 图中符 号说明 ： 11- 心电采集子系统采集到的心电特征信号， 12- 血压采集子系统采集到的血压特 征信号， 13- 处理器对心电采集子系统的控制信号， 14- 处理器对血压采集子系统的控制信 号， 15- 综合心电、 血压等生理参数智能诊断的结果， 16-GPS 定位子系统返回给处理器的方 位信息， 17- 处理器发送给液晶屏的刷屏信息， 18- 触屏驱动器返回给处理器的触屏按键位 说 明 书 CN 103371814 A 6 4/7 页。
21、 7 置， 19- 处理器发送给无线模块的指令和诊断信息。 具体实施方式 0037 如图 1 所示， 本发明的基于智能诊断的远程无线心电监护系统包括 ： 心电采集子 系统、 血压采集子系统、 智能诊断单元、 GPS 定位子系统、 人机交互子系统、 无线传输子系统、 存储器单元、 系统处理器单元。通过以上八个子系统来完成对病人的 24 小时智能监控。处 理器通过心电采集子系统的控制信号13实现导联的切换和A/D转换功能的实现 ； 处理器通 过血压采集子系统的控制信号 14 来实现气泵的充气、 放气以及血压值的读取 ； 心电采集子 系统提取到的心电特征信号 11 送给处理器中的智能诊断单元做进一步。
22、分析 ； 血压采集子 系统采集到的血压值信号 12 送给处理器作为心电特征诊断的辅助信息， 智能诊断单元得 出的最终诊断结果 15 由处理器进行存储 ； GPS 定位子系统通过串口以中断的形式向处理器 发送时间、 坐标、 速度等信息 16， 处理器根据格式匹配的方式从中提取出所需的坐标信息 ； 人机交互子系统有触屏和显示屏组成， 通过信号 17、 18 实现了用户和处理器之间的各种交 互需求和功能显示 ； 无线传输子系统与处理器之间通过串口方式通信， 处理器得到的诊断 信息和位置坐标 19 通过无线模块发送给监护端。 0038 如图 2(a) 所示， 心电采集子系统主要包括心电电极、 导联控制。
23、器和信号预处理部 分 ； 心电电极片与导联控制器相连， 导联控制器的输出经由放大、 滤波电路后进行 A/D 转 换， A/D 转换模块的输出端连接小波滤波电路， 小波滤波电路的输出结果送给智能诊断单 元。 0039 本实施方案中导联控制器由模拟开关 MAX397 实现， 放大电路和滤波电路由运算 放大器 AD620 和 TLV2254 搭建， A/D 转换芯片采用 TLC1549， 处理器采用 Altera Cyclonell 2C70FPGA， coif5 小波在 FPGA 上以硬件电路的方式设计实现。 0040 该心电采集子系统通过电极片从人体表面采集心电信号， 导联控制器以时分复用 的方。
24、式分时段连通不同的导联， 实现了十二导联心电信号的同时采集。采集到的心电信号 非常微弱， 通常在 0.5 2mv 之间， 该差分信号叠加在由电极和皮肤接触所产生的 300mv 左右的直流电压分量之上， 其中还包含着由电极和地之间电势所产生的 1.5v 的共模电压， 因此本发明采用了仪用放大器 AD620 搭建放大电路， 该运放在接近 1kHz 共模抑制比高达 100dB， 50uV 的最大输入失调电压， 1nA 的最大输入偏置电流， 性能能够满足这里心电信号 发达的要求。由于心电信号的频率大约在 0.05 100Hz 之间， 所以需要一个低通滤波器 来消除高于 100Hz 的高频噪声， 一个高。
25、通滤波器用于滤除低于 0.05Hz 的低频噪声， 同时还 需要一个陷波电路用于消除 50Hz 的工频干扰。滤波器的截止频率为 f 1/2RC， 分别将 0.05Hz、 50Hz、 100Hz 带入上式， 可得出各滤波器的参数， 滤波器利用运放 TLV2254 搭建实 现。本设计采用了 TLC 1549 作为模数转换器， 该芯片是一款串行 10 位 A/D 转换芯片， 数据 端口少， 而且可以满足本发明对数据精度的要求。从身体表面采集到的心电信号经过以上 预处理后以数字信号形式交给 FPGA 中的小波滤波电路， 本发明采用 coif5 滤波器， 经过三 层分解、 重构， 该滤波过程在 FPGA 。
26、中的 FIFO( 先进先出 ) 单元完成， 每次对 512 个数据进行 滤波处理。 0041 血压采集子系统实施方案 ： 很多心脏病患者往往伴有血压异常的问题， 因此本发 明添加了血压采集功能。如图 3(a) 所示， 血压采集子系统主要包括气泵控制器、 气泵、 压力 说 明 书 CN 103371814 A 7 5/7 页 8 传感器、 A/D 转换模块。处理器连接气泵控制器， 气泵控制器连接气泵， 气泵控制器可控制 气泵充气、 放气， 并可以通过压力传感器读取袖带中的气压值。压力传感器安装在袖带中， 处理器通过 TLC1549 模数转换器读取压力传感器测得的当前压力值， 通过压力的变化可以 。
27、确定脉搏情况。 0042 如图 3(b) 所示， 当处理器需要采集当前血压时， 气泵控制器控制气泵对袖带充 气， 压力传感器时刻检测当前脉搏情况， 当动脉搏动消失后， 气泵控制器控制气泵缓慢放 气， 当压力传感器再次检测到脉搏时， 记下袖带中的气压值 P1， 气泵继续放气， 当脉搏再次 消失时， 记录袖带中的气压值 P2， 两次得到的气压值经过数据矫正后即可得到最终包含收 缩压和舒张压的血压值。 0043 如图 2(b) 所示， 智能诊断单元是本发明的核心， 该单元包括特征提取和智能诊断 两个主要功能， 该单元通过FPGA中设计的硬件电路实现。 一个标准的心电图主要包含P波、 QRS 波群、 。
28、T 波、 PR 间期、 ST 段、 QT 间期、 HRV 间期等， 心脏病变往往会引起这些波的幅值和 间期的变化， 这是本发明智能诊断的依据。 0044 智能诊断单元中的特征提取功能单元， 主要检测 R 波、 平均心率 Rate 与 HRV 间期。 QRS 波群是一个心电周期特征最显著的部分， 它具有高幅度和高斜率、 具有一定的宽度的时 域波形特征和频谱分布在 ECG 信号的中、 高频区域的频域特征 ( 峰值频率在 10 20Hz)。 0045 如图4(a)所示， GPS模块采用RS232电平逻辑， 因此在与处理器连接前需要经过串 行接口芯片进行电平转换。GPS 定位子系统工作流程参见图 4(。
29、b)， GPS 模块以串口中断形 式向处理器发送包含时间、 坐标、 高度、 速度等的信息， 处理器在中断处理程序中， 通过判断 开始位、 结束位来定位完整的数据包， 并从中提取出坐标信息， 来确定当前患者所处位置。 0046 如图 5(a) 所示， 人机交互子系统主要包括液晶显示和触屏输入功能。液晶显示部 分本发明采用 Y70-4024-65K 彩色液晶模块， 该屏幕是广泛地应用于工业控制等设备上的 彩色 TFT 液晶显示屏， 具有 7 寸超大显示面积的同时兼有 400 240 的分辨率， 大大降低了 系统数据传输和存储的压力， 采用 16 位标准 8080 总线接口方式， 色彩支持 6553。
30、6 色使图像 更加细腻， 独有 2 页显存， 单独操作一页不影响其它页， 便于满足本项目中需要实时部分刷 屏的需求。本发明基于该屏开发了一套友好的人机交互界面， 可以很方便的显示各项检测 信息和功能菜单。触屏显示部分主要包括触屏、 触屏驱动器两部分， 本发明采用 7 寸电阻式 触屏， 在触屏上电后， 按压触屏的不同位置可引起输出电位的相应变化， 根据输出电压的情 况可以精确的定位出用户的触摸位置和功能选项， 本发明采用 ADS7843 作为触屏控制器， ADS7843 是 TI 公司生产的 4 线电阻触摸屏转换接口芯片， 它是一款具有同步串行接口的 12 位取样模数转换器， 通过两次 A/D 。
31、转换可以定位出与之连接的触屏的触摸位置的横坐标和 纵坐标。该子系统可完成设备运行情况的显示， 和用户命令输入， 即人机交互功能。 0047 如图 6(a) 所示， 本发明将心电诊断设备以移动终端的形式挂载在通信网络上， 这样就可以非常方便的借助公共无线网络实现终端和监护端之间的通信。本发明采用 MC323CDMA 模块， 该模块采用 RS232 电平逻辑， 与处理器间通过串口转换芯片连接， 处理器 通过 AT 指令的形式通过串口给无线模块发送指令设计采用短信的方式向监护端发送患者 诊断信息， 短信发送流程参见图 6(b) 示， 处理器对无线模块完成波特率等的初始化后进入 短信发送过程， 该过程。
32、主要括监护端 ( 信息接收端 ) 号码录入、 发送信息录入、 结束符插入 等， 完成以上操作后， 无线模块即可以通过电信 CDMA 网络向监护端发出诊断信息。 说 明 书 CN 103371814 A 8 6/7 页 9 0048 基于智能诊断的远程无线心电监护系统的特征提取方法， 使用 R 波检测算法， 采 用基于阈值和奇点检测相结合的方法来实现对心电信号中R波和S-T段的定位以及其他心 电特征的提取， 其步骤为 ： 0049 (1) 设定一个时间窗， 在该时间段内采用二次差分方法， 寻找波形奇异点 ； 对长度 为 N( 取 N 4096) 的数据进行奇点检测 ； 0050 奇点检测采用的算。
33、法为对一个采样周期的心电信号数据， 按照公式 (1) 至 (3) 进 行 diff(sign(diff(N) 运算 ： 0051 0052 0053 0054 其中diff为信号差分， sign为符号函数， N为一个采样周期的心电信号， 运算结果 为 -2 的点即为心电信号的极大值点 ； 0055 保留得到的所有极大值在一个数组 A 内， 并记录其对应波形位置 ； 0056 取阈值Rth， 与上述数组A的值比较， 保留大于Rth的极大值在数组B内， 并记录其 对应波形位置 ； 0057 认为 B 中数据为检测到的 R 波幅值， 其对应位置为 R 波位置 ； 0058 (2) 自适应阈值 Rth。
34、 设置 ： 0059 选取的N点ECG滤波之后的数据， 做奇异点检测之后， 寻找奇异点中的极大值信号 幅值范围， 按公式 (4)， 均分为 15 份， 并将每份幅值记录在 Th 数组中 ； 0060 1i15且iZ (4) 0061 设积分投影函数为 ： 0062 1 i 15 且 i Z (5) 0063 其中 . 0064 1 i 15 且 i Z (6) 0065 说 明 书 CN 103371814 A 9 7/7 页 10 0066 即做极大值点在划分的各幅值段 Th(i) 上的积分投影 ； 0067 选取零分布与非零分布的交接点处 i 值， 计算 Th(i)， 确定 Rth 为 T。
35、h(i) ； 0068 阈值确定出 R 波位置和幅值之后， 计算平均心率 Rate， 单位 ： 次 / 分， 0069 Rate 60*Nr/(nr(end)-nr(1)/fs) (8) 0070 其中， Nr 是固定时间窗中测得的 R 波个数， nr(end) 是固定时间窗内最后一个 R 波 的位置， nr(1) 是固定时间窗内第一个 R 波的位置， fs是数据采样率 ； 0071 根据平均心率加入相关误检、 漏检更正， 最终修正出正确 R 波位置和平均心率 Rate ； 0072 然后用公式 (9)， 计算 HRV 间期。 0073 HRV RR(i+1)-RR(i) (9) 0074 其。
36、中 RR 为相邻两个 R 波间期。 0075 综上算法， 可得到后续诊断所需特征参数 ： Rate 和 H RV ； 0076 (3) 根据此指标通过查表的方式将提取到的心电特征与 MIT-BI H 数据库分析总 结得出诊断结果。 0077 本发明可以诊断的心脏病包括房性早搏、 室性早搏、 室性心动过速、 停搏、 房性颤 动、 室性颤动、 房扑、 室颤、 室上速和交接逸搏。 说 明 书 CN 103371814 A 10 1/7 页 11 图 1 说 明 书 附 图 CN 103371814 A 11 2/7 页 12 图 2(a) 说 明 书 附 图 CN 103371814 A 12 3/7 页 13 图 2(b) 说 明 书 附 图 CN 103371814 A 13 4/7 页 14 图 3(a) 图 3(b) 图 4(a) 说 明 书 附 图 CN 103371814 A 14 5/7 页 15 图 4(b) 说 明 书 附 图 CN 103371814 A 15 6/7 页 16 图 5(a) 图 5(b) 说 明 书 附 图 CN 103371814 A 16 7/7 页 17 图 6(a) 图 6(b) 说 明 书 附 图 CN 103371814 A 17 。