US 11,793,449 B2
Heart rate variability analysis method, device and use thereof
Luming Li, Beijing (CN); Hongyun Liu, Beijing (CN); and Zhao Yang, Beijing (CN)
Assigned to TSINGHUA UNIVERSITY, Beijing (CN); and BEIJING PINS MEDICAL CO., LTD, Beijing (CN)
Filed by TSINGHUA UNIVERSITY, Beijing (CN); and BEIJING PINS MEDICAL CO., LTD, Beijing (CN)
Filed on Apr. 1, 2021, as Appl. No. 17/301,390.
Application 17/301,390 is a continuation of application No. 16/086,399, granted, now 10,973,455, previously published as PCT/CN2017/095953, filed on Aug. 4, 2017.
Claims priority of application No. 201611094149.6 (CN), filed on Dec. 2, 2016.
Prior Publication US 2021/0219906 A1, Jul. 22, 2021
This patent is subject to a terminal disclaimer.
Int. Cl. A61B 5/024 (2006.01); A61B 5/00 (2006.01); A61B 5/0245 (2006.01); A61B 5/316 (2021.01); A61B 5/339 (2021.01); A61B 5/341 (2021.01); A61B 5/349 (2021.01); A61B 5/352 (2021.01)
CPC A61B 5/4094 (2013.01) [A61B 5/02405 (2013.01); A61B 5/4035 (2013.01); A61B 5/4848 (2013.01); A61B 5/0245 (2013.01); A61B 5/316 (2021.01); A61B 5/339 (2021.01); A61B 5/341 (2021.01); A61B 5/349 (2021.01); A61B 5/352 (2021.01)] 6 Claims
OG exemplary drawing
 
1. An apparatus for analyzing heart rate variability, comprising:
an acquiring module configured to acquire electrocardiography data externally;
a processing module configured to perform a digitization process and a denoising process on the electrocardiography data;
a sequence forming module configured to form a sinus NN interval sequence with the processed electrocardiography data;
a selecting module configured to select sinus NN interval data of a subject during a long term of more than 20 hours and during short terms of 2 to 10 minutes in an awake state and in a sleep state; and
a calculating module to calculate at least one of a time domain index, a frequency domain index and a nonlinear index for the sinus NN interval data of the subject during the long term of more than 20 hours and during the short terms of 2 to 10 minutes in the awake state and in the sleep state, the calculation result are used to analyze the efficacy of vagus nerve stimulation (VNS) treatment for patients with medically intractable epilepsy, wherein the calculation formulas for time domain indexes are expressed as follows, and one or more of the following indexes are selected to be calculated during the analysis;
 
Name of index Definition of index Unit
 
mean value (MEAN)

OG Complex Work Unit Math
 ms
 
overall standard deviation (SDNN)

OG Complex Work Unit Math
 ms
 
root-mean-square of difference of adjacent interval (RMSSD)

OG Complex Work Unit Math
 ms
 
percentage of adjacent pNN50 = NN50/TotalNN × 100% %
intervals with difference
 
 
greater than 50 ms
 
 
(pNN50)
 
the calculation formulas for frequency domain indexes are expressed as follows, and one or more of the following indexes are selected to be calculated during the analysis;
 
Name of index Definition of index Unit Frequency range
 
TP Total signal power, sum of ms2 ≤0.4 Hz
 
 VLF, LF and HF
 
 
VLF Power in very low ms2 0.0033 to 0.04 Hz
 
 frequency band
 
 
LF Power in low frequency ms2 0.04 to 0.15 Hz
 
 band
 
 
HF Power in high frequency ms2 0.15 to 0.4 Hz
 
 band
 
 
LF/HF Ratio of power in low
 
 
 
 frequency band to power
 
 
 
 in high frequency band
 
the nonlinear index comprises SD1 and SD2 indexes in a Poincare Plot, and one or more of the indexes are selected to be calculated during the analysis;
wherein the threshold of RMSSD is 37 ms, the threshold of pNN50 is 27%, the threshold of SD1 is 35 ms, and the threshold of HF is 798 ms.