Physiological state control system, physiological state control method, and vibration wave calculation program product

Abstract

The invention provides a physiological state control system, a physiological state control method, and a vibration wave calculation program product capable of reducing burden on a living body and making the autonomic nervous system close to normal. A physiological state control system according to an embodiment of the present invention is provided with: a waveform measurement device that acquires waveform information relating to a living body; a physiological state control device that calculates a vibration wave applied to the living body on the basis of the waveform information acquired by the waveform measurement device; and a drive device that applies the calculated vibration wave from outside the living body to at least one of the carotid bifurcation and the vertebral artery of the living body.

Description

Technical Field

[0001] The present disclosure relates to a physiological state control system, a physiological state control method, and a vibration (oscillation) wave calculation program. Background Art

[0002] Patent Document 1 discloses a technique for normalizing the autonomic nervous system by sensing a person's activity and applying appropriate nerve stimulation based on the activity.

[0003] Prior Art Documents

[0004] Patent Document 1: Japanese Patent Application Laid-Open No. 2010-508969

[0005] Non-Patent Document 1: W Zong, T Heldt, GB Moody and RG Mark, "An Open-source Algorithm to Detect Onset of Arterial Blood Pressure Pulses", Computers in Cardiology 2003; 30: 259-262. [online], [retrieved on August 12, 2021], Internet <https: / / lcp.mit.edu / pdf / Zong03a.pdf> Summary of the Invention

[0006] In the method of Patent Document 1, a nerve stimulation signal is applied by an implantable medical device such as a pacemaker, and the burden on the living body caused by the installation of the device in the body and the like increases.

[0007] The present disclosure has been completed to solve such problems, and provides a physiological state control system, a physiological state control method, and a vibration wave calculation program that can reduce the burden on the living body and bring the autonomic nervous system closer to normal.

[0008] The physiological state control system according to the present embodiment includes: a waveform measuring device that acquires waveform information of pulse waves of the carotid artery bifurcation and the vertebral artery related to a living body and brain waves (electroencephalograms) of each part measured by the 10-20 method (system) or the like; a physiological state control device that calculates a vibration wave to be applied to the living body based on the waveform information acquired by the waveform measuring device; and a driving device that applies the calculated vibration wave from the outside of the living body to at least one of the carotid artery bifurcation and the vertebral artery of the living body. According to this configuration, the burden on the living body can be reduced and the autonomic nervous system can be brought closer to normal.

[0009] In order to regulate the activities of the circulatory organs, digestive organs, respiratory organs, etc., the autonomic nerves are constantly working 7×24 hours regardless of one's will. Since they are nerves that react automatically, they unconsciously regulate functions such as breathing, blood circulation, body temperature regulation, digestion, excretion, reproduction, and immunity, and are indispensable for maintaining life. When the autonomic nerves are hypertonic, symptoms of autonomic nerve disorders such as physical weakness, constipation or diarrhea, headache, hot flashes, palpitations, and paralysis will appear. It is believed that the control of these physiological states is a manifestation of brain function centered on the brainstem.

[0010] After investigation, the path of cerebral blood flow that provides the energy source for brain activity is as follows: The blood that inhales oxygen in the lungs is transported from the left atrium to the left ventricle, and then the blood is transported from the aortic arch to the whole body. The left and right common carotid arteries branched from the aortic arch branch into the internal carotid arteries at the carotid bifurcation and supply blood to the frontal lobe, temporal lobe, and parietal lobe. The left and right subclavian arteries branched from the aortic arch branch into the left and right vertebral arteries, and supply blood to the brainstem, cerebellum, and occipital lobe along the cervical vertebrae.

[0011] By measuring the waveforms of these parts using the waveform measuring device, information on the blood flow supplied by these arteries can be obtained. In addition, by observing the brain waves of each part of the brain, the activity state of the brain can be observed from both aspects of cerebral blood flow and brain waves.

[0012] In the above physiological state control system, the physiological state control device filters the obtained waveform information with at least one frequency band (frequency range), performs a Hilbert transform on the waveform information of the filtered frequency band, extracts the real part and the imaginary part of the complex waveform equation obtained by performing a Hilbert transform on the waveform information. If the real part is set as the instantaneous amplitude and the imaginary part is set as the instantaneous phase, the time differential of the instantaneous amplitude and the instantaneous phase, that is, the instantaneous frequency, is calculated during the period that can be regarded as a stable state physiologically, and the distributions of the instantaneous amplitude and the instantaneous frequency are calculated, and they are respectively distributed like a Gaussian distribution.

[0013] In the above physiological state control system, the physiological state control device may include at least any one of the frequency bands of 0.004 - 0.015 Hz (VLF2), 0.015 - 0.04 Hz (VLF1), 0.04 - 0.15 Hz (LF), 0.15 - 0.4 Hz (HF), 0.4 - 1.5 Hz (δ1), 1.5 - 4 Hz (δ2), 4 - 8 Hz (θ), 8 - 13 Hz (α), 13 - 30 Hz (β), and 30 - 100 Hz (γ) as the frequency band. In addition, the δ1, δ2, θ, α, β, γ frequency bands above 1.5 Hz (hertz) can be further divided into multiple divisions. With this configuration, the physiological state can be obtained from the brain waves and the pulse waves.

[0014] In the above physiological state control system, when autonomic nerve disorder symptoms occur, it is observed that the average value and variance value of the distribution of the amplitude and frequency in at least one of the frequency bands VLF2, VLF1, LF, HF, δ1, δ2, θ, α, β, and γ deviate from the normal distribution. The physiological state control device calculates to make the average value and variance value of the amplitude and frequency during autonomic nerve abnormality become the instantaneous amplitude and instantaneous frequency of the distribution of the amplitude and frequency during normal autonomic nerve. By applying their vibration or pressure to the driving device, the blood flow supplied to the brainstem and the brain is made close to the normal value, thereby being able to improve autonomic nerve disorder. According to this configuration, the physiological state can be made close to normal.

[0015] In the above physiological state control system, the physiological state control device may also calculate the vibration wave applied to the living body based on at least any one of the electroencephalogram, pulse wave, and pulse interval wave as the waveform information. In addition, it is also possible to confirm whether the physiological state is close to the normal value by examining the distribution shape of the instantaneous amplitude and instantaneous frequency of the electroencephalogram band waveform, and adjust the amplitude of driving the driving device. According to this configuration, the physiological state can be made close to the normal state.

[0016] The physiological state control method according to this embodiment includes: a step of obtaining waveform information of at least any one of the electroencephalogram and pulse wave related to the living body; a step of calculating a vibration wave applied to the living body based on the obtained waveform information; and a step of applying the calculated vibration wave from the outside of the living body to at least any one of the left and right carotid artery bifurcations and the left and right vertebral arteries (left and right cervical vertebrae) of the living body. According to this configuration, the burden on the living body can be reduced and the autonomic nervous system can be made close to normal.

[0017] In the above physiological state control method, the step of calculating the vibration wave includes: a step of filtering the obtained waveform information with at least one frequency band; a step of performing a Hilbert transform on the waveform information of the filtered frequency band; a step of calculating the instantaneous amplitude corresponding to the real part of the complex waveform equation obtained by performing a Hilbert transform on the waveform information, and the instantaneous frequency corresponding to the time differential value of the instantaneous phase corresponding to the imaginary part of the complex waveform equation; a step of calculating the distribution of the instantaneous value similar to the Gaussian distribution; and a step of calculating the vibration wave in such a way that the calculated distribution of the instantaneous value similar to the Gaussian distribution becomes the distribution of the instantaneous value similar to the Gaussian distribution in the normal physiological state, and having a unit for storing the distribution shape of the instantaneous amplitude and instantaneous frequency in the normal physiological state and the distribution shape of the instantaneous amplitude and instantaneous frequency in the state where the autonomic nerve is abnormal.

[0018] In the above physiological state control method, in the filtering step, at least any one of the frequency bands of 0.004 to 0.015 Hz (VLF2), 0.015 to 0.04 Hz (VLF1), 0.04 to 0.15 Hz (LF), 0.15 to 0.4 Hz (HF), 0.4 to 1.5 Hz (δ1), 1.5 to 4 Hz (δ2), 4 to 8 Hz (θ), 8 to 13 Hz (α), 13 to 30 Hz (β), and 30 to 100 Hz (γ) may be used as the frequency band.

[0019] In the above physiological state control method, as the waveform information, based on at least any one of electroencephalogram, pulse wave, and pulse interval waveform, the distribution of instantaneous amplitude and instantaneous frequency of each frequency band and each part is measured, the frequency band and part deviating from the distribution shape of the normal physiological state are determined, and a vibration waveform having the distribution of the instantaneous amplitude and instantaneous frequency of the frequency band in the normal physiological state is input to a driving device of any one of the left and right carotid bifurcations and the left and right vertebral arteries related to the frequency band and part, so that the blood flow in this part approaches the normal state, thereby being able to improve the abnormality of the autonomic nervous system.

[0020] When there are multiple frequency bands deviating from the distribution shape of the normal physiological state in a specific part, the vibration waveform of the normal physiological state input to the driving device may also be a waveform formed by overlapping the vibration waveforms of multiple frequency bands.

[0021] According to the present embodiment, it is possible to provide a physiological state control system, a physiological state control method, and a vibration wave calculation program that can reduce the burden on a living body and make the autonomic nervous system approach the normal physiological state.

[0022] The above and other objects, features, and advantages of the present disclosure will be more fully understood from the following detailed description and the accompanying drawings. Description of the Drawings

