CN114469019B - 脉搏波信号的滤波方法、装置和计算机设备 - Google Patents
脉搏波信号的滤波方法、装置和计算机设备 Download PDFInfo
- Publication number
- CN114469019B CN114469019B CN202210391331.7A CN202210391331A CN114469019B CN 114469019 B CN114469019 B CN 114469019B CN 202210391331 A CN202210391331 A CN 202210391331A CN 114469019 B CN114469019 B CN 114469019B
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- wave signal
- filtering
- signal
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Signal Processing (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本申请涉及医疗器械技术领域，公开了一种脉搏波信号的滤波方法、装置和计算机设备，首先对第一脉搏波信号进行第一带宽滤波处理，得到第二脉搏波信号，再计算第二脉搏波信号的信噪比，若信噪比不满足预期阈值，对第二脉搏波信号进行第二带宽滤波，将进行了第二带宽滤波的第二脉搏波信号与预设脉搏波信号进行对比，若第二脉搏波信号不是预设脉搏波信号，可根据预设的重复步骤，重复上述步骤，这样通过对第二脉搏波信号进行多次滤波，能够进一步滤除掉杂波以及不必要的噪声，从而使得得到的第三脉搏波信号、第四脉搏波信号的质量被提高，便于后期使用第三脉搏波信号、第四脉搏波信号对人体数据进行测算时，提高测算的准确率。
Description
技术领域
本申请涉及医疗器械技术领域，特别涉及一种脉搏波信号的滤波方法、装置和计算机设备。
背景技术
目前如智能手环、智能手表等可穿戴设备，都是通过传感器采集脉搏波信号(Photo PlethysmoGraphy，简称PPG)来实现对于心率值、血氧值等人体数据的测算。
但是脉搏波信号本质上是一种光信号，所以随着佩戴者的运动、环境光的变化等情况，会导致脉搏波信号的采集比较不稳定，噪声信号明显。为了准确的实现相关测算，对脉搏波信号进行滤波处理便至关重要。而现有技术中，往往只能够针对静态场景进行比较基本的滤波，难以满足使用需求。
发明内容
本申请的主要目的为提供一种脉搏波信号的滤波方法，旨在解决现有技术中脉搏波信号在动态场景中进行滤波后，仍存在大量杂波的技术问题。
本申请提出一种脉搏波信号的滤波方法，包括：
对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
根据所述第二脉搏波信号计算信噪比；
判断所述信噪比是否满足预期阈值；
若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号。
同时，本申请还提供一种脉搏波信号的滤波装置，包括：
第一带宽滤波模块，用于对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
计算模块，用于根据所述第二脉搏波信号计算信噪比；
第一判断模块，用于判断所述信噪比是否满足预期阈值；
第二带宽滤波模块，用于若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
判断模块，用于判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
重复模块，用于若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
滤波运算处理模块，用于对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
带通滤波模块，若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号。
同时，本申请还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述脉搏波信号的滤波方法的步骤。
同时，本申请还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述脉搏波信号的滤波方法的步骤。
有益效果：本申请首先对第一脉搏波信号进行第一带宽滤波处理，从而得到第二脉搏波信号，再计算第二脉搏波信号的信噪比，这样能够根据信噪比了解到第二脉搏波信号的质量情况，若信噪比不满足预期阈值，则代表当前的第二脉搏波信号的质量不佳，仍然还存在部分噪声与杂波，此时可对第二脉搏波信号进行第二带宽滤波，从而将多余的噪声与杂波滤除，之后，将进行了第二带宽滤波的第二脉搏波信号与预设脉搏波信号进行对比，若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，可根据预设的重复步骤，重复上述步骤，这样能够得到多个第二脉搏波信号，再对多个第二脉搏波信号进行滤波运算处理，最后得到的第三脉搏信号 ；若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，那么可对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号；这样通过对第二脉搏波信号进行多次滤波，能够进一步滤除掉杂波以及不必要的噪声，从而使得得到的第三脉搏波信号、第四脉搏波信号的质量被提高，便于后期使用第三脉搏波信号、第四脉搏波信号对人体数据（例如心率值、血氧值等）进行测算时，提高测算的准确率。
附图说明
图1为本申请一实施例的脉搏波信号的滤波方法流程示意图。
图2为本申请一实施例的脉搏波信号的滤波装置结构示意图。
图3为本申请一实施例的计算机设备内部结构示意图。
本申请目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
如图1-图3所示，本申请提出一种脉搏波信号的滤波方法，包括：
S1、对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
S2、根据所述第二脉搏波信号计算信噪比；
S3、判断所述信噪比是否满足预期阈值；
S4、若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
S5、判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
S6、若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
S7、对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
S8、若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号。
如上述步骤S1-S8所述，首先可对第一脉搏波信号按照常规的带宽（0.5hz-4hz）进行第一带宽滤波处理，从而得到第二脉搏波信号，再计算第二脉搏波信号的信噪比，这样能够根据信噪比了解到第二脉搏波信号的质量情况，若信噪比不满足预期阈值，则代表当前的第二脉搏波信号的质量不佳，仍然还存在部分噪声与杂波，此时可对第二脉搏波信号进行第二带宽滤波，从而将多余的噪声与杂波滤除，之后，将进行了第二带宽滤波的第二脉搏波信号与预设脉搏波信号进行对比，若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，可根据预设的重复步骤，重复步骤S2- S4，这样能够得到多个第二脉搏波信号，再对多个第二脉搏波信号进行滤波运算处理，最后得到的第三脉搏信号 ；若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，那么可对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号；这样通过对第二脉搏波信号进行多次滤波，能够进一步滤除掉杂波以及不必要的噪声，从而使得得到的第三脉搏波信号、第四脉搏波信号的质量被提高，便于后期使用第三脉搏波信号、第四脉搏波信号对人体数据（例如心率值、血氧值等）进行测算时，提高测算的准确率。
在一个实施例中，所述根据所述第二脉搏波信号计算信噪比的步骤S2，包括：
S21、获取所述第二脉搏波信号的热噪声电流值以及散粒噪声值；
S22、根据所述热噪声电流值、所述散粒噪声值计算总噪声输出电流值，其中，计算公式为：
其中，Iz表示总噪声输出电流值，
表示热噪声电流值，b表示散粒噪声值；
S23、获取所述第二脉搏波信号的电流平均值；
S24、根据所述电流平均值与所述总噪声输出电流值计算信噪比，其中，计算公式为：
其中，S表示信噪比，
表示电流平均值。
如上述步骤S21-S24所述，由于第二脉搏波信号存在的噪声主要是热噪声与散粒噪声，因此，可通过热噪声的电压值以及传感器的电阻值来计算热噪声电流值，从而得到热噪声电流值，再通过热噪声电流值与散粒噪声值，从而可计算出总噪声输出电流值，在根据电流平均值以及总噪声输出电流值计算信噪比，这样使得计算出来的信噪比较为客观、准确。
在一个实施例中，所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤S4，包括：
S41、获取第一带宽滤波的第一带通滤波值；
S42、根据所述第一带通滤波值设定第二带宽滤波的第二带通滤波值；
S43、根据所述第二带通滤波值对所述第二脉搏波信号进行偏离修正处理与降噪处理，得到第二带宽滤波后的所述第二脉搏波信号。
如上述步骤S41-S43所述，通过获取第一带宽滤波的第一带通滤波值，这样能够将第二带通滤波值设定为窄于第一带通滤波值的值，再根据第二带通滤波值对第二脉搏波信号进行偏离修正处理与降噪处理，这样能够防止第二脉搏波信号出现基线漂移的问题，便于第二脉搏波信号进行后续的滤波运算处理。
在一个实施例中，所述按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号的步骤S6，包括：
S61、获取预设重复次数；
S62、根据所述预设重复次数以及所述信噪比，设定第二带宽滤波的第二带通滤波值，得到与预设重复次数对应的多个第二带通滤波值；
S63、将每一个所述第二带通滤波值依次代入所述对所述第二脉搏波信号进行第二带宽滤波的重复步骤中，得到多个第二脉搏波信号。
如上述步骤S61-S63所述，通过获取预设重复次数，从而能够根据重复次数设定第二脉搏波信号每一次进行第二带宽滤波的第二带通滤波值，例如重复次数为三次，那么可将第一次重复时的第二带通滤波值设置在3.5 Hz -4 Hz之间，将第二次重复时的第二带通滤波值设置在1.5 Hz -2.5 Hz之间，将第三次重复时的第二带通滤波值设置在0.5 Hz -1.5Hz之间，这样使得第二脉搏波信号能够在不同的带宽频率下进行滤波，从而使得进行了第二带宽滤波所得到的第二脉搏波信号也不尽相同，这样便于对多个第二脉搏波信号进行整体滤波运算处理，从而得到质量相对较好的第三脉搏波信号。
在一个实施例中，所述对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号的步骤S7，包括：
S71、获取多个所述第二脉搏波信号的波峰信息；
S72、根据所述波峰信息对多个所述第二脉搏波信号进行波形特征提取，得到多个所述第二脉搏波信号的波形特征数据；
S73、根据所述波形特征数据对多个第二脉搏波信号进行分数计算，得到与多个所述第二脉搏波信号对应的多个分数值；
S74、从多个分数值中筛选出最高的分数值，并将最高的分数值所对应的第二脉搏波信号作为第三脉搏波信号。
如上述步骤S71-S74所述，通过获取每一个第二脉搏波信号的波峰信息，这样可根据波峰信息对每一个第二脉搏波信号的波形特征进行提取，从而得到每个第二脉搏波信号的波形特征数据，再根据波形特征数据计算每一个第二脉搏波信号的分数值，分数值越高，则该第二脉搏波信号的质量越好，这样能够从多个第二脉搏波信号中筛选出质量最好的第二脉搏波信号作为第三脉搏波信号，从而使得第三脉搏波信号作为基本数据对人体数据进行测算时，一稿人体数据测算的准确率。
在一个实施例中，波形特征数据包括：脉搏波的起始位置、波峰位置、波峰振幅、相邻波峰间隔、收缩时间和舒张时间。
本申请还提出一种脉搏波信号的滤波装置，其特征在于，包括：
第一带宽滤波模块1，用于对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
计算模块2，用于根据所述第二脉搏波信号计算信噪比；
第一判断模块3，用于判断所述信噪比是否满足预期阈值；
第二带宽滤波模块4，用于若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
判断模块5，用于判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
重复模块6，用于若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
滤波运算处理模块7，用于对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
带通滤波模块8，若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号。
在一个实施例中，所述计算模块2包括：
第一获取单元，用于获取所述第二脉搏波信号的热噪声电流值以及散粒噪声值；
第一计算单元，用于根据所述热噪声电流值、所述散粒噪声值计算总噪声输出电流值，其中，计算公式为：
其中，Iz表示总噪声输出电流值，
表示热噪声电流值，b表示散粒噪声值；
第二获取单元，用于获取所述第二脉搏波信号的电流平均值；
第二计算单元，用于根据所述电流平均值与所述总噪声输出电流值计算信噪比，其中，计算公式为：
其中，S表示信噪比，
表示电流平均值。
在一个实施例中，所述第二带宽滤波模块4，包括：
第三获取单元，用于获取第一带宽滤波的第一带通滤波值；
第一设定单元，用于根据所述第一带通滤波值设定第二带宽滤波的第二带通滤波值；
处理单元，用于根据所述第二带通滤波值对所述第二脉搏波信号进行偏离修正处理与降噪处理，得到第二带宽滤波后的所述第二脉搏波信号。
在一个实施例中，所述重复模块6，包括：
第四获取单元，用于获取预设重复次数；
第二设定单元，用于根据所述预设重复次数以及所述信噪比，设定第二带宽滤波的第二带通滤波值，得到与预设重复次数对应的多个第二带通滤波值；
代入单元，用于将每一个所述第二带通滤波值依次代入所述对所述第二脉搏波信号进行第二带宽滤波的重复步骤中，得到多个第二脉搏波信号。
在一个实施例中，所述滤波运算处理模块，包括：
第五获取单元，用于获取多个所述第二脉搏波信号的波峰信息；
波形特征提取单元，用于根据所述波峰信息对多个所述第二脉搏波信号进行波形特征提取，得到多个所述第二脉搏波信号的波形特征数据；
分数计算单元，用于根据所述波形特征数据对多个第二脉搏波信号进行分数计算，得到与多个所述第二脉搏波信号对应的多个分数值；
筛选单元，用于从多个分数值中筛选出最高的分数值，并将最高的分数值所对应的第二脉搏波信号作为第三脉搏波信号。
在一个实施例中，所述波形特征数据包括：脉搏波的起始位置、波峰位置、波峰振幅、相邻波峰间隔、收缩时间和舒张时间。
如图3所示，本申请还提供了一种计算机设备，该计算机设备可以是服务器，其内部结构可以如图3所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设计的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***、计算机程序和数据库。该内存器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的数据库用于存储脉搏波信号的滤波方法的过程需要的所有数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现脉搏波信号的滤波方法。
本领域技术人员可以理解，图3中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定。
本申请一实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述任意一个脉搏波信号的滤波方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储与一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的和实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可以包括只读存储器（ROM）、可编程ROM（PROM）、电可编程ROM（EPROM）、电可擦除可编程ROM（EEPROM）或闪存。易失性存储器可包括随机存取存储器（RAM）或者外部高速缓冲存储器。作为说明而非局限，RAM通过多种形式可得，诸如静态RAM（SRAM）、动态RAM（DRAM）、同步DRAM（SDRAM）、双速据率SDRAM（SSRSDRAM）、增强型SDRAM（ESDRAM）、同步链路（Synchlink）DRAM（SLDRAM）、存储器总线（Rambus）直接RAM（RDRAM）、直接存储器总线动态RAM（DRDRAM）、以及存储器总线动态RAM（RDRAM）等。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其它变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、装置、物品或者方法不仅包括那些要素，而且还包括没有明确列出的其它要素，或者是还包括为这种过程、装置、物品或者方法所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、装置、物品或者方法中还存在另外的相同要素。
以上所述仅为本申请的优选实施例，并非因此限制本申请的专利范围，凡是利用本申请说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本申请的专利保护范围内。
Claims (8)
1.一种脉搏波信号的滤波方法，其特征在于，包括：
对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
根据所述第二脉搏波信号计算信噪比；
判断所述信噪比是否满足预期阈值；
若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号；
其中，所述按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号的步骤，包括：
获取预设重复次数；
根据所述预设重复次数以及所述信噪比，设定第二带宽滤波的第二带通滤波值，得到与预设重复次数对应的多个第二带通滤波值；
将每一个所述第二带通滤波值依次代入所述对所述第二脉搏波信号进行第二带宽滤波的重复步骤中，得到多个第二脉搏波信号；
其中，所述对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号的步骤，包括：
获取多个所述第二脉搏波信号的波峰信息；
根据所述波峰信息对多个所述第二脉搏波信号进行波形特征提取，得到多个所述第二脉搏波信号的波形特征数据；
根据所述波形特征数据对多个第二脉搏波信号进行分数计算，得到与多个所述第二脉搏波信号对应的多个分数值；
从多个分数值中筛选出最高的分数值，并将最高的分数值所对应的第二脉搏波信号作为第三脉搏波信号。
2.根据权利要求1所述的脉搏波信号的滤波方法，其特征在于，所述根据所述第二脉搏波信号计算信噪比的步骤，包括：
获取所述第二脉搏波信号的热噪声电流值以及散粒噪声值；
根据所述热噪声电流值、所述散粒噪声值计算总噪声输出电流值，其中，计算公式为：
其中，Iz表示总噪声输出电流值，
表示热噪声电流值，b表示散粒噪声值；
获取所述第二脉搏波信号的电流平均值；
根据所述电流平均值与所述总噪声输出电流值计算信噪比，其中，计算公式为：
其中，S表示信噪比，
表示电流平均值。
3.根据权利要求1所述的脉搏波信号的滤波方法，其特征在于，所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，包括：
获取第一带宽滤波的第一带通滤波值；
根据所述第一带通滤波值设定第二带宽滤波的第二带通滤波值；
根据所述第二带通滤波值对所述第二脉搏波信号进行偏离修正处理与降噪处理，得到第二带宽滤波后的所述第二脉搏波信号。
4.根据权利要求1所述的脉搏波信号的滤波方法，其特征在于，所述波形特征数据包括：脉搏波的起始位置、波峰位置、波峰振幅、相邻波峰间隔、收缩时间和舒张时间。
5.一种脉搏波信号的滤波装置，其特征在于，包括：
第一带宽滤波模块，用于对第一脉搏波信号进行第一带宽滤波，得到第二脉搏波信号；
计算模块，用于根据所述第二脉搏波信号计算信噪比；
第一判断模块，用于判断所述信噪比是否满足预期阈值；
第二带宽滤波模块，用于若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波；
判断模块，用于判断第二带宽滤波后的所述第二脉搏波信号是否为预设脉搏波信号；
重复模块，用于若进行了第二带宽滤波的所述第二脉搏波信号不是预设脉搏波信号，按照预设重复次数，重复所述根据所述第二脉搏波信号计算信噪比的步骤至所述若所述信噪比不满足预期阈值，对所述第二脉搏波信号进行第二带宽滤波的步骤，得到多个滤波后的第二脉搏波信号；
滤波运算处理模块，用于对多个第二脉搏波信号进行滤波运算处理，得到第三脉搏波信号；
带通滤波模块，若第二带宽滤波后的所述第二脉搏波信号是预设脉搏波信号，对所述第二脉搏波信号进行带通滤波，得到第四脉搏波信号；
其中，所述重复模块，包括：
第四获取单元，用于获取预设重复次数；
第二设定单元，用于根据所述预设重复次数以及所述信噪比，设定第二带宽滤波的第二带通滤波值，得到与预设重复次数对应的多个第二带通滤波值；
代入单元，用于将每一个所述第二带通滤波值依次代入所述对所述第二脉搏波信号进行第二带宽滤波的重复步骤中，得到多个第二脉搏波信号；
其中，所述滤波运算处理模块，包括：
第五获取单元，用于获取多个所述第二脉搏波信号的波峰信息；
波形特征提取单元，用于根据所述波峰信息对多个所述第二脉搏波信号进行波形特征提取，得到多个所述第二脉搏波信号的波形特征数据；
分数计算单元，用于根据所述波形特征数据对多个第二脉搏波信号进行分数计算，得到与多个所述第二脉搏波信号对应的多个分数值；
筛选单元，用于从多个分数值中筛选出最高的分数值，并将最高的分数值所对应的第二脉搏波信号作为第三脉搏波信号。
6.根据权利要求5所述的一种脉搏波信号的滤波装置，其特征在于，所述计算模块包括：
第一获取单元，用于获取所述第二脉搏波信号的热噪声电流值以及散粒噪声值；
第一计算单元，用于根据所述热噪声电流值、所述散粒噪声值计算总噪声输出电流值，其中，计算公式为：
Iz= 
；
其中，Iz表示总噪声输出电流值，
表示热噪声电流值，b表示散粒噪声值；
第二获取单元，用于获取所述第二脉搏波信号的电流平均值；
第二计算单元，用于根据所述电流平均值与所述总噪声输出电流值计算信噪比，其中，计算公式为：
其中，S表示信噪比，
表示电流平均值。
7.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述脉搏波信号的滤波方法的步骤。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4中任一项所述脉搏波信号的滤波方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210391331.7A CN114469019B (zh) | 2022-04-14 | 2022-04-14 | 脉搏波信号的滤波方法、装置和计算机设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210391331.7A CN114469019B (zh) | 2022-04-14 | 2022-04-14 | 脉搏波信号的滤波方法、装置和计算机设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114469019A CN114469019A (zh) | 2022-05-13 |
| CN114469019B true CN114469019B (zh) | 2022-06-21 |
Family
ID=81489053
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210391331.7A Active CN114469019B (zh) | 2022-04-14 | 2022-04-14 | 脉搏波信号的滤波方法、装置和计算机设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114469019B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11276448A (ja) * | 1998-03-31 | 1999-10-12 | Seiko Epson Corp | 信号抽出装置および信号抽出方法 |
| CN101911704A (zh) * | 2008-01-14 | 2010-12-08 | 汤姆森许可贸易公司 | 使用基于多格稀疏性的滤波来进行去伪像滤波的方法和设备 |
| CN103823910A (zh) * | 2012-11-16 | 2014-05-28 | 镇江市顶智微电子科技有限公司 | 一种基于智能算法的空间日照分析平台实现方法 |
| CN108801131A (zh) * | 2018-06-11 | 2018-11-13 | 华中师范大学 | 北斗高频变形监测数据的滤波方法及*** |
| CN109316172A (zh) * | 2018-10-18 | 2019-02-12 | 大连理工大学 | 基于脉搏波特征参数对人体健康状况的检测方法 |
| CN109565098A (zh) * | 2016-08-18 | 2019-04-02 | 华为技术有限公司 | 一种滤波器 |
| CN110801210A (zh) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | 一种脉搏波信号的滤波方法、装置、可读介质及电子设备 |
| CN111513752A (zh) * | 2020-04-01 | 2020-08-11 | 中国海洋大学 | 一种基于脉搏声信号的脉诊仪 |
| CN112152725A (zh) * | 2019-06-28 | 2020-12-29 | 阿尔卡塔尔海底网络公司 | 抑制传输信号不稳定性所产生噪声的方法和装置 |
| CN112262512A (zh) * | 2018-04-10 | 2021-01-22 | 伊顿智能动力有限公司 | 使用组合式断路器和继电器在电动移动应用中进行电源分配的***、方法和装置 |
| CN212588311U (zh) * | 2020-04-23 | 2021-02-23 | 深圳芯森微电子有限公司 | 一种可配置低频高通滤波器 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2003222001A1 (en) * | 2002-03-18 | 2003-10-08 | Sonomedica, Llc | Method and system for generating a likelihood of cardiovascular disease from analyzing cardiovascular sound signals. |
| US7924936B2 (en) * | 2006-06-19 | 2011-04-12 | Xg Technology, Inc. | System and method for wave damping |
| US8212941B2 (en) * | 2008-04-30 | 2012-07-03 | Mediatek Inc. | Digitized analog TV signal processing system |
| US8372013B2 (en) * | 2008-09-19 | 2013-02-12 | Medtronic, Inc. | Method and apparatus for determining a respiration parameter in a medical device |
| CN201333055Y (zh) * | 2008-12-05 | 2009-10-28 | 中国医学科学院生物医学工程研究所 | 基于分布式算法实现的高频超声信号动态滤波器 |
| US8483628B2 (en) * | 2010-06-03 | 2013-07-09 | Broadcom Corporation | Multiple-phase frequency translated filter |
| CN102004264B (zh) * | 2010-10-18 | 2015-09-23 | 中国石油化工股份有限公司 | 一种地震采集资料质量定量分析与评价方法 |
| CN107595249B (zh) * | 2017-09-30 | 2021-03-23 | 深圳前海全民健康科技有限公司 | 基于脉搏波的怀孕女性筛查方法 |
-
2022
- 2022-04-14 CN CN202210391331.7A patent/CN114469019B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11276448A (ja) * | 1998-03-31 | 1999-10-12 | Seiko Epson Corp | 信号抽出装置および信号抽出方法 |
| CN101911704A (zh) * | 2008-01-14 | 2010-12-08 | 汤姆森许可贸易公司 | 使用基于多格稀疏性的滤波来进行去伪像滤波的方法和设备 |
| CN103823910A (zh) * | 2012-11-16 | 2014-05-28 | 镇江市顶智微电子科技有限公司 | 一种基于智能算法的空间日照分析平台实现方法 |
| CN109565098A (zh) * | 2016-08-18 | 2019-04-02 | 华为技术有限公司 | 一种滤波器 |
| CN112262512A (zh) * | 2018-04-10 | 2021-01-22 | 伊顿智能动力有限公司 | 使用组合式断路器和继电器在电动移动应用中进行电源分配的***、方法和装置 |
| CN108801131A (zh) * | 2018-06-11 | 2018-11-13 | 华中师范大学 | 北斗高频变形监测数据的滤波方法及*** |
| CN109316172A (zh) * | 2018-10-18 | 2019-02-12 | 大连理工大学 | 基于脉搏波特征参数对人体健康状况的检测方法 |
| CN112152725A (zh) * | 2019-06-28 | 2020-12-29 | 阿尔卡塔尔海底网络公司 | 抑制传输信号不稳定性所产生噪声的方法和装置 |
| CN110801210A (zh) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | 一种脉搏波信号的滤波方法、装置、可读介质及电子设备 |
| CN111513752A (zh) * | 2020-04-01 | 2020-08-11 | 中国海洋大学 | 一种基于脉搏声信号的脉诊仪 |
| CN212588311U (zh) * | 2020-04-23 | 2021-02-23 | 深圳芯森微电子有限公司 | 一种可配置低频高通滤波器 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114469019A (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20190295729A1 (en) | Universal non-invasive blood glucose estimation method based on time series analysis | |
| DE60207717T2 (de) | Signalkomponenten-prozessor | |
| US8099250B2 (en) | Impedance parameter values | |
| TWI583355B (zh) | Heart rate detection method and heart rate detection device | |
| JP5283381B2 (ja) | 心調律信号のシリーズ（ｒｒ）を処理するための方法、及び心調律の変動性を分析するための、特に生物の痛み又はストレスを評価するためのその使用 | |
| US20160270699A1 (en) | Personal authentication apparatus, personal authentication method, and recording medium | |
| CN108652640B (zh) | 一种基于心电信号的无创血糖检测方法及*** | |
| CN113288131B (zh) | 基于图卷积网络的无创血糖检测方法、处理器及装置 | |
| CN112043258A (zh) | 动态心率预测方法、装置、设备及可读存储介质 | |
| EP1639497A1 (en) | Method and apparatus for extracting causal information from a chaotic time series | |
| CN110840401A (zh) | 一种心率值测定方法、装置、设备及可读存储介质 | |
| CN114469019B (zh) | 脉搏波信号的滤波方法、装置和计算机设备 | |
| WO2020169719A1 (en) | Methods and apparatus for analyte concentration monitoring using harmonic relationships | |
| CN112842362A (zh) | 用于植入式电刺激装置的自响应检测参数优化方法及*** | |
| DE112018001399T5 (de) | Blutdruckdatenverarbeitungsvorrichtung, blutdruckdatenverarbeitungsverfahren und -programm | |
| CN110881972A (zh) | St检测方法、装置、计算机设备及存储介质 | |
| CN108261196B (zh) | 心电电极脱落检测方法、装置、计算机设备及存储介质 | |
| CN117883082A (zh) | 一种异常情绪识别方法、***、设备及介质 | |
| DE112016002969T5 (de) | Blutdruckverhältnisberechnungseinrichtung, blutdruckverhältnisberechnungsverfahren, blutdruckverhältnisberechnungsprogramm und aufzeichungsmedium, auf dem das programm aufgezeichnet ist. | |
| CN113647905A (zh) | 深度神经网络预测血压的训练方法、计算机装置和存储介质 | |
| CN117357080A (zh) | 近红外光谱信号去噪方法及装置、终端设备、存储介质 | |
| CN109063652B (zh) | 一种信号处理方法、***及计算机存储介质 | |
| DE112016005420T5 (de) | Blutdruckinformations-berechnungsvorrichtung, blutdruckinformations-berechnungsverfahren, blutdruckinformations- berechnungsprogramm und aufzeichnungsmedium zum aufzeichnendes programms | |
| US20080208017A1 (en) | Method and Apparatus for Monitoring a Sedated Patient | |
| CN113100735A (zh) | 一种基于静脉血流图的血栓检测装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |