一种心率测量的数字信号处理方法、系统及设备的制造方法
【技术领域】
[0001] 本发明涉及数字信号处理技术领域，尤其涉及一种心率测量的数字信号处理方 法、系统及设备。
【背景技术】
[0002] 心率是一个反应心脏功能的指标，经常测试心率可以预知心脏疾病以便及时治疗 疾病。在运动时使用心率监测可以快速清楚地反映运动时身体的信息，以使运动强度保持 在适当的水平，使健身效果更好。随着科技的发展进步，如今越来越多的监测设备开始配备 心率传感器，包括心率监测带、智能手表、智能手环甚至是耳机。心率监测设备的多样化，鉴 于传统对心率信号处理的复杂性，使得设备监测的功能总是不够完善，甚至很多设备需要 安静的状态下测量，不适合运动中持续测量心率。
【发明内容】
[0003] 本发明实施例的目的在于提出一种心率测量的数字信号处理方法，旨在解决现有 技术心率监测设备的功能不够完善，需要安静的状态下测量，不适合运动中持续测量心率 问题。
[0004] 本发明实施例是这样实现的，一种心率测量的数字信号处理方法，所述方法包括 以下步骤：
[0005] 获取数字信号A，并得到数字信号A的采样频率f，采样时长t，信号采样容量η ;
[0006] 取数字信号A的前M个样本，评估样本方差；
[0007] 对数字信号A进行低通滤波处理，得到数字信号Α' ；
[0008] 对所述数字信号Α'做差分操作，得到差分数字信号B ;
[0009] 从所述数字信号B中提取满足设定阈值的所有点，并记录每个点的位置，形成点 集T ;
[0010] 对所述点集Τ，按距离d进行点集分组，并提取每组中点对应位置在数字信号A中 的对应点的最大值，记录数字信号A中每组最大值的点位置，形成新的点集S ;
[0011] 计算点集S的平均距离根据点集S的平均距离€和采样频率f计算心率q ;
[0012] 输出心率q。
[0013] 本发明实施例的另一目的在于提出一种心率测量的数字信号处理系统，所述系统 包括：
[0014] 数字信号获取模块，用于获取数字信号A，并得到数字信号A的采样频率f，采样时 长t，信号采样容量η;
[0015] 样本方差评估模块，用于取数字信号A的前M个样本，评估样本方差；
[0016] 具体为：分割M个样本为N组，计算每组的方差σ,:，其中i = 1，2, "·，Ν，Ν为正整 数；评估样本方差：
[0017] 低通滤波处理模块，用于对数字信号A进行低通滤波处理，得到数字信号A' ；
[0018] 具体为：
[0019] 滤波模板半径
[0020] 滤波模板为
[0021] 用此滤波模板Ll对数字信号A做卷积操作，得到数字信号A'，即?表 示卷积算子，e为自然底数；
[0022] 差分操作模块，用于对所述数字信号A'做差分操作，得到差分数字信号B ;
[0023] 具体为：
[0024] 滤波模板为：滤波模板为：L2 = [1 -2 1];
[0025] 用滤波模板L2与信号A'卷积操作，得到差分数字信号B ;即及表示卷 积算子，e为自然底数；
[0026] 点集T获取模块，用于从所述数字信号B中提取满足设定阈值的所有点，并记录每 个点的位置，形成点集T。；
[0027] 所述设定阈值为b' λ i，由经验获取的。，
[0028] 其中bj为数字信号B中信号值，j = 1，2,…，k，k为正整数； A1G (-0.01,-0.055)；
[0029] 点集S获取装置，用于对所述点集T，按距离d进行点集分组，并提取每组中点对应 位置在数字信号A中的对应点的最大值，记录数字信号A中每组最大值的点位置，形成新的 点集S。；
[0030] 其中
[0031] 心率计算装置，用于计算点集S的平均距离石，根据点集S的平均距离I和采样 频率f计算心率q ;
[0032] 计算方法为：
[0033] 心率输出模块，用于输出心率q。
[0034] 本发明实施例的另一目的在于提出一种包含所述心率测量的数字信号处理系统 的设备。
[0035] 本发明的有益效果
[0036] 本发明提出一种心率测量的数字信号处理方法、系统及设备，本发明方法可在廉 价单片机上完全实时处理，在时域信号中处理计算，简单，易于实现，可移植性强，优化后的 流程可处理采样率500hz的数字信号，在主频72M、RAM 64k的廉价单片机上完全实时处理。 在低采样率下，测量精度高。抗噪能力极强，不会出现心跳尖峰漏检的情况，可测环境适应 性强，甚至在多变的信号样本中可实现零概率心跳尖峰漏检。
【附图说明】
[0037] 图1是本发明优选实施例一种心率测量的数字信号处理方法流程图；
[0038] 图2是数字信号A的部分波形图；
[0039] 图3是低通滤波处理后数字信号A' ；的波形图；
[0040] 图4是点集T在数字信号B中的位置图；
[0041] 图5是点集S在数字信号A中的位置图；
[0042] 图6是本发明优选实施例一种心率测量的数字信号处理系统结构图；
[0043] 图7是图6系统中的点集S获取装置的结构图；
[0044] 图8是图6系统中的心率计算装置的结构图。
【具体实施方式】
[0045] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对 本发明进行进一步详细说明，为了便于说明，仅示出了与本发明实施例相关的部分。应当理 解，此处所描写的具体实施例，仅仅用于解释本发明，并不用以限制本发明。
[0046] 本发明实施例提出一种心率测量的数字信号处理方法、系统及设备，本发明实施 例方法先对数字信号A进行数据预处理，然后进行数据过滤，形成点集T，对所述点集T按距 离d进行点集分组，并提取每组中点对应位置在数字信号A中的对应点的最大值，并记录数 字信号A中每组最大值的点位置，形成新的点集S ;计算点集S的平均距离?根据点集S ，+' 的平均距离&和采样频率f计算心率q ;输出心率q。本发明方法可在廉价单片机上完全实 时处理，在时域信号中处理计算，简单，易于实现，可移植性强；在低采样率下，测量精度高。 抗噪能力极强，不会出现心跳尖峰漏检的情况，可测环境适应性强，甚至在多变的信号样本 中可实现零概率心跳尖峰漏检。
[0047] 实施例一
[0048] 图1是本发明优选实施例一种心率测量的数字信号处理方法流程图，所述方法包 括以下步骤：
[0049] S101、获取数字信号A，并得到数字信号A的采样频率f，采样时长t，信号采样容量 η ;图2是数字信号A的部分波形图，横坐标表示样点序号，纵坐标表示信号幅值。
[0050] S102、取数字信号A的前M个样本，评估样本方差；
[0051] 具体为：
[0052] 分割M个样本为N组，计算每组的方差其中i = 1，2,…，N，N为正整数；
[0053] 评估样本方差：
[0054] S103、对数字信号A进行低通滤波处理，得到数字信号A' ；
[0055] 图3是低通滤波处理后数字信号A'的波形图；
[0056] 具体为：
[0057] 滤波模板半径
[0058] 滤波模板为
[0059] 用此滤波模板LI对数字信号A做卷积操作，得到数字信号A'，即，=J ? LI，?表 示卷积算子，e为自然底数；
[0060] S104、对所述数字信号A'做差分操作，得到差分数字信号B ;
[0061] 具体为：
[0062] 滤波模板为 L2 = [1 -2 1];
[0063] 用滤波模板L2与信号A'卷积操作，得到差分数字信号B ;即B = fiaEX，#表示卷 积算子，e为自然底数；
[0064] S105、从所述数字信号B中提取满足设定阈值的所有点，并记录每个点的位置，形 成点集T。图4是点集T在数字信号B中的位置图；
[0065] 所述设定阈值为b' λ i，由经验获取的。
[0066] 其中bj为数字信号B中信号值，j = 1，2,…，k，k为正整数； A1G (-0.01,-0.055)；
[0067] S106、对所述点集T，按距离d进行点集分组，并提取每组中点对应位置在数字信 号A中的对应点的最大值，记录数字信号A中每组最大值的点位置，形成新的点集S。图5 是点集S在数字信号A中的位置图；
[0068] 其中
[0069] S107、计算点集S的平均距离ζ，根据点集S的平均距离ξ和采样频率f
[0070] 计算心率q ;
[0071] 计算心率q的方法为：
[0072] 步骤S107也可以替换为：
[0073] 先对点集S的进行离群过滤，然后再计算离群过滤处理后的点集S'的平均距离 f、根据点集S'的平均距离'和采样频率f计算心率q ;
[0074] 计算心率q的方法为：
[0075] 离群过滤可以采用高斯正态分布方法，估计平均值.与标准差〇，只保留满足
范围内的点，超出的部分舍去。离群过滤中的高斯正态分布方法为本技 术领域的通用方法，在此不再赘述。
[0076] S108、输出心率 q。
[0077] 实施例二
[0078] 图6是本发明优选实施例一种心率测量的数字信号处理系统结构图；所述系统包 括：
[0079] 数字信号获取模块，用于获取数字信号A，并得到数字信号A的采样频率f，采样时 长t，信号采样容量η;
[0080] 样本方差评估模块，用于取数字信号A的前M个样本，评估样本方差；
[0081] 具体为：分割M个样本为N组，计算每组的方差G?其中i = 1，2, "·，Ν，N为正整 f- 数；评估样本方差：
[0082] 低通滤波处理模块，用于对数字信号A进行低通滤波处理，得到数字信号A' ；
[0083] 具体为：
[0084] 滤波模板半径