[0023] Figure 1 It is a diagram illustrating the configuration of the physiological state control system according to Embodiment 1.

[0024] Figure 2 It is a block diagram illustrating the physiological state control device according to Embodiment 1.

[0025] Figure 3 It is a flowchart illustrating the physiological state control method according to Embodiment 1.

[0026] Figure 4 It is a graph of waveform information related to a living body according to Embodiment 1, where the horizontal axis represents time and the vertical axis represents intensity.

[0027] Figure 5It is a flowchart of a method for calculating a vibration wave by the physiological state control device according to Exemplary Embodiment 1.

[0028] Figure 6 It is a chart of a signal waveform obtained by filtering a pulse wave as waveform information related to a living body in each frequency band. The horizontal axis represents time, and the vertical axis represents intensity.

[0029] Figure 7 It is a chart of a signal waveform obtained by filtering an electroencephalogram as waveform information related to a living body in each frequency band. The horizontal axis represents time, and the vertical axis represents intensity.

[0030] Figure 8 It is a diagram of each frequency band in the case of filtering waveform information related to a living body.

[0031] Figure 9 It is a chart of a signal waveform obtained by performing a Fourier transform on waveform information related to a living body. The horizontal axis represents frequency, and the vertical axis represents spectral power.

[0032] Figure 10 It is a diagram showing a complex waveform equation obtained by performing a Hilbert transform on waveform information related to a living body in the complex plane according to Exemplary Embodiment 1.

[0033] Figure 11A It is a diagram of the probability density distribution of instantaneous values and a Gaussian distribution according to Exemplary Embodiment 1.

[0034] Figure 11B It is a diagram of the probability density distribution of instantaneous values and a Gaussian distribution according to Exemplary Embodiment 1.

[0035] Figure 11C It is a diagram of the probability density distribution of instantaneous values and a Gaussian distribution according to Exemplary Embodiment 1.

[0036] Figure 11D It is a diagram of the probability density distribution of instantaneous values and a Gaussian distribution according to Exemplary Embodiment 1.

[0037] Figure 12A It is a diagram of the probability density distribution of the instantaneous phase difference in the HF frequency band according to Exemplary Embodiment 1.

[0038] Figure 12B It is a diagram of the probability density distribution of the instantaneous phase difference in the LF frequency band according to Exemplary Embodiment 1.

[0039] Figure 13A It is a diagram of the probability density distribution of the logarithmic amplitude and frequency in the θ frequency band according to Exemplary Embodiment 1.

[0040] Figure 13B It is a graph showing the probability density distribution of the logarithmic amplitude and frequency in the α band related to Example Embodiment 1.

[0041] Figure 13C It is a graph showing the probability density distribution of the logarithmic amplitude and frequency in the θ band related to Example Embodiment 1.

[0042] Figure 13D It is a graph showing the probability density distribution of the logarithmic amplitude and frequency in the α band related to Example Embodiment 1.

[0043] Figure 14A It is a graph showing the probability density distribution of the phase difference between the brain wave and the pulse wave related to Example Embodiment 1.

[0044] Figure 14B It is a graph showing the probability density distribution of the phase difference between the brain wave and the pulse wave related to Example Embodiment 1.

[0045] Figure 15A It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest (quiet) related to Example Embodiment 1.

[0046] Figure 15B It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest related to Example Embodiment 1.

[0047] Figure 15C It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest related to Example Embodiment 1.

[0048] Figure 15D It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest related to Example Embodiment 1.

[0049] Figure 16A It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest related to Example Embodiment 1.

[0050] Figure 16B It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the left brain during eyes-open rest related to Example Embodiment 1.

[0051] Figure 17A It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the right brain during eyes-open rest related to Example Embodiment 1.

[0052] Figure 17B It is a graph showing the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain wave of the right brain during eyes-open rest related to Example Embodiment 1.

[0053] Figure 17C It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-open rest related to Exemplary Embodiment 1.

[0054] Figure 17D It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-open rest related to Exemplary Embodiment 1.

[0055] Figure 18A It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-open rest related to Exemplary Embodiment 1.

[0056] Figure 18B It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-open rest related to Exemplary Embodiment 1.

[0057] Figure 19A It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0058] Figure 19B It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0059] Figure 19C It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0060] Figure 19D It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0061] Figure 20A It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0062] Figure 20B It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the left brain during eyes-closed rest related to Exemplary Embodiment 1.

[0063] Figure 21A It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest related to Exemplary Embodiment 1.

[0064] Figure 21B It is a graph showing the probability density distribution and Gaussian distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest related to Exemplary Embodiment 1.

[0065] Figure 21C It is a graph showing the probability density distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest in Example Embodiment 1 and the Gaussian distribution.

[0066] Figure 21D It is a graph showing the probability density distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest in Example Embodiment 1 and the Gaussian distribution.

[0067] Figure 22A It is a graph showing the probability density distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest in Example Embodiment 1 and the Gaussian distribution.

[0068] Figure 22B It is a graph showing the probability density distribution of the instantaneous frequency in the brain wave of the right brain during eyes-closed rest in Example Embodiment 1 and the Gaussian distribution.

[0069] Figure 23 It is a graph illustrating the complex transform expressions in each frequency band related to Example Embodiment 1 in a three-dimensional manner. The horizontal axis and the vertical axis in the complex plane respectively represent the and The axis orthogonal to the complex plane represents the frequency.

[0070] Figure 24 It is a graph illustrating the complex transform expressions in each frequency band related to Example Embodiment 1 in a three-dimensional manner. The horizontal axis and the vertical axis in the complex plane respectively represent the and The axis orthogonal to the complex plane represents the frequency.

[0071] Figure 25 It is a graph showing the relationship between the total hemoglobin concentration and the frequency. The horizontal axis represents the frequency, and the vertical axis represents the total hemoglobin concentration.

[0072] Figure 26 It is a graph showing the definition of the pulse wave waveform and the pulse interval related to Example Embodiment 2.

[0073] Figure 27 It is a graph showing the definition of the pulse interval wave related to Example Embodiment 2.

[0074] Figure 28 It is a graph showing the time series data of the instantaneous amplitude of the pulse wave in each frequency band in the normal physiological state.

[0075] Figure 29 It is a graph showing the time series data of the instantaneous frequency of the pulse wave in each frequency band in the normal physiological state. Specific Embodiments

[0076] Hereinafter, embodiments of the present disclosure will be described. However, the scope of protection of the present disclosure is not limited to the following embodiments. In addition, not all structures described in the embodiments are essential as means for solving the problems. For clarity, the following descriptions and drawings are appropriately omitted and simplified. In each drawing, the same reference numerals are assigned to the same elements, and repeated descriptions are omitted as needed.

[0077] (Embodiment 1)

[0078] A physiological state control system according to Embodiment 1 will be described. The physiological state control system of this embodiment, for example, acquires waveform information related to a living body such as an electroencephalogram waveform and a pulse wave waveform, and calculates a vibration wave applied to the living body based on the acquired waveform information. The pulse wave waveform can be a photoplethysmogram waveform, or a piezoelectric waveform obtained by detecting pulse pressure. In addition, it can also be a Doppler effect processing waveform of ultrasonic waves reflected on the arterial wall, or an arterial diameter time series data calculated based on a dynamic image of an arterial cross-section obtained by ultrasonic tomography. The electroencephalogram is obtained by simultaneously measuring multiple parts of the head using the 10-20 method or the like. The living body is, for example, a human being as a subject. For example, the physiological state control system acquires waveform information of a normal physiological state related to the living body, calculates characteristic quantities of a predetermined vibration wave with a set center frequency and amplitude, and records the characteristic quantities of the normal physiological state (the average values of the amplitudes and frequencies in each frequency band and the variance values of the amplitudes and frequencies). Moreover, the physiological state control system calculates the characteristic quantities of the vibration wave in the current physiological state, and applies a vibration wave from the outside of the living body to compensate for the characteristic quantities that are different from the characteristic quantities of the vibration wave in the normal physiological state (the average values of the amplitudes and frequencies in each frequency band and the variance values of the amplitudes and frequencies), thereby controlling the physiological state of the living body. Hereinafter, first, the <configuration of the physiological state control system> will be described. Then, the <physiological state control method> using the physiological state control system will be described.

[0079] <Configuration of the Physiological State Control System>

[0080] Figure 1 is a diagram illustrating the configuration of the physiological state control system according to Embodiment 1. As Figure 1 shown, the physiological state control system 1 includes an electroencephalogram detector 10, a pulse wave detector 20, a driving device 30, and a physiological state control device 50. Devices that acquire waveform information related to a living body, such as the electroencephalogram detector 10 and the pulse wave detector 20, are called waveform detectors.

[0081] The electroencephalogram detector 10 acquires the electroencephalogram of a living body as waveform information related to the living body. The electroencephalogram detector 10 includes a sensor 11 and a main body 12. The sensors 11 are arranged at a plurality of positions according to the 10-20 method or the like, and measurement is performed simultaneously in the main body 12. The sensor 11 is attached to the scalp of a person's head, for example, and senses the electroencephalogram information of a person from the outside of the living body. The electroencephalogram information of a person is voltage, for example. In addition, the sensor 11 can sense current, magnetic field, etc. in addition to voltage as the electroencephalogram information of a person. The sensor 11 is attached to the living body in a non-invasive manner.

