CN113349752A - 一种基于传感融合的可穿戴设备实时心率监测方法 - Google Patents
一种基于传感融合的可穿戴设备实时心率监测方法 Download PDFInfo
- Publication number
- CN113349752A CN113349752A CN202110500563.7A CN202110500563A CN113349752A CN 113349752 A CN113349752 A CN 113349752A CN 202110500563 A CN202110500563 A CN 202110500563A CN 113349752 A CN113349752 A CN 113349752A
- Authority
- CN
- China
- Prior art keywords
- heart rate
- signal
- value
- frequency
- acceleration
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012544 monitoring process Methods 0.000 title claims abstract description 26
- 230000004927 fusion Effects 0.000 title claims abstract description 20
- 230000001133 acceleration Effects 0.000 claims abstract description 76
- 238000001228 spectrum Methods 0.000 claims abstract description 43
- 230000003044 adaptive effect Effects 0.000 claims abstract description 32
- 238000001914 filtration Methods 0.000 claims abstract description 23
- 238000003066 decision tree Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000013145 classification model Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000013186 photoplethysmography Methods 0.000 claims description 45
- 230000008859 change Effects 0.000 claims description 32
- 101100166455 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) ccg-4 gene Proteins 0.000 claims description 23
- 101150114608 PPG1 gene Proteins 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000036541 health Effects 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 7
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 208000024172 Cardiovascular disease Diseases 0.000 description 4
- 208000017667 Chronic Disease Diseases 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- 210000000707 wrist Anatomy 0.000 description 3
- 230000005189 cardiac health Effects 0.000 description 2
- 230000036996 cardiovascular health Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 210000000624 ear auricle Anatomy 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004630 mental health Effects 0.000 description 1
- 230000017311 musculoskeletal movement, spinal reflex action Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009323 psychological health Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000276 sedentary effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 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
- 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
- A61B5/02444—Details of sensor
-
- 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
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
Abstract
本发明属于人体心率健康监测技术领域，具体为一种基于传感融合的可穿戴设备实时心率监测方法。通过构建多传感融合的最小二乘自适应滤波去噪算法模型以去除原始光电脉搏信号和加速度信号中的干扰因素，通过提取三轴加速度信号的峰峰值、均方根特征建立决策树分类模型，同时根据该模型定义出心率分类区间，然后将滤波器输出的PPG信号输入到分类模型中，以确定出当前心率值所在频谱区间，并在该区间查找最大频谱峰值确定心率频率。最后，融合三路去噪后的光电脉搏波信号计算得到的心率值，并将其作为最终心率值。与现有技术相比，本发明有效去除了计算过程中的噪声干扰，计算过程简单，且获得的心率值准确度更高。
Description
技术领域
本发明属于人体心率健康监测技术领域，具体为一种基于传感融合的可穿戴设备实时心率监测方法。
背景技术
随着社会城市化的快速发展，人们的工作节奏加快，各种生活压力日益增大。根据《中国心血管健康疾病越报告2019概要》显示可知，目前，由心血管健康问题引起的疾病风险最高，据推算我国现有心脑血管疾病患者人数大约为3.3亿，心血管患病率及死亡率仍处于上升阶段，且患病趋于年轻化。尽管中国心血管病防治工作取得初步成效，但仍然面临着严峻挑战。2019年新冠疫情爆发，远程也成为热门话题，合理的慢性疾病管理，能够避免看急诊和住院医疗，减少就医次数，节省费用和人力成本。心血管疾病是慢性病之首，慢性疾病往往患病时间长，需要长期的监测与管理。心率作为重要的体征信号，对慢性疾病的管理和控制、社交网络和身心健康水平的评估以及运动健康管理都起着至关重要的作用，如何实时管理和控制慢性心血管疾病，有效防范自身心脏健康风险，实时心率监测及预防心脏健康疾病成为了当前业内的一个重要研究课题。
可穿戴设备即直接穿在身上,或是整合到用户的衣服或配件的一种便携式设备。它不仅仅是一种硬件设备，更是通过软件支持以及数据交互、云端交互来实现强大的功能，利用可穿戴设备监测体征信号能够客观分析个体的身心健康状况。在心率监测方面，目前医疗中多采用单导联或多导联设备监测心电图来计算心率，受制于检测流程，无法实时获得监测结果,这在日常生活中存在不便。为改善这一情况，人们提出了光电容积脉搏波描述法(Photoplethysmography,PPG)。
光电容积脉搏波描述法(Photoplethysmography,PPG)是通过PPG传感器将光线反射到皮肤上，测量光线强度的变化，光线通过皮肤反射或透射，可置于耳垂、指尖或手腕处进行长时间无创心率监测。在实际使用中，由于受试者的自主或不自主运动，可穿戴设备采集的PPG脉搏波信号存在运动伪影(Motion Artifact,MA)干扰，而这些运动伪影将严重影响心率计算的速率和准确性，使得监测结果的实时性和准确度不够。目前针对这些运动伪影的去除，普遍采用信号分解与重构的方法，通过将PPG信号分解为若干部分后选择有效部分作为去噪后的PPG信号。然而这些方法在对信号的分解与重构过程中，均涉及到大量矩阵运算，去噪算法的复杂度较高，可穿戴设备的低算力难以实现。
发明内容
本发明的目的在于:提供一种基于传感融合的可穿戴设备实时心率监测方法，以克服运动伪影对心率计算速率和准确性的影响，提升心率监测的时效性和准确度，计算过程简单，且无需对PPG信号进行分解与重构。
为实现上述发明目的，本发明采取如下技术方案：
一种基于多传感融合的实时心率监测方法，包括以下步骤：
步骤1、获取原始光电脉搏波信号和加速度信号
利用可穿戴设备上的光电脉搏波传感器和三轴加速度传感器分别获取光电脉搏波信号PPG和加速度信号ACC；
步骤2、构建多传感融合的最小二乘自适应滤波去噪算法模型
2.1、对原始光电脉搏波信号和加速度信号预处理，以去除信号中的工频干扰、高频噪声和环境光影响；
2.2、将步骤2.1所得预处理后的光电脉搏PPG波信号和加速度ACC信号输入基于最小二乘法自适应滤波器去噪框架中：
以X轴加速度信号为最小二乘自适应滤波器输入信号，PPG信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG1信号；
以Y轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG2信号；
以Z轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG3信号；
按照上述方式融合得到三路去噪后的光电脉搏波信号，分别记为PPG1、PPG2和PPG3；
步骤3、构建基于决策树频谱选择算法
3.1、三轴加速度特征提取：从步骤2.1所得预处理后的加速度信号ACC中，取4-12s的时间窗处理，每次2s-4s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算得到时间窗内每一个信号的峰峰值Ppval＝|ymax-ymin|、均方根
其中ymax为当前时间窗内的信号最大值，ymin为当前时间窗内的信号最小值，N为时间窗内的信号数组长度，i为时间窗内对应指定信号的长度，其取值范围由时间窗的时间长度和采样频率确定；
3.2、根据步骤3.1提取的三轴加速度的特征构建决策树分类模型，并依据构建的决策树分类模型定义出分类区间；分类区间的定义为：依据静止情况下人体心率值不会大于120bpm的特点，将心率值小于等于120心跳每分钟(bpm)的为A类，大于120心跳每分钟(bpm)的为B类；根据心率计算方法，将分类区间按照频谱进行划分则A类为频谱中频率小于或等于2Hz的区间：B类为频谱中频率大于2Hz的区间；
3.3、通过离散傅立叶变换将滤波器输出的信号PPG1、PPG2和PPG3从时域数据转换到频域数据；
3.4、将步骤3.3所得频域数据输入步骤3.2构建的决策树分类模型中，预测出当前心率可查找频谱区间，并在对应频谱区间中查找最大谱峰值，若为A类，则在频谱中频率小于或等于2Hz的区间进行查找，若为B类，则在频谱中频率高于2Hz的区间进行查找。确定当前心率频率f1，并转换成当前心率值Hr'；
Hr'＝f1×60S
3.5、计算步骤3.4所得的当前心率值Hr'与上一次计算心率值Hrlast的误差，若|Hr'-Hlrast|＞5则返回步骤3.4；若|Hr'-Hrlast|≤5,则分别计算出PPG1的心率值Hr1、PPG2的心率值Hr2和PPG3的心率值Hr3，并通过
的方式进行融合得到最终心率值。
进一步的，所述构建基于决策树频谱选择算法还可以通过验证当前心率变化合理性来实现；验证当前心率变化合理性的详细过程过：
从步骤2所得加速度信号ACC中，取4-12s的时间窗处理，每次2s-4s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算每一个窗内信号的绝对平均值
计算当前时间窗内三轴加速度的绝对均值之和，确定出心率增减区间，得到预测心率变化趋势；
trend＝(mean_x+mean_y+mean_z)
Arr_Trend.append(trend)
其中，trend为预测的心率趋势变化值，mean_x,mean_y,mean_z分别表示加速度x轴，y轴，z轴的绝对平均值，计算出的当前时间窗内三轴加速度的绝对均值之和添加入心率趋势变化数组Arr_Trend中，设数组元素个数为20～30，当大于所设元素组个数时则移除首项，小于等于所设元素组个数时时依次从数组尾端添加，Arr_Trend数组中数据变化率表示为心率趋势变化；
根据预测心率变化趋势确定出心率增减区间的极值点来验证心率变化的合理性；若出现实际计算心率变化与预测心率变化趋势不一致时，则取当前心率值与上一次计算的心率值绝对误差大于5、且同时大于当前心率前4s计算的心率值数组的均值作为当前心率值。
进一步的，所述光电脉搏波信号和加速度信号预处理采取方式为：先将光电脉搏波信号进行四阶IIR带通滤波，带通频率为(0.8Hz-4Hz)，再对光电脉搏波信号进行均值滤波；将加速度信号进行二阶IIR带通滤波，带通频率为(0.5Hz-15Hz)，再对加速度信号进行均值滤波。
本发明提供的一种基于多传感融合的可穿戴实时心率监测方法，针对实际的原始光电脉搏波信号中存在运动伪影，工频干扰，高频噪声，环境光影响等严重影响了心率计算的准确性的问题；先采用了带通滤波、均值滤波进行预处理，消除工频干扰和高频噪声；而对于影响最为严重且难以消除的运动伪影；则根据加速度信号可被模拟运动伪影的特点，利用自适应相消的原理，对原始PPG信号中的运动伪影进行去除或衰减；之后将预处理后的光电脉搏波信号PPG和加速度信号输入基于最小二乘自适应滤波器去噪框架中得到三路去噪后的光电脉搏波信号，记为PPG1、PPG2和PPG3；由于不同轴向的运动伪影权重不一，因此在此过程中，采用了三种不同的融合方式以提升计算准确。
由于心率的估计常用方法是将干净的PPG信号，如PPG1、PPG2和PPG3从时域信号转换到频域信号，通过选择频谱中的最大峰值对应的频率点计算心率值。
为计算得到正确的心率值，将三路去噪后的光电脉搏波信号转化到频域，实际过程中PPG1、PPG2和PPG3信号对应的频谱最大谱峰值可能并非每次都对应实际心率频率点，通过选择合适的频谱查找区间有助于减少查找复杂度和减小计算误差。因此，本发明通过提取预处理后加速度信号特征构建决策树分类模型，并基于决策树分类模型定义出心率分类区间，然后将滤波器输出的信号输入到分类模型中，确定出当前心率值所在频谱区间，并在该区间查找最大频谱峰值确定心率频率。最后，融合三路去噪后的光电脉搏波信号计算得到的心率值。为使心率计算结果更加准确，本发明还从加速度特征提取了绝对平均值，通过计算当前时间窗内三轴加速度的绝对均值之和，确定趋势函数增减区间，得到心率变化趋势；通过验证心率变化的合理性来提升心率计算的准确性。
综上所述，本发明由于采用了上述技术方案，因此具有了以下有益效果：
(1)本发明适用于对被测者长时间实时监控，自动在可穿戴设备上进行心率计算，不用在电脑上进行这些复杂的计算，而且还具有保护个人隐私的功能。
(2)本发明计算方式简单，可嵌入实时系统中，实现了数据的采集、特征提取、噪声去除和数据的保存同步进行。
附图说明
图1是本发明基于多传感融合的可穿戴实时心率监测流程图；
图2是实施例基于最小二乘自适应滤波去噪的算法框图；
图3是实施例中静止(站立)下PPG信号及滤波后时频图；
图4是实施例中静止(静坐)下PPG信号及滤波后时频图；
图5是实施例中走路时PPG信号及滤波后时频图；
图6是实施例中跑步时PPG信号及滤波后时频图；
图7是实施例中决策树的频谱选择树状图；
图8是实施例中基于加速度特征的心率趋势预测图；
图9是可穿戴设备进行手腕处光电脉搏波信号去噪和Python仿真计算对比图；
图10为使用本发明方法与现有技术计算的心率计算结果对比图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。
在本实时例中，被测对象主要是18-40岁的人群，以可穿戴设备为载体，实时计算监测佩戴者的腕部光电脉搏波信号和加速度信号，最主要的是通过最小二乘自适应滤波器融合光电脉搏波信号和加速度信号对光电脉搏波信号的运动伪影进行去除，为准确计算到心率值，通过决策树模型对心率查找区间进行预测，在预测的心率区间内查找光电脉搏波频谱峰值，并通过预测心率趋势的变化对心率估计值进行验证。图1是本发明基于多传感融合的可穿戴实时心率监测流程图。如图1所示，本发明的可穿戴设备实时心率监测方法，包括以下步骤：
步骤1、利用可穿戴设备上的光电脉搏波传感器和三轴加速度传感器分别获取光电脉搏波(PPG)信号和加速度(ACC)信号。本实施例设计了时长为16min的自然生活状态下的连续心率监测实验，共采集数据12组。为获取真实心率值，实验者佩戴ECG设备，通过蓝牙连接手机APP获得真实心率数据。
步骤2、构建多传感融合的最小二乘自适应滤波去噪算法模型
2.1对采集的原始光电脉搏波信号和加速度信号预处理，以去除工频干、高频噪声和环境光影响；实施时，先将光电脉搏波信号进行四阶IIR带通滤波，带通频率为(0.8Hz-4Hz)，再对光电脉搏波信号进行均值滤波；将加速度信号进行二阶IIR带通滤波，带通频率为(0.5Hz-15Hz)，再对加速度信号进行均值滤波。
2.2、由于加速度信号具有被模拟运动伪影的特点，基于该特点，本身实施例将步骤2.1所得预处理后的光电脉搏PPG波信号和加速度ACC信号输入基于最小二乘发自适应滤波器去噪框架中，利用自适应相消的原理构建多路最小二乘自适应滤波去噪算法框架以消除运动伪影；如图2所示，
以X轴加速度信号为最小二乘自适应滤波器输入信号，PPG信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG1信号；
以Y轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG2信号；
以Z轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG3信号。
如图3所示，在静止下站立的PPG信号分别经本发明提出的去噪算法后分别得到PPG1、PPG2、PPG3，其中频谱图中REF为真实参考心率值的心率频率点，可知在PPG2与PPG3得到的结果中最大谱峰值的频率点与实际参考心率频率点相近；如图4所示，在静止下静坐的PPG信号滤波后得到的结果，PPG频谱图中最大谱峰值横坐标与实际心率频率点相距较远，而PPG1与PPG2的频谱图中最大谱峰值对应频率点与实际心率参考相近，可知，在静止下，通过使用三轴加速度信号采用RLS自适应滤波对于原始PPG信号滤波有效。
图5为走路下PPG信号及滤波后时频图，走路情况下的PPG信号形态出现部分畸变，且在走路下PPG信号频谱图中干扰频率点较多，左图分别为三轴加速度信号的时频图，可知在部分峰值频率点处，加速度信号与PPG信号存在相同频率点的谱峰值，存在运动伪影干扰，经三路滤波器滤波可知，如右图所示，频谱图中干扰频率点减少，并且PPG1中最大谱峰值频率点与实际参考心率点相近，可近似作为当前计算值。
如图6所示，在跑步下PPG信号受到运动伪影干扰较大，在实际心率参考值对应频率点处，由于加速度信号存在与实际心率频率分量，在经过RLS自适应滤波器后，频谱图中干扰频率点减少，实际心率点虽然在在PPG1，PPG2信号中频谱中被淹没，但在PPG3中最大谱峰值频率点与实际参考心率点相近，可近似作为当前计算值。
图9为可穿戴设备和电脑端应用本发明提出的去噪算法的结果对比图，其中图9(a)为原始信号图，图9(b)为去噪后的信号图，其中Filtered PPG-Computer表示电脑端处理的结果，Filtered PPG-Wearable表示可穿戴设备处理的结果；对比图9(a)和图9(b)可以看出，利用上述方式后，本发明能够有效去除运动伪影。
步骤3、构建基于决策树频谱选择算法，图7是本发明基于基于决策树的频谱选择树状图，从图7可以看出，本发明计算更简单。具体构建过程如下：
3.1、三轴加速度特征提取：从步骤2所得加速度信号ACC中，取8s的时间窗处理，每次2s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算得到时间窗内每一个信号的峰峰值Ppval＝|ymax-ymin|、均方根
其中ymax为当前时间窗内的信号最大值，ymin为当前时间窗内的信号最小值；N为时间窗内的信号数组长度，i为时间窗内对应指定信号的长度，i的取值范围由时间窗的时间长度和采样频率确定。
3.2、根据步骤3.1提取的三轴加速度的特征作为输入，选用CART算法构建出决策树分类模型，并依据构建的决策树分类模型定义出分类区间；分类区间的定义为：依据静止情况下人体心率值不会大于120bpm的特点，将心率值小于等于120心跳每分钟(bpm)的为A类，大于120心跳每分钟(bpm)的为B类；根据心率计算方法，将分类区间按照频谱进行划分则A类为频谱中频率小于或等于2Hz的区间：B类为频谱中频率大于2Hz的区间。具体实施时，为避免模型过拟合，设置决策树层数最大为4层，分割内部节点所需的最小样本数设置为5。
3.3、通过离散傅立叶变换将滤波器输出的信号PPG1、PPG2、PPG3从时域数据转换到频域数据。
3.4、将步骤3.3所得频域数据输入步骤3.2的决策树分类模型中，基于分类模型预测出当前心率可查找频谱区间，并在对应频谱区间中查找最大谱峰值，若为A类，则在频谱中频率小于或等于2Hz的区间进行查找，若为B类，则在频谱中频率高于2Hz的区间进行查找。确定心率频率f1，并转换成当前心率值Hr'；
Hr'＝f1×60S
3.5、计算步骤3.2.3所得的当前心率值Hr'与上一次计算心率值Hrlast的误差，若|Hr'-Hrlast|＞5则返回步骤3.4；若|Hr'-Hrlast|≤5,则分别计算PPG1、PPG2和PPG3各自对应的心率值Hr1、Hr2和Hr3，并通过
的方式进行融合得到最终心率值。
图10为使用本发明方法与现有技术计算的心率计算结果对比图。从图10可以看出，依据本发明所提供的方法计算出的心率值准确性更高。
鉴于心率对人体健康的重要性，为进一步提升心率计算的准确性，所构建的基于决策树频谱选择算法还可以通过对当前心率变化的合理性验证来获得更加准确的心率值，详述如下：
从步骤2所得加速度信号ACC中，取8s的时间窗处理，每次2s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算每一个窗内信号的绝对平均值
计算当前时间窗内三轴加速度的绝对均值之和，确定趋势函数增减区间，得到预测的心率变化趋势；
trend＝(mean_x+mean_y+mean_z)
Arr_Trend.append(trend)
其中，trend为预测的心率趋势变化值，mean_x,mean_y,mean_z分别表示加速度x轴，y轴，z轴的绝对平均值；
将计算出的trend添加入心率趋势变化数组Arr_Trend中，数组元素个数为30，当大于30则移除首项，不满30时依次从数组尾端添加；
Arr_Trend数组中数据变化率表示为心率趋势变化。
根据预测的心率变化趋势判定心率增减区间的极值点来验证心率变化的合理性。由于心率短时间内不可突变性，因此，可以结合前4s的心率变化数组平滑计算其均值作为最终心率值。也就是说，当实际计算心率增减趋势变化情况与预测得到心率增减区间趋势不一致时,则取当前心率值与上一次计算的心率值绝对误差大于5且小于当前心率前4s计算的心率值数组的均值作为最终的心率值。
为更清楚的说明通过验证当前心率合理性方法的优点，下面结合图8来进行阐述。图8是本实施例基于加速度特征的心率趋势预测图。其中True Heart Trend表示真实心率归一化后的趋势变化，Fitted Heart Trend表示预测的心率趋势变化。12组数据按上述本发明的算法得到结果如下表所示,其中AAE表示绝对误差均值，AAEP表示绝对均值误差百分比，ASTPW表示选取的8s PPG数据窗口所需的心率计算时间。Pearson correlation表示皮尔逊相关系数。
平均绝对误差(Average Absolute Err，AAE)表示，对计算的所有真实心率与估计心率之间绝对误差之和求平均：
考虑到不同实验个体的身体因素和心率变化范围有所差异，本实验引用平均绝对误差百分比对(Average Absolute Error Percentage,AAEP)，求取所有实验的预测结果和真实心率值之间的绝对误差占真实心率值的百分比的平均值作为衡量结果鲁棒性的指标：
皮尔逊相关系数主要针对有序的离散或连续性变量，考察变量之间的线性相关程度；斯皮尔曼相关系数和肯德尔相关系数则分别是针对变量的单调趋势和分类变量。
其中，BPMest表示预测心率所构成的集合，BPMtrue表示真实心率值的结合
Cov(BPMest,BPMtrue)表示俩者之间的协方差，
表示预测心率的标准差
由此可见，通过验证心率变化的合理性所计算出的心率值，其准确性更高。
以上所述仅为本申请的实施例而已，并不限制于本申请。在不脱离本发明的范围的情况下对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围内。
Claims (3)
1.一种基于多传感融合的实时心率监测方法，其特征在于：包括以下步骤：
步骤1、获取原始光电脉搏波信号和加速度信号
利用可穿戴设备上的光电脉搏波传感器和三轴加速度传感器分别获取光电脉搏波信号PPG和加速度信号ACC；
步骤2、构建多传感融合的最小二乘自适应滤波去噪算法模型
2.1、对原始光电脉搏波信号和加速度信号预处理，以去除信号中的工频干扰、高频噪声和环境光影响；
2.2、将步骤2.1所得预处理后的光电脉搏PPG波信号和加速度ACC信号输入基于最小二乘发自适应滤波器去噪框架中:
以X轴加速度信号为最小二乘自适应滤波器输入信号，PPG信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG1信号；
以Y轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG2信号；
以Z轴加速度信号为最小二乘自适应滤波器输入信号，PPG1信号为最小二乘自适应滤波器期望信号，滤波器输出得到PPG3信号；
步骤3、构建基于决策树频谱选择算法
3.1、三轴加速度特征提取：从步骤2.1预处理后的加速度信号ACC中，取4-12s的时间窗处理，每次2s-4s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算得到时间窗内每一个信号的峰峰值Ppval＝|ymax-ymin|、均方根
其中ymax为当前时间窗内的信号最大值，ymin为当前时间窗内的信号最小值,N为时间窗内的信号数组长度，i为时间窗内对应指定信号的长度，i的取值范围由时间窗的时间长度和采样频率确定；
3.2、根据步骤3.1提取的三轴加速度的特征构建决策树分类模型，并依据构建的决策树分类模型定义出分类区间；分类区间的定义为：依据静止情况下人体心率值不会大于120bpm的特点，将心率值小于等于120心跳每分钟(bpm)的为A类，大于120心跳每分钟(bpm)的为B类；根据心率计算方法，将分类区间按照频谱进行划分则A类为频谱中频率小于或等于2Hz的区间：B类为频谱中频率大于2Hz的区间；
3.3、通过离散傅立叶变换将滤波器输出的信号PPG1、PPG2、PPG3从时域数据转换到频域数据；
3.4、将步骤3.3所得频域数据输入步骤3.2的决策树分类模型中，预测出当前心率可查找频谱区间，并在对应频谱区间中查找最大谱峰值，若为A类，则在频谱中频率小于或等于2Hz的区间进行查找，若为B类，则在频谱中频率高于2Hz的区间进行查找。确定心率频率f1，并转换成当前心率值Hr'；
Hr'＝f1×60S
3.5、计算步骤3.4所得的当前心率值Hr'与上一次计算心率值Hrlast的误差，若|Hr'-Hrlast|＞5则返回步骤3.4，若|Hr'-Hrlast|≤5,则分别计算PPG1的心率值Hr1、PPG2的心率值Hr2和PPG3的心率值Hr3，并通过
的方式进行融合得到最终心率值。
2.根据权利要求1所述的一种基于多传感融合的实时心率监测方法，其特征在于：所述构建基于决策树频谱选择算法还可以通过验证当前心率变化合理性来实现；验证当前心率变化合理性的详细过程为：
从步骤2.1所得加速度信号ACC中，取4-12s的时间窗处理，每次2s-4s的数据滑窗，得到当前时间窗数据流yi，通过当前时间窗数据流yi计算每一个窗内信号的绝对平均值
计算当前时间窗内三轴加速度的绝对均值之和，确定心率增减区间，得到预测心率变化趋势；
trend＝(mean_x+mean_y+mean_z)
Arr_Trend.append(trend)
其中，trend为预测的心率趋势变化值，mean_x,mean_y,mean_z分别表示加速度x轴，y轴，z轴的绝对平均值；
计算出当的前时间窗内三轴加速度的绝对均值之和添加入心率趋势变化数组Arr_Trend中，数组元素个数为20～30，当大于所设元素组个数时则移除首项，小于等于所设元素组个数时时依次从数组尾端添加；Arr_Trend数组中数据变化率表示为心率趋势变化；
根据预测心率变化趋势确定出心率增减区间的极值点来验证心率变化的合理性；若出现实际计算心率变化与预测心率变化趋势不一致时，则取当前心率值与上一次计算的心率值绝对误差大于5、且同时大于当前心率前4s计算的心率值数组的均值作为当前心率值。
3.根据权利要求1所述的一种基于多传感融合的实时心率监测方法，其特征在于：所述光电脉搏波信号和加速度信号预处理采取方式为：
先将光电脉搏波信号进行四阶IIR带通滤波，带通频率为(0.8Hz-4Hz)，再对光电脉搏波信号进行均值滤波；将加速度信号进行二阶IIR带通滤波，带通频率为(0.5Hz-15Hz)，再对加速度信号进行均值滤波。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110500563.7A CN113349752B (zh) | 2021-05-08 | 2021-05-08 | 一种基于传感融合的可穿戴设备实时心率监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110500563.7A CN113349752B (zh) | 2021-05-08 | 2021-05-08 | 一种基于传感融合的可穿戴设备实时心率监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113349752A true CN113349752A (zh) | 2021-09-07 |
| CN113349752B CN113349752B (zh) | 2022-10-14 |
Family
ID=77525890
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110500563.7A Expired - Fee Related CN113349752B (zh) | 2021-05-08 | 2021-05-08 | 一种基于传感融合的可穿戴设备实时心率监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113349752B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113397497A (zh) * | 2021-06-10 | 2021-09-17 | 维沃移动通信有限公司 | 信号处理方法、装置和电子设备 |
| CN113854990A (zh) * | 2021-10-27 | 2021-12-31 | 青岛海信日立空调系统有限公司 | 一种心跳检测方法及装置 |
| CN114010174A (zh) * | 2021-11-05 | 2022-02-08 | 广东电网有限责任公司 | 一种心率监测方法、装置及系统 |
| CN115120217A (zh) * | 2022-08-25 | 2022-09-30 | 首都医科大学附属北京同仁医院 | 一种运动伪影去除方法、装置、存储介质及电子设备 |
| CN116269285A (zh) * | 2022-11-28 | 2023-06-23 | 电子科技大学 | 一种非接触式常态化心率变异性估计系统 |
| CN116919373A (zh) * | 2023-09-15 | 2023-10-24 | 中国地质大学（武汉） | 一种基于双通道ppg的非麻醉动物心率监测系统及方法 |
| CN114010174B (zh) * | 2021-11-05 | 2024-05-14 | 广东电网有限责任公司 | 一种心率监测方法、装置及系统 |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1354553A1 (en) * | 2002-04-15 | 2003-10-22 | Samsung Electronics Co., Ltd. | Apparatus and method for detecting heartbeat using PPG |
| US20060229519A1 (en) * | 2005-04-06 | 2006-10-12 | Konica Minolta Sensing, Inc. | Biological information processing apparatus and operation program product for the same |
| JP2011229767A (ja) * | 2010-04-28 | 2011-11-17 | Nintendo Co Ltd | 生体情報処理プログラム、生体情報処理装置、生体情報処理方法および生体情報処理システム |
| US20160089086A1 (en) * | 2014-09-26 | 2016-03-31 | Pixart Imaging Inc. | Heart rate detection module, and detection and denoising method thereof |
| CN105919584A (zh) * | 2016-06-23 | 2016-09-07 | 电子科技大学 | 用于可穿戴心率监测设备的心率估计方法及装置 |
| CN105943015A (zh) * | 2016-06-04 | 2016-09-21 | 浙江大学 | 一种主动降噪的穿戴式心率变异性监测设备 |
| CN106137167A (zh) * | 2016-07-21 | 2016-11-23 | 浙江师范大学 | 一种基于光电容积脉搏波信号的运动噪声检测方法 |
| CN106236041A (zh) * | 2016-08-23 | 2016-12-21 | 电子科技大学 | 一种实时且准确的测量心率及呼吸率的算法及系统 |
| CN106264505A (zh) * | 2016-07-21 | 2017-01-04 | 浙江师范大学 | 一种基于支持向量机的心率谱峰选择方法 |
| CN106691425A (zh) * | 2016-12-30 | 2017-05-24 | 北京工业大学 | 一种运动手环的腕部心率监测方法 |
| CN206400473U (zh) * | 2016-08-25 | 2017-08-11 | 成都正盛源科技有限公司 | 一种可穿戴社交感知系统 |
| CN107438398A (zh) * | 2015-01-06 | 2017-12-05 | 大卫·伯顿 | 移动式可穿戴的监控系统 |
| EP3406189A1 (en) * | 2017-05-25 | 2018-11-28 | Tata Consultancy Services Limited | System and method for heart rate estimation |
| CN109864713A (zh) * | 2019-04-04 | 2019-06-11 | 北京邮电大学 | 基于多通道并行滤波和谱峰加权选择算法的心率监测方法 |
| CN110169764A (zh) * | 2019-05-06 | 2019-08-27 | 上海理工大学 | 一种lms自适应滤波ppg信号心率提取方法 |
| CN110664390A (zh) * | 2019-10-11 | 2020-01-10 | 南京航空航天大学 | 基于腕带式ppg和深度学习的心率监测系统及方法 |
| CN110801214A (zh) * | 2019-11-27 | 2020-02-18 | 青岛歌尔智能传感器有限公司 | 一种心率实时检测方法及系统 |
| CN111481190A (zh) * | 2020-04-02 | 2020-08-04 | 南京润楠医疗电子研究院有限公司 | 基于双路ppg信号的稳健心率测量方法和装置 |
| US20200359919A1 (en) * | 2014-08-07 | 2020-11-19 | Apple Inc. | Motion Artifact Removal by Time Domain Projection |
-
2021
- 2021-05-08 CN CN202110500563.7A patent/CN113349752B/zh not_active Expired - Fee Related
Patent Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1354553A1 (en) * | 2002-04-15 | 2003-10-22 | Samsung Electronics Co., Ltd. | Apparatus and method for detecting heartbeat using PPG |
| US20060229519A1 (en) * | 2005-04-06 | 2006-10-12 | Konica Minolta Sensing, Inc. | Biological information processing apparatus and operation program product for the same |
| JP2011229767A (ja) * | 2010-04-28 | 2011-11-17 | Nintendo Co Ltd | 生体情報処理プログラム、生体情報処理装置、生体情報処理方法および生体情報処理システム |
| US20200359919A1 (en) * | 2014-08-07 | 2020-11-19 | Apple Inc. | Motion Artifact Removal by Time Domain Projection |
| US20160089086A1 (en) * | 2014-09-26 | 2016-03-31 | Pixart Imaging Inc. | Heart rate detection module, and detection and denoising method thereof |
| CN107438398A (zh) * | 2015-01-06 | 2017-12-05 | 大卫·伯顿 | 移动式可穿戴的监控系统 |
| CN105943015A (zh) * | 2016-06-04 | 2016-09-21 | 浙江大学 | 一种主动降噪的穿戴式心率变异性监测设备 |
| CN105919584A (zh) * | 2016-06-23 | 2016-09-07 | 电子科技大学 | 用于可穿戴心率监测设备的心率估计方法及装置 |
| CN106137167A (zh) * | 2016-07-21 | 2016-11-23 | 浙江师范大学 | 一种基于光电容积脉搏波信号的运动噪声检测方法 |
| CN106264505A (zh) * | 2016-07-21 | 2017-01-04 | 浙江师范大学 | 一种基于支持向量机的心率谱峰选择方法 |
| CN106236041A (zh) * | 2016-08-23 | 2016-12-21 | 电子科技大学 | 一种实时且准确的测量心率及呼吸率的算法及系统 |
| CN206400473U (zh) * | 2016-08-25 | 2017-08-11 | 成都正盛源科技有限公司 | 一种可穿戴社交感知系统 |
| CN106691425A (zh) * | 2016-12-30 | 2017-05-24 | 北京工业大学 | 一种运动手环的腕部心率监测方法 |
| EP3406189A1 (en) * | 2017-05-25 | 2018-11-28 | Tata Consultancy Services Limited | System and method for heart rate estimation |
| CN108926333A (zh) * | 2017-05-25 | 2018-12-04 | 塔塔咨询服务有限公司 | 用于心率估计的系统和方法 |
| CN109864713A (zh) * | 2019-04-04 | 2019-06-11 | 北京邮电大学 | 基于多通道并行滤波和谱峰加权选择算法的心率监测方法 |
| CN110169764A (zh) * | 2019-05-06 | 2019-08-27 | 上海理工大学 | 一种lms自适应滤波ppg信号心率提取方法 |
| CN110664390A (zh) * | 2019-10-11 | 2020-01-10 | 南京航空航天大学 | 基于腕带式ppg和深度学习的心率监测系统及方法 |
| CN110801214A (zh) * | 2019-11-27 | 2020-02-18 | 青岛歌尔智能传感器有限公司 | 一种心率实时检测方法及系统 |
| CN111481190A (zh) * | 2020-04-02 | 2020-08-04 | 南京润楠医疗电子研究院有限公司 | 基于双路ppg信号的稳健心率测量方法和装置 |
Non-Patent Citations (14)
| Title |
|---|
| AYDEMIR, TUGBA；SAHIN, MEHMET；AYDEMIR, ONDER: "A New Method for Activity Monitoring Using Photoplethysmography Signals Recorded by Wireless Sensor", 《JOURNAL OF MEDICAL AND BIOLOGICAL ENGINEERING》 * |
| GU, KAIYUN；GUAN, ZHIYUAN；CHANG, YU: "Hemodynamic effects of pulsatile unloading of left ventricular assist devices (LVAD) on intraventricular flow and ventricular stress", 《JOURNAL OF BIOMECHANICS》 * |
| ZHANG, ZHILIN: "Photoplethysmography-Based Heart Rate Monitoring in Physical Activities via Joint Sparse Spectrum Reconstruction", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
| 付光明,穆平安,蒋睿杰.: "基于级联型LMS自适应滤波PPG信号心率估算", 《电子测量技术》 * |
| 庄海洋,左熹,陈国兴.: "软弱地基上地铁区间隧道周围场地位移和加速度反应分析", 《第九届全国岩土力学数值分析与解析方法讨论会论文集》 * |
| 徐礼胜,靳雁冰,王琦文,李锡勇,印重.: "多传感器融合的穿戴式心率监测系统", 《哈尔滨工业大学学报》 * |
| 徐菁,朱敏: "基于三轴加速度信号的心率修正算法", 《计算机应用与软件》 * |
| 杨刚,薛挺,高伟: "一种抗运动干扰的实时心率提取方法", 《计算机工程与应用》 * |
| 沈傲立: "基于带状分布的多个加速度传感器对人体基本运动状态的区分", 《上海体育学院》 * |
| 董金陇,杨录,刘康驰,高斌.: "脉搏信号伪基线去除方法研究", 《中国测试》 * |
| 陈暘: "基于PPG的去运动伪影及心率估计方法研究", 《万方》 * |
| 陈炎,李丹,李彦海,马皛源,魏建明.: "基于加速度传感器的心率信号处理及检测方法", 《科学技术与工程》 * |
| 陈龙聪,刘改琴,刘亚涛: "基于STM32和USB接口的开放式高性能新型医用信号仪的研制", 《中国生物医学工程学报》 * |
| 马梦瑶: "心率监测系统中噪声干扰消除技术的研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113397497A (zh) * | 2021-06-10 | 2021-09-17 | 维沃移动通信有限公司 | 信号处理方法、装置和电子设备 |
| WO2022257909A1 (zh) * | 2021-06-10 | 2022-12-15 | 维沃移动通信有限公司 | 信号处理方法、装置和电子设备 |
| CN113854990A (zh) * | 2021-10-27 | 2021-12-31 | 青岛海信日立空调系统有限公司 | 一种心跳检测方法及装置 |
| CN114010174A (zh) * | 2021-11-05 | 2022-02-08 | 广东电网有限责任公司 | 一种心率监测方法、装置及系统 |
| CN114010174B (zh) * | 2021-11-05 | 2024-05-14 | 广东电网有限责任公司 | 一种心率监测方法、装置及系统 |
| CN115120217A (zh) * | 2022-08-25 | 2022-09-30 | 首都医科大学附属北京同仁医院 | 一种运动伪影去除方法、装置、存储介质及电子设备 |
| CN115120217B (zh) * | 2022-08-25 | 2022-11-25 | 首都医科大学附属北京同仁医院 | 一种运动伪影去除方法、装置、存储介质及电子设备 |
| CN116269285A (zh) * | 2022-11-28 | 2023-06-23 | 电子科技大学 | 一种非接触式常态化心率变异性估计系统 |
| CN116919373A (zh) * | 2023-09-15 | 2023-10-24 | 中国地质大学（武汉） | 一种基于双通道ppg的非麻醉动物心率监测系统及方法 |
| CN116919373B (zh) * | 2023-09-15 | 2023-12-19 | 中国地质大学（武汉） | 一种基于双通道ppg的非麻醉动物心率监测系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113349752B (zh) | 2022-10-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113349752B (zh) | 一种基于传感融合的可穿戴设备实时心率监测方法 | |
| Venkatesan et al. | ECG signal preprocessing and SVM classifier-based abnormality detection in remote healthcare applications | |
| Sun et al. | Photoplethysmography-based heart rate monitoring using asymmetric least squares spectrum subtraction and bayesian decision theory | |
| CN108478206B (zh) | 运动状态下基于脉搏波的心率监测方法 | |
| Shyam et al. | PPGnet: Deep network for device independent heart rate estimation from photoplethysmogram | |
| Yang et al. | A novel adaptive spectrum noise cancellation approach for enhancing heartbeat rate monitoring in a wearable device | |
| Chung et al. | Deep learning for heart rate estimation from reflectance photoplethysmography with acceleration power spectrum and acceleration intensity | |
| Song et al. | Nas-ppg: Ppg-based heart rate estimation using neural architecture search | |
| CN105249925B (zh) | 一种中医脉象采集装置及降噪系统和降噪方法 | |
| CN105997043B (zh) | 一种基于腕式可穿戴设备的脉率提取方法 | |
| KR102043376B1 (ko) | 심층 신경망 알고리즘을 이용하는 실시간 스트레스 분석 방법 | |
| Azar et al. | Deep recurrent neural network-based autoencoder for photoplethysmogram artifacts filtering | |
| Gupta et al. | Higher order derivative-based integrated model for cuff-less blood pressure estimation and stratification using PPG signals | |
| CN112022126A (zh) | 一种基于CNN-BiLSTM模型和PPG的智能血压预测方法 | |
| Chang et al. | Deepheart: accurate heart rate estimation from PPG signals based on deep learning | |
| Kumar et al. | STSR: spectro-temporal super-resolution analysis of a reference signal less photoplethysmogram for heart rate estimation during physical activity | |
| Wan et al. | Combining parallel adaptive filtering and wavelet threshold denoising for photoplethysmography-based pulse rate monitoring during intensive physical exercise | |
| CN114569096A (zh) | 一种基于视频流的非接触式连续血压测量方法及系统 | |
| CN105816165A (zh) | 一种实时动态心率监测装置及监测方法 | |
| Liu et al. | Power and area efficient wavelet-based on-chip ECG processor for WBAN | |
| Guo et al. | An effective photoplethysmography heart rate estimation framework integrating two-level denoising method and heart rate tracking algorithm guided by finite state machine | |
| CN112274145A (zh) | 近红外脑功能成像数据的处理方法、装置和存储介质 | |
| Sivanjaneyulu et al. | Cnn based ppg signal quality assessment using raw ppg signal for energy-efficient ppg analysis devices in internet of medical things | |
| Tu et al. | Fitbeat: a lightweight system for accurate heart rate measurement during exercise | |
| KR101831064B1 (ko) | Ppg 신호를 이용한 동잡음 제거 장치 및 그 방법 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20221014 |