CN105433931A - Processing device and method for describing waveform by light volume change - Google Patents
Processing device and method for describing waveform by light volume change Download PDFInfo
- Publication number
- CN105433931A CN105433931A CN201510565554.0A CN201510565554A CN105433931A CN 105433931 A CN105433931 A CN 105433931A CN 201510565554 A CN201510565554 A CN 201510565554A CN 105433931 A CN105433931 A CN 105433931A
- Authority
- CN
- China
- Prior art keywords
- volume
- signal
- digital filtering
- describe
- smooth change
- 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.)
- Pending
Links
- 230000008859 change Effects 0.000 title claims abstract description 77
- 238000012545 processing Methods 0.000 title abstract description 16
- 238000000034 method Methods 0.000 title abstract description 8
- 238000001914 filtration Methods 0.000 claims abstract description 80
- 230000003044 adaptive effect Effects 0.000 claims abstract description 62
- 230000004044 response Effects 0.000 claims description 19
- 239000008280 blood Substances 0.000 claims description 17
- 210000004369 blood Anatomy 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000003672 processing method Methods 0.000 claims description 6
- 230000002452 interceptive effect Effects 0.000 claims description 5
- 230000005520 electrodynamics Effects 0.000 claims 2
- 230000006698 induction Effects 0.000 claims 2
- 230000035479 physiological effects, processes and functions Effects 0.000 claims 2
- 238000010606 normalization Methods 0.000 claims 1
- 230000000737 periodic effect Effects 0.000 abstract description 21
- 238000001514 detection method Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 21
- 230000001121 heart beat frequency Effects 0.000 description 9
- 238000000605 extraction Methods 0.000 description 7
- 210000000707 wrist Anatomy 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开了一种光体积变化描述波形的处理装置及其方法,该处理装置包含由一光体积变化描述波形(PPG)信号撷取单元撷取PPG信号,并于撷取PPG信号期间,以一运算信息撷取单元撷取运动信息,将该PPG信号及运动信息一并输出至一第一可适性数字滤波单元,将运动伪影的干扰信号自该PPG信号中滤除,再将滤除后的PPG信号输出至一第二可适性数字滤波单元,以撷取该第二PPG信号中的周期性信号并提高周期性信号的信噪比,作为读出精确的心跳信号的依据。因此,本发明的PPG信号处理装置可应用于穿戴式心跳检测装置上,提供正确的心跳信息。
The invention discloses a processing device and method for a light volume change description waveform. The processing device includes a light volume change description waveform (PPG) signal acquisition unit that acquires a PPG signal, and during the acquisition of the PPG signal, An operation information acquisition unit acquires motion information, outputs the PPG signal and motion information to a first adaptive digital filtering unit, filters out interference signals of motion artifacts from the PPG signal, and then filters the The divided PPG signal is output to a second adaptive digital filtering unit to capture the periodic signal in the second PPG signal and improve the signal-to-noise ratio of the periodic signal as a basis for reading an accurate heartbeat signal. Therefore, the PPG signal processing device of the present invention can be applied to a wearable heartbeat detection device to provide correct heartbeat information.
Description
技术领域technical field
本发明关于一种光体积变化描述波形的处理装置及其方法,尤指一种应用于穿戴式心跳检测装置的光体积变化描述波形的处理装置及其方法。The present invention relates to a processing device and method for describing waveforms of light volume changes, in particular to a processing device and method for describing waveforms of light volume changes applied to wearable heartbeat detection devices.
背景技术Background technique
穿戴式电子装置通常提供有监测使用者运动功能(如计步功能),或提供检测及记录人体生理特征(如心跳、血压或血氧)等功能,因此颇受消费者喜爱。然而,穿戴式电子装置的心跳检测结果仍较一般固定式血压计检测心跳结果差,主要原因在于该穿戴式电子装置于撷取人体的生理特征信号时,使用者的肢体活动会对该生理特征信号产生干扰信号。Wearable electronic devices usually provide the function of monitoring the user's movement (such as the step counting function), or provide the function of detecting and recording the physiological characteristics of the human body (such as heartbeat, blood pressure or blood oxygen), so they are very popular among consumers. However, the heartbeat detection results of wearable electronic devices are still worse than those of general fixed sphygmomanometers. The signal generates an interfering signal.
一般会将撷取到的生理特征信号自时域波形(如图9A所示)转换至频域波形(如图9B所示),再以频域中的周期信号作为量测心跳频率的依据,如量测到正确心跳频率即如图9D所示。然而,但在频域信号中的除了心跳信号Sheart外,如图9B所示,还包含有因手部摆动与身体晃动造成的低频信号,或包含因手部快速摆动与身体快速晃动所造成的高频信号,当其中部份的低频或高频信号的能量会高于心跳信号的能量。由于频率信号中包含有运动干扰信号Smotion(或称运动伪影;MotionArtifact),会被误测为心跳信号,因此检测出的心跳信号转Sheart换至时域信号后,即会出现非心跳信号Smotion,如图9C所示。Generally, the extracted physiological characteristic signal is converted from the time-domain waveform (as shown in FIG. 9A ) to the frequency-domain waveform (as shown in FIG. 9B ), and then the periodic signal in the frequency domain is used as the basis for measuring the heartbeat frequency. If the correct heartbeat frequency is measured, it will be as shown in FIG. 9D . However, in addition to the heartbeat signal S heart in the frequency domain signal, as shown in Figure 9B, it also contains low-frequency signals caused by hand swings and body shaking, or includes low-frequency signals caused by rapid hand swings and body shaking The high-frequency signal, when part of the energy of the low-frequency or high-frequency signal will be higher than the energy of the heartbeat signal. Since the frequency signal contains motion interference signal S motion (or motion artifact; MotionArtifact ), it will be misdetected as a heartbeat signal, so after the detected heartbeat signal is transferred from Shear to time domain signal, non-heartbeat will appear The signal S motion is shown in FIG. 9C .
因此,目前穿戴式电子装置检测心跳的准确度有必要进一步改善之。Therefore, it is necessary to further improve the accuracy of heartbeat detection by current wearable electronic devices.
发明内容Contents of the invention
有鉴于穿戴式电子装置检测生理信号易受到运动伪影的干扰降低其检测准确度,本发明主要目的系提供一种光体积变化描述波形(PPG)的处理装置及其方法,有助于提高生理信号的量测准确度。In view of the fact that the detection of physiological signals by wearable electronic devices is susceptible to the interference of motion artifacts and reduces its detection accuracy, the main purpose of the present invention is to provide a processing device and method for describing a photovolume change description waveform (PPG), which helps to improve physiological Signal measurement accuracy.
欲达上述目的所使用的主要技术手段系令该光体积变化描述波形(PPG)的处理装置包含有:The main technical means used to achieve the above purpose is to make the processing device of the light volume change description waveform (PPG) include:
一光体积变化描述波形(PPG)信号撷取单元,系用以撷取第一PPG信号;A light volume change description waveform (PPG) signal acquisition unit is used to extract the first PPG signal;
一运动信息撷取单元,系用以撷取一运动信息;A motion information extraction unit is used to retrieve a motion information;
一第一可适性数字滤波单元,系连接至该PPG信号撷取单元及该运动信息撷取单元,以接收该第一PPG信号及该运动信息,并依据该运动信息,将该第一PPG信号所包含的运动伪影干扰信号予以消除,并输出第二光体积变化描述波形(PPG)信号;A first adaptive digital filtering unit is connected to the PPG signal acquisition unit and the motion information acquisition unit to receive the first PPG signal and the motion information, and according to the motion information, the first PPG The motion artifact interference signal contained in the signal is eliminated, and the second photovolume change description waveform (PPG) signal is output;
一第二可适性数字滤波单元,系连接至该第一可适性数字滤波单元,以接收该第二PPG信号,并撷取该第二PPG信号中的周期性信号,并提高该周期性信号的信噪比后输出一生理特征信号。A second adaptive digital filtering unit is connected to the first adaptive digital filtering unit to receive the second PPG signal, extract the periodic signal in the second PPG signal, and improve the periodicity output a physiological characteristic signal after the signal-to-noise ratio of the signal.
欲达上述目的所使用的主要技术手段系令该光体积变化描述波形(PPG)的处理方法包含有:The main technical means used to achieve the above-mentioned purpose is to make the processing method of the light volume change description waveform (PPG) include:
(a)获得一第一光体积变化描述波形(PPG)信号以及运动信息;(a) Obtaining a first photovolume description waveform (PPG) signal and motion information;
(b)依据该运动信息,将该第一PPG信号所包含运动伪影干扰信号予以消除,并输出一第二PPG信号;以及(b) canceling the motion artifact interference signal included in the first PPG signal according to the motion information, and outputting a second PPG signal; and
(c)撷取该第二PPG信号中的一周期性信号,并提高该周期性信号的信噪比后输出一该生理特征信号。(c) extracting a periodic signal in the second PPG signal, and outputting a physiological characteristic signal after increasing the signal-to-noise ratio of the periodic signal.
由上述说明可知,本发明主要于撷取PPG信号期间,以一运算信息撷取单元撷取运动信息,由于自人体撷取到的PPG信号会包含有运动伪影的干扰信号,故该第一可适性数字滤波单元依据运动信息,将大部份的运动伪影干扰信号自该PPG信号中滤除,再将滤除后的PPG信号输出至一第二可适性数字滤波单元,以撷取该第二PPG信号中的周期性信号并提高周期性信号的信噪比,作为读出精确的心跳信号的依据。因此,本发明的PPG信号处理装置可应用于穿戴式心跳检测装置上,提供正确的心跳信息。It can be seen from the above description that the present invention mainly uses an operation information extraction unit to extract motion information during the acquisition of PPG signals. Since the PPG signals extracted from the human body will contain interference signals with motion artifacts, the first The adaptive digital filtering unit filters most of the motion artifact interference signals from the PPG signal according to the motion information, and then outputs the filtered PPG signal to a second adaptive digital filtering unit to extract Taking the periodic signal in the second PPG signal and improving the signal-to-noise ratio of the periodic signal is used as the basis for reading out the accurate heartbeat signal. Therefore, the PPG signal processing device of the present invention can be applied to a wearable heartbeat detection device to provide correct heartbeat information.
附图说明Description of drawings
图1为本发明PPG信号处理装置的第一较佳实施例的功能方块图。FIG. 1 is a functional block diagram of the first preferred embodiment of the PPG signal processing device of the present invention.
图2为本发明应用于一穿戴式电子装置的使用示意图。FIG. 2 is a schematic diagram of the application of the present invention to a wearable electronic device.
图3A至图3G为图1部份功能方块的输出信号波形图。3A to 3G are output signal waveform diagrams of some functional blocks in FIG. 1 .
图4A为本发明PPG信号处理装置的第二较佳实施例的一部份功能方块图。FIG. 4A is a partial functional block diagram of the second preferred embodiment of the PPG signal processing device of the present invention.
图4B为本发明PPG信号处理装置的第二较佳实施例的另一部份功能方块图。FIG. 4B is another partial functional block diagram of the second preferred embodiment of the PPG signal processing device of the present invention.
图5A为本发明图1的第一可适性数字滤波单元的一功能方块图。FIG. 5A is a functional block diagram of the first adaptive digital filtering unit in FIG. 1 of the present invention.
图5B为本发明图4A的第一可适性数字滤波单元的一功能方块图。FIG. 5B is a functional block diagram of the first adaptive digital filtering unit in FIG. 4A of the present invention.
图6为本发明图1的第二可适性数字滤波单元的一功能方块图。FIG. 6 is a functional block diagram of the second adaptive digital filtering unit in FIG. 1 of the present invention.
图7A为图3A及图3D转换为一频域波形图。FIG. 7A is a waveform diagram converted from FIG. 3A and FIG. 3D into a frequency domain.
图7B为图3A及图3G转换为一频域波形图。FIG. 7B is a frequency-domain waveform diagram converted from FIG. 3A and FIG. 3G .
图7C为图3A、图3D及图3G转换为一频域波形图。FIG. 7C is a frequency-domain waveform diagram converted from FIG. 3A , FIG. 3D and FIG. 3G .
图8为依据图7C中转换图3G的频域波形图,获得心跳频率的时域波形图。FIG. 8 is a time domain waveform diagram of the heartbeat frequency obtained by converting the frequency domain waveform diagram of FIG. 3G in FIG. 7C .
图9A为量测人体之一生理特征信号的时域波形图。FIG. 9A is a time-domain waveform diagram of measuring a physiological characteristic signal of a human body.
图9B为图9A转换为一频域波形图。FIG. 9B is a waveform diagram converted from FIG. 9A into a frequency domain.
图9C为依据图9B获得心跳频率的时域波形图。FIG. 9C is a time-domain waveform diagram of the heartbeat frequency obtained according to FIG. 9B .
图9D为正常心跳频率的一时域波图。FIG. 9D is a time-domain wave diagram of a normal heartbeat frequency.
其中,附图标记:Among them, reference signs:
10PPG信号撷取单元11第一带通滤波器10PPG signal extraction unit 11 first bandpass filter
12第一正规化器20运动信息撷取单元12 first normalizer 20 motion information acquisition unit
20a三轴动力感应器21第二直流准位调整单元20a Three-axis power sensor 21 Second DC level adjustment unit
21a第二正规化器22延迟电路21a second normalizer 22 delay circuit
30第一可适性数字滤波单元30 The first adaptive digital filter unit
301第一有限脉冲响应数字滤波单元301 first finite impulse response digital filter unit
302减法器303权重调整单元302 subtractor 303 weight adjustment unit
31第一误差收敛系数自动调整单元32去极端值单元31 First error convergence coefficient automatic adjustment unit 32 Extreme value removal unit
33第二带通滤波器34第三正规化器33 second bandpass filter 34 third normalizer
40第二可适性数字滤波单元40 second adaptive digital filtering unit
401第二有限脉冲响应数字滤波单元401 second finite impulse response digital filter unit
402减法器403权重调整单元402 subtractor 403 weight adjustment unit
41第二误差收敛系数自动调整单元41 The second error convergence coefficient automatic adjustment unit
50手腕51血管50 wrist 51 blood vessels
60穿戴式电子装置60 Wearable Electronic Devices
具体实施方式detailed description
本发明系主要针对一种反应人体生理信号的光体积变化描述波形(PPG)信号,并对其进行的处理装置及其方法,以下以实施例详加说明本发明PPG信号的处理装置及方法的技术内容。The present invention is mainly aimed at a light volume change description waveform (PPG) signal that reflects the physiological signal of a human body, and a processing device and method thereof. The following examples illustrate the PPG signal processing device and method of the present invention in detail technical content.
首先请参阅图1所示,为本发明PPG信号处理装置的一较佳实施例的功能方块图,其包含有一PPG信号撷取单元10、一运动信息撷取单元20、一第一可适性数字滤波单元30及一第二可适性数字滤波单元40;其中该第一可适性数字滤波单元30电连接至该PPG信号撷取单元10及该运动信息撷取单元20,该第二可适性数字滤波单元40电连接至该第一可适性数字滤波单元30。First please refer to Fig. 1, which is a functional block diagram of a preferred embodiment of the PPG signal processing device of the present invention, which includes a PPG signal acquisition unit 10, a motion information acquisition unit 20, a first adaptability Digital filter unit 30 and a second adaptive digital filter unit 40; wherein the first adaptive digital filter unit 30 is electrically connected to the PPG signal acquisition unit 10 and the motion information acquisition unit 20, the second adaptive The adaptive digital filtering unit 40 is electrically connected to the first adaptive digital filtering unit 30 .
请配合参阅图2,上述PPG信号撷取单元10若应用于一种如电子手表的穿戴式电子装置60,可采用红外光元件、绿光元件、红光元件或激光元件等发光元件,藉由照射手腕50的血管51获得对应该血管51收缩的一第一PPG信号PPG(如图3A所示)。但,于检测手腕同时人体活动非为静态,故该第一PPG信号PPG包含有运动伪影的干扰信号,即如图7A所示,当图3A的第一PPG信号PPG经快速傅立叶转换(FFT)至一频域波形图,频域波形中包含二个明显波峰,即分别为心跳信号Sheart以及运动伪影的干扰信号Smotion。Please refer to FIG. 2, if the above-mentioned PPG signal acquisition unit 10 is applied to a wearable electronic device 60 such as an electronic watch, light-emitting elements such as infrared light elements, green light elements, red light elements or laser elements can be used. The blood vessel 51 of the wrist 50 is irradiated to obtain a first PPG signal PPG corresponding to the contraction of the blood vessel 51 (as shown in FIG. 3A ). However, the human body activity is not static when the wrist is detected, so the first PPG signal PPG contains interference signals with motion artifacts, that is, as shown in FIG. 7A, when the first PPG signal PPG in FIG. 3A is fast Fourier transformed (FFT ) to a frequency-domain waveform diagram, the frequency-domain waveform contains two obvious peaks, which are respectively the heartbeat signal S heart and the interference signal S motion of motion artifacts.
如图1所示,当该第一PPG信号PPG应用于检测人体的特定生理信号(如心跳)时,由于人类的心跳信号有特定的频率范围,因此可先将该PPG信号撷取单元10所输出的第一PPG信号PPG输出至一第一带通滤波器11,该第一带通滤波器11会自该第一PPG信号PPG中滤除心跳信号的频率范围以外的噪声,再将经带通滤波过的第一PPG信号PPG_F(如图3B所示)再输出至一第一正规化器12,调整经带通滤波过的该第一PPG信号PPG_F的尺度范围后的第一PPG信号PPG_C(如图3C所示),再输出至该第一可适性数字滤波单元30。此外,该第一正规化器12亦可使用一第一直流准位调整单元等信号调整元件,以调整经带通滤波过的该第一PPG信号PPG_F的直流准位。As shown in Figure 1, when the first PPG signal PPG is used to detect specific physiological signals (such as heartbeat) of the human body, since the human heartbeat signal has a specific frequency range, the PPG signal acquisition unit 10 can first The output first PPG signal PPG is output to a first band-pass filter 11, and the first band-pass filter 11 will filter out the noise outside the frequency range of the heartbeat signal from the first PPG signal PPG, and then pass the band-pass filter The filtered first PPG signal PPG_F (as shown in FIG. 3B ) is then output to a first normalizer 12, and the first PPG signal PPG_C after adjusting the scale range of the band-pass filtered first PPG signal PPG_F (as shown in FIG. 3C ), and then output to the first adaptive digital filtering unit 30 . In addition, the first normalizer 12 can also use a signal adjustment element such as a first DC level adjustment unit to adjust the DC level of the band-pass filtered first PPG signal PPG_F.
如图1示,上述运动信息撷取单元20用以撷取该手腕的运动信息,特别是于撷取上述第一PPG信号PPG的期间内同时撷取该运动信息M,该运动信息撷取单元20输出至一第二直流准位调整单元21及一延迟电路22,以调整直流准位并延迟一段时间后,再输出至该第一可适性数字滤波单元30。在本实施例中,再配合图4A所示,该运动信息撷取单元20采用一个三轴动力感应器20a,故该运动信息M包含有第一轴、第二轴及第三轴运动信息X、Y、Z。此外,该三轴动力感应器20a所输出的第一轴、第二轴及第三轴运算信息X、Y、Z可较佳地依序分别输出至一第二直流准位调整单元21及一延迟电路22,以调整其尺度范围直流准位并延迟一段时间后,再输出至该第一可适性数字滤波单元30。此外,该第二直流准位调整单元21及一延迟电路22可以一第二正规化器21a取代之,即如图4B所示,该三轴动力感应器20a所输出的第一轴、第二轴及第三轴运算信息X、Y、Z可较佳地依序分别输出至该第二正规化器21a,以调整其尺度范围,再输出至该第一可适性数字滤波单元30。As shown in FIG. 1, the above-mentioned motion information extraction unit 20 is used to capture the motion information of the wrist, especially to simultaneously capture the motion information M during the period of capturing the first PPG signal PPG. The motion information extraction unit 20 is output to a second DC level adjustment unit 21 and a delay circuit 22 to adjust the DC level and delay for a period of time before outputting to the first adaptive digital filter unit 30 . In this embodiment, as shown in FIG. 4A , the motion information acquisition unit 20 uses a three-axis dynamic sensor 20a, so the motion information M includes the motion information X of the first axis, the second axis, and the third axis. , Y, Z. In addition, the first-axis, second-axis, and third-axis calculation information X, Y, and Z output by the three-axis power sensor 20a can be preferably sequentially output to a second DC level adjustment unit 21 and a The delay circuit 22 adjusts the DC level of the scale range and delays for a period of time before outputting to the first adaptive digital filter unit 30 . In addition, the second DC level adjustment unit 21 and a delay circuit 22 can be replaced by a second normalizer 21a, that is, as shown in FIG. 4B, the first axis, second The operation information X, Y, and Z of the axis and the third axis can preferably be respectively output to the second normalizer 21 a in order to adjust its scale range, and then output to the first adaptive digital filter unit 30 .
如图1所示的第一可适性数字滤波单元30获得经带通滤波过与直流准位或尺度范围调整过的该第一PPG信号PPG_C,以及经过直流准位或尺度范围调整后的运动信息M_C,该第一可适性数字滤波单元可依据该运动信息M_C,将该第一PPG信号PPG_C所包含的大部份运动伪影干扰信号予以消除,并输出一第二PPG信号Error_C。在本实施例,该第一可适性数字滤波单元30可较佳地使用最小平方演算法,其成本函数为:Error_C(k)=PPG_C(k)-m(k);其中Error_C(k)为本次的误差值;PPG_C(k)为目前第一PPG信号;m(k)为目前运动信息经数字滤波后的数值。该最小平方演算法可将该第一可适性数字滤波单元的误差值最小化,以消除该第一PPG信号中的大部份运动伪影干扰信号。The first adaptive digital filtering unit 30 shown in FIG. 1 obtains the first PPG signal PPG_C that has been band-pass filtered and adjusted with the DC level or scale range, and the motion after the DC level or scale range adjustment. Information M_C, the first adaptive digital filter unit can eliminate most of the motion artifact interference signals contained in the first PPG signal PPG_C according to the motion information M_C, and output a second PPG signal Error_C. In this embodiment, the first adaptive digital filtering unit 30 may preferably use the least square algorithm, and its cost function is: Error_C(k)=PPG_C(k)-m(k); where Error_C(k) is the current error value; PPG_C(k) is the current first PPG signal; m(k) is the value of the current motion information after digital filtering. The least square algorithm can minimize the error value of the first adaptive digital filter unit, so as to eliminate most of the motion artifact interference signals in the first PPG signal.
又该第一可适性数字滤波单元30实现上述最小平方演算法的架构,如图1及图5A所示,其包含有一第一有限脉冲响应数字滤波单元301、一减法器302及一权重调整单元303。其中该第一有限脉冲响应数字滤波单元301连接至该运动信息撷取单元20,将依序所接收的该运动信息m予以数字滤波后输出至该减法器302,由该减法器302将该第一PPG信号PPG_C与经数字滤波后的该运动信息M_C予以相减,以获得该误差值Error_C。由于该权重调整单元303连接至该第一有限脉冲响应数字滤波单元301及该减法器302,故可依据该第一PPG信号PPG_C与经数字滤波后的该运动信息M_C予以相减后的误差值Error_C,来调整该第一有限脉冲响应数字滤波单元301于每次数字滤波的权重值W,直到误差值最小化后输出该第二PPG信号Error_C,即如图3D所示。再请配合参阅图7A所示,为图3D的第二PPG信号Error_C转换至频域波形,与图3A同样转换至频域波形相比较后可知,运动伪影的干扰信号Smotion(1)已有效被抑制,波形中的最高峰即为心跳信号Sheart(1)。And this first adaptability digital filtering unit 30 realizes the framework of above-mentioned least square algorithm, as shown in Figure 1 and Figure 5A, it comprises a first finite impulse response digital filtering unit 301, a subtractor 302 and a weight adjustment Unit 303. Wherein the first finite impulse response digital filter unit 301 is connected to the motion information extraction unit 20, digitally filters the motion information m received in sequence and outputs it to the subtractor 302, and the subtractor 302 outputs the first A PPG signal PPG_C is subtracted from the digitally filtered motion information M_C to obtain the error value Error_C. Since the weight adjustment unit 303 is connected to the first finite impulse response digital filter unit 301 and the subtractor 302, the error value after subtraction between the first PPG signal PPG_C and the digitally filtered motion information M_C can be obtained. Error_C to adjust the weight W of the first FIR digital filtering unit 301 for each digital filtering until the error value is minimized and then output the second PPG signal Error_C, as shown in FIG. 3D . Please also refer to FIG. 7A, which is the conversion of the second PPG signal Error_C in FIG. 3D to the frequency domain waveform. Compared with the same conversion to the frequency domain waveform in FIG. 3A, it can be seen that the interference signal S motion (1) of the motion artifact has been Effectively suppressed, the highest peak in the waveform is the heartbeat signal S heart (1).
上述第一有限脉冲响应数字滤波器301的更新函数为:The updating function of the above-mentioned first finite impulse response digital filter 301 is:
m(k+1)=m(k)+w(k)×M_C(k);其中m(k+1)为下一次撷取的该运动信息经数字滤波后的数值,w(k)为该第一有限脉冲响应数字滤波单元的本次权重值,M_C(k)为目前运动信息。该权重值由以下公式计算:m(k+1)=m(k)+w(k)×M_C(k); where m(k+1) is the digitally filtered value of the motion information captured next time, and w(k) is The current weight value of the first FIR digital filter unit, M_C(k), is the current motion information. This weight value is calculated by the following formula:
w(k+1)=w(k)+2×μ×Error_C(k)×X_C(k);其中w(k+1)为该第一有限脉冲响应数字滤波单元的下一次权重值,μ为误差收敛系数。w(k+1)=w(k)+2×μ×Error_C(k)×X_C(k); wherein w(k+1) is the next weight value of the first finite impulse response digital filter unit, μ is the error convergence coefficient.
再请参阅图4A及图5B所示,本实施例为配合使用三轴动力感应器20a的运动信息撷取单元,该三轴动力感应器20a连接有三个第一有限脉冲响应数字滤波单元301,且该第一可适性数字滤波单元30a的本次误差值为:Please refer to FIG. 4A and FIG. 5B again, the present embodiment is a motion information acquisition unit using a three-axis dynamic sensor 20a, and the three-axis dynamic sensor 20a is connected with three first finite impulse response digital filter units 301, And the current error value of the first adaptive digital filtering unit 30a is:
Error_C(k)=PPG_C(k)-y_0(k)-y_1(k)-y_2(k),而各该第一有限脉冲响应数字滤波单元301的更新函数分别如下:Error_C(k)=PPG_C(k)-y_0(k)-y_1(k)-y_2(k), and the update functions of each of the first finite impulse response digital filter units 301 are as follows:
y_0(k+1)=y_0(k)+w_0(k)×X_C(k);y_0(k+1)=y_0(k)+w_0(k)×X_C(k);
y_1(k+1)=y_1(k)+w_1(k)×Y_C(k);y_1(k+1)=y_1(k)+w_1(k)×Y_C(k);
y_2(k+1)=y_2(k)+w_2(k)×Z_C(k);其中:y_2(k+1)=y_2(k)+w_2(k)×Z_C(k); where:
X_C(k)为目前第一轴运动信息;X_C(k) is the current movement information of the first axis;
y_0(k+1)为下一次该第一轴运动信息经数字滤波后的数值;y_0(k+1) is the value after digital filtering of the first axis motion information next time;
Y_C(k)为目前第二轴运动信息;Y_C(k) is the current movement information of the second axis;
y_1(k+1)为下一次该第二轴运动信息经数字滤波后的数值;y_1(k+1) is the digitally filtered value of the second axis motion information next time;
Z_C(k)为目前第三轴运动信息;Z_C(k) is the current movement information of the third axis;
y_2(k+1)为下一笔该第三轴运动信息经数字滤波后的数值。y_2(k+1) is the digitally filtered value of the next piece of motion information of the third axis.
此外,为使得该第一可适性数字滤波单元30具的误差收敛系数可动态调整,加快误差最小化的速度,如图1所示,可进一步包含有一第一误差收敛系数自动调整单元31,其包含以下动态产生该误差收敛系数的步骤:In addition, in order to enable the error convergence coefficient of the first adaptive digital filter unit 30 to be dynamically adjusted and speed up the speed of error minimization, as shown in FIG. 1 , a first error convergence coefficient automatic adjustment unit 31 may be further included, It includes the following steps to dynamically generate the error convergence coefficient:
(a)由该第一PPG信号获得其向量值;(a) obtaining its vector value from the first PPG signal;
(b)将该第一PPG信号转换为矩阵,并取得该矩阵与转置矩阵的内积,以获得该第一PPG信号的强度;(b) converting the first PPG signal into a matrix, and obtaining the inner product of the matrix and the transposed matrix to obtain the strength of the first PPG signal;
(c)设定一先进先出型缓冲区;(c) setting a first-in-first-out buffer;
(d)自该先进先出型缓冲区中连续获得步骤(b)的乘积;以及(d) continuously obtaining the product of step (b) from the FIFO buffer; and
(e)选择大于0且于缓冲区内的最大值的倒数的数值为该误差收敛系数;其中该最大值为该先进先出型缓冲区内的最大乘积。(e) Selecting the value greater than 0 and the reciprocal of the maximum value in the buffer as the error convergence coefficient; wherein the maximum value is the maximum product in the first-in-first-out buffer.
再如图1所示,上述第一可适性数字滤波单元20所输出一第二PPG信号Error_C虽然可滤除大部份运动伪影的干扰信号,但仍有其它的干扰信号,故可进一步将该第二PPG信号Error_C依序输出至一去极端值单元32、一第二带通滤波器33及一第三正规化器34,由该去极端值单元32设定一临界值,将时域下的第二PPG信号Error_C中能量大于临界值的异常值消除,再由该第二带通滤波器33滤除心跳信号的频率范围以外的噪声(如图3E所示),最后经第三正规化器34将该第二PPG信号的能量予以正规化,使其能量分布于1及-1之间,如图3E所示,以增加系统的工作动态范围,如图3F所示。As shown in FIG. 1 again, although a second PPG signal Error_C output by the above-mentioned first adaptive digital filtering unit 20 can filter out most of the interference signals of motion artifacts, there are still other interference signals, so further The second PPG signal Error_C is sequentially output to an extreme value removal unit 32, a second bandpass filter 33, and a third normalizer 34, and a critical value is set by the extreme value removal unit 32, and the time In the second PPG signal Error_C under the domain, the abnormal value whose energy is greater than the critical value is eliminated, and then the noise (as shown in Figure 3E) outside the frequency range of the heartbeat signal is filtered by the second bandpass filter 33, and finally through the third The normalizer 34 normalizes the energy of the second PPG signal so that its energy is distributed between 1 and −1, as shown in FIG. 3E , so as to increase the working dynamic range of the system, as shown in FIG. 3F .
接着,将正规化后的第二PPG信号PPG_N输出至该第二可适性数字滤波单元40,将其中的周期性信号的信噪比(SNR)提高之后输出一生理特征信号。请配合参阅图1及图6所示,在本实施例中该第二可适性数字滤波单元40包含有一第二有限脉冲响应数字滤波单元401、一减法器402及一权重调整单元403;其中该第二有限脉冲响应数字滤波单元401透过一时间延迟单元404接收正规化后的第二PPG信号PPG_N,以接收到延迟后的第二PPG信号,并将该第二PPG信号予以数字滤波后输出至该减法器402,该减法器402则将本次接收到正规化后的第二PPG信号PPG_N,与在本次之前Δ次的正规化第二PPG信号PPG_N相减,以输出一误差值Error_P。该误差值Error_P输入至该权重调整单元403,以依据该误差值调整该第二有限脉冲响应数字滤波单元401下次数字滤波的权重值W_P;当Error_P趋近于0时,经该第二有限脉冲响应数字滤波单元401数字滤波后输出信号,即为该周期性信号;即如图7B所示,图3G的周期性信号转换为频域波形,波形最高峰仍是心跳信号Sheart(2),再参阅图7C所示,将图7A的图3D转换的频域波形两相比较后发现,图3G周期性信号Sheart(2)于频域波形的峰值较图3D第二PPG信号Error_C于频域波形的峰值高出Sd;是以,正规化后的第二PPG信号Error_C经过第二可适性数字滤波单元40处理后,确实可将周期性信号Sheart(2)的信噪比(SNR)提高。Next, output the normalized second PPG signal PPG_N to the second adaptive digital filter unit 40 , and output a physiological characteristic signal after improving the signal-to-noise ratio (SNR) of the periodic signal therein. Please refer to FIG. 1 and FIG. 6, in this embodiment, the second adaptive digital filtering unit 40 includes a second finite impulse response digital filtering unit 401, a subtractor 402 and a weight adjustment unit 403; The second finite impulse response digital filtering unit 401 receives the normalized second PPG signal PPG_N through a time delay unit 404 to receive the delayed second PPG signal, and digitally filters the second PPG signal output to the subtractor 402, and the subtractor 402 subtracts the normalized second PPG signal PPG_N received this time from the normalized second PPG signal PPG_N Δ times before this time to output an error value Error_P. The error value Error_P is input to the weight adjustment unit 403, so as to adjust the weight value W_P of the next digital filtering of the second finite impulse response digital filter unit 401 according to the error value; The output signal after digital filtering by the impulse response digital filtering unit 401 is the periodic signal; that is, as shown in FIG. 7B, the periodic signal in FIG. 3G is converted into a frequency domain waveform, and the highest peak of the waveform is still the heartbeat signal S heart (2) Referring to Fig. 7C again, after comparing the frequency domain waveforms of Fig. 3D conversion in Fig. 7A, it is found that the peak value of the periodic signal S heart (2) in Fig. 3G in the frequency domain waveform is higher than that in Fig. The peak value of the frequency domain waveform is higher than S d ; therefore, after the normalized second PPG signal Error_C is processed by the second adaptive digital filter unit 40, the signal-to-noise ratio of the periodic signal S heart (2) can indeed be reduced to (SNR) improved.
在本实施例中,上述该第二有限脉冲响应数字滤波单元401输出的该周期性信号为:PPG_P(k+1)=PPG_P(k)+w_P(k)×PPG_N(k-Δ);其中:PPG_P(k+1)为下一次该周期性信号经数字滤波后的数值;PPG_P(k)为目前周期性信号经数字滤波后的数值;w_P(k)为本次该有限脉冲响应数字滤波单元的权重值;PPG_N(k-Δ)为本次之前Δ次的第二PPG信号。而该第二可适性数字滤波单元的权重值为:w_P(k+1)=w_P(k)+2×μ×Error_P(k)×PPG_N(k-Δ);其中w_P(k+1)为下次该有限脉冲响应数字滤波单元的权重值;μ为该误差收敛系数;及Error_P(k)为本次误差值。该第二可适性数字滤波单元的误差值为:Error_P(k)=PPG_N(k)-PPG_P(k);其中PPG_N(k)为目前第二PPG信号。是以,当误差接近为0,则该第二有限脉冲响应数字滤波单元输出的该周期性信号即接近实际的生理特征信号(心跳次数),故可依据周期性信号计算心跳次数。In this embodiment, the periodic signal output by the second finite impulse response digital filter unit 401 is: PPG_P(k+1)=PPG_P(k)+w_P(k)×PPG_N(k-Δ); : PPG_P(k+1) is the value of the periodic signal after digital filtering next time; PPG_P(k) is the value of the current periodic signal after digital filtering; w_P(k) is the digital filtering of the finite impulse response this time The weight value of the unit; PPG_N(k-Δ) is the second PPG signal Δ times before this time. And the weight value of this second adaptive digital filtering unit is: w_P(k+1)=w_P(k)+2×μ×Error_P(k)×PPG_N(k-Δ); wherein w_P(k+1) is the weight value of the finite impulse response digital filter unit next time; μ is the error convergence coefficient; and Error_P(k) is the current error value. The error value of the second adaptive digital filtering unit is: Error_P(k)=PPG_N(k)−PPG_P(k); where PPG_N(k) is the current second PPG signal. Therefore, when the error is close to 0, the periodic signal output by the second finite impulse response digital filter unit is close to the actual physiological characteristic signal (heartbeat frequency), so the heartbeat frequency can be calculated according to the periodic signal.
如图1所示,为使得该第二可适性数字滤波单元30具的误差收敛系数可动态调整,以加快误差最小化的速度,同样可进一步包含有一第二误差收敛系数自动调整单元41,其包含以下步骤:As shown in Figure 1, in order to make the error convergence coefficient of the second adaptive digital filtering unit 30 dynamically adjustable, to speed up the speed of error minimization, it can also further include a second error convergence coefficient automatic adjustment unit 41, It consists of the following steps:
(a)由该第二PPG信号获得其向量值;(a) obtaining its vector value from the second PPG signal;
(b)将该第二PPG信号转换为矩阵,并对该矩阵与一转置矩阵的内积,以获得该第二PPG信号的强度;(b) converting the second PPG signal into a matrix, and inner producting the matrix with a transposed matrix to obtain the strength of the second PPG signal;
(c)设定一先进先出型缓冲区;(c) setting a first-in-first-out buffer;
(d)自该先进先出型缓冲区中连续获得步骤(b)的乘积;(d) continuously obtaining the product of step (b) from the FIFO buffer;
(e)选择大于0且于缓冲区内最大值的倒数数值为该误差收敛系数;其中该最大值为该先进先出型缓冲区内的最大乘积。(e) Selecting the reciprocal value greater than 0 and the maximum value in the buffer as the error convergence coefficient; wherein the maximum value is the maximum product in the first-in-first-out buffer.
欲使用图3G的周期性信号计算心跳次数,则经快速传立叶转换(FFT)后即如图7B所示的频域波形(X轴以指数(index)表示之,Y轴则是经过FFT后的能量频谱),依据频域波形计算每分钟的心跳次数Fheart=index×(fs/n)×60,其中fs为取样频率,n为FFT点数。请配合图8所示,即为所计算出的该心跳频率于时域下波形图,相较图9C受运动伪影干扰的心跳波形图可知,本发明的心跳频率量测即不再受到运动伪影的干扰。If you want to use the periodic signal in Figure 3G to calculate the number of heartbeats, the frequency-domain waveform shown in Figure 7B will be obtained after Fast FFT Transformation (FFT) (the X-axis is represented by an index, and the Y-axis is represented by an FFT After energy spectrum), calculate the number of heartbeats per minute F heart = index × (fs/n) × 60 according to the frequency domain waveform, where fs is the sampling frequency, and n is the number of FFT points. Please refer to FIG. 8, which is the waveform diagram of the calculated heartbeat frequency in the time domain. Compared with the waveform diagram of the heartbeat disturbed by motion artifacts in FIG. 9C, it can be seen that the heartbeat frequency measurement of the present invention is no longer affected by motion. Artifact interference.
由上述说明可知,本发明主要于撷取PPG信号期间,以一运算信息撷取单元撷取运动信息,由于自人体撷取到的PPG信号会包含有运动伪影的干扰信号,故该第一可适性数字滤波单元依据运动信息,将大部份的运动伪影干扰信号自该PPG信号中滤除,再将滤除后的PPG信号输出至一第二可适性数字滤波单元,以撷取该第二PPG信号中的周期性信号并提高周期性信号的信噪比,作为读出精确的心跳信号的依据。因此,本发明的PPG信号处理装置可应用于穿戴式心跳检测装置上,提供正确的心跳信息。It can be seen from the above description that the present invention mainly uses an operation information extraction unit to extract motion information during the acquisition of PPG signals. Since the PPG signals extracted from the human body will contain interference signals with motion artifacts, the first The adaptive digital filtering unit filters most of the motion artifact interference signals from the PPG signal according to the motion information, and then outputs the filtered PPG signal to a second adaptive digital filtering unit to extract Taking the periodic signal in the second PPG signal and improving the signal-to-noise ratio of the periodic signal is used as the basis for reading out the accurate heartbeat signal. Therefore, the PPG signal processing device of the present invention can be applied to a wearable heartbeat detection device to provide correct heartbeat information.
以上所述仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例公开如上,然而并非用以限定本发明,任何所属技术领域的技术人员,在不脱离本发明技术方案的范围内,当可利用上述公开的技术内容作出些许更动或修改为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修改,均仍属于本发明技术方案的范围内。The above descriptions are only preferred embodiments of the present invention, and do not limit the present invention in any form. Although the present invention has been disclosed as above with preferred embodiments, it is not intended to limit the present invention, and any technology in the technical field Personnel, without departing from the scope of the technical solution of the present invention, when the technical content disclosed above can be used to make some changes or modifications to equivalent embodiments with equivalent changes, but any content that does not depart from the technical solution of the present invention, according to the present invention Any simple modifications, equivalent changes and modifications made to the above embodiments by the technical essence still belong to the scope of the technical solution of the present invention.
Claims (21)
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462051957P | 2014-09-18 | 2014-09-18 | |
US62/051,957 | 2014-09-18 | ||
TW104125755A TWI558375B (en) | 2014-09-18 | 2015-08-07 | Ppg signal processing device and method thereof |
TW104125755 | 2015-08-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105433931A true CN105433931A (en) | 2016-03-30 |
Family
ID=55544898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510565554.0A Pending CN105433931A (en) | 2014-09-18 | 2015-09-08 | Processing device and method for describing waveform by light volume change |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105433931A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106691425A (en) * | 2016-12-30 | 2017-05-24 | 北京工业大学 | Wrist heart rate monitoring method of sports bracelet |
CN107223036A (en) * | 2017-04-10 | 2017-09-29 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for adaptive filtering motion interference |
CN108778109A (en) * | 2016-03-01 | 2018-11-09 | 皇家飞利浦有限公司 | Equipment, system and method for the vital sign for determining object |
CN110101372A (en) * | 2019-04-24 | 2019-08-09 | 上海工程技术大学 | A kind of municipal rail train driver physiological status monitoring system |
US10466783B2 (en) | 2018-03-15 | 2019-11-05 | Sanmina Corporation | System and method for motion detection using a PPG sensor |
CN111374657A (en) * | 2018-12-28 | 2020-07-07 | 宏碁股份有限公司 | Method and device for detecting arrhythmia based on photoplethysmography signal |
CN114224309A (en) * | 2017-11-16 | 2022-03-25 | 原相科技股份有限公司 | Physiological detection system and physiological detection device |
CN116269260A (en) * | 2023-03-01 | 2023-06-23 | 亿慧云智能科技(深圳)股份有限公司 | A smart watch heart rate monitoring method and system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1989895A (en) * | 2003-05-09 | 2007-07-04 | 三星电子株式会社 | Ear device for measuring biological signal and its measuring method |
WO2012134395A1 (en) * | 2011-03-25 | 2012-10-04 | Nitto Denko Corporation | A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system |
CN103099615A (en) * | 2013-01-23 | 2013-05-15 | 深圳市理邦精密仪器股份有限公司 | Method and device for eliminating exercise electrocardiosignal interference |
CN103596491A (en) * | 2011-01-21 | 2014-02-19 | 伍斯特理工学院 | Physiological parameter monitoring with a mobile communication device |
US20140058217A1 (en) * | 2012-08-25 | 2014-02-27 | The Board Of Trustees Of The Leland Stanford Junior University | Motion artifact mitigation methods and devices for pulse photoplethysmography |
CN103781414A (en) * | 2011-09-16 | 2014-05-07 | 皇家飞利浦有限公司 | Device and method for estimating the heart rate during motion |
US20140288436A1 (en) * | 2012-06-22 | 2014-09-25 | Fitbit, Inc. | Wearable heart rate monitor |
-
2015
- 2015-09-08 CN CN201510565554.0A patent/CN105433931A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1989895A (en) * | 2003-05-09 | 2007-07-04 | 三星电子株式会社 | Ear device for measuring biological signal and its measuring method |
CN103596491A (en) * | 2011-01-21 | 2014-02-19 | 伍斯特理工学院 | Physiological parameter monitoring with a mobile communication device |
WO2012134395A1 (en) * | 2011-03-25 | 2012-10-04 | Nitto Denko Corporation | A method of measuring an artefact removed photoplethysmographic (ppg) signal and a measurement system |
CN103781414A (en) * | 2011-09-16 | 2014-05-07 | 皇家飞利浦有限公司 | Device and method for estimating the heart rate during motion |
US20140288436A1 (en) * | 2012-06-22 | 2014-09-25 | Fitbit, Inc. | Wearable heart rate monitor |
US20140058217A1 (en) * | 2012-08-25 | 2014-02-27 | The Board Of Trustees Of The Leland Stanford Junior University | Motion artifact mitigation methods and devices for pulse photoplethysmography |
CN103099615A (en) * | 2013-01-23 | 2013-05-15 | 深圳市理邦精密仪器股份有限公司 | Method and device for eliminating exercise electrocardiosignal interference |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108778109A (en) * | 2016-03-01 | 2018-11-09 | 皇家飞利浦有限公司 | Equipment, system and method for the vital sign for determining object |
CN106691425A (en) * | 2016-12-30 | 2017-05-24 | 北京工业大学 | Wrist heart rate monitoring method of sports bracelet |
CN106691425B (en) * | 2016-12-30 | 2019-06-21 | 北京工业大学 | A wrist heart rate monitoring method for a sports bracelet |
CN107223036A (en) * | 2017-04-10 | 2017-09-29 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for adaptive filtering motion interference |
WO2018187895A1 (en) * | 2017-04-10 | 2018-10-18 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for self-adaptive filtering of motion interference |
CN107223036B (en) * | 2017-04-10 | 2019-03-29 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for adaptively filtering out motion interference |
CN114224309B (en) * | 2017-11-16 | 2023-11-17 | 原相科技股份有限公司 | Physiological detection system and physiological detection device |
CN114224309A (en) * | 2017-11-16 | 2022-03-25 | 原相科技股份有限公司 | Physiological detection system and physiological detection device |
WO2019177700A3 (en) * | 2018-03-15 | 2020-02-20 | Sanmina Corporation | System method for motion detection using a ppg sensor |
US10466783B2 (en) | 2018-03-15 | 2019-11-05 | Sanmina Corporation | System and method for motion detection using a PPG sensor |
US11675434B2 (en) | 2018-03-15 | 2023-06-13 | Trilinear Bioventures, Llc | System and method for motion detection using a PPG sensor |
CN111374657A (en) * | 2018-12-28 | 2020-07-07 | 宏碁股份有限公司 | Method and device for detecting arrhythmia based on photoplethysmography signal |
CN111374657B (en) * | 2018-12-28 | 2023-08-08 | 宏碁股份有限公司 | Method and device for detecting arrhythmia based on photoplethysmograph signals |
CN110101372A (en) * | 2019-04-24 | 2019-08-09 | 上海工程技术大学 | A kind of municipal rail train driver physiological status monitoring system |
CN116269260A (en) * | 2023-03-01 | 2023-06-23 | 亿慧云智能科技(深圳)股份有限公司 | A smart watch heart rate monitoring method and system |
CN116269260B (en) * | 2023-03-01 | 2023-07-25 | 亿慧云智能科技(深圳)股份有限公司 | Smart watch heart rate monitoring method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105433931A (en) | Processing device and method for describing waveform by light volume change | |
Han et al. | Electrocardiogram signal denoising based on a new improved wavelet thresholding | |
EP2451359B1 (en) | Noise reduction of breathing signals | |
CN113349752A (en) | Wearable device real-time heart rate monitoring method based on sensing fusion | |
WO2019036908A1 (en) | Heart rate measuring method and apparatus, and electronic terminal | |
JP7177443B2 (en) | Biological vibration signal detector | |
CN107949321A (en) | Temporal Interference Removal and Improved Heart Rate Measurement Tracking Mechanism | |
CN102551726A (en) | Respiratory signal processing apparatus, respiratory signal processing method, and program | |
CN105286845A (en) | Movement noise elimination method suitable for wearable heart rate measurement device | |
CN103099615A (en) | Method and device for eliminating exercise electrocardiosignal interference | |
EP2501277A1 (en) | Apparatus and methods for monitoring heart rate and respiration | |
CN105816165B (en) | A real-time dynamic heart rate monitoring device and monitoring method | |
CN104244127A (en) | Heart rate detection method applied to ear phone and ear phone capable of detecting heart rate | |
CN112587133A (en) | Method for measuring blood oxygen saturation | |
CN104244125A (en) | Heart rate detection method applied to ear phone and ear phone capable of detecting heart rate | |
CN108567419B (en) | Signal processing device and method based on respiratory heartbeat characteristic wave | |
CN110292372B (en) | Detection device | |
TWI558375B (en) | Ppg signal processing device and method thereof | |
Bai et al. | Using an adaptive filter to remove ECG motion artifact interference | |
CN108294736A (en) | Continuous BP measurement system and measurement method | |
CN106197480B (en) | A kind of processing system of Low SNR signal | |
JP2013106837A (en) | Heart rate detection method, heart rate detector, and mental stress measuring apparatus | |
JP6464004B2 (en) | Measuring apparatus, measuring method and program | |
JP6388448B2 (en) | Biological information signal source identification apparatus and biological information signal source identification method | |
Altay et al. | Method for improvement of biosignal measurement accuracy in the IoB system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160330 |
|
WD01 | Invention patent application deemed withdrawn after publication |