[0082] The sensor 11 outputs the sensed electroencephalogram information to the main body 12 of the electroencephalogram detector 10. The main body 12 of the electroencephalogram detector 10 measures the time change of voltage or the like output from the sensor 11. The sensor 11 is connected to the main body 12 through a wired or wireless communication line. In addition, the main body 12 is connected to the physiological state control device 50 through a wired or wireless communication line. The main body 12 outputs the measured electroencephalogram to the physiological state control device 50.

[0083] The pulse wave detector 20 measures the pulse wave of a living body as waveform information related to the living body. The pulse wave is waveform information of a living body formed by the pulse interval, blood output, and physical characteristics of blood vessels. The pulse wave detector 20 includes a sensor 21 and a main body 22. The pulse wave detector 20 can be either a photoelectric type or a piezoelectric type. In the case of a photoelectric detector, near-infrared light with a light wavelength of 800 to 1000 nm can be used. The sensor 21 is attached to the skin on a person's neck (carotid artery bifurcation) or cervical vertebra (vertebral artery), for example, and measures the pulse wave information flowing into the human brain from the outside of the living body. Specifically, the sensor 21 can also be arranged at least at either the vicinity of the left and right carotid artery bifurcations or the vicinity of the left and right vertebral arteries (cervical vertebrae). The pulse wave information of a person is pulse pressure, for example. In addition, the sensor 21 can sense blood flow rate, etc. in addition to pulse pressure as the pulse wave information of a person. The sensor 21 is attached to the living body in a non-invasive manner.

[0084] The sensor 21 outputs the sensed pulse wave information to the main body 22 of the pulse wave detector 20. The main body 22 of the pulse wave detector 20 measures the time change of pulse pressure or the like output from the sensor 21. The sensor 21 is connected to the main body 22 through a wired or wireless communication line. In addition, the main body 22 is connected to the physiological state control device 50 through a wired or wireless communication line. The main body 22 outputs the measured pulse wave to the physiological state control device 50.

[0085] The driving device 30 applies the vibration wave calculated by the physiological state control device 50 from the outside of the living body to at least one of the vicinity of the left and right internal carotid arteries (arteries flowing into the brain from the carotid bifurcation) that supply blood to the brain and the vicinity of the left and right vertebral arteries (arteries flowing along the cervical vertebrae into the brainstem). The driving device 30 includes an oscillator 31 and a main body 32. The oscillator 31 is attached, for example, to the skin on a person's neck, collarbone, or cervical vertebra (through which the vertebral artery passes), and applies a vibration wave to the person. Specifically, the oscillator 31 may also be arranged at least at any one of the vicinity of the left and right carotid bifurcations and the vicinity of the vertebral artery (cervical vertebra). In this case, the driving device 30 applies the calculated vibration wave from the outside of the living body to at least any one of the left and right carotid bifurcations and the left and right vertebral arteries (cervical vertebra) of the living body. Thereby, the pulse wave of the living body is changed. One or more oscillators 31 may be attached to the living body. The oscillator 31 is attached to the living body in a non-invasive manner.

[0086] The oscillator 31 is connected to the main body 32 through a wired or wireless communication line. In addition, the main body 32 is connected to the driving device 30 through a wired or wireless communication line. The main body 32 receives the information of the vibration wave calculated by the physiological state control device 50 from the physiological state control device 50. The information of the vibration wave input to the driving device is, for example, the average value of the amplitude and frequency of the vibration wave and the variance value of the amplitude and frequency. In the vibration wave, the amplitude and frequency in at least one of the frequency bands VLF2, VLF1, LF, HF, δ1, δ2, θ, α, β, and γ are distributed like a Gaussian distribution. The main body 32 amplifies the information of the vibration wave received from the physiological state control device 50 to drive the oscillator 31. The oscillator 31 can apply vibrations of an arbitrary waveform overlapping multiple frequency bands to the living body and apply the vibration wave from the outside of the living body. In this way, the driving device 30 applies the vibration wave from the outside of the living body to at least any one of the left and right carotid bifurcations and the left and right vertebral arteries (cervical vertebra) of the living body.

[0087] The physiological state control device 50 stores the characteristic quantities of the vibration wave of the living body in the normal physiological state based on the waveform information obtained by the waveform measuring device. The characteristic quantities are the average value of the amplitude and frequency of the vibration wave and the variance value of the amplitude and frequency. The vibration wave is a waveform that has passed through a band-pass filter in the VLF2, VLF1, LF, HF, δ1, δ2, θ, α, β, and γ frequency bands. If the probability density distribution of its amplitude and its frequency is calculated during a period that can be regarded as a stable state physiologically, they will be distributed like a Gaussian distribution respectively. Moreover, when autonomic nervous system disorder symptoms occur, the same characteristic quantities are calculated, the differences from the characteristic quantities in the normal physiological state are detected, and a vibration waveform that compensates for the differences is generated. By inputting this vibration waveform into the driving device and driving the oscillator 31, the blood flow into the brain is made close to the normal state, and the autonomic nervous system disorder symptoms are improved.

[0088] Specifically, assume that due to autonomic neuropathy, the average of the amplitude distribution of a certain frequency band of the pulse wave of the right internal carotid artery changes from M0 to M1, the variance changes from S0 to S1, the average of the frequency distribution changes from m0 to m1, and the variance changes from s0 to s1. At this time, when an abnormal (abnormal) site with the same amplitude and frequency distribution in this frequency band is found in the electroencephalogram in the right temporal lobe, parietal lobe, and frontal lobe, while monitoring the degree of recovery of this abnormality, the oscillator 31 is placed on the right internal carotid artery to give vibrations that make the average of the amplitude distribution approach M0 from M1, the variance approach S0 from S1, and make the average of the frequency approach m0 from m1, and the variance approach s0 from s1. Accordingly, autonomic neuropathy symptoms can be improved. The reason for not suddenly giving the vibration distribution in the normal physiological state is to avoid burdening the organism. The same applies even if the arterial system that flows into the brain is different or the abnormal frequency band is different. In addition, when abnormalities occur simultaneously in multiple frequency bands, it is only necessary to overlap the vibration distribution waveforms that improve the abnormal frequency bands and give them to the desired oscillator.

[0089] Figure 2 is a block diagram illustrating the physiological state control device 50 according to Embodiment 1. As Figure 2 shown, the physiological state control device 50 includes a control unit 50a, a communication unit 50b, a storage unit 50c, an interface unit 50d, a waveform information acquisition unit 51, a filtering unit 52, a conversion unit 53, an instantaneous value calculation unit 54, a distribution calculation unit 55, and a vibration wave calculation unit 56. The control unit 50a, the communication unit 50b, the storage unit 50c, the interface unit 50d, the waveform information acquisition unit 51, the filtering unit 52, the conversion unit 53, the instantaneous value calculation unit 54, the distribution calculation unit 55, and the vibration wave calculation unit 56 respectively have functions as a control unit, a communication unit, a storage unit, an interface unit, a waveform information acquisition unit, a filtering unit, a conversion unit, an instantaneous value calculation unit, a distribution calculation unit, and a vibration wave calculation unit.

[0090] The physiological state control device 50 is an information processing device including a computer. The control unit 50a includes, for example, processors such as a CPU (Central Processing Unit), an MPU (Micro Processing Unit), an ECU (Electronic Control Unit), an FPGA (Field-Programmable Gate Array), and an ASIC (Application Specific Integrated Circuit). The control unit 50a functions as an arithmetic unit that performs control processing, arithmetic processing, and the like. In addition, the control unit 50a controls the operations of various components such as the communication unit 50b, the storage unit 50c, the interface unit 50d, the waveform information acquisition unit 51, the filtering unit 52, the conversion unit 53, the instantaneous value calculation unit 54, the distribution calculation unit 55, and the vibration wave calculation unit 56.

[0091] Each component of the physiological state control device 50 can be implemented, for example, by executing a program under the control of the control unit 50a. More specifically, each component can be implemented by the control unit 50a executing a program stored in the storage unit 50c. In addition, each component can also be implemented by recording the required program on an arbitrary non-volatile recording medium and installing it as needed. In addition, each component is not limited to being implemented based on software of a program, and can also be implemented by any combination of hardware, firmware, and software.

[0092] The communication unit 50b receives waveform information measured by a waveform measuring device such as an electroencephalogram measuring device 10 and a pulse wave measuring device 20 from the waveform measuring device. The waveform information is, for example, time series data of an electroencephalogram and a pulse wave. The communication unit 50b sends the information of the vibration wave calculated by the vibration wave calculation unit 56 to the driving device 30.

[0093] The storage unit 50c can have, for example, a storage device such as a memory or a hard disk. The storage device is, for example, a ROM (ReadOnly Memory) or a RAM (Random Access Memory). The storage unit 50c has a function of storing control programs, arithmetic programs, and the like executed by the control unit 50a. In addition, the storage unit 50c has a function of temporarily storing processing data and the like.

[0094] The storage unit 50c can also store waveform information such as an electroencephalogram and a pulse wave received by the communication unit 50b. The storage unit 50c can also store the vibration wave calculated by the vibration wave calculation unit 56. The storage unit 50c can also store various parameters used by the vibration wave calculation unit 56 to calculate the vibration wave.

[0095] The interface unit 50d is, for example, a User Interface. The interface unit 50d has input devices such as a keyboard, a touch panel, or a mouse, and output devices such as a display or a speaker. The interface unit 50d accepts operations for inputting data by a user (such as an operator) and outputs information to the user.

