WO2022111203A1 - 心率检测方法和装置 - Google Patents
心率检测方法和装置 Download PDFInfo
- Publication number
- WO2022111203A1 WO2022111203A1 PCT/CN2021/126987 CN2021126987W WO2022111203A1 WO 2022111203 A1 WO2022111203 A1 WO 2022111203A1 CN 2021126987 W CN2021126987 W CN 2021126987W WO 2022111203 A1 WO2022111203 A1 WO 2022111203A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- acceleration data
- axis
- data
- target acceleration
- heart rate
- Prior art date
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 112
- 230000001133 acceleration Effects 0.000 claims abstract description 292
- 238000012545 processing Methods 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000004927 fusion Effects 0.000 claims description 64
- 238000001914 filtration Methods 0.000 claims description 23
- 238000004422 calculation algorithm Methods 0.000 claims description 22
- 230000000737 periodic effect Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 11
- 238000000513 principal component analysis Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 abstract description 11
- 238000009532 heart rate measurement Methods 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 27
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 239000000284 extract Substances 0.000 description 4
- 238000013186 photoplethysmography Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000010432 diamond Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000004984 smart glass Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
- A61B5/6803—Head-worn items, e.g. helmets, masks, headphones or goggles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
- A61B5/681—Wristwatch-type devices
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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
-
- 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/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本申请提出一种心率检测方法和装置,其中,方法包括:采集用户的多轴原始加速度数据,根据多轴原始加速度数据判断用户是否满足预设检测条件;若获知用户满足预设检测条件,则对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据;根据融合频域加速度数据中的峰值数据确定用户的心率值。由此,基于多轴加速度来确定心率,降低了测量功耗,并且在频域中提取心率值,提高了心率的测量准确性。
Description
相关申请的交叉引用
本申请要求安徽华米健康科技有限公司于2020年11月25日提交的、发明名称为“心率检测方法和装置”的、中国专利申请号“202011339438.4”的优先权。
本申请涉及数字信号处理技术领域,尤其涉及一种心率检测方法和装置
随着智能设备的逐渐普及,基于智能设备来进行心率测量,由于可以随时监测用户的身体健康状态,因此,逐渐成为流行。
相关技术中,利用光电容积描记(PhotoPlethySmogram,PPG)技术进行人体运动心率的检测,PPG是红外无损检测技术,需要频繁的开关PPG传感器,导致测量的功耗较大,大大缩短了测量设备的续航时间。
发明内容
本申请旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本申请的第一个目的在于提出一种心率检测方法,以基于多轴加速度来确定心率,降低了测量功耗,并且在频域中提取心率,提高了心率的测量准确性。
本申请的第二个目的在于提出一种心率检测装置。
本申请的第三个目的在于提出一种计算机设备。
本申请的第四个目的在于提出一种非临时性计算机可读存储介质。
为达上述目的,本申请第一方面实施例提出了一种心率检测方法,包括:包括:采集用户的多轴原始加速度数据,根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件;若获知所述用户满足所述预设检测条件,则对所述多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;对所述多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据;根据所述融合频域加速度数据中的峰值数据确定所述用户的心率值。
为达上述目的,本申请第二方面实施例提出了一种心率检测装置,包括:判断模块,用于采集用户的多轴原始加速度数据,根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件;滤波处理模块,用于在获知所述用户满足所述预设检测条件时,对所述多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;获取模块,用于对所述 多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据;确定模块,用于根据所述融合频域加速度数据中的峰值数据确定所述用户的心率值。
为达上述目的,本申请第三方面实施例提出了一种计算机设备,包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时,实现如上述实施例所描述的心率检测方法。
为了实现上述目的,本申请第四方面实施例提出了一种非临时性计算机可读存储介质,当所述存储介质中的指令由处理器被执行时,使得处理器能够实现如上述实施例所描述的心率检测方法。
本申请实施例至少包括如下有益技术效果:
采集用户的多轴原始加速度数据,根据多轴原始加速度数据判断用户是否满足预设检测条件,进而,若获知用户满足预设检测条件,则对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,最后,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,并根据融合频域加速度数据中的峰值数据确定用户的心率值。由此,基于多轴加速度来确定心率,降低了测量功耗,并且在频域中提取心率值,提高了心率的测量准确性。
本申请附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为本申请实施例所提供的一种心率检测方法的流程示意图;
图2为本申请实施例所提供的另一种心率检测方法的流程示意图;
图3为本申请实施例所提供的又一种心率检测方法的流程示意图;
图4是根据本申请实施例所提供的一种心率检测场景图;
图5是根据本申请实施例所提供的一种心率检测模块场景示意图;
图6是根据本申请实施例所提供的另一种心率检测场景图;
图7是根据本申请实施例所提供的又一种心率检测场景图;
图8是根据本申请实施例所提供的还一种心率检测场景图;
图9是根据本申请实施例所提供的还一种心率检测方法的流程示意图;
图10是根据本申请实施例所提供的再一种心率检测场景图;
图11为本申请实施例提供的一种心率检测装置的结构示意图;以及
图12为本申请实施例提供的另一种心率检测装置的结构示意图。
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参考附图描述本申请实施例的心率检测方法和装置。其中,本申请心率检测方法和装置的执行主体,可以是任意便携式终端设备,该终端设备可以是手机、平板电脑、个人数字助理、穿戴式设备等具有各种操作系统的硬件设备,该穿戴式设备可以是智能手环、智能手表、智能眼镜等。
图1为本申请实施例所提供的一种心率检测方法的流程示意图。如图1所示,该心率检测方法包括以下步骤:
步骤101,采集用户的多轴原始加速度数据,根据多轴原始加速度数据判断用户是否满足预设检测条件。
在本实施例中,可以根据加速度计等设备来采集用户的多轴原始加速度数据,该加速计可以设置在如上述提到的便携式终端设备中,另外,本实施例中提到的多轴原始加速度数据可以为x,y,z三轴加速度数据中的至少两种。
不难理解的是,因为活动导致的加速度计测量的多轴原始加速度数据变化远远超过心率导致的变化,在佩戴者活动的情况下,采集到的多轴原始加速度数据无法用以准确计算心率,因此,在本实施例中,仅在判断出佩戴者处于安静状态等活动量较小的情况下进行心率检测。
其中,本实施例中的预设检测条件对应于佩戴者处于安静状态等活动量较小的情况,下面示例说明如何根据多轴原始加速度数据判断所述用户是否满足预设检测条件:
示例一:
在本示例中,如图2所示,根据多轴原始加速度数据判断用户是否满足预设检测条件,包括:
步骤201,按照第一预设算法对每轴原始加速度数据进行计算,以获取与多轴原始加速度数据对应的多轴特征数据。
其中,第一预设算法可以是对每轴原始加速度数据取方差计算得到方差值,将方差值作为对应轴的原始加速度数据的轴特征数据,第一预设算法可以是对每轴原始加速度数据进行标准差计算得到标准差,将标准差作为将方差值作为对应轴的原始加速度数据的轴特征数据,第一预设算法也可以为对每轴原始加速度数据进行幅度值计算,将得到的幅度值等反应每轴原始加速度大小的特征数据作为将方差值作为对应轴的原始加速度数据的轴特 征数据。
当然,也可以根据实验数据确定每轴原始加速度数据的预设门限值,从而,第一预设算法为统计原始加速度数据超过对应预设门限值的次数,将次数作为对应轴的原始加速度数据的轴特征数据。
步骤202,将每轴特征数据与对应轴的第一预设阈值进行比较,获取多轴特征数据的比较结果。
应当理解的是,预先根据大量实验数据设置每轴特征数据所在轴的预设阈值,其中,该预设阈值可以与加速度计所在的设备硬件有关,进而,将每轴特征数据与对应轴的预设阈值进行比较,获取多轴特征数据的比较结果,其中,该比较结果可以为每轴特征数据与对应轴的预设阈值差值。
步骤203,若根据比较结果获知多轴特征数据均满足预设的第一检测范围,则获知用户满足预设检测条件。
在本实施例中,预先根据大量实验数据设置第一检测范围,该第一检测范围可以对应于上述差值的取值范围,当差值位于该第一检测范围内,则认为用户处于相对安安静的状态,认为满足预设检测条件。
步骤204,若根据比较结果获知至少一轴特征数据不满足预设的检测范围,则获知用户不满足预设检测条件。
在本实施例中,多轴特征数据均满足预设的检测范围才认为用户处于安静的状态,否则若根据比较结果获知至少一轴特征数据不满足预设的第一检测范围,则获知用户不满足预设检测条件。
示例二:
在本示例中,如图3所示,根据多轴原始加速度数据判断用户是否满足预设检测条件,包括:
步骤301,对多轴原始加速度数据求和处理获取融合原始加速度数据。
在本实施例中,对多轴原始加速度数据求和处理获取融合原始加速度数据,将多轴原始加速度数据作为一个整体进行判断。
其中,对多轴原始加速度数据求和处理可以理解为将相同时间点采集到的多轴原始加速度值求和得到对应的原始加速度数据,该原始加速数据整体反应了多轴原始加速度数据的大小。
步骤302,按照第二预设算法对融合原始加速度数据进行计算,获取融合特征数据。
第二预设算法可以为融合原始加速度数据取方差计算得到的方差值,将方差值作为融 合特征数据,第二预设算法可以为融合原始加速度数据进行标准差计算得到的标准差,将标准差作为融合特征数据,第二预设算法也可以为对融合原始加速度数据进行幅度值计算得到的幅度值等反应融合原始加速度数据大小的特征数据,将特征数据作为融合特征数据。
当然,也可以根据实验数据确定融合原始加速度数据的预设门限值,从而,第二预设算法为统计融合原始加速度数据超过对应预设门限值的次数,将次数作为融合特征数据。
步骤303,将融合特征数据与对应的第二预设阈值进行比较,获取融合特征数据的比较结果。
应当理解的是,预先根据大量实验数据设置融合特征数据的第二预设阈值,其中,该第二预设阈值可以与加速度计所在的设备硬件有关,进而,将融合特征数据与对应轴的第二预设阈值进行比较,获取融合特征数据的比较结果,其中,该比较结果可以为融合特征数据与对应的第二预设阈值差值。
步骤304,若根据比较结果获知融合特征数据满足预设的第二检测范围,则获知用户满足预设检测条件。
在本实施例中,预先根据大量实验数据设置第二检测范围,该第二检测范围可以对应于上述差值的取值范围,当差值位于该第二检测范围内,则认为用户处于相对安安静的状态,认为满足预设检测条件。
步骤305,若根据比较结果获知融合特征数据不满足预设的检测范围,则获知用户不满足预设检测条件。
在本实施例中,如根据比较结果获知融合特征数据不满足预设的检测范围,则获知用户不满足预设检测条件,用户可能处于运动状态等,此时测量的心率是不准确的。
步骤102,若获知用户满足预设检测条件,则对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据。
在本实施例中,若是获知用户满足预设检测条件,则认为此时获取心率比较准确,从而,对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,这里的高通滤波处理可以理解为预处理操作,以去除基线漂移以及呼吸率对于心率检测的影响。其中,高通滤波处理用于将较低频值的多轴原始加速度筛选掉,其中,高通滤波处理的截止频率可以根据实验数据标定。
在本申请的一个实施例中,对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,可以为对每轴原始加速度数据进行N阶的滑动平均获取对应的轴平均值,其中,N为大于1的整数,从而保证对应的多轴原始加速度数据进行高通滤波处理后不失真,上述N可以根据信号的特征以及采样率适配得到,得到方式可以为基于预选训练的模型得到 等现有技术实现,在此不再赘述。其中,滑动平均法又称移动平均法。在简单平均数法基础上,通过顺序逐期增减新旧数据求算移动平均值,借以消除偶然变动因素,找出事物发展趋势,并据此进行预测的方法。上述阶数可以理解为滑动平均算法中的窗口的宽度。
进一步的,分别用每轴原始加速度数据减去对应的轴平均值,以获取多轴原始加速度数据对应的多轴目标加速度数据,以达到高通滤波的效果。
在本实施例中,考虑到高通滤波后的数据仅仅包含很强的振动信息,因此,为了方便观测,对处理后的每轴目标加速度数据分别上下移动M个单位以方便观测,其中,M可以为任意根据场景需要设置的整数。比如,M可以为40。
步骤103,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据。
应当理解的是,在本实施例中,对多轴目标加速度数据进行傅里叶变换处理,将多轴目标加速度数据由时域转化为频域,由此,即使在时域中的心率峰值不明显,比如,如图4所示,在时域上,每轴目标加速度数据的峰值(幅度值)并不十分明显,导致提取心率困难。
本申请中基于频域中的融合频域加速度数据来提取对应的心率,对即使在时域中峰值不明显的数据,也能够提取对应的心率。
在本申请的一个实施例中,为了能够更好的提取心率,在对多轴目标加速度数据进行傅里叶变换处理之前,还可以按照第三预设算法对每轴目标加速度数据进行增强处理,以强化每轴目标加速度数据的包络特征等,其中,在一些可能的示例中,上述第三预设算法可以为平方算法,即对每轴目标加速度数据求取平方值作为新的每轴目标加速度数据等。在另一些可能的示例中,上述第三预设算法可以为对每轴目标加速度数据加上对应轴的预设值,以实现对每轴目标加速度数据的加强。
步骤104,根据融合频域加速度数据中的峰值数据确定用户的心率值。
具体的,在获取融合频域加速度数据后,根据融合频域加速度数据中的峰值确定用户心率值,也可以理解,选取频率响应中能量最强的频点当作心率值输出,其中,由于此时的峰值对应的心率值是在频域的心率值,因此,可以将该峰值数据转换为时域。
为了使得本领域的技术人员,对本申请实施例的心率检测方法更加清楚,下面结合具体的示例来说明,其中,该示例中,多轴原始加速度数据对应X、Y、Z原始加速度数据。
在本实施例中,如图5所示的结合执行模块的是心率检测流程,采集用户的三轴原始加速度数据后,活动量模块,通过对三轴原始加速度数据进行方差或标准差或幅度值的提取,来确定用户此前的活动量,以保证仅在判断出佩戴者处于安静状态下进行心率检测。
进而,若是活动量对应与用户处于安静状态,则将三轴原始加速度数据发送至预处理 (高通滤波)模块,对三轴原始加速度数据进行高通滤波,以去除基线漂移以及呼吸率对于心率检测的影响,举例而言,当三轴原始加速度数据为如图6左图所示时,分别对三轴原始加速度数据进行N阶的滑动平均,再分别用三轴原始加速度数据减去滑动平均后的数据以达到高通滤波的效果。
高通滤波后的加速度数据仅仅包含很强的振动信息,如图6右图,对处理后的y轴,z轴数据分别上下移动40个单位以方便观测,以获取得到三轴目标加速度数据。此时可以获知,y轴基本不包含明显的周期性信号,x轴与z轴中能够标记处代表心率振动时刻的特征点(以菱形标识)。
进而,将三轴目标加速度数据输入到数据增强模块,对三轴目标加速度数据分别进行平方处理,以增强包络信息,为后续的心率提取提供便利,主要寻找峰值的包络等信息进行处理,直接从频域进行心率结果的检测,克服了从时域寻找峰值点的不足。
在获取到数据增强模块输出的三轴目标加速度数据后,将该三轴目标加速度数据输入至数据融合模块,数据融合模块通过主成分分析(PCA,Principle Component Analysis)等技术对增强后的三轴目标加速度数据进行分离,提取出特征最明显的分量当作心率信号。如图7的左图到右图所示(图7中左图在图中没有清晰的区分出每轴的目标加速度数据,但是并不影响本申请合并过程的表达),示出了将三轴目标加速度数据合并成一轴的过程。
时频转换模块对融合后的三轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,即将时域数据转换到频域获取融合频域加速度数据来进行心率的提取。最后,心率计算模块对融合频域加速度数据响应中能量最强的频点当作心率值输出。如图8所示的融合频域加速度数据中以菱形标识标注出对应的能量最强的频点为心率值。
由此,本实施例中的不但降低了心率值提取的成本,而且基于频域中的加速度数据来提取心率值,对传感器的摆放和佩戴位置等的要求较低,因为加速数据较为敏感,因此,实用性高。
综上,本申请实施例的心率检测方法,采集用户的多轴原始加速度数据,根据多轴原始加速度数据判断用户是否满足预设检测条件,进而,若获知用户满足预设检测条件,则对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,最后,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,并根据融合频域加速度数据中的峰值数据确定用户的心率值。由此,基于多轴加速度来确定心率,降低了测量功耗,并且在频域中提取心率,提高了心率的测量准确性。
基于上述实施例,在不同的应用场景中,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据的方式不同。
在本申请的一个实施例中,如图9所示,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,包括:
步骤401,对多轴目标加速度数据进行数据处理,生成融合时域加速度数据。
在本实施例中,对多轴目标加速度数据进行数据处理,生成融合时域加速度数据,即将多轴目标加速度数据合并为一轴进行分析。
作为一种可能的实现方式,检测多轴目标加速度数据中是否至少一轴目标加速度数据具有周期性信息,比如,检测每轴目标加速度数据中每个采集点的采集目标加速度数据,确定与该采集目标加速度数据的差小于预设阈值参考采集点出现的时间点,根据该时间点确定时间间隔,如大于预设个数的采集点和参考采集点的时间间隔均一致,则认为对应轴目标加速度数据具有周期性信息。
又比如,根据每轴目标加速度数据中各个采集点绘制形状信息,若是该形状信息与预设的包络形状匹配,则认为对应轴目标加速度数据具有周期性信息。
进而,若获知至少一轴目标加速度数据具有周期性信息,则认为该目标加速度数据频域下可以提取出心率值,则分别对每轴目标加速度数据进行平方处理,进而对平方处理后的多轴目标加速度数据求和开根号处理,生成融合时域加速度数据。
在本示例中,如获知多轴目标加速度数据都不具有周期性信息,此时如采用上述方法来分别对每轴目标加速度数据进行平方处理,则反而会进一步弱化周期性信息,导致心率提取困难,因此,在本实施例中,为了保证心率提取的准确性,根据主成分分析技术提取多轴目标加速度数据的特征分量数据,生成融合时域加速度数据。其中,主成分分析技术通过正交变换将一组可能存在相关性的多轴目标加速度数据转换为一组线性不相关的变量,转换后的这组变量叫主成分(特征分量数据)。当然,在获知多轴目标加速度数据具有周期性信息时,也可以直接使用对应的主成分分析技术提取多轴目标加速度数据的特征分量数据,生成融合时域加速度数据。
步骤402,对融合时域加速度数据进行傅里叶变换处理,获取融合频域加速度数据。
具体的,对融合时域加速度数据进行傅里叶变换处理,获取融合频域加速度数据,该融合频域加速度数据基于频域反应了心率值。
在本申请的一个实施例中,可以首先对每轴目标加速度数据进行傅里叶变换处理,获取多轴频域加速度数据,进而,按照第四预设算法对多轴频域加速度数据进行数据处理,获取融合频域加速度数据。比如,第四预设算法为直接对多轴频域加速度数据求和,将求和结果作为融合频域加速度数据,又比如,第四预设算法为对多轴频域加速度数据求和后去平方值等得到对应的融合频域加速度数据。
在本示例中,如图10所示,当多轴目标加速度数据分别为X、Y、Z轴时,则分别对X、Y、Z轴目标加速度数据进行傅里叶变换处理,获取多轴频域加速度数据,进而对X、Y、Z轴频域加速度数据求和等获取到融合频域加速度数据。
在本申请的另一个实施例中,也可以对多目标加速度数据求和,对求和后的目标加速度数据来进行整体的傅里叶变化处理得到融合频域加速度数据。
综上,本申请实施例的心率检测方法,可以根据场景需要灵活的将加速度数据转换为频域,克服了时域检测不准确的问题,提高了心率值检测的精确度。
为了实现上述实施例,本申请还提出一种心率检测装置。
图11为本申请实施例提供的一种心率检测装置的结构示意图。
如图11所示,该心率检测装置包括:判断模块10、滤波处理模块20、获取模块30和确定模块40。
其中,判断模块10,用于采集用户的多轴原始加速度数据,根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件;
滤波处理模块20,用于在获知用户满足预设检测条件时,对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;
获取模块30,用于对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据;
确定模块40,用于根据融合频域加速度数据中的峰值数据确定用户的心率值。
需要说明的是,前述对心率检测方法实施例的解释说明也适用于该实施例的心率检测装置,此处不再赘述。
由此,本实施例中的不但降低了心率值提取的成本,而且基于频域中的加速度数据来提取心率值,对传感器的摆放和佩戴位置等的要求较低,因为加速数据较为敏感,因此,实用性高。
综上,本申请实施例的心率检测装置,采集用户的多轴原始加速度数据,根据多轴原始加速度数据判断用户是否满足预设检测条件,进而,若获知用户满足预设检测条件,则对多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,最后,对多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,并根据融合频域加速度数据中的峰值数据确定用户的心率值。由此,基于多轴加速度来确定心率,降低了测量功耗,并且在频域中提取心率,提高了心率的测量准确性。
基于上述实施例,在不同的应用场景中,对多轴目标加速度数据进行傅里叶变换处理, 获取融合频域加速度数据的方式不同。
在本申请的一个实施例中,如图12所示,在如图11所示的基础上,该获取模块30包括:生成单元31和获取单元32,
其中,生成单元31,用于对多轴目标加速度数据进行数据处理,生成融合时域加速度数据;
获取单元32,用于对融合时域加速度数据进行傅里叶变换处理,获取融合频域加速度数据。
在一些可能的示例中,生成单元31,具体用于:
检测多轴目标加速度数据中是否至少一轴目标加速度数据具有周期性信息;
若获知至少一轴目标加速度数据具有周期性信息,则分别对每轴目标加速度数据进行平方处理,进而对平方处理后的多轴目标加速度数据求和开根号处理,生成融合时域加速度数据;
若获知多轴目标加速度数据都不具有周期性信息,则根据主成分分析技术提取多轴目标加速度数据的特征分量数据,生成融合时域加速度数据。
在一些可能的示例中,获取单元32,具体用于:
对每轴目标加速度数据进行傅里叶变换处理,获取多轴频域加速度数据;
按照预设算法对多轴频域加速度数据进行数据处理,获取融合频域加速度数据。
需要说明的是,前述对心率检测方法实施例的解释说明也适用于该实施例的心率检测装置,此处不再赘述。
综上,本申请实施例的心率检测装置,可以根据场景需要灵活的将加速度数据转换为频域,克服了时域检测不准确的问题,提高了心率值检测的精确度。
为了实现上述实施例,本申请还提出一种计算机设备,,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时,实现如上述实施例所述的心率检测方法。
为了实现上述实施例,本申请还提出一种非临时性计算机可读存储介质,当所述存储介质中的指令由处理器被执行时,使得能够执行上述实施例所述的心率检测方法。
为了实现上述实施例,本申请还提出一种计算机程序产品,当所述计算机程序产品中的指令处理器执行时,执行上述实施例所述的心率检测方法。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本申请的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本申请的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。
应当理解,本申请的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实 施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如,如果用硬件来实现和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本申请各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (14)
- 一种心率检测方法,其特征在于,包括:采集用户的多轴原始加速度数据,根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件;若获知所述用户满足所述预设检测条件,则对所述多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;对所述多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据;根据所述融合频域加速度数据中的峰值数据确定所述用户的心率值。
- 根据权利要求1所述的方法,其特征在于,所述根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件,包括:按照第一预设算法对每轴所述原始加速度数据进行计算,以获取与所述多轴原始加速度数据对应的多轴特征数据;将每轴所述特征数据与对应轴的第一预设阈值进行比较,获取所述多轴特征数据的比较结果;若根据所述比较结果获知所述多轴特征数据均满足预设的第一检测范围,则获知所述用户满足所述预设检测条件;若根据所述比较结果获知至少一轴特征数据不满足预设的检测范围,则获知所述用户不满足所述预设检测条件。
- 根据权利要求1所述的方法,其特征在于,所述根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件,包括:对所述多轴原始加速度数据求和处理获取融合原始加速度数据;按照第二预设算法对所述融合原始加速度数据进行计算,获取融合特征数据;将所述融合特征数据与对应的第二预设阈值进行比较,获取所述融合特征数据的比较结果;若根据所述比较结果获知所述融合特征数据满足预设的第二检测范围,则获知所述用户满足所述预设检测条件;若根据所述比较结果获知所述融合特征数据不满足预设的检测范围,则获知所述用户不满足所述预设检测条件。
- 根据权利要求1至3中任一项所述的方法,其特征在于,所述对所述多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据,包括:对每轴所述原始加速度数据进行N阶的滑动平均获取对应的轴平均值,其中,N为大于 1的整数;分别用每轴所述原始加速度数据减去对应的所述轴平均值,以获取与多轴所述原始加速度数据对应的多轴目标加速度数据。
- 根据权利要求1至4中任一项所述的方法,其特征在于,在对所述多轴目标加速度数据进行傅里叶变换处理之前,还包括:按照第三预设算法对每轴所述目标加速度数据进行增强处理。
- 根据权利要求1至5中任一项所述的方法,其特征在于,所述对所述多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,包括:对所述多轴目标加速度数据进行数据处理,生成融合时域加速度数据;对所述融合时域加速度数据进行傅里叶变换处理,获取融合频域加速度数据。
- 根据权利要求6所述的方法,其特征在于,所述对所述多轴目标加速度数据进行数据处理,生成融合时域加速度数据,包括:检测所述多轴目标加速度数据中是否至少一轴目标加速度数据具有周期性信息;若获知至少一轴目标加速度数据具有周期性信息,则分别对每轴所述目标加速度数据进行平方处理,进而对平方处理后的所述多轴目标加速度数据求和开根号处理,生成融合时域加速度数据;若获知所述多轴目标加速度数据都不具有周期性信息,则根据主成分分析技术提取所述多轴目标加速度数据的特征分量数据,生成融合时域加速度数据。
- 根据权利要求6所述的方法,其特征在于,所述对所述多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速度数据,包括:对每轴所述目标加速度数据进行傅里叶变换处理,获取多轴频域加速度数据;按照第四预设算法对所述多轴频域加速度数据进行数据处理,获取融合频域加速度数据。
- 一种心率检测装置,其特征在于,包括:判断模块,用于采集用户的多轴原始加速度数据,根据所述多轴原始加速度数据判断所述用户是否满足预设检测条件;滤波处理模块,用于在获知所述用户满足所述预设检测条件时,对所述多轴原始加速度数据进行高通滤波处理获取多轴目标加速度数据;获取模块,用于对所述多轴目标加速度数据进行傅里叶变换处理,获取融合频域加速 度数据;确定模块,用于根据所述融合频域加速度数据中的峰值数据确定所述用户的心率值。
- 根据权利要求9所述的装置,其特征在于,所述获取模块,包括:生成单元,用于对所述多轴目标加速度数据进行数据处理,生成融合时域加速度数据;获取单元,用于对所述融合时域加速度数据进行傅里叶变换处理,获取融合频域加速度数据。
- 根据权利要求10所述的装置,其特征在于,所述生成单元,具体用于:检测所述多轴目标加速度数据中是否至少一轴目标加速度数据具有周期性信息;若获知至少一轴目标加速度数据具有周期性信息,则分别对每轴所述目标加速度数据进行平方处理,进而对平方处理后的所述多轴目标加速度数据求和开根号处理,生成融合时域加速度数据;若获知所述多轴目标加速度数据都不具有周期性信息,则根据主成分分析技术提取所述多轴目标加速度数据的特征分量数据,生成融合时域加速度数据。
- 根据权利要求10或11所述的装置,其特征在于,所述获取单元,具体用于:对每轴所述目标加速度数据进行傅里叶变换处理,获取多轴频域加速度数据;按照第四预设算法对所述多轴频域加速度数据进行数据处理,获取融合频域加速度数据。
- 一种计算机设备,其特征在于,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时,实现如权利要求1-8中任一所述的心率检测方法。
- 一种非临时性计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1-8中任一所述的心率检测方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011339438.4A CN114533010A (zh) | 2020-11-25 | 2020-11-25 | 心率检测方法和装置 |
CN202011339438.4 | 2020-11-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022111203A1 true WO2022111203A1 (zh) | 2022-06-02 |
Family
ID=81659325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2021/126987 WO2022111203A1 (zh) | 2020-11-25 | 2021-10-28 | 心率检测方法和装置 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114533010A (zh) |
WO (1) | WO2022111203A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114732383A (zh) * | 2022-06-13 | 2022-07-12 | 深圳市华屹医疗科技有限公司 | 体征指标监测方法、装置、设备、存储介质和程序产品 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040039420A1 (en) * | 2002-08-26 | 2004-02-26 | Medtronic Physio-Control Manufacturing Corp. | Apparatus, software, and methods for cardiac pulse detection using accelerometer data |
CN102458246A (zh) * | 2009-06-05 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 运动确定设备 |
CN106441536A (zh) * | 2016-09-07 | 2017-02-22 | 广州视源电子科技股份有限公司 | 一种箱包的自称重方法、装置与具有自称重功能的箱包 |
CN106491138A (zh) * | 2016-10-26 | 2017-03-15 | 歌尔科技有限公司 | 一种运动状态检测方法及装置 |
CN108056769A (zh) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | 一种生命体征信号分析处理方法、装置和生命体征监测设备 |
CN109381168A (zh) * | 2017-08-09 | 2019-02-26 | 三星电子株式会社 | 用户使用的电子装置和用于用户佩戴的电子装置的方法 |
CN111616695A (zh) * | 2020-06-29 | 2020-09-04 | 歌尔科技有限公司 | 一种心率获取方法、装置、系统和介质 |
CN111930230A (zh) * | 2020-07-27 | 2020-11-13 | 歌尔光学科技有限公司 | 姿态检测方法、可穿戴设备及计算机可读存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10201286B2 (en) * | 2014-08-22 | 2019-02-12 | Apple Inc. | Frequency domain projection algorithm |
US10765332B2 (en) * | 2015-04-29 | 2020-09-08 | Brainlab Ag | Detection of the heartbeat in cranial accelerometer data using independent component analysis |
CN106237604A (zh) * | 2016-08-31 | 2016-12-21 | 歌尔股份有限公司 | 可穿戴设备及利用其监测运动状态的方法 |
CN111568433A (zh) * | 2020-05-18 | 2020-08-25 | 复旦大学附属中山医院 | 一种基于三轴加速度传感器的生理和行为监测方法 |
-
2020
- 2020-11-25 CN CN202011339438.4A patent/CN114533010A/zh active Pending
-
2021
- 2021-10-28 WO PCT/CN2021/126987 patent/WO2022111203A1/zh active Application Filing
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040039420A1 (en) * | 2002-08-26 | 2004-02-26 | Medtronic Physio-Control Manufacturing Corp. | Apparatus, software, and methods for cardiac pulse detection using accelerometer data |
CN102458246A (zh) * | 2009-06-05 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 运动确定设备 |
CN106441536A (zh) * | 2016-09-07 | 2017-02-22 | 广州视源电子科技股份有限公司 | 一种箱包的自称重方法、装置与具有自称重功能的箱包 |
CN106491138A (zh) * | 2016-10-26 | 2017-03-15 | 歌尔科技有限公司 | 一种运动状态检测方法及装置 |
CN109381168A (zh) * | 2017-08-09 | 2019-02-26 | 三星电子株式会社 | 用户使用的电子装置和用于用户佩戴的电子装置的方法 |
CN108056769A (zh) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | 一种生命体征信号分析处理方法、装置和生命体征监测设备 |
CN111616695A (zh) * | 2020-06-29 | 2020-09-04 | 歌尔科技有限公司 | 一种心率获取方法、装置、系统和介质 |
CN111930230A (zh) * | 2020-07-27 | 2020-11-13 | 歌尔光学科技有限公司 | 姿态检测方法、可穿戴设备及计算机可读存储介质 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114732383A (zh) * | 2022-06-13 | 2022-07-12 | 深圳市华屹医疗科技有限公司 | 体征指标监测方法、装置、设备、存储介质和程序产品 |
Also Published As
Publication number | Publication date |
---|---|
CN114533010A (zh) | 2022-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3369375A1 (en) | Method and apparatus for identifying state of motion of human body | |
US20210267551A1 (en) | Noise detection method and apparatus | |
CN111091116A (zh) | 一种用于判断心律失常的信号处理方法及系统 | |
US9936889B2 (en) | Apparatus and method of controlling threshold for detecting peaks of physiological signals | |
US10588574B2 (en) | System and methods for adaptive noise quantification in dynamic biosignal analysis | |
CN104546007B (zh) | 胎动检测抗干扰处理方法及装置 | |
US11540748B2 (en) | Method and system for gait detection of a person | |
CN109645979A (zh) | 动态心电信号伪差识别方法及装置 | |
CN107550484B (zh) | 一种心磁信号质量评价方法及系统 | |
WO2022111203A1 (zh) | 心率检测方法和装置 | |
Salvi et al. | An optimised algorithm for accurate steps counting from smart-phone accelerometry | |
KR101510522B1 (ko) | 심전도 신호 분류 시스템 | |
WO2023061565A1 (en) | Apparatus, system and method for determining whether a person is asleep | |
TW201811261A (zh) | 訊號偵測方法 | |
KR101941172B1 (ko) | 생체 신호의 피크를 검출하는 임계값 제어 방법 및 장치. | |
CN106419884A (zh) | 一种基于小波分析的心率计算方法与系统 | |
EP4198933A1 (en) | Fall risk prevention method and device for carrying out same | |
Zhao et al. | Periodicity-based swimming performance feature extraction and parameter estimation | |
CN106326672A (zh) | 入睡检测方法与系统 | |
CN111743668B (zh) | 假肢控制方法、装置、电子设备和存储介质 | |
CN115336981A (zh) | 一种可穿戴设备佩戴状态的检测方法及相关组件 | |
CN113331824B (zh) | 一种基于激光测量的眨眼判断方法及装置 | |
CN112890830A (zh) | 一种基于睡眠脑网络的抑郁症患者数据分类方法及装置 | |
KR20180128634A (ko) | 휴대용 생체정보 측정 단말기를 이용한 R-peak 검출 방법 및 시스템 | |
CN110368020B (zh) | 一种心音信号预处理方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 21896705 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 21896705 Country of ref document: EP Kind code of ref document: A1 |