The invention relates to methods, systems and computer programs to measure the phase of a pseudo-periodic physiological signal of a user.
More precisely, the invention relates to methods to measure the phase of a pseudo-periodic physiological signal of a user. A physiological signal of a user can be measured in different ways. One typical way is to measure a physiological signal of a user using a suitable sensor, as part of a medical examination procedure. This is generally performed using specific material by professional personal (such as nurses or doctors) which are experienced with this material.
There is a recent trend for measuring physiological signals of someone by untrained personal, very often by the person itself. This would enable to obtain signal measurements over a longer period of time than what is typically done during a medical examination procedure. However, the measurement conditions are adverse because the measurement will be set by the untrained user and may be subjected to aggressive every-day life conditions.
Such physiological signals are typically pseudo-periodic. This reflects the cycles of human living. By “pseudo-periodic”, it is meant that the signal is not exactly periodic according to the strict mathematical definition, but roughly follows a shape similar to that of a trigonometric function. Taking as an example the breathing cycle, it comprises an inspiration phase and an expiration phase which together toughly form a period of a trigonometric signal. For one cycle, the period is the time duration between the start of two successive inspirations. Due to biological processes, the period is not exactly constant over time, but varies with time. The phase is defined as the instantaneous position within a given cycle. Another typical example of a pseudo-periodic physiological signal which can be used within the frame of the invention is an electro-encephalogram system.
It is possible, based on a previously measured signal, to determine the exact period of a full oscillation, and, a posteriori, the exact phase within the period at any given time within the period. However, this requires, at said given time, to know the shape of the signal after this given time. In other words, it is not possible to determine the instantaneous phase of the signal.
One strives to determine the instantaneous phase of a pseudo-periodic biolological signal in real time. This would be useful to successfully drive a process based on the instantaneous phase of the biological signal.
According to a first object, the invention relates to a method to measure in real-time the phase of a pseudo-periodic physiological signal of a user, wherein :
Thanks to these features, the instantaneous phase of the signal can be measured in real time.
According to some embodiments, one may use one or more of the following features:
the method further comprises:
the method further comprises measuring said physiological signal with a sensor;
the method further comprises emitting an acoustic stimulation signal based on the phase determined by the phase determination module.
According to another aspect, the invention relates to a computer program comprising instructions adapted to perform the steps of the above methods wherein the computer program is run on a processor.
According to another aspect, the invention relates to a system to measure in real-time the phase of a pseudo-periodic physiological signal of a user, comprising, a physiological signal of the user being provided, which was measured continuously during a time interval ]ta; tb[ extending between times ta and tb, and of duration Dt, where the physiological signal does not extend after tb:
According to some embodiments, the system may comprise one or more features adapted to perform the above methods.
The list of drawings hereby follows:
On the figures, identical or similar elements are designated using the same reference sign.
Thereafter, one or more embodiments of the invention will be described.
Referring firstly to
The device 1 is adapted to be worn by a person P, in particular during the person's sleep period.
The device is adapted in particular to be worn on the head of the person P.
To this end, the device 1 comprises a supporting member 2. The supporting member 2 is adapted to surround the head of the person P at least partially so as to be held thereon. In one embodiment of the invention illustrated in
In the embodiment illustrated in
For example, a first arm 2a surrounds a back of the head, and a second arm 2b surrounds the top of the head. The first and second arms 2a, 2b are respectively connected at their respective ends at a left lateral arm connection point 2e and a right lateral arm connection point 2f, respectively located near the left and right temples of the person P. Finally, the third and fourth arms 2c, 2d respectively extend from the left lateral 2e and right lateral 2f arm connection points, towards the front of the person P.
The device 1 further comprises a plurality of electrodes 3, at least one acoustic transducer 4, and embedded conditioning and control electronics 5.
The electrodes 3, acoustic transducer 4, and electronics 5 are operatively connected to each other. Thus, the embedded conditioning and control electronics 5 are particularly suitable for controlling and for receiving information from the plurality of electrodes 3, and are also able to command and control the emission of an acoustic signal A by the acoustic transducer 4.
To this end, the electrodes 3, the acoustic transducer 4, and the electronics 5 are mounted on the supporting member 2. In this manner the electrodes 3, the acoustic transducer 4, and the electronics 5 are close to each other so that communication between these members 3, 4, 5 is particularly fast and high speed. In the example of
This allows implementing an operation of stimulating the brain waves of a person P in soft real-time.
Thus, in particular, the electronics 5 are capable, in soft real-time, of receiving a measurement signal S from the plurality of electrodes 3 and controlling the emission by the acoustic transducer of an acoustic signal A synchronized with a predefined temporal pattern T of a brain wave of the person P.
“Synchronized with a predefined temporal pattern of a brain wave” is understood to mean that the acoustic signal emitted by the device is temporally synchronized with a brain wave of the person. More precisely, it means that the acoustic signal emitted by the device is temporally synchronized with an instantaneous phase of a brain wave of the person as detailed below.
“Soft real-time” is understood to mean an implementation of the stimulation operation such that the time constraints on this operation, in particular the duration of the operation or the frequency at which it is repeated, are satisfied on the average over a predefined total implementation duration, for example a few hours. It is understood that the implementation of said operation may at certain times exceed said time constraints as long as the average operation of the device 1 and the average implementation of the method satisfies these constraints over the predefined total implementation duration. Time limits may be predefined, beyond which the implementation of the stimulation operation is to be stopped or paused.
To enable such an implementation in soft real-time, a maximum distance between the electrodes 3, the acoustic transducer 4, and the electronics 5 may be less than approximately one meter and preferably less than a few decimeters, enabling them to be connected through a wire embedded in the wearable. In this manner, sufficiently rapid communication between the elements of the device 1 can be guaranteed.
The electrodes 3, the acoustic transducer 4, and the electronics 5 may for example be housed in cavities of the supporting member 2, snapped onto the supporting member 2, or attached to the supporting member 2 for example by gluing, screwing, or other suitable means of attachment. In one embodiment of the invention, the electrodes 3, the acoustic transducer 4, and the electronics 5 may be detachably mounted on the supporting member 2.
Referring now to
The acoustic transducer 4 is adapted to emit an acoustic signal A stimulating at least one inner ear of the person P.
In a first embodiment illustrated in particular in
This osteophonic device 4 may for example be adapted for placement near the ear, for example above it as illustrated in
In a second embodiment, the acoustic transducer 4 is a speaker stimulating the inner ear of the person P via an ear canal leading to said inner ear.
This speaker may be placed outside the ear of the person P or in the ear canal.
The acoustic signal A is a modulated signal that at least partially lies within a frequency range audible to a person P, for example the range of 20 Hz to 30 kHz.
The electrodes 3 are adapted to be in contact with the person P, and in particular with the skin of the person P, in order to capture at least one measurement signal S representative of a physiologic electrical signal E of the person P.
The physiological electrical signal E may in particular be an electroencephalogram (EEG), electro-myogram (EMG), electrooculogram (EOG), photoplethysmogram, pulse-oxygram and accelerometer or any other biosignal measurable in a person P, in particular those with a waveform close to a sinus function.
In particular, the physiological electrical signal E advantageously is an electroencephalogram (EEG) of the person P.
To this end, in one embodiment of the invention, the device 1 comprises at least two electrodes 3 of which at least one is a reference electrode 3a and at least one is an EEG measurement electrode 3b.
The device 1 may further comprise a ground electrode 3c.
In one particular embodiment, the device 1 comprises at least three EEG measurement electrodes 3, so as to capture physiological electrical signals E comprising at least three electroencephalogram measurement channels.
The EEG measurement electrodes 3 are for example arranged on the surface of the scalp of the person P.
In other embodiments, the device 1 may further comprise an EMG measurement electrode, and possibly an EOG measurement electrode.
The measurement electrodes 3 may be reusable electrodes or disposable electrodes. Advantageously, the measurement electrodes 3 are reusable electrodes in order to simplify the everyday use of the device.
The measurement electrodes 3 may be dry electrodes or electrodes coated with contact gel. The electrodes 3 may also be textile or silicone electrodes.
In one embodiment of the invention, the measurement electrodes 3 are active electrodes adapted to preprocess the measurement signal S, for example to perform at least one of the following preprocessing operations:
Such preprocessing of the measurement signal S may for example be implemented by an analog module of the measurement electrode 3 or by an analog module located near the measurement electrode 3.
The embedded conditioning and control electronics 5 receive the measurement signals S from the electrodes 3, possibly preprocessed as detailed above.
Alternatively, one may use other kinds of sensors to measure physiological signals of the user. These may include one or more of a pulse-oxymeter and/or inertial sensors. Physiological signals may thus include a movement signal, such as a respiratory movement signal, obtained from an accelerometer, or a signal representative of respiration, such as obtained from a pulse-oxymeter and/or a signal representative of the cardiac rythm. Other considered signals may include body temperature, body sound and/or body vibrations signals.
If the measurement signals S received by the electronics 5 are not preprocessed, the electronics 5 may apply one and/or more preprocessing operations as detailed above.
The embedded conditioning electronics 5 include one or more microchips, for example at least one microprocessor.
As detailed above, the embedded conditioning and control electronics 5 are adapted to implement an operation of stimulating brain waves of the person P, an operation which will later be described in more detail.
Said means of the embedded conditioning and control electronics 5 are for example microchips, microprocessors, and/or electronic memories, where appropriate mounted and interconnected on flexible or rigid printed circuit boards and operatively connected to the electrodes 3 and to the transducer 4 via wired connections 10.
The device 1 may further comprise a memory 6 as illustrated in
The memory 6 is operatively connected to the electronics 5. The memory 6 may be controlled by the embedded conditioning and control electronics 5 so as to store the measurement signals S.
In one advantageous embodiment of the invention, the memory 6 is capable of storing measurement signals S for a duration of several hours, for example at least eight hours so as to cover an average sleep period of a person P.
The device 1 may further comprise a communication module 7 for communicating with an external server 100. The communication module 7 may be mounted on the supporting member 1 as described above for the electrodes 3, the acoustic transducer 4, and electronics 5. The communication module 7 may be controlled by the embedded conditioning and control electronics 5.
The electronics 5 may in particular be adapted to control the communication module 7 to transfer the measurement signals S stored in memory 6 to an external server 100. The transfer operation may be implemented after a sleep period of the person P.
The communication module 7 may advantageously be a wireless communication module, for example a module implementing a protocol such as Bluetooth or Wi-Fi.
In this manner, when the P person is in a sleep period, he or she is not disturbed by cables, in particular if it is necessary to transmit data during the sleep period.
The device 1 may also comprise a battery 8. The battery 8 may be mounted on the supporting member 1 as described above for the electrodes 3, the acoustic transducer 4, and the electronics 5. The battery 8 may be capable of supplying power to the plurality of electrodes 3, the acoustic transducer 4, and the electronics 5, and where appropriate the memory 6 and the communication module 7. The battery 8 is preferably adapted to supply power for several hours without recharging, more preferably for at least eight hours so as to cover an average sleep period of a person P.
The device 1 can thus operate autonomously during a sleep period of the person P. In this manner in particular, the device 1 is self-contained and adapted to implement one or more operations of stimulating slow brain waves without communicating with an external server 100, in particular without communicating with an external server 100 for several minutes, more preferably several hours, more preferably at least eight hours. This reduces the exposure of the person P to electromagnetic radiation. In particular, the device 1 may also be used to assist the person with falling asleep.
“Self-contained” is thus understood to mean that the device can operate for an extended period of several minutes, preferably several hours, in particular at least eight hours, without needing to be recharged with electrical energy, communicate with external elements such as an external server, or be structurally connected to an external device such as a securing member such as an arm or a bracket.
In this manner the device is suitable for use in the everyday life of a person P without imposing particular constraints.
Furthermore, the supporting member 2 advantageously comprises a device 9 for adjusting to the diameter of the head of the person P. This allows adjusting device 1 to the person P and therefore enables particularly good contact between the electrodes 3 and the skin of the person P.
The adjustment device 9 allows changing a dimension of the supporting member 2 according to a diameter of the head of a person P, to allow fine-tuned adjustment to said diameter.
In one embodiment illustrated in particular in
In a variant of this embodiment, the adjustment device 9 may also include a lock adapted to prevent or allow a relative movement of said two parts 9a, 9b. The lock may be an integral part of one of parts 9a, 9b or may be an element independent of the two parts 9a, 9b.
In another embodiment of the invention, the adjustment device 9 is a soft and flexible portion of the supporting member 2. This portion may be a portion of fabric or elastomer, for example of stretch fabric.
The physiological signal is a pseudo-periodic signal. For the example of slow brain waves, as shown on
The electronics 5 comprise a signal characterization module 11. The signal characterization module 11 may comprise a pulsation determination module 12 and a phase determination module 13.
The pulsation determination module 12 determines the pulsation wps of the physiological signal during time interval ]ta ; tb[. This may for example be based on a trigonometric base function having a pulsation wbf. According to such an example, the pulsation wps of the physiological signal can thus be estimated as equal to wbf.
The phase determination module 13 determines the phase fps of the physiological signal at time tb. This may for example be based on a trigonometric base function having a phase fbf at time t1. According to such an example, the phase fps of the physiological signal at tb can thus be estimated as equal to fbf.
According to some embodiments, the pulsation determination module 12 and the phase determination module 13 can be inter-mixed, so that the pulsation wps of the physiological signal during time interval ]t0; t1[ and the phase fps of the physiological signal at time tb can be determined together.
According to a first step, one provides with signal data y=[y1, . . . , yn] on time window frame t=[t1, . . . , tn] (t1=ta; tn=tb). This signal can be modelised using a trigonometric function, such as a cosinus function:
One defines the value y′ at point k of the base function as follows: y′[k]=Y0 cos(w*t[k]+phi[k]), where w is the pulsation of the base function, phi[k] is the phase of the base function at point k, and Y0 its amplitude.
Alternately, y′[k] can be written as y′[n]=A cos(w*t[k])+B sin(w*t[k]), where A and B are unknown parameters.
One defines a plurality of candidates wj for the pulsation of the signal. In particular, the candidates are selected within a pre-determined pulsation range. The pre-determined pulsation range can be pre-determined in the system, and will depend of the nature of the physiological signal. For example, the pre-determined pulsation range can be a maximal pulsation range designed to cover every likely possibilities of pulsation for the physiological signal under study. For example, the maximal pulsation range will be between 50% and 150% of the average pulsation for this kind of physiological signal.
If the pulsation determination module has already been used recently (in particular since the present process is supposed to be used in real-time at a high frequency), the pre-determined pulsation range can be determined based on the estimated pulsation at the previous use. This is likely to enable a much narrower pre-determined pulsation range. For example, the pre-determined pulsation range could be between 90% and 110% of the estimated pulsation at the previous use of the pulsation determination module. This is for example the case when the phase of the input signal is needed at any time. This is for example the case when a stimulation is performed based on the instantaneous phase of the input signal, such as, for example, when the stimulation is performed based on the respiratory rhythm of the user.
The number of candidates wj is for example 10 to 50. These can be selected uniformly, under a pre-determined repartition key (for example providing more candidates close to the middle of the pre-determined pulsation range than on the edges), or randomly in the pre-determined pulsation range.
A trigonometric function determination module 14 determines the trigonometric function y′ which best approximates the signal. Determining the pulsation w for interval ]ta; tb[ and the phase at time tb can be performed by minimizing a distance between the functions y and y′ over interval ]ta; tb[. A distance determination module 15 can be used to determine this distance.
The distance can be any suitable distance.
In particular, the distance determination module may affect a greater weight to points of the physiological signal closer to tb than far away from tb.
For example, for a given pulsation candidate wj, the distance determination module may solve the following least-square problem:
e[theta]=1/Nsumk((y[k]−y′[k])2*lambda[k])
where lambda is a forgetting factor vector, theta is the vector [A, B] that we are trying to optimize, and N is the number of points of vector y on the time interval ]ta; tb[ under study.
One defines D=[cos(wjt[i]), sin(wjt[i]] as a matrix of size (n, 2).
The distance determination module attempts at determining y′=D*theta by minimizing the error vector e[theta]=1/N*(y−D*theta)T*lambda*(y−D*theta), where “T” designates the transpose of a matrix.
The expression of theta which minimizes the error is theta=(DT*lambda0.5D)−1*lambda0.5y. In this expression “lambda0.5” designates the vector of coordinates the square roots of the coordinates of the vector lambda.
Practically speaking, the above minimization is performed for all pulsation candidate wj. The method first determines the best phase which brings the trigonometric function of pulsation wj closest to the measured signal. Then the error between the measured signal and the determined trigonometric functions is assessed for each trigonometric function, and the best trigonometric function is selected. The phase of the physiological signal at time tb is determined based on the phase of the best trigonometric function at time tb. In more details, this is performed according to the following scheme:
for time tb et for a range of pulsations wj:
Selection, among the various functions y′ of that for which the error with respect to y is minimum
The phase at time tb is the phase of the elected function y′ at tb.
The phase at time tb is thus determined not taking into account any signal after tb. This allows to determine the phase in real time.
At a later time td, the process can be performed again, to determine the phase of the measured signal at td. The calculation is performed for a time interval ]tc; td[. For example, the duration of the time interval is ]tc; td[ is the same as that of ]ta; tb[. The process is repeated rapidly, so that tb will be between tc and td. For example, the process is performed at a frequency of the order of the sampling frequency of the input signal. For example, it is performed at a frequency of at least 25 Hz. For an input EEG signal, it can be performed at a frequency over 200 Hz.
According to another example, one applies a recursive least square method to determine the unknowns to be determined.
As discussed above, the expression of theta which minimizes the error is theta=(DT*lambda0.5D)−1*lambda0.5y.
A recursive expression of this problem can be written as:
Theta(n+1)=(D(n+1)T*lambda0.5D(n+1)−1*lambda0.5y(n+1), where n is an integer, n+1 is the next integer for which the problem is to be solved, D(n+1) is the expression of matrix D at integer (n+1), and y(n+1) is the expression of vector y at integer (n+1).
In more details, this is performed according to the following scheme:
For time t et for a range of pulsations wj (in the present case, we know and use the solution to the problem at the previous time t−1):
The pulsation wj which minimizes the error will be selected.
The phase at time t will be selected as the phase for the selected signal.
In some embodiments, and for example when the signal is of poor quality, one may rely on the calculations performed for previous times when performing the calculations for the present time. This may be true notably for the respiratory signal. In particular, if the respiratory signal is obtained from acceleration measurements, the occurrence of artifact movements of the person from the acceleration measurements may impair the ability to properly detect the respiratory movement. Detection, for example from the accelerometer, of such artifact movements may classify the respiratory signal as being of poor quality.
According to some embodiments, a stimulation of the user may take into account the phase of the physiological signal determined in real-time.
According to an aspect, a computer program comprises instructions to execute the steps as described above when executed on a computer.
Number | Date | Country | Kind |
---|---|---|---|
17305703.5 | Jun 2017 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2018/065390 | 6/11/2018 | WO | 00 |