[0096] The waveform information acquisition unit 51 acquires waveform information obtained by the communication unit 50b from a waveform measuring device. The filtering unit 52 filters the acquired waveform information with at least one frequency band. The conversion unit 53 performs a conversion that complexifies the waveform information in the filtered frequency band. The conversion that complexifies the waveform information is, for example, a Hilbert transform. The instantaneous value calculation unit 54 calculates an instantaneous value including at least any one of an instantaneous amplitude, an instantaneous frequency, and an instantaneous phase difference. The instantaneous amplitude is an amplitude term corresponding to the real part of the complex waveform obtained by complexifying the waveform information of a pulse wave and an electroencephalogram. The instantaneous frequency is the time differential value of the phase term corresponding to the imaginary part. The instantaneous phase difference corresponds to the difference between the phase terms of the electroencephalogram and the pulse wave. The distribution calculation unit 55 calculates the probability density distribution of the instantaneous value. In addition, the distribution calculation unit 55 calculates a Gaussian distribution that fits the probability density distribution. The vibration wave calculation unit 56 calculates information on a vibration wave in such a manner that the calculated Gaussian distribution becomes a predetermined Gaussian distribution. The vibration wave calculation unit 56 outputs the calculated information on the vibration wave to the communication unit 50b. The information on the vibration wave includes an amplitude and a center frequency. The communication unit 50b outputs the information on the vibration wave to the drive device 30, the drive oscillator 31.

[0097] <Physiological state control method>

[0098] Next, a physiological state control method will be described. Figure 3 is a flowchart illustrating the physiological state control method according to Embodiment 1. As Figure 3 shown, the physiological state control method includes a step of acquiring waveform information related to a living body (step S11), a step of calculating a vibration wave applied to the living body (step S12), and a step of applying the vibration wave from the outside of the living body (step S13).

[0099] As Figure 3 shown in step S11, first, waveform information related to a living body is acquired. For example, waveform measuring devices such as an electroencephalogram measuring device 10 and a pulse wave measuring device 20 acquire an electroencephalogram and a pulse wave as waveform information related to a living body.

[0100] Figure 4 is a graph of waveform information related to a living body according to Embodiment 1, where the horizontal axis represents time and the vertical axis represents intensity. In Figure 4Among them, as an example of waveform information related to an organism, time series data of a pulse wave is shown. In addition, the waveform information related to the organism is not limited to the time series data of the pulse wave, but can also be the time series data of the brain wave, or the time series data of the pulse interval (resampled at a desired frequency). The pulse wave detector 20 acquires the time series data of the pulse wave via, for example, the sensor 21 disposed at the bifurcation of the left and right carotid arteries or the left and right vertebral arteries of the subject. For example, for a subject in a stable state, the pulse wave is measured for several minutes to several tens of minutes. The sampling frequency is, for example, 500 Hz. In addition, any frequency between 10 and 1000 Hz can be adopted as the sampling frequency. The pulse wave data obtained in this way becomes Figure 4 such a time series waveform. The waveform detector outputs the acquired waveform information to the physiological state control device 50.

[0101] Next, as shown in step S12, the vibration wave applied to the organism is calculated. The physiological state control device 50 calculates the vibration wave applied to the organism based on the acquired waveform information. The calculation method of the vibration wave will be described later. The physiological state control device 50 outputs the information of the calculated vibration wave to the drive device 30.

[0102] Next, as shown in step S13, the calculated vibration wave is applied from the outside of the organism. For example, the drive device 30 causes the oscillator 31 to generate a vibration wave based on the information of the vibration wave output from the physiological state control device 50. As a result, the oscillator 31 applies the calculated vibration wave from the outside of the organism to at least any one of the left carotid artery bifurcation on the left side of the organism, the right carotid artery bifurcation on the right side, the left vertebral artery on the left side, and the right vertebral artery on the right side.

[0103] <Calculation method of vibration wave>

[0104] Next, the calculation method of the vibration wave in step S12 described above will be described. Figure 5 is a flowchart illustrating the method of calculating the vibration wave performed by the physiological state control device 50 according to Embodiment 1. As Figure 5As shown, the method for calculating a vibration wave includes: a waveform information acquisition step of acquiring waveform information related to a living body (step S21); a filtering step of filtering the acquired waveform information using at least one frequency band (step S22); a transformation step of transforming the waveform information of the filtered frequency band into a complex form (step S23); an instantaneous value calculation step of calculating an instantaneous value including at least any one of an instantaneous amplitude, an instantaneous frequency, and an instantaneous phase (step S24); a distribution calculation step of calculating a distribution of the instantaneous value (step S25); and a vibration wave calculation step of calculating a vibration wave applied to the living body in such a manner that the calculated distribution becomes a predetermined distribution (step S26). Calculations are performed in advance in a normal physiological state through steps S21 to S25, and the numerical values are stored in the storage unit 50c. Hereinafter, each step will be described.

[0105] <Waveform information acquisition step>

[0106] The waveform information acquisition unit 51 of the physiological state control device 50 acquires waveform information received by the communication unit 50b from waveform measuring devices such as the electroencephalogram measuring device 10 and the pulse wave measuring device 20. The waveform information acquisition unit 51 acquires an electroencephalogram waveform and a pulse wave waveform based on, for example, time series data of the electroencephalogram and the pulse wave.

[0107] <Filtering step>

[0108] Figure 6 is a chart showing a signal waveform obtained by filtering a pulse wave as waveform information related to a living body in each frequency band according to Embodiment 1. The horizontal axis represents time, and the vertical axis represents intensity. In Figure 6 HF, LF, VLF1, and VLF2 are shown as each frequency band. In addition, the pulse interval waveform described later is also shown in Figure 6 Figure 7 is a chart showing a signal waveform obtained by filtering an electroencephalogram as waveform information related to a living body in each frequency band according to Embodiment 1. The horizontal axis represents time, and the vertical axis represents intensity. In Figure 7 the frequency band waveforms of the electroencephalogram (δ1, δ2, θ, α, β, γ) are shown as frequency bands shorter than the heartbeat interval. In addition, BVP (Blood Volume Pulse) is also shown in Figure 7 Figure 8 is a diagram of each frequency band in the case of filtering waveform information related to a living body according to Embodiment 1. In Figure 8 the approximate center frequencies in several frequency bands are also shown.

[0109] As Figure 6 and Figure 7 ​​As shown, the filtering unit 52 of the physiological state control device 50 filters the acquired waveform information using at least one frequency band. Specifically, the filtering unit 52 of the physiological state control device 50 filters the waveform information acquired by the waveform information acquisition unit 51 for each frequency band.

[0110] As Figures 6 to 8 shown, as the frequency bands for filtering, for example, VLF2 (0.004 - 0.015 Hz), VLF1 (0.015 - 0.04 Hz), LF (0.04 - 0.15 Hz), HF (0.15 - 0.4 Hz), δ1 (0.4 - 1.5 Hz), δ2 (1.5 - 4 Hz), θ (4 - 8 Hz), α (8 - 13 Hz), β (13 - 30 Hz), and γ (30 - 100 Hz) etc. are selected. Here, VLF1 and VLF2 are frequency bands associated with autonomic nerve functions such as body temperature regulation, digestion, excretion, reproduction, and immunity, HF is the respiratory variation band, and LF is the blood pressure variation band. Signals in frequency bands higher than them are loaded on the pulse wave flowing into the brain from the internal carotid artery at the carotid bifurcation and the vertebral artery branched from the subclavian artery, including the frequency bands of δ waves (0.4 - 4.0 Hz), θ waves (4.0 - 8.0 Hz), α waves (8.0 - 13.0 Hz), β waves (13.0 - 30.0), and γ waves (30 - 100 Hz) in brain waves.

[0111] Figure 9 is a chart illustrating the signal waveform obtained by performing Fourier transform on the waveform information measured by a commonly commercially available pulse meter. The horizontal axis represents frequency, and the vertical axis represents spectral power. Since a high-pass filter is added near 0.1 Hz, the signals of frequencies below 0.1 Hz are lost. In addition, since green light is irradiated onto the blood vessels and the amount of reflected light is used as the pulse wave signal, only the information of the capillaries on the skin surface can be obtained and it cannot be used in this embodiment. For use in this embodiment, it is necessary to observe the pulse wave signal of the artery at a depth of 2 - 3 cm from the skin with a frequency band width of at least 0.001 - 100 Hz using near-infrared light (wavelength 800 - 1000 nm). Furthermore, in order to observe the pulse wave of the vertebral artery that covers a wide frequency band of the living body and is deeper than the carotid artery, it is necessary to observe the pulse wave signal of the artery at a depth of 2 - 5 cm from the skin with a frequency band width of at least 0.0001 - 400 Hz using near-infrared light (wavelength 800 - 1600 nm).

[0112] Regarding the brain wave waveform, the display of the figure is also omitted, but filtering processing can be performed in the same frequency band as the pulse wave waveform.

[0113] <Transformation step>

[0114] Next, the physiological state control device 50 performs a Hilbert transform that complexifies the waveform information of each filtered frequency band. Specifically, the transform unit 53 of the physiological state control device 50 performs a Hilbert transform on the time series data of the waveforms of each filtered frequency band shown in Figure 6 and Figure 7 and transforms the waveform information such as the waveforms shown in into a complex waveform equation as shown in the following formula (1).

[0115]

[0116] where a k (t) is the instantaneous value of the logarithmic amplitude of the k-band waveform, and ψ k (t) is the instantaneous value of the phase of the k-band waveform.

