The invention relates to a device and to a method for detecting and reporting of a stress condition of a person.
Using the pulse rate or the heart rate variability for the detection of stress conditions has been known for a long time. Thereby the interval between two heart beats is defined—in the sense of the present invention—as the time between the onsets of two contractions of the cardiac chambers. This onset of the chamber contraction shows up in the electrocardiogram (ECG) as the R wave. The distance between two R waves is usually denoted RR interval. After having averaged over a defined number of RR intervals, the heart rate can be determined by calculation. The individual values of the RR intervals vary around the mean value thus obtained. Thereby the variations can change from beat to beat. The variation is usually termed as heart rate variability (HRV). In principle, the heart rate can also be determined by a pressure measurement carried out on an artery.
Physiologically, the heart rate variability (HRV) depends on the ability of the human organism to adapt the rate of the cardiac rhythm. Variations of the heart rate, i.e. variations of the temporal interval between two heart beats, can occur in a resting state, in which case they are mostly spontaneous, but also in the course of specific variations of the surrounding conditions, e.g. under stress. A healthy organism continuously adapts the heart beat rate to the current conditions via physiological regulation pathways of the vegetative nervous system. Therefore, physical or psychological stress usually results in an increase of the heart frequency which ordinarily decreases again upon relief and relaxation. Thereby, a good adaptability to stress results in a higher variability of the heart rate. Under chronic stress burden, the adaptability is reduced. In this respect, it is known that the heart rate variability taken by itself already provides a certain—albeit still very unreliable—indicator for the current stress burden and the ability of a person to cope with stress.
Several methods for determining the stress condition of a person have been proposed in the prior art, including the proposal to use further measurement parameters in addition to the pulse rate. Thus, DE 103 19 361 A1 proposes to use the pulse wave latency in addition to the heart rate variability.
Regarding the analysis of the heart rate variability, reference is made to DE 100 06 154 A1, DE 10 2006 039 957 A1, and also to DE 10 2008 030 956 A1 and EP 1 156 851 B1, wherein the person skilled in the art can find various determination methods.
It is an object of the invention to provide a device and a corresponding method for detecting and reporting of a stress condition of a person having a higher reliability as compared to the prior art.
The object of the invention is achieved according to a first aspect by a method for detecting and reporting of a stress condition of a person, the method comprising the following steps:
Z=A*P+B*HRV
P
Hist
=ΣE(t)(P−G*F1(P0)).
On the one hand, the features of the invention imply that the values of the current pulse frequency and of the current heart rate variability, each one weighted with respective weighting factors, are added together. This reflects the findings already known from prior art that at a very high pulse frequency the heart rate variability temporarily decreases sharply and thus has a correspondingly reduced significance. On the other hand, in the resting state the pulse frequency has only small significance regarding the stress tolerance, so that then the heart rate variability becomes more important. By providing the above mentioned status function with at least one, preferably additive correcting value that includes the history of the person at least within the past 2 hours, preferably within an interval of about 5 hours to maximally 72 hours, the status value of the stress condition becomes—according to the present invention—considerably more informative.
Thereby, the history is taken into account by means of a sum PHist of the measured pulse frequencies P minus a function F1 at least of the resting heart rate P0, the sum being weighted with a weighting factor E according to
P
Hist
=E*Σ(P−G*F1(P0))
following the method of a moving window (“moving average”). Thereby, E can be a constant, but also a quantity that diminishes linearly with time and which has—assessed from the current point of time—the full value at the beginning of the summation and a negligible value at the end of time. Alternatively, the history can also be implemented by a filter, preferably by a digital lowpass filter that includes the entire history with regard to stress and recovery, which predominantly is a recent history.
In this context, it is reasonable if the function F1 mentioned above still depends on the age of the person and/or the maximal pulse of the person, preferably as determined by the Conconi test.
According to a second aspect a lowpass filter for the pulse rate history is used instead of the sum function of the pulse history.
According to a third aspect as history a sum of a function F2 of the ratios LFtot/HFtot, which sum is weighted with a weighing factor H, according to
HRVHist=H*Σ(F2(LFtot/HFtot))
is taken into account instead of the sum function of the pulse rate history, again according to the method of a moving window (“moving average”).
Thereby, H can again be a constant, but also a quantity that diminishes linearly with time and which has—assessed from the current point of time—the full value at the beginning of the summation and a negligible value at the end of time. Again, the history can alternatively also be implemented by a filter, preferably by a digital lowpass filter that includes the entire history with regard to stress and recovery, but predominantly a recent history.
It is particularly simple and advantageous if the function F2(LFtot/HFtot) is a norm function having the values 1 at a current HRV smaller than a first threshold value of a predetermined standardized HRVnorm, 0 at a current HRV larger than the first threshold value of the predetermined standardized HRVnorm but smaller than a second threshold value of the predetermined standardized HRVnorm, and −1 at a current HRV larger than the second threshold value of the predetermined standardized HRVnorm.
According to a further aspect of the present invention the object is achieved by means of a device for detecting and reporting of a stress condition of a person, the device comprising:
Z=A*P+B*HRV
P
Hist
=ΣE(t)(P−G*F1(P0))
According to a further aspect of the present invention the object is achieved by means of a device for detecting and reporting of a stress condition of a person, the device comprising:
Z=A*P+B*HRV
HRVHist=ΣH(t)F2(LFtot/HFtot)
The features of the invention according to the two last-mentioned aspects provide a particularly informative device that is relatively simple and reliable as compared to corresponding devices according to the prior art, in particular because it does not determine any superfluous measurement parameters. The device advantageously comprises a display device that is configured in such manner that it can display, preferably graphically, at least the pulse rate, the HRV value during a predetermined or preselected time and also the status function.
According to an aspect of the present invention it is advantageous if the device is configured in such manner that the history includes a sum PHist of the measured pulse frequencies P minus a function F1 at least of the resting heart rate P0, the sum being weighted with a weighting factor E, according to
P
Hist
=E*Σ(P−G*F1(P0))
Thereby, E can be a constant, but also a quantity that diminishes linearly with time and which has—assessed from the current point of time—the full value at the beginning of the summation and a negligible value at the end of time.
Also in this case it is appropriate if the above mentioned function F1 further depends on the age of the person and/or the maximal pulse of the person, preferably as determined by the Conconi test.
For the function F in the device the following functions are proposed alternatively:
F=(c1−(Age−20)·c2)·P0
F=P
0+(P−P0)·c3
F=Pmean
It is particularly advantageous if the device is configured in such manner that the history also takes into account a sum of a function F2 of the ratios LFtot/HFtot, which sum is weighted with a weighting factor H, according to
HRVHist=H*Σ(F2(LFtot/HFtot))
Thereby, H can be a constant, but also a quantity that diminishes linearly with time and which has—assessed from the current point of time—the full value at the beginning of the summation and a negligible value at the end of time.
It is particularly simple and advantageous if the function F2(LFtot/HFtot) is a norm function having the values 1 at a current HRV smaller than a first threshold value of a predetermined standardized HRVnorm, 0 at a current HRV larger than the first threshold value of the predetermined standardized HRVnorm, but smaller than a second threshold value of the predetermined standardized HRVnorm, and −1 at a current HRV larger than the second threshold value of the predetermined standardized HRVnorm. Again, the history can also be implemented alternatively as a filter, preferably a digital lowpass filter that takes into account the entire history regarding stress and recovery, which predominantly is a recent history.
The elements mentioned above and claimed, which shall be used according to the invention, as well as those described in the following exemplary embodiments, are not subject to any particular limitations in terms of their size, shape, use of material and technical design, so that the selection criteria known in the respective application field can be adopted without any restrictions. In particular, the method of the present invention is not intended to assess the health status or the pathological status of the person.
The device according to the invention comprises, according to a preferred exemplary embodiment of the invention, a measuring device for detecting the pulse rate and the values that are necessary for calculating heart rate variability. In the present case this is a pulse measuring sensor, but alternatively it can also be an electrical sensor for measuring electric cardiographic measurement values, as well as a display device. Moreover, the device comprises an interface for the input of person-related parameters, which are particularly needed for the detection of the history to be used according to the invention. A key component of the device is a computing device that controls the necessary acquisition of the measurement data, processes the measurement data in their necessary digital form, carries out the data processing and controls the display.
The proposed device is generally used according to the following scheme:
Firstly, a parametrization is performed, which obviously is required only once for every person. The parameterization comprises data on age, sex, and if necessary any other correcting factors, if such shall be used for fine tuning.
Moreover, a step shall be performed that will henceforth be called calibration and that has to be performed periodically, e.g. annually on the occasion of an aptitude test. In this process are determined e.g. the resting heart rate p0—e.g. in a relaxed state in the morning before breakfast or stress-relieved after 5′ lying down—as well as—if this option is activated—the maximal heart rate pmax—e.g. by means of the Conconi test or of a similar stress test. Moreover, a normalized value of the heart rate variability HRVnorm is determined as base value in a resting state. If moreover e.g. the anaerobic threshold shall also be included in the function, this can also be determined and entered.
It is intended that the device shall be worn in all actions. This results in an operation with a turning on about at least 2 min before the action, although for a complete detection of the history it is obviously better to do the activation several hours before the action. Turning off should also not be done immediately after the effort—if possible—so that relaxation data can be detected (relaxation curve).
The analysis is carried out by the computing device. In the specified exemplary embodiment the display comprises a display of the continuously collected data that depends on the use, e.g. with a temporal display of 10 values per second for the value HRV, but also an averaged data display as well as an analysis of the value HRV by means of a fast Fourier transform FFT with a window of about 2 minutes.
The calculation of a status value (stress value) in the present exemplary embodiment is a function of the pulse rate P and the heart rate variability HRV, and also of the history of the two values (phist, HRVhist). Moreover, in the present exemplary embodiment it is contemplated to carry out adaptations to the data of the person, which adaptations are made on the occasion of cyclic calibrations. The history serves the purpose of taking into account previous actions and should include whether a person was already previously exposed to stress (reduction of the performance capability and of the duration of the action, respectively). The status value being used is calculated according to
Stress=A·P+B·Phist+C·HRV+D·HRVhist
p: current pulse
HRV: current value of the heart rate variability
Phist: pulse history
HRVhist: HRV history
Thereby,
is used with
P0: resting heart rate, determined e.g. after getting up, before breakfast
Pmax: maximal pulse, determined e.g. by the Conconi test
As a simple function F, the device allows selection between different functions, namely
F=(c1−(Alter−20)·c2)·P0
F=P
0+(Pmax−P0)·c3
F=Pmean
However, other functions can be selected which include further history- and person-related values, e.g. sex, body size values such as body mass index, etc.
The interval between two heart beats is usually defined as the time between the onsets of two contractions of the cardiac chambers (R wave). Therefore, the distance between two R waves is denoted as RR interval. The real RR intervals vary around the mean heart rate, wherein these variations can even change from beat to beat; this is termed as variability of the cardiac frequency, or heart rate variability (HRV).
There are various characteristic values that are useful for the analysis in the stress sensor, but in the present exemplary embodiment the quotient of the integrals of different frequency ranges obtained from a Fourier transform is used
Thereby the value HF is mostly attributed to a parasympathetic influence whereas the value LF is attributed to a sympathetic influence.
Alternatively, also a statistical time analysis, namely the variation between the RR intervals (e.g. RMSSD or pNN50) can be used. The device is configured in such manner that the above option can be chosen.
In both cases lower values are worse than higher values because larger variations are indicative of a healthy heart and indicate a normal interplay between the sympathetic and parasympathetic nervous system.
According to the present exemplary embodiment the Fourier transform is carried out in the device by means of the RR data over an interval of about 2 minutes. However, this can be adapted depending on the particular application with the object to obtain a not too large latency but nonetheless sufficient datapoints for a significant value. The calculation of the HRV is carried out according to known algorithms for a fast Fourier transform.
The calculation of the heart rate variability is carried out in the present exemplary embodiment as follows: The norm value for the heart rate variability HRVnorm for the person is determined on the occasion of the calibration. In general the value lies between 1.5 and 2.0.
Example for Function F
1: HRV<d1·HRVnorm
F=0: d1·HRVnorm<HRV<d2·HRVnorm
−1: HRV>d2·HRVnorm
Examples are shown hereinbelow for subjects