[0117] Figure 10 is a diagram showing the complex waveform equation obtained by performing a Hilbert transform on the waveform information related to the living body in Embodiment 1 in the complex plane. As shown in Figure 10 , the waveform information (vibration) represented by the complex waveform equation can be represented on the complex plane. When the horizontal axis in the complex plane is set to and the vertical axis is set to , a k (t) represents the radius centered on the origin, and ψ k (t) represents the angle with the horizontal axis. Here, k = 1, 2, 3, 4... represents each frequency band.

[0118] Perform a Hilbert transform on each frequency band waveform of the pulse wave waveform and the brain wave waveform.

[0119] <Instantaneous value calculation step>

[0120] Next, the instantaneous value calculation unit 54 of the physiological state control device 50 calculates instantaneous values. The instantaneous values include at least any one of the instantaneous logarithmic amplitude, the time derivative of the instantaneous phase, i.e., the instantaneous frequency, and the instantaneous phase difference. Here, the instantaneous phase difference is the difference between the instantaneous phases of the pulse wave and the brain wave, or the difference between the instantaneous phases of the pulse wave and the heart rate interval. The instantaneous logarithmic amplitude corresponds to the logarithmic real part of the complex waveform equation obtained by performing a Hilbert transform on the waveform information. The instantaneous frequency corresponds to the time derivative of the logarithmic imaginary part of the complex waveform equation. The instantaneous phase difference corresponds to the difference between the imaginary terms of the logarithms of the complex transform expressions of the brain wave and the pulse wave. Specifically, the instantaneous logarithmic amplitude, the instantaneous frequency, and the instantaneous phase difference are defined by the following formulas (2), (3), and (4), respectively.

[0121] Performing the processing of formula (1) on the pulse wave data can obtain the instantaneous logarithmic amplitude and the instantaneous frequency of formulas (2) and (3).

[0122] ak (t)(2)

[0123] ω k (t) = dψ k (t) / dt (3)

[0124] By processing the brain wave data according to equation (1), the instantaneous logarithmic amplitude and instantaneous frequency of equations (4) and (5) can be obtained.

[0125] A k (t)(4)

[0126] Ω k (t) = dΨ k (t) / dt (5)

[0127] The instantaneous phase difference between the time series data of the brain wave and the pulse wave is shown in formula (6).

[0128] θ k (t) = Ψ k (t)-ψ k (t) (6)

[0129] This is a difference that represents the phase relationship of how the blood flow supplied to the brain is consumed by each center in the brain, and is therefore an important feature quantity for understanding the physiological state of each center.

[0130] In θ k When the distribution of (t) follows a random process, it can be approximated by the Von Mises distribution shown in equation (7). j (κ) is the deformed Bessel function of the first kind of order j.

[0131] f(θ) = exp(κ cos(θ-μ)) / 2πI0(κ) (7)

[0132] I j (κ) = (κ / 2) j Σ(κ 2 / 4) i / i!Γ(j+i+1) (8)

[0133] Thus, the instantaneous value calculation unit 54 calculates the instantaneous logarithmic amplitude a including the real part of the exponential part of the complex waveform equation obtained by Hilbert transforming the waveform information. k (t) and A k (t), the instantaneous frequency ω equivalent to the time differential value of its imaginary part k (t) and Ω k (t) and the instantaneous phase difference θ corresponding to the phase difference between the frequency band waveforms of the brain wave and the pulse wave k(t).

[0134] <Distributed calculation>

[0135] Next, the distributed calculation unit 55 of the physiological state control device 50 calculates the distribution of the instantaneous values of the pulse wave and the brain wave. Specifically, the distributed calculation unit 55 calculates the instantaneous logarithmic amplitudes a k (t) and A k (t), the instantaneous frequencies ω k (t) and Ω k (t), and the probability density distributions of the instantaneous phase differences θ k (t). In addition, the distributed calculation unit 55 obtains, in the normal physiological state and the state with autonomic nervous disorder symptoms, the instantaneous logarithmic amplitudes a k (t) and A k (t), the instantaneous frequencies ω k (t) and Ω k (t), and the instantaneous phase differences θ k (t) that are fitted to a unimodal distribution, such as a Gaussian probability density distribution, within a desired time interval that can be regarded as a stable physiological state. Here, since the value of the instantaneous phase difference is limited to the range of -π to +π, for example, the von Mises distribution is applied.

[0136] Figures 11A to 11D is a diagram showing the probability density distribution of the instantaneous value and the Gaussian distribution according to Embodiment 1. The bar graph represents the probability density distribution, and the solid line represents the Gaussian distribution. Figure 11A shows the probability density distribution of the instantaneous logarithmic amplitude of HF and the Gaussian distribution, Figure 11B shows the probability density distribution of the instantaneous logarithmic amplitude of LF and the Gaussian distribution. Figure 11C shows the probability density distribution of the instantaneous logarithmic frequency of HF and the Gaussian distribution, Figure 11D shows the probability density distribution of the instantaneous frequency of LF and the Gaussian distribution. The mean and variance values of each Gaussian distribution are also shown. The mean and variance values are also shown in the following Gaussian distribution diagrams. It can be seen that the logarithmic amplitude and frequency of one oscillator are distributed like a Gaussian distribution for each frequency band. Figures 12A to 12B represents the probability density distribution of the instantaneous phase difference between the HF and LF frequency bands according to Embodiment 1, and the solid line represents the von Mises distribution. Figures 13A to 13D and Figures 14A to 14BShows the probability density distributions of the logarithmic amplitudes, frequencies, and phase differences between brain waves and pulse waves in the θ and α frequency bands related to Embodiment 1. In the figure, the solid lines representing the logarithmic amplitudes and frequencies are Gaussian distributions, and the solid line representing the phase difference represents the von Mises distribution. For frequencies higher than 12 Hz, such as the β and γ frequency bands, although not shown in the figure, the logarithmic amplitudes and frequencies of each signal waveform are distributed like Gaussian distributions respectively.

[0137] Figures 15A to 15D and Figures 16A to 16B Is a graph exemplifying the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain waves (frontal lobe) of the left brain during eyes-open rest related to Embodiment 1. Figure 15A Shows HF, Figure 15B Shows LF, Figure 15C Shows the θ wave, Figure 15D Shows the α wave. Figure 16A Shows the β wave, Figure 16B Shows the γ wave. Figures 17A to 17D and Figures 18A to 18B Is a graph exemplifying the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain waves of the right brain during eyes-open rest related to Embodiment 1. Figure 17A Shows HF, Figure 17B Shows LF, Figure 17C Shows the θ wave, Figure 17D Shows the α wave. Figure 18A Shows the β wave, Figure 18B Shows the γ wave.

[0138] Figures 19A to 19D and Figures 20A to 20B Is a graph exemplifying the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain waves (frontal lobe) of the left brain during eyes-closed rest related to Embodiment 1. Figure 19A Shows HF, Figure 19B Shows LF, Figure 19C Shows the θ wave, Figure 19D Shows the α wave. Figure 20A Shows the β wave, Figure 20B Shows the γ wave. Figures 21A to 21D and Figures 22A to 22B Is a graph exemplifying the probability density distribution of the instantaneous frequency and the Gaussian distribution in the brain waves of the right brain during eyes-closed rest related to Embodiment 1. Figure 21A Shows HF, Figure 21B Shows LF, Figure 21C Shows the θ wave, Figure 21D Shows the α wave. Figure 22A Shows the β wave, Figure 22B Shows the γ wave.

[0139] Such as Figures 15A to 15D 、 Figures 16A to 16B 、 Figures 17A to 17D andFigures 18A to 18B As shown, when the eyes are open and at rest, it is generally a waking state and a state of reduced parasympathetic nerve activity. Thus, when the eyes are open and at rest, it is a so-called energetic state. For example, as Figure 15C and Figure 17C shown, the central frequency of the theta wave when the eyes are open and at rest is shifted towards the lower frequency side compared to the central frequency when the eyes are closed and at rest as shown in Figure 19C and Figure 21C shown. Additionally, as Figure 15D and Figure 17D shown, the central frequency of the alpha wave when the eyes are open and at rest is shifted towards the higher frequency side compared to the central frequency when the eyes are closed and at rest as shown in Figure 19D and Figure 21D shown.

[0140] On the other hand, as Figures 19A to 19D , Figures 20A to 20B , Figures 21A to 21D as well as Figures 22A to 22B shown, when the eyes are closed and at rest, it is generally a resting state and a state of increased parasympathetic nerve activity. Thus, when the eyes are closed and at rest, it is a so-called relaxed state. For example, as Figure 19C and Figure 21C shown, the central frequency of the theta wave when the eyes are closed and at rest is shifted towards the higher frequency side compared to the central frequency when the eyes are open and at rest as shown in Figure 15C and Figure 17C shown. Additionally, as Figure 19D and Figure 21D shown, the central frequency of the alpha wave when the eyes are closed and at rest is shifted towards the lower frequency side compared to the central frequency when the eyes are open and at rest as shown in Figure 15D and Figure 17D shown.

[0141] <Vibration Wave Calculation>

[0142] Next, the physiological state control device 50 checks the difference in the characteristic quantities of the distribution of each frequency band between when autonomic imbalance occurs and when in a normal physiological state, and calculates a vibration wave that makes this distribution approach a normal distribution. For example, assume that due to autonomic imbalance, the average of the amplitude distribution of a certain frequency band of the pulse wave in the right internal carotid artery changes from M0 to M1, the variance changes from S0 to S1, the average of the frequency distribution changes from m0 to m1, and the variance changes from s0 to s1. At this time, when an abnormal part with the same amplitude and frequency distribution in this frequency band is found in the brain waves in the right temporal lobe, parietal lobe, and frontal lobe, while monitoring the degree of recovery of this abnormality, the oscillator 31 is placed on the right internal carotid artery part to give a vibration that makes the average of the amplitude distribution approach M0 from M1, the variance approach S0 from S1, and makes the average of the frequency approach m0 from m1, and the variance approach s0 from s1. By doing so, autonomic imbalance symptoms can be improved. The reason for not suddenly giving the vibration distribution in the normal physiological state is to avoid burdening the living body. The same applies even if the arterial system through which blood flows into the brain is different or the frequency band where the abnormality occurs is different. The vibration wave can also be a wave formed by overlapping for each frequency band (k = 1, 2, 3, 4...). The vibration wave calculation unit 56 can also calculate the amplitude and center frequency of the vibration wave of a predetermined Gaussian distribution as the amplitude and center frequency of the vibration wave applied to the living body.

[0143] In Figure 28 , Figure 29 represent samples of the time series data of the instantaneous amplitude and instantaneous frequency of the pulse wave of each frequency band in the normal physiological state. The instantaneous amplitude and instantaneous frequency data of the pulse wave are obtained at time intervals such that the distribution becomes Gaussian-like in a physiologically stable state for each frequency band, and time series data in which the instantaneous amplitude and instantaneous frequency change continuously are generated using this waveform. At this time, at least the second-order time differential values are made consistent and connected at the start point and end point of the time series data, and a vibration waveform is circularly generated using this sample data. When the average value and variance value of this sample data are different from the desired values, the difference in the average value is added and the ratio of the variance value is multiplied so that the average value and variance value become the desired values. When applying a vibration waveform to multiple frequency bands, the vibration waveform is calculated for each frequency band according to the above steps, a waveform formed by overlapping them is generated, and input to the drive device 30.

[0144] The vibration waveform in the normal state can also be approximated by a mathematical function so as to have the desired average value and variance value and input to the drive device 30. In this case, at least the second-order time differential values need to be continuous.

[0145] In the case where the autonomic imbalance is relatively mild and the difference from the normal state is small, the vibration waveform of the frequency band with a difference can also be used as the waveform in the normal physiological state and applied to the oscillator 31 from the beginning.

[0146] If the calculated vibration wave is applied to the left and right carotid bifurcations or the vertebral arteries branching from the left and right subclavian arteries, the pulse pressure vibration wave can be loaded onto the cerebral blood flow via the carotid arteries or the vertebral arteries. Brain activities include the activities of the autonomic nervous system centers in the brainstem in the frequency bands of VLF2 (0.004 - 0.015 Hz), VLF1 (0.015 - 0.04 Hz), LF (0.04 - 0.15 Hz), and HF (0.15 - 0.4 Hz), and the activities of various centers in the brain (cortex, limbic system) in the frequency bands of delta wave (0.4 - 4.0 Hz), theta wave (4.0 - 8.0 Hz), alpha wave (8.0 - 13.0 Hz), beta wave (13.0 - 30.0), and gamma wave (30 - 100 Hz). Therefore, the physiological state control system 1 measures the brain waves and generates a vibration wave that becomes a pulse pressure signal for assisting the brain wave activities in each frequency band. Moreover, the physiological state control system 1 applies the generated vibration wave from the carotid bifurcation to the frontal lobe, temporal lobe, and parietal lobe via the internal carotid artery, and applies it to the cerebellum and occipital lobe via the brainstem from the vertebral artery. Thus, the physiological state control system 1 can assist brain activities and assist and activate the brain activities in each frequency band.

[0147] In addition, as an example of calculating the vibration wave, it can also be calculated based on the equal - energy surface of 1 / f (also called the 1 / f spectrum). Hereinafter, the method of calculating based on the equal - energy surface of 1 / f will be described.

[0148] <1 / f equal - energy surface>

[0149] Figure 23 、 Figure 24 is a diagram illustrating the complex transformation formula in each frequency band related to Embodiment 1 in a three - dimensional manner. The horizontal axis and the vertical axis in the complex plane respectively represent the and The axis orthogonal to the complex plane represents the frequency (energy axis). For example, k = 4 in the frequency band longer than the heartbeat interval is represented as HF (respiration), k = 3 is represented as LF (blood pressure), k = 2 is represented as VLF1 (possibly related to the autonomic nerve), and k = 1 is represented as VLF2 (possibly related to the autonomic nerve). For the brain wave frequency band, k = 5 is set as the delta1 wave, k = 6 is set as the delta2 wave, k = 7 is set as the theta wave, k = 8 is set as the alpha wave, etc. In this way, the moving radius in the complex plane becomes the logarithmic amplitude of each frequency - band waveform, the energy axis becomes the frequency axis, and the rotation direction around the energy axis is Figure 23 in the case of Figure 24In this case, it is equivalent to the phase angle between the pulse wave and the brain wave. Moreover, it can be seen that on each axis, it is distributed like a Gaussian distribution and rotates around the energy axis at the same time. In this way, a depiction can be obtained as if electrons in each energy level rotate in orbits like electron clouds represented by the probability of existence around the atomic nucleus. Moreover, the Gaussian-like distribution of each frequency band similar to the electron cloud is formed along the constant energy surface of 1 / F.

[0150] Therefore, the vibration wave calculation unit 56 of the physiological state control device 50 can also use the Gaussian distribution along the constant energy surface of 1 / F as the predetermined Gaussian distribution. That is, the vibration wave calculation unit 56 can also calculate the vibration wave in such a way that the Gaussian distribution of the acquired waveform information becomes the Gaussian distribution along the constant energy surface of 1 / F.

[0151] Figure 25 It shows the measurement result of cerebral blood flow measured by optical topography technology (NIRS) measured in the frontal lobe. This is a graph illustrating the relationship between the total hemoglobin concentration (o2hb+hhb) and the frequency. The horizontal axis represents the frequency, and the vertical axis represents the total hemoglobin concentration. As Figure 25 shown, the spectrum of the total hemoglobin concentration (o2hb+hhb) is carried on the 1 / f line. Therefore, it can be considered that the variation spectrum of blood flow is carried on 1 / f.

[0152] When vibrating the blood vessel at a frequency below HF (below 0.4 Hz), the smaller the frequency, the larger the amplitude of the vibration that must be imparted. Therefore, if such a large-amplitude vibration is imparted to the blood vessel, it may cause blood flow to stop, so such an action cannot be performed. To avoid this, by performing AM modulation on the low-frequency vibration waveform and then loading it on the vibration waveform of about several Hz to 10 Hz, the low-frequency signal can be loaded on the arterial blood flow without stopping the arterial blood flow.

[0153] In addition, the physiological state control system 1 can normalize the autonomic nervous system for abnormalities in physiological states such as intra-day variations and inter-day variations by expanding and processing waveform information in the same way as HF, LF, VLF1, and VLF2 to lower frequency bands such as ULF1 (1.5 to 4 mHz), ULF2 (0.4 to 1.5 mHz), etc.

[0154] The physiological state control system 1, for example, sets the amplitude and the center frequency for each frequency band of the pulse wave based on the heartbeat variation, and calculates the vibration wave with a Gaussian-like fluctuation loaded on the 1 / f spectrum. Moreover, by applying the calculated vibration wave to the carotid artery bifurcation or the vertebral artery, the pulse pressure vibration of the 1 / f spectrum can be imparted to the cerebral artery flow. Moreover, by optimizing the vibration wave while feeding back the pulse wave and the brain wave, the autonomic nerve or brain activity can be nursed or assisted to activate the adjustment of the autonomic nerve or the brain wave activity.

[0155] (Embodiment 2)

[0156] Next, the physiological state control system according to Embodiment 2 will be described. The physiological state control system of this embodiment only considers the fluctuations of the autonomic nervous system that are longer than the heartbeat cycle as waveform information related to the living body. Therefore, it does not process the brain wave frequency band above 0.4 Hz, but only uses the pulse wave waveform. The pulse wave waveform can be a photoplethysmogram waveform, a piezoelectric waveform based on pulse pressure, a Doppler effect processing waveform of ultrasonic waves reflected on the artery surface, or arterial diameter time series data calculated from the dynamic image of the arterial cross-section based on ultrasonic tomography. Since the signals in the frequency band below 0.4 Hz are generated by the heartbeat fluctuations of the heart under the control of the brainstem, the pulse wave can be obtained at any position on the body. For the photoplethysmogram, it is not necessary to use near-infrared light, and ordinary green light can be used. However, it is necessary to use a signal waveform near 0.1 Hz without a high-pass filter like a commercially available pulse wave sensor.

[0157] <Waveform information acquisition step>

[0158] The waveform information acquisition unit 51 of the physiological state control device 50 acquires the waveform information received by the communication unit 50b from a waveform measuring device such as the pulse wave detector 20. The waveform information acquisition unit 51 obtains the pulse wave waveform based on the time series data of the pulse wave.

[0159] Figure 26 It is a diagram illustrating the definitions of the pulse wave waveform and the pulse interval according to Embodiment 2. Figure 27 It is a diagram illustrating the definition of the pulse interval wave according to Embodiment 2. As Figure 26 shown, by using the pulse wave rising position detection algorithm proposed in Non-Patent Document 1 to calculate the pulse interval (PPI: Peak-Peak Interval) from the pulse wave waveform data B(t), and plotting the rising time of the pulse wave and the pulse interval for spline interpolation, the Figure 27 shown pulse interval wave R(t) can be obtained. In view of this, the Figure 27 shown pulse interval wave can also be used as waveform information for the calculation of the vibration wave. That is, the physiological state control device 50 can also calculate the vibration wave applied to the living body based on the pulse interval wave as waveform information.

[0160] <Filtering step>

[0161] As frequency bands for filtering the pulse wave waveform and the pulse interval waveform, for example, VLF2 (0.004 to 0.015 Hz), VLF1 (0.015 to 0.04 Hz), LF (0.04 to 0.15 Hz), and HF (0.15 to 0.4 Hz) are selected. Here, VLF1 and VLF2 are variation bands associated with autonomic nerve functions such as body temperature regulation, digestion, excretion, reproduction, and immunity, HF is a respiratory variation band, and LF is a blood pressure variation band.

[0162] If the same processing as in Embodiment 1 is performed, the instantaneous logarithmic amplitude of Equation (9) and the instantaneous frequency of Equation (10) can be obtained from the pulse wave data. In addition, the symbols in the equations of Embodiment 2 may have meanings different from those of the symbols in the equations of Embodiment 1.

[0163] a k (t)(9)

[0164] ω k (t) = dψ k (t) / dt (10)

[0165] The instantaneous logarithmic amplitude of Equation (11) and the instantaneous frequency of Equation (12) can be obtained from the pulse interval data.

[0166] A k (t)(11)

[0167] Ω k (t) = dΨ k (t) / dt (12)

[0168] As the instantaneous phase difference between the pulse interval waveform and the time series data of the pulse wave, Equation (13) can be obtained.

[0169] θ k (t) = Ψ k (t)-ψ k (t) (13)

[0170] In this way, the instantaneous value calculation unit 54 calculates the instantaneous values including the instantaneous logarithmic amplitude a k (t) and A k (t), the instantaneous frequencies ω k (t) and Ω k (t), which correspond to the time differential value of the imaginary part of the exponential part of the complex waveform equation obtained by performing the Hilbert transform on the waveform information, and the phase difference θ k (t) between the pulse interval and the pulse wave.

[0171] <Distribution calculation>

[0172] Next, the distribution calculation unit 55 of the physiological state control device 50 calculates the distribution of the instantaneous values of the pulse wave and the pulse interval waveform. Specifically, the distribution calculation unit 55 calculates the instantaneous logarithmic amplitudes a k (t) and A k (t), the instantaneous frequencies ω k (t) and Ω k (t), and the probability density distributions of the instantaneous phase difference θ k (t) between the pulse wave and the pulse interval, fit these probability density distributions to, for example, a Gaussian distribution, obtain the distribution shapes of the respective instantaneous logarithmic amplitudes, instantaneous frequencies, and instantaneous phase differences, and save them in the storage unit 50c. Here, since the value of the instantaneous phase difference is limited to the range of -π to +π, for example, the von Mises distribution is applied.

[0173] <Vibration wave calculation>

[0174] Next, the physiological state control device 50 checks the differences in the characteristic quantities of the distributions in each frequency band between the occurrence of autonomic neuropathy and the normal physiological state, and calculates a vibration wave that makes the distribution close to the normal distribution. For example, it is assumed that due to autonomic neuropathy, the mean of the amplitude distribution of a certain frequency band of the pulse wave of VLF2 changes from m0 to m1, the variance changes from S0 to S1, the mean of the frequency distribution changes from m0 to m1, and the variance changes from s0 to s1. Although the heartbeat interval variation in the frequency band below 0.4 Hz is caused by the brainstem function, since the brain wave is not monitored, the physiological characteristic quantities of each frequency band can be calculated by measuring the pulse wave at the fingertips or the like. While monitoring these values, the oscillator 31 is placed at the left and right vertebral arteries supplying blood to the brainstem or the carotid artery bifurcation supplying blood to the frontal lobe, temporal lobe, and parietal lobe, and vibrations are given to make the mean of the amplitude distribution close from M1 to M0, the variance close from S1 to S0, and the mean of the frequency close from m1 to m0 and the variance close from s1 to s0. Thus, autonomic neuropathy symptoms can be improved. The reason for not suddenly giving the vibration distribution in the normal physiological state is to avoid burdening the organism. The same applies even if the arterial system supplying blood to the brain is different or the abnormal frequency band is different. The vibration wave can also be a wave formed by overlapping for each frequency band (k = 1, 2, 3, 4). The vibration wave calculation unit 56 can also calculate the amplitude and center frequency of a predetermined Gaussian distribution vibration wave as the amplitude and center frequency of the vibration wave applied to the organism.

[0175] When applying vibrations of less than HF (less than 0.4 Hz) to blood vessels, it is possible to apply vibrations with a large amplitude to the blood vessels, resulting in blood flow stoppage. To avoid this, a signal obtained by AM modulating a vibration waveform of a low frequency below HF with a vibration waveform of several Hz to about 10 Hz is input to a driving device, and vibrations are applied from the outside to the vicinity of the vertebral arteries of the left and right cervical vertebrae or the vicinity of the bifurcation of the left and right carotid arteries via an oscillator. Here, regarding the vertebral artery, since the vertebral artery supplies blood to the brainstem, it is expected to have a greater effect on autonomic nerve disorders.

[0176] In the case where the autonomic nerve disorder is relatively mild and the difference from the normal state is small, the vibration waveform of the frequency band with the difference can also be used as the vibration waveform of the normal physiological state and applied to the oscillator 31 from the beginning.

[0177] In addition, the pulse wave measurement sensor 21 needs to detect information on blood vessels several centimeters deep from the surface of the living body. Considering the absorption spectra of hemoglobin and water, which are the main light-absorbing substances present in the living body, a wavelength of 700 to 2000 nm, which is called the optical window, is used. At the same time, it is necessary to be able to observe 0.004 Hz to 200 Hz, where 0.004 Hz is the lower limit of the VLF2 frequency, and 200 Hz is obtained by doubling the upper limit of 100 Hz of the brain wave called gamma wave in consideration of the Nyquist condition.

[0178] In addition, the physiological state control system 1 can normalize the autonomic nervous system for abnormalities in physiological states such as diurnal variations and daily variations by expanding and processing waveform information in the same way as LF and VLF to a lower frequency band lower than VLF2, such as ULF1 (1.5 to 4 mHz), ULF2 (0.4 to 1.5 mHz), etc.

[0179] Furthermore, the present disclosure is not limited to the above-described embodiments and can be appropriately changed without departing from the gist. For example, combinations of the respective configurations of Embodiments 1 and 2 are also included in the scope of the technical idea of this embodiment. In addition, the following vibration wave calculation program that causes a computer to execute the calculation of the vibration wave performed by the physiological state control device 50 is also included in the scope of the technical idea.

[0180] (Supplementary Note 1)

[0181] A vibration wave calculation program that causes a computer to execute:

[0182] A step of acquiring waveform information related to a living body;

[0183] A step of filtering the acquired waveform information with at least one frequency band;

[0184] A step of performing a Hilbert transform on the waveform information of the filtered frequency band;

[0185] A step of calculating an instantaneous value including at least any one of an instantaneous logarithmic amplitude, an instantaneous frequency, and an instantaneous phase, where the instantaneous logarithmic amplitude is equivalent to the absolute value of the logarithmic amplitude term of a complex waveform equation obtained by performing a Hilbert transform on the waveform information, the instantaneous frequency is equivalent to the time differential value of the phase term of the logarithm of the complex waveform equation, and the instantaneous phase is equivalent to the phase term of the logarithm of the complex waveform equation;

[0186] A step of calculating the Gaussian distribution of the instantaneous value; and

[0187] A step of calculating the vibration wave in such a manner that the calculated Gaussian distribution becomes the predetermined Gaussian distribution.

[0188] (Supplementary Note 2)

[0189] According to the vibration wave calculation program described in Supplementary Note 1,

[0190] In the filtering step, at least any one frequency band among 0.004 - 0.015 Hz, 0.015 - 0.04 Hz, 0.04 - 0.15 Hz, 0.15 - 0.4 Hz, 0.4 - 1.5 Hz, 1.5 - 4 Hz, 4 - 15 Hz, and 15 - 40 Hz is included as the frequency band.

[0191] (Supplementary Note 3)

[0192] According to the vibration wave calculation program described in Supplementary Note 1,

[0193] In the step of calculating the vibration wave,

[0194] Calculate the vibration wave such that the center frequency of the theta wave band during eyes - open rest shifts to a lower - frequency side than the center frequency of the theta wave band during eyes - closed rest,

[0195] Calculate the vibration wave such that the center frequency of the alpha wave band during eyes - open rest shifts to a higher - frequency side than the center frequency of the alpha wave band during eyes - closed rest.

[0196] (Supplementary Note 4)

[0197] According to the vibration wave calculation program described in any one of Supplementary Notes 1 - 3,

[0198] In the step of calculating the vibration wave, as the waveform information, based on at least any one of electroencephalogram, pulse wave, and pulse interval wave, calculate the vibration wave applied to the organism.

[0199] Any type of non-transitory computer-readable medium can be used to store a program and provide it to a computer. Non-transitory computer-readable media include any type of tangible storage medium. Examples of non-transitory computer-readable media include magnetic storage media (e.g., floppy disks, magnetic tapes, hard disk drives, etc.), magneto-optical storage media (e.g., magneto-optical disks), CD-ROM (Compact Disc Read-Only Memory), CD-R (Recordable Compact Disc), CD-R / W (Rewritable Compact Disc), and semiconductor memories (e.g., mask ROM, PROM (Programmable ROM), EPROM (Erasable PROM), flash ROM, RAM (Random Access Memory), etc.). Any type of transitory computer-readable medium can be used to provide a program to a computer. Examples of transitory computer-readable media include electrical signals, optical signals, and electromagnetic waves. A transitory computer-readable medium can provide a program to a computer via a wired communication line (e.g., wires and optical fibers) or a wireless communication line.

[0200] According to the present disclosure as thus described, it is apparent that the embodiments of the present disclosure can be varied in many ways. Such variations should not be regarded as departing from the spirit and scope of the present disclosure, and all such modifications that are readily conceivable by those skilled in the art are intended to be included within the scope of the appended claims.

Claims

1. A physiological state control system, comprising: A waveform measuring device for obtaining waveform information related to a living body; a physiological state control device that calculates a vibration wave to be applied to the biological body based on the waveform information obtained by the waveform measuring device; and The driving device applies the calculated vibration wave from the outside of the biological body to at least any one of the left carotid bifurcation, the right carotid bifurcation, the vertebral artery branched from the left subclavian artery, and the vertebral artery branched from the right subclavian artery.

2. The physiological state control system according to claim 1, The waveform information includes first waveform information and second waveform information, The physiological state control device, filtering the obtained first waveform information and the obtained second waveform information using at least one frequency band, performing Hilbert transform on the first waveform information and the second waveform information of the frequency band obtained by filtering, calculating an instantaneous value of at least one of an instantaneous logarithmic amplitude corresponding to the real part of the logarithm of the complex waveform equation obtained by Hilbert transforming the first waveform information, an instantaneous logarithmic amplitude corresponding to the real part of the logarithm of the complex waveform equation obtained by Hilbert transforming the second waveform information, an instantaneous frequency corresponding to the time differential value of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the first waveform information, an instantaneous frequency corresponding to the time differential value of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the second waveform information, and an instantaneous phase difference corresponding to the difference between the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the first waveform information and the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the second waveform information, calculating a probability density distribution of the instantaneous value during a predetermined period that can be regarded as a physiologically stable state, and approximating the calculated probability density distribution to a predetermined distribution, The approximated mean value and variance value of the predetermined distribution are stored as feature quantities of the physiological state.

3. The physiological state control system according to claim 2, The first waveform information is a pulse wave, and the second waveform information is a brain wave, or the first waveform information is a pulse wave, and the second waveform information is a waveform obtained by interpolating and resampling pulse interval data calculated from the pulse wave at a desired frequency.

4. The physiological state control system according to claim 2, The physiological state control device includes at least any one of 0.004-0.015 Hz, 0.015-0.04 Hz, 0.04-0.15 Hz, 0.15-0.4 Hz, 0.4-1.5 Hz, 1.5-4 Hz, 4-8 Hz, 8-13 Hz, 13-30 Hz and 30-100 Hz as the frequency band.

5. The physiological state control system according to claim 4, In a case where the probability density distribution of the instantaneous frequency becomes a multimodal distribution, the frequency band is divided into a plurality of bands so that the probability density distribution becomes a unimodal distribution.

6. The physiological state control system according to claim 2, The physiological state control device calculates the characteristic quantity of the physiological state in a physiologically normal state and a physiologically abnormal state different from the normal state, compares the two to detect the differences, and calculates the vibration wave to reduce the differences between the two and make the characteristic quantity of the latter gradually approach the characteristic quantity of the former.

7. The physiological state control system according to claim 6, The physiological state control device detects the difference between a physiologically normal state and a physiologically abnormal state different from the normal state based on at least any one of the brain wave, the pulse wave, the phase difference between the brain wave and the pulse wave, and the phase difference between the pulse wave and the pulse interval waveform as the waveform information.

8. The physiological state control system according to claim 1, The vibration wave is input to the driving device, and the calculated vibration wave is applied to at least one of the left carotid bifurcation, the right carotid bifurcation, the left vertebral artery, and the right vertebral artery via a vibrator.

9. The physiological state control system according to claim 7, While the vibration wave is being applied, the characteristic quantity is monitored to approach the physiologically normal state.

10. The physiological state control system according to claim 8, In a frequency band below 0.4 Hz, a vibration wave obtained by AM modulation, i.e. amplitude modulation, using a vibration wave in the range of 2 Hz to 13 Hz is input into the driving device, and applied via the vibrator to at least any one of the left carotid bifurcation, the right carotid bifurcation, the left vertebral artery, and the right vertebral artery.

11. A physiological state control method, comprising: Steps for obtaining waveform information related to a biological body; A step of calculating a vibration wave to be applied to the biological body based on the acquired waveform information; and The step of applying the calculated vibration wave from the outside of the biological body to at least any one of the left carotid bifurcation, the right carotid bifurcation, the vertebral artery branched from the left subclavian artery, and the vertebral artery branched from the right subclavian artery.

12. The physiological state control method according to claim 11, The step of calculating the vibration wave comprises: The step of filtering the acquired waveform information using at least one frequency band; A step of performing Hilbert transform on the waveform information of the frequency band obtained by filtering; The step of calculating the instantaneous value of at least one of an instantaneous logarithmic amplitude corresponding to the real part of the logarithm of the complex waveform equation obtained by Hilbert transforming the waveform information, and an instantaneous frequency corresponding to the time differential value of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the waveform information; The step of calculating the probability density distribution of the instantaneous values ​​during a predetermined period which can be regarded as a physiologically stable state; The step of approximating the calculated probability density distribution to a Gaussian distribution; and The step of calculating the vibration wave in such a manner that the probability density distribution becomes the predetermined Gaussian distribution.

13. The physiological state control method according to claim 12, The probability density distribution is a unimodal distribution including a Gaussian distribution.

14. The physiological state control method according to claim 12, The step of obtaining waveform information related to the living body obtains pulse waves and brain waves as the waveform information, The step of calculating the vibration wave comprises: The step of filtering the acquired pulse wave, brain wave and pulse interval waveform using at least one frequency band; A step of performing Hilbert transform on the pulse wave, the brain wave and the pulse interval waveform of the frequency band obtained by filtering; a step of calculating an instantaneous value, the instantaneous value including at least any one of an instantaneous logarithmic amplitude which is a real part of the logarithm of the complex waveform equation obtained by Hilbert transforming the pulse wave, the brain wave and the pulse interval waveform, an instantaneous frequency which is a time differential value of the phase which is an imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the pulse wave, the brain wave and the pulse interval waveform, an instantaneous phase difference which is a difference between the phase of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the brain wave and the phase of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the pulse wave, and an instantaneous phase difference which is a difference between the phase of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the pulse interval waveform and the phase of the imaginary part of the logarithm of the complex waveform equation obtained by Hilbert transforming the pulse wave; The step of calculating the probability density distribution of the instantaneous values ​​during a predetermined period which can be regarded as a physiologically stable state; and A step of calculating the vibration wave so that the calculated probability density distribution becomes a predetermined distribution.

15. The physiological state control method according to claim 12, The step of obtaining waveform information related to the biological body obtains the pulse wave by using a photoplethysmometer that measures the pulse wave flowing into the brain, and the wavelength of the light used is 700 to 2000 nm and the bandwidth is in the range of 0.004 Hz to 200 Hz.

16. A vibration wave calculation program product, comprising a vibration wave calculation program, the vibration wave calculation program causing a computer to execute: Steps for obtaining waveform information related to a biological body; The step of filtering the acquired waveform information using at least one frequency band; A step of performing Hilbert transform on the waveform information of the frequency band obtained by filtering; The step of calculating the instantaneous value of at least one of an instantaneous logarithmic amplitude corresponding to the absolute value of the logarithmic amplitude term of the complex waveform equation obtained by Hilbert transforming the waveform information, an instantaneous frequency corresponding to the time differential value of the logarithmic phase term of the complex waveform equation, and an instantaneous phase corresponding to the logarithmic phase term of the complex waveform equation; a step of calculating a Gaussian distribution of said instantaneous value; and The step of calculating the vibration wave so that the calculated Gaussian distribution becomes the predetermined Gaussian distribution.

17. The vibration wave calculation program product according to claim 16, In the filtering step, The frequency band includes at least any one of 0.004 to 0.015 Hz, 0.015 to 0.04 Hz, 0.04 to 0.15 Hz, 0.15 to 0.4 Hz, 0.4 to 1.5 Hz, 1.5 to 4 Hz, 4 to 15 Hz, and 15 to 40 Hz.

18. The vibration wave calculation program product according to claim 16, In the step of calculating the vibration wave, calculating the vibration wave so as to shift the center frequency of the frequency band of the theta wave when the eyes are open and at rest to a frequency lower than the center frequency of the frequency band of the theta wave when the eyes are closed and at rest, The vibration wave is calculated so as to shift the center frequency of the frequency band of the α wave when the eyes are open and at rest to a frequency higher than the center frequency of the frequency band of the α wave when the eyes are closed and at rest.

19. The vibration wave calculation program product according to claim 16, In the step of calculating the vibration wave, the vibration wave applied to the living body is calculated based on at least any one of a brain wave, a pulse wave, and a pulse interval wave as the waveform information.

Images