US20040230252A1 - Method and apparatus for affecting the autonomic nervous system - Google Patents

Method and apparatus for affecting the autonomic nervous system Download PDF

Info

Publication number
US20040230252A1
US20040230252A1 US10/778,306 US77830604A US2004230252A1 US 20040230252 A1 US20040230252 A1 US 20040230252A1 US 77830604 A US77830604 A US 77830604A US 2004230252 A1 US2004230252 A1 US 2004230252A1
Authority
US
United States
Prior art keywords
stimulus
activity
stimuli
nervous system
human subject
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US10/778,306
Inventor
Saul Kullok
Yossi Kullok
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AUTO BALANCE Ltd
Original Assignee
AUTO BALANCE Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by AUTO BALANCE Ltd filed Critical AUTO BALANCE Ltd
Priority to US10/778,306 priority Critical patent/US20040230252A1/en
Publication of US20040230252A1 publication Critical patent/US20040230252A1/en
Assigned to AUTO BALANCE LTD. reassignment AUTO BALANCE LTD. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: KULLOK, JOSE R., KULLOK, SAUL
Priority to US12/169,914 priority patent/US8442632B2/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • A61M2021/0005Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
    • A61M2021/0016Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus by the smell sense
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • A61M2021/0005Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
    • A61M2021/0022Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus by the tactile sense, e.g. vibrations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • A61M2021/0005Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
    • A61M2021/0027Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus by the hearing sense
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • A61M2021/0005Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
    • A61M2021/0044Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus by the sight sense
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M21/00Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis
    • A61M2021/0005Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus
    • A61M2021/0066Other devices or methods to cause a change in the state of consciousness; Devices for producing or ending sleep by mechanical, optical, or acoustical means, e.g. for hypnosis by the use of a particular sense, or stimulus with heating or cooling
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/06Radiation therapy using light
    • A61N2005/0658Radiation therapy using light characterised by the wavelength of light used
    • A61N2005/0662Visible light
    • A61N2005/0663Coloured light

Definitions

  • the present invention relates generally to a method and apparatus for affecting a subject's health or condition by using information regarding the sympathetic and/or parasympathetic branches of the subject's autonomic nervous system to apply and/or modulate stimuli to the subject.
  • the present invention also relates to a method and apparatus for affecting the autonomic nervous system, wherein stimuli is applied in coordination with cyclical activities of the subject's body.
  • ANS autonomic nervous system
  • pathology The autonomic nervous system (ANS) and its role in health and pathology is a field of medicine that has been explored and written about at great length.
  • ANS autonomic nervous system
  • prior art methods and devices that use the concept of applying sensory stimuli to a subject's body to affect the subject's health or condition.
  • the prior art includes methods and devices for applying stimuli to a patient's body.
  • U.S. Pat. No. 5,577,990 to Widjaja et al. discloses a device that directs light and sound toward a patient, apparently eliciting a relaxation response from the patient. It is also known in the art to combine stimuli with feedback from the body.
  • U.S. Pat. No. 5,562,719 to Lopez-Claros discloses a method and apparatus for treating disorders such as Seasonal Affective Disorder by preferentially directing light therapy to the non-dominant hemisphere of the brain.
  • 4,289,121 to Kupriyanovich discloses a method and device for controlling the functional state of the central nervous system using audio and light signals applied according to the body's biorhythms that correspond to a stable state of the central nervous system.
  • This reference describes modulating frequencies depending on the electroencephalogram (EEG), electrocardiogram (ECG), or the measured respiration rate of the patient, wherein the amplitude or rhythmic signals correspond to the volume of sound and brightness of light.
  • EEG electroencephalogram
  • ECG electrocardiogram
  • Kupriyanovich also suggests that the device include a feedback system to automatically vary the illumination in response to certain changes in the patient's vital signs, including the patient's pulse rate, temperature, and respiratory rate.
  • a patient's vital signs do not provide complete information about the autonomic nervous system. The same vital sign reading for one patient may represent different levels of sympathetic and parasympathetic activity. Therefore, the complexity or the various dimensions of the autonomic nervous system activity is not fully reflected in a patient's vital sign measurements or in the patient
  • the prior art also includes stimulation treatments and biofeedback treatments that involve the patient's cognitive awareness and involvement in the treatment. These treatments are also known to include adjustment of the stimuli through trial and error.
  • the patient chooses the initial light and audio signals.
  • a patient's subjective feelings do not accurately reflect the complex interactions of the patient's autonomic nervous system.
  • the nervous system comprises the central nervous system and the peripheral nervous system.
  • the central nervous system comprises the brain and spinal cord
  • the peripheral nervous system comprises a network of nerves that connects the brain and spinal cord to the rest of the body.
  • the brain which is the site of cognitive awareness and the control center for the rest of the body, comprises the cerebrum, the brain stem, and the cerebellum.
  • the brain coordinates the ability to move, touch, taste, smell, hear, and see.
  • the cerebrum regulates a variety of voluntary activities of the body, including speech, thought, planning, and initiating communication or action.
  • a variety of critical body functions are automatically regulated by the brain stem. These functions include adjusting posture, regulating breathing, swallowing, and heart rate, controlling the rate at which the body burns food, and increasing alertness when needed.
  • the autonomic nervous system is a part of the peripheral nervous system and comprises the nerves that communicate between the brain stem and the body's internal organs.
  • the autonomic nervous system comprises the sympathetic and parasympathetic branches or systems, and it functions below the conscious level through complex interactions between its two branches to respond quickly and continuously to perturbations that threaten the stability of the body's internal environment.
  • Responses to sympathetic and parasympathetic stimulation are frequently antagonistic. For example, they have opposing or antagonistic effects on heart rate. Whereas stimulation of the sympathetic branch increases heart rate, stimulation of the parasympathetic branch decreases heart rate. In addition, the body's response to activity in one branch depends on the level of activity in the other branch.
  • a useful, albeit simplistic, analogy for the parasympathetic and sympathetic branches is that the sympathetic branch functions as the body's gas pedal and the parasympathetic branch functions as the body's brakes.
  • Sympathetic and parasympathetic activity make up a complex, dynamic system that is continuously adjusting to changing conditions in the body and in the outside environment.
  • the autonomic nervous system strives to optimize activity in each branch and to balance the two branches at every passing moment, depending on both internal and external conditions.
  • SA sinoatrial
  • the heart rate i.e., the time interval between heartbeats, is determined by how long it takes the membranes of these pacemaker cells to spontaneously depolarize to the threshold level.
  • the heart beats at a spontaneous or intrinsic rate, which is approximately 100 beats per minute, in the absence of outside influences. Outside influences are required to increase or decrease the heart rate from its intrinsic rate.
  • the two most important outside influences on heart rate come from the autonomic nervous system. Fibers from both the sympathetic branch and parasympathetic branch of the autonomic nervous system terminate on cells in the SA node and both can modify the intrinsic heart rate. Activating the cardiac sympathetic nerves increases cardiac sympathetic tone, thereby increasing heart rate. Increasing cardiac parasympathetic tone, on the other hand, slows the heart rate. Both sympathetic and parasympathetic nerves influence heart rate by altering the course of spontaneous depolarization of the resting potential in SA pacemaker cells.
  • Heart rate variability is the amount of heart rate fluctuation around a mean heart rate. Such fluctuations reflect oscillations in sympathetic-parasympathetic balance associated with a variety of factors, including respiration, baroreceptor reflexes, vasomotor control, and thermoregulatory processes. The main periodic fluctuations found are respiratory sinus arrhythmia and baroreflex-related and thermoregulation-related heart rate variability.
  • Heart rate variability is demonstrated by every normal person's heart, regardless of that person's state of health and regardless of the presence of stress or disturbances. Even a sleeping person displays heart rate variability. Each person has a measurable baseline heart rate variability even in the absence of external stressors, such as traffic, screaming babies, and looming deadlines.
  • Heart rate variability can be used as a mirror of the cardiorespiratory control system, and it is a valuable tool to investigate the sympathetic and parasympathetic function of the autonomic nervous system.
  • Heart rate variability provides information about sympathetic-parasympathetic interplay and balance, which includes other valuable information about the nervous system, including, for example, the risk for sudden cardiac death in patients after myocardial infarction.
  • Heart rate variability measurements are easy to perform, are noninvasive, and are easily and accurately reproducible. In addition, heart rate variability has been found to be largely unaffected by placebos.
  • Heart rate variability can be influenced by physiologic and maturational factors. Maturation of the autonomic nervous system results in an increase in heart rate variability with gestational age and during early post-natal life. Heart rate variability decreases with age, and this decline begins in childhood. In addition, heart rate variability is influenced by provocation and physical disorders.
  • the present invention relates generally to a method and apparatus for affecting a subject's health or condition by using information regarding the sympathetic and/or parasympathetic branch of the autonomic nervous system to modulate and/or apply stimuli to the patient.
  • This invention may be used to treat patients with various pathologies as well as to treat healthy subjects in order to improve or refine the functioning of their autonomic nervous system.
  • One aspect of the present invention involves stimulation of the autonomic nervous system based on separate analysis of the sympathetic branch and/or the parasympathetic branch of the autonomic nervous system.
  • this invention does not require any physical or mental efforts on the part of the patient, nor does it involve the patient's cognitive awareness or involvement.
  • the present invention may even be used to treat a sleeping patient.
  • the autonomic nervous system is monitored and the stimuli is modulated by continuous assessment of heart rate variability.
  • Another aspect of the present invention is to provide a real-time feedback loop or system, wherein information about the autonomic nervous system is continuously monitored and then fed back or conveyed to the patient through afferent neural pathways, i.e., nerves that conduct impulses from the periphery of the body inward to the spinal cord.
  • this feedback loop is carried out by application of sensory stimuli.
  • Another aspect of the present invention is to coordinate the stimuli with specific phases of cyclical activities of the body, such as systole or diastole in the cardiac cycle, and inspiration or expiration in the breathing cycle.
  • FIG. 1 is a block diagram, generally illustrating an embodiment of an apparatus of the present invention
  • FIG. 2 is an electrocardiogram (ECG) of a typical series of heartbeats, each showing the QRS complex, from which a fiducial point is identified in order to measure heart rate and heart rate variability;
  • ECG electrocardiogram
  • FIG. 3 is a schematic of an amplifier, which may be used in the apparatus of the present invention.
  • FIG. 4 is a block diagram, generally illustrating an embodiment of the method of the present invention.
  • FIG. 5 illustrates a typical power spectrum
  • FIG. 6 is a flowchart depicting the function of the signal processor in one embodiment of the present invention.
  • FIG. 7 is a flowchart depicting the function of the stimuli modulation software in one embodiment of the present invention.
  • FIG. 8 is a chart summarizing the stimuli parameters in one embodiment of the present invention as shown in FIG. 4.
  • the present invention is generally directed toward affecting a patient's autonomic nervous system, whether to treat a condition or simply to enhance performance of the autonomic nervous system.
  • the invention involves applying one or more default sensory stimuli, such as a visual stimulus and/or an audio stimulus to the patient, obtaining separate information about the parasympathetic and/or sympathetic branches of the autonomic nervous system by monitoring the patient, which information is preferably used to continuously alter the default stimuli according to the information obtained.
  • information about the autonomic nervous system is obtained by transforming the patient's ECG electrical signal into a measurement of heart rate variability parameters as a function of frequency.
  • FIG. 1 An apparatus according to one embodiment of the present invention is depicted schematically in FIG. 1.
  • the apparatus generally comprises patient monitor 100 , stimuli modulating unit 200 , and stimuli applicator 300 .
  • Patient or subject monitor 100 is utilized to monitor a condition of the patient or subject that may be analyzed to separately assess sympathetic and/or parasympathetic activity of the patient's autonomic nervous system.
  • the monitored patient condition or parameter may be any patient activity, including physiological, cognitive, and behavioral activity of the patient or subject.
  • patient monitor 100 comprises commercially available ECG machine.
  • Other suitable patient monitors include digital cameras to quantify the amount of the body's agitation, and a skin resistance galvanometer.
  • stimuli modulating unit 200 Signals representative of the patient condition monitored are transmitted to stimuli modulating unit 200 by patient monitor 100 .
  • the apparatus is controlled by stimuli modulating unit 200 , which includes central processing unit (CPU) 210 , memory 240 , in which data may be stored, signal processor 250 , and stimuli modulation software 260 , which contains the expert system comprising the algorithm to control the stimulus or stimuli.
  • CPU 210 is coupled with monitor 230 and/or keyboard 220 to provide an interface with the operator of the apparatus.
  • these interface elements are optional and may comprise alternative elements known in the art.
  • these components i.e., CPU 210 , keyboard 220 , monitor 230 , memory 240 , signal processor 250
  • these components are well known in the art. In fact, these components are included in most existing computers.
  • the apparatus includes an electrocardiograph connected to a standard Intel based personal computer (PC) or workstation with a color monitor, a soundboard, and an audio headset.
  • the electrocardiograph is used to measure electrical currents associated with heart muscle activity, from which the subject's HRV may be analyzed and the subject's sympathetic and/or parasympathetic activity assessed.
  • the algorithm is executed by the PC, the visual output is displayed on the color monitor, and the audio output is transmitted by the headset.
  • step 510 signal processor 250 amplifies the signal received from patient monitor 100 .
  • step 520 the analog signal is converted to a digital signal. It should be noted that some of the steps shown in FIG. 6 may not be required. For example, if the signals received by signal processor 250 are digital signals, then converting the signal to a digital signal, i.e., step 520 , is obviously not necessary.
  • step 530 a time domain analysis is performed on the signal. This is followed by step 540 comprising a frequency domain analysis of the signal.
  • step 540 comprising a frequency domain analysis of the signal.
  • signal processor 250 may be incorporated in software as part of stimuli modulation software 260 or separately in memory 240 , as is known in the art.
  • signal processor 250 produces output parameters which are indicative of sympathetic and/or parasympathetic activities of patient 400 .
  • the stimulation modulation software 260 determines the appropriate stimuli or the appropriate modulation of the stimuli.
  • the algorithm by which this is accomplished is depicted in the flowchart of FIG. 7, and is described in detail herein.
  • Stimuli applicator 300 applies appropriate stimulus to the patient 400 based on information received from modulating unit 200 .
  • applicator 300 may include headphones for sound stimulus and a television or computer monitor for light stimulus.
  • Light stimuli may be applied in any suitable manner, such as by LCD or CRT monitor, incandescent, fluorescent, or neon lighting.
  • a variety of different types of stimuli also may be used, including pressure applied to the patient and thermal radiation.
  • virtual reality helmets may be used as applicator 300 , including an open or closed helmet with glasses and earphones attached.
  • the apparatus may be completely wireless, incorporating infrared technology.
  • a subject provided with ECG electrodes and an appropriate interface to his or her own home computer could access modulating unit 200 or both patient monitor 100 and modulating unit 200 at a remote location via the Internet or a direct dial-in connection.
  • the multimedia features of current home computers make them ideal stimuli applicators 300 , providing both sound and light stimuli.
  • patient monitor 100 and stimuli applicator 300 may be incorporated together into a home unit which would communicate with modulating unit 200 via a standard telephone connection.
  • subjects without home computers or computer skills. would simply dial in to the care provider for treatment without having to make regular visits to a clinic or other treatment center.
  • signal processor 250 and stimuli modulation software 260 may be downloaded from the Internet to home computers, such as personal computers, Web TVs, and Digital TVs.
  • the patient is monitored in order to ascertain the activity of the sympathetic and/or parasympathetic branch of the patient's autonomic nervous system.
  • this activity is assessed through analysis of the patient's heart rate variability (HRV).
  • HRV heart rate variability
  • the measurement of heart rate variability contains information about several physical parameters of autonomic nervous system activity, including heart rate. Therefore, this measurement provides a multi-dimensional picture of the patient's condition that could not otherwise be obtained from simply measuring a patient's heart rate, which is simply the average number of heartbeats in a certain time interval.
  • Information about the autonomic nervous system provided by heart rate variability as compared to heart rate is analogous to information about a moving body provided by acceleration and velocity measurements as compared to only velocity measurements. Because heart rate variability measurements provide a more complete picture, better diagnosis and treatment may be obtained based on analysis of heart rate variability as opposed to simply heart rate.
  • any method of monitoring or assessing the activity of the autonomic nervous system for example breathing, may be incorporated in the present invention as long as sympathetic and/or parasympathetic branch of the autonomic nervous system may be separately analyzed.
  • step 610 the signal processor output parameters are mathematically manipulated.
  • autonomic nervous system activity is measured through statistical and spectral analysis of a series of RR intervals from a patient's ECG.
  • a resting baseline ECG of the patient is first obtained by a standard electrocardiograph.
  • FIG. 2 Such an ECG is illustrated in FIG. 2.
  • each heartbeat includes what is referred to as a QRS complex.
  • the average QRS complex spans approximately 75 msec.
  • Errors resulting from noise effects, missing data, arrhythmic or ectopic beats will distort the HRV analysis. These errors may be eliminated or reduced through interpolation of previous or successive RR intervals in the ECG signal. Other algorithms known in the art such as linear regression or auto-correlation functions may also be used to reduce these errors.
  • the subject is instructed not to engage in certain activities for approximately four hours prior to taking the ECG, which may affect his heart rate and heart rate variability, such as taking stimulants, including chocolate.
  • the ECG is taken while the subject is at rest in a semi-supine position. While the heart rate is being monitored, the subject is instructed to keep still and to avoid physical and mental activities such as chewing gum making mathematical calculations, or memorizing information.
  • the ECG signal which is in the millivolt range, is then amplified into the volt range (see step 510 in FIG. 6).
  • Several amplification methods are known in the art. Typically, amplifiers are manufactured to order, and special hardware may be designed for this particular application.
  • a type of ECG amplifier schematic that is appropriate for the present invention is shown in FIG. 3. This circuit diagram employs known symbols for circuit elements.
  • This ECG amplifier has two input bootstrapped buffers as an input stage and a differential amplifier as an output stage. The bootstrapped buffer makes the input impedance very high both at d.c. and at 50 Hz. The total gain for this amplifier is 180.
  • the analog signal is then sampled and transformed into a digital signal by a standard digital signal processing card preferably having 12-bit resolution (see step 520 in FIG. 6).
  • the electrocardiograph itself may be provided with a digital output.
  • it is sufficient to sample the ECG signal at a rate of 1000 times per second. Obviously, however, sampling at a higher rate up until the Nyquist rate will improve the accuracy of the digital representation, whereas sampling at a lower rate reduces the signal to noise ratio. In any case, a sampling rate lower than 500 times per second is not recommended.
  • the digital signal is then processed using the pattern recognition method described herein, or any other suitable mathematical technique known in the art (see step 530 in FIG. 6).
  • the pattern recognition method is preferred because it renders accurate readings even for erratic signal patterns, such as may be obtained, for example, from hyperactive kids. Through more sophisticated pattern recognition methods, increasingly stable and noise independent fiducial points may be obtained based on fewer data points.
  • RR values resulting from sinoatrial (SA) node depolarization are referred to as normal to normal (NN) intervals, and in the present invention, only the NN values are considered.
  • SA sinoatrial
  • NN normal to normal
  • the method of transforming an RR series into an NN series is well known in the art.
  • This process is repeated continuously until the end of the treatment over successive segments or time windows, wherein the starting points of each time window are separated by a predetermined interval.
  • this predetermined interval is determined by M heartbeats.
  • the preferred values for M are 4, 10, and 23, and 4 is the default value for M.
  • the subject's ECG is obtained both to measure the subject's baseline HRV and to monitor the subject during treatment.
  • the type of stimuli to be used and the parameters of the stimuli to be modulated is initially determined based on evaluation of the patient's baseline HRV.
  • This baseline HRV may be obtained from a single ECG segment of at least 1500 RR intervals, which typically will require no more than 24 minutes of recording.
  • the subject's ECG is continuously obtained until the end of the treatment, wherein several successive segments of heartbeats are identified, and wherein each segment comprises N heartbeats, as explained above. Stimuli is applied and modulated based on information about the autonomic nervous system obtained from these segments.
  • Heart rate variability may be determined in at least the following two ways: (1) by calculation of indices based on statistical operations on NN intervals (time domain analysis) or (2) by spectral (frequency domain) analysis of an array of NN intervals.
  • time domain parameters of HRV may be used in the present method:
  • SDNN standard deviations of the NN intervals
  • these parameters and other time domain parameters of ANS activity can provide information to apply or modulate sensory stimuli to the patient.
  • Frequency domain analysis or power spectral density (PSD) analysis is a preferred method for measuring heart rate variability according to the present invention (see step 540 in FIG. 6).
  • PSD analysis provides a measurement of power as a function of frequency.
  • PSD analysis enables one to evaluate the contribution to variance of frequency-specific oscillations. Therefore, one can measure not only the amount of variability, but also its distribution in relation to oscillation frequency (i.e., number of heart rate or NN value fluctuations per second).
  • the point of PSD analysis is to pass from a continuous time function, i.e., the ECG tracing or curve, to a discrete-time function, which is representative of the instantaneous heart rate.
  • Other methods of performing frequency domain analysis includes coarse-graining spectral analysis, useful in extracting the harmonic components from a broad band noise spectrum and wavelet analysis, capable of processing data at different scales or resolutions, a method well suited for approximating data containing sharp discontinuities.
  • Fourier analysis is a simple, widely used technique that involves decomposing the series of sequential RR intervals into a sum of sinusoidal functions of different amplitudes and frequencies.
  • the resulting power spectrum can be graphed as the magnitude of variability as a function of frequency.
  • the power spectrum therefore reflects the amplitude of the heart rate fluctuations present at different oscillation frequencies.
  • Fourier analysis may be performed on a short ECG recordings of 0.5 minute to several minutes to much longer recordings.
  • the autoregressive model is a linear prediction formula that attempts to predict an output y[n] of a system based on previous outputs (y[n ⁇ 1], y[n ⁇ 2], etc.) and inputs (x[n], x[n ⁇ 1], x[n ⁇ 2], etc.). Deriving the linear prediction formula involves determining the coefficients a1, a2, and b0, b1, b2, etc. in the following equation:
  • y[n ](estimated) a 1 *y[n ⁇ 1]+ a 2 *y[n ⁇ 2] . . . + b 0 *x[n]+b 1 *x[n ⁇ 1]+. . .
  • the system is determined by calculating a set of coefficients that yield an accurate prediction y[n].
  • the model will differ depending on the order of the model (n) chosen, higher orders corresponding to larger n values. Higher order models will produce higher resolution and more noise, whereas lower order models will produce less noise, or a smoother result, comprising less information.
  • geometrical methods can also provide information about HRV.
  • the most common geometrical methods include: a) sample density distribution of NN interval duration; b) sample density distribution of differences between adjacent NN intervals; and c) Lorenz plot of NN or RR intervals (i.e., XY scattergram).
  • These analyses yield shapes or patterns (e.g., the Lorenz plot characteristically yield linear, triangular, or elliptical shapes), and the geometrical and/or graphical properties of these resulting patterns provide information about various aspects of HRV.
  • any method that provides separate information about the sympathetic and/or parasympathetic branches of the ANS may be incorporated in the present invention.
  • the method of the present invention generally focuses on short-term ECG segments. For these short term segments, the following spectral bands may be identified: a) a very low frequency (VLF) band from 0.003 Hz to 0.04 Hz; b) a low frequency (LF) band from 0.04 Hz to 0.15 Hz; and c) a high frequency (HF) band from 0.15 Hz to 0.45 Hz in adults and to 0.5 Hz in children. Because the VLF band may contain many non-periodic components and may be significantly affected by certain algorithms used in the HRV analysis, it is not considered in a preferred embodiment of the present invention.
  • VLF very low frequency
  • LF low frequency
  • HF high frequency
  • PSD analysis is used to measure heart rate variability.
  • the 16th order transformation of the AR model is used.
  • An example of a power spectrum obtained by such an analysis is shown in FIG. 5, depicting a graph of variability (i.e., power) as a function of frequency (Hz).
  • Variability measurements below 0.15 Hz on the variability curve is generally a reflection of sympathetic activity
  • variability measurements above 0.15 Hz is generally a reflection of parasympathetic activity.
  • PSD analysis Separate evaluation of the parasympathetic and sympathetic branches of the ANS is possible through PSD analysis.
  • the power spectrum may also be divided into several other frequency bands or regions, in addition to or other than these two classical frequency regions. For example, the power spectrum above 0.04 Hz may be divided into three regions, wherein the lower region is below approximately 0.08 Hz, the middle region is located between approximately 0.08 Hz and 0.15 Hz, and the high region is above 0.15 Hz.
  • Characterizing each frequency region are two parameters of particular interest: the area under the curve of the power spectrum (S), which represents variance and is measured in msec 2 ; and the height of the curve at its peak (A), representing the maximum value of variability within the frequency region and is measured in msec 2 /Hz.
  • the methods of the present invention also includes other parameters which are derived from these parameters, including: (a) absolute value in the change in S over consecutive time windows (
  • the stimuli applied to the patient or subject may comprise any stimuli, including the sensory stimuli of taste, smell, and touch.
  • a background stimuli is first applied to the patient.
  • the background stimuli includes a visual background stimulus and/or an audio background stimulus.
  • the stimuli is modulated according to the activity of the sympathetic and/or parasympathetic branches of the subject's autonomic nervous system (see step 620 ).
  • the background stimuli may have default values for certain parameters, depending on the type of stimuli. For light and sound stimuli, frequency default values may be given. Background values for color brightness and sound loudness will preferably be adjusted to a minimum, and according to each individual subject, and his level of perception.
  • the visible color spectrum frequency range is taken to be approximately 380 nm to approximately 770 nm.
  • This visible frequency range is divided into a low frequency color spectrum of approximately 380 nm (red) to approximately 530 nm (green) and a high frequency color spectrum of approximately 530 nm to approximately 770 nm (blue).
  • the visual background stimulus generally consists of one of two colors. One default color has a hue within the high frequency color spectrum (i.e., approximately 490 nm), and the other default color has a hue within the low frequency color spectrum (i.e., approximately 600 nm).
  • the audible frequency range is taken to be approximately 25 Hz to approximately 1320 Hz.
  • This audible frequency range is divided into a low frequency range from approximately 25 Hz to approximately 440 Hz, and a high frequency range from approximately 440 Hz to approximately 1320 Hz.
  • One default sound has a frequency of approximately 880 Hz, which falls within the high frequency range, and the other default sound has a frequency of approximately 234 Hz, which falls within the low frequency range.
  • the audio background stimulus comprises any type of sound, such as pink noise or violin sounds and has a given loudness.
  • the background values for the sound stimulus is total silence (i.e., frequencies below 16 Hz) and total darkness.
  • the stimuli is applied during a predefined time interval, which lasts at least 60 msec, referred to as an “action time” (T A ).
  • T A time interval
  • the stimuli consists of light
  • T A should not be shorter than 60 msec to avoid fundamental frequencies above 16 Hz because these may trigger epileptic seizures.
  • Other types of stimuli may have similar limitations.
  • the action time T A repeats itself after an interval designated as time off or T OFF. during which the stimulus will return to its background values.
  • the stimuli modulation follows a cycle time equal to T A +T OFF .
  • T W waiting time
  • T W is the interval between the R N fiducial point and the beginning of T A .
  • T A may end before the occurrence of the next R N+1 fiducial point, resulting in an additional “gap time” (T G ) between the end of an action time and the beginning of the next waiting time, starting at the next fiducial point.
  • T G gap time
  • the stimuli are altered away from their default values (see step 620 in FIG. 7), according to two parameters referred to as the Stimulus Control Index (SCI) and the Directional Index (D).
  • the Stimulus Control Index (SCI) determines the degree of change from the default for any stimulus parameter, and the Direction Index determines the direction of change from the default.
  • there are four SCIs including an SCI for color hue ( ⁇ hu), an SCI for color brightness ( ⁇ br), an SCI for sound volume ( ⁇ lo), and an SCI for sound frequency ( ⁇ sf).
  • the SCI is either the absolute fractional change of S or the absolute fractional change of A, and will be a number between 0 and 1.
  • the Directional Index is either D+ or D ⁇ , depending on the increase or decrease of S or A.
  • D(s) is D+ if S increases
  • D(s) is D ⁇ if S decreases
  • D(a) is D+ if A increases
  • D(a) is D ⁇ if A decreases.
  • are first calculated.
  • the absolute fractional change of S is calculated, and this value determines the degree of the change from the default value of color hue and sound loudness.
  • the SCI for color hue ( ⁇ hu) and sound loudness ( ⁇ lo) is
  • the direction of this change is determined by the Directional Index for S.
  • the Directional Index is D(s)+ if S increases from one time window to the next, and the Directional Index is D(s) ⁇ if S decreases from one time window to the next. If the light stimulus default value is between green and blue, D(s)+ indicates a change in color toward the blue region of the color spectrum, and D(s) ⁇ indicates a change in color toward the green region of the color spectrum. In addition, D(s)+ indicates an louder audio signal, and D(s) ⁇ indicates a quieter audio signal.
  • the absolute fractional change of A is also calculated, and, as summarized in FIG. 8, this value determines the value for the change from the default value of color brightness and sound frequency.
  • the SCI for color brightness ( ⁇ br) and sound frequency ( ⁇ sf) is
  • the Directional Index is D(a)+ if A increases from one time window to the next, and the Directional Index is D(a) ⁇ if A decreases from one time window to the next.
  • D(a)+ indicates an increase in brightness
  • D(a) ⁇ indicates a decrease in brightness.
  • D(a)+ indicates an increase in sound frequency
  • D(a) ⁇ indicates a decrease in sound frequency.
  • one default audio stimuli is a sound having a frequency of 234 Hz.
  • the SCI for sound frequency, or ⁇ sf is used to determine the change in frequency required. This is determined by multiplying ⁇ sf with the appropriate sound frequency scale range.
  • the appropriate sound frequency scale range represents the maximum possible range in which the sound frequency may shift.
  • the change in frequency would be a negative shift of 0.5 (209 Hz) or approximately 104 Hz.
  • the sound frequency will shift down by 104 Hz to a sound frequency of 130 Hz, which is applied during T A .
  • the SCI for color hue ( ⁇ hu) and sound loudness ( ⁇ lo) are determined by the absolute fractional change in S
  • the SCI for color brightness ( ⁇ br) and sound frequency ( ⁇ sf) are determined by the absolute fractional change in A.
  • the four SCIs may also be determined by any combination of the absolute fractional changes in S and A.
  • the Directional Index for any of the stimuli parameters may also be determined by any combination of changes in S and A. Appropriate scale ranges must be determined for each stimuli parameter to be modulated, and these scale ranges are stored within stimuli modulation software 260 .
  • stimuli modulation is determined by further multiplying the SCI with a feedback index value (FIV) between 0 and 1 (0 ⁇ FIV ⁇ 1) depending on the effect intended.
  • FIV feedback index value
  • SCI*FIV for each D+ that appears.
  • An FIV of 1 indicates equal shifts in both directions.
  • FIV values will depend on the effect desired on sympathetic or parasympathetic activity. Initial assessment of a patient, based on baseline HRV measurement, together with accumulated experience and empirical data, will assist in determining FIV values at the expert system in stimuli modulation software 260 of the apparatus. Other alternative means of introducing a negative or positive feedback component in the applied and/or modulated stimuli may be incorporated in the present invention.
  • the stimuli is applied in coordination with cyclical activity of the body, such as with systole or diastole of the cardiac cycle and/or with inspiration or expiration of the breathing cycle.
  • an alternative method of modulating the stimuli comprises establishing a veto for predefined sections or phases of a selected biological cycle. Halting or withdrawing the stimulus during predefined phases of a biological cycle optimizes the chances of selectively influencing the sympathetic and/or the parasympathetic branches of the ANS.
  • Stimuli may therefore be applied in coordination with systole to stimulate the sympathetic branch, while stimuli may be applied in coordination with diastole to stimulate the parasympathetic branch.
  • T S The time period between the fiducial point, calculated by the method shown above, and the end of the mechanical systole, is the systolic interval or T S .
  • This systolic interval can be calculated to a precision of a few milliseconds by the following equation:
  • T W is the waiting time in msec, between the fiducial point and the 10 beginning of T A . Any combinations for T W and T A values are possible.
  • the value of Z may differ for different stimuli acting on the same biological cycle.
  • the value of Z is at least 40 msec.
  • stimuli may be applied in coordination with inspiration to stimulate the sympathetic branch, while stimuli is applied in coordination with expiration to stimulate the parasympathetic branch.
  • Stimuli coordination with the breathing cycle may require monitoring the lung movements in order to determine the beginning and the end of the inspiration period, as well as the total duration of the respiratory cycle.
  • the breathing cycle T B is four times longer than the cardiac cycle and generally lasts more than 3000 msec.
  • the inspiration period generally lasts around 1000 msec. Therefore, even if the breathing monitoring device may not be highly accurate, synchronizing stimuli application and/or modulation lasting around 200 or 300 msc, for example, with relatively long inspiration or expiration periods, will not present technical difficulties.
  • Devices for monitoring the breathing cycle include chest belts with, for example, piezo-electric sensors.
  • Another type of monitor includes finger photopletismographic devices which are sensitive to venous pressure oscillations coupled with the mechanical wave generated by the lung's periodical volume changes during the breathing cycle. After eliminating noise by proper filtering, the first derivative of the breathing wave form curve will identify the beginning and end of the inspiration slope. Once the respiratory fiducial point at the beginning of the inspiration phase is obtained, average values for the breathing cycle (mT B ) together with its two main phases, can then be calculated. Stimuli coordination with breathing can be achieved with the following algorithm, with values given in msec. To coordinate T A with the inspiration period T 1 the following equations are used:
  • T W is the time interval between the respiratory fiducial point and the beginning of the action time T A.
  • mT B is the average T B value of at least 5 breathing cycles, in msec.
  • the value of Z may change according to the case and in the preferred method Z is at least 100 msec.
  • One aspect of the invention relates to the timing of application of the stimuli.
  • the stimuli may be applied at certain times that are significant to sympathetic or parasympathetic activities.
  • the sympathetic branch of the autonomic nervous system is more active during the day, and that the parasympathetic branch of the autonomic nervous system is more active at night. Therefore, one aspect of the invention is to apply sensory stimuli either during the day or during the night depending on which branch of the autonomic nervous system to be stimulated or the particular health situation to be addressed.
  • HRV measurements for clinical evaluation may require short time ECG recordings, including at least nine hundred RR intervals. Although these recordings may be taken at any time, they should be renormalized for a predefined hour, such as 12 am. To compare HRV changes in one patient before and after treatment, more accurate results will be obtained if the ECG recordings of that patient are taken having the same length and at the same hour in the day.
  • a subject is treated for insomnia and/or stress by stimulating or increasing parasympathetic activity of the autonomic nervous system.
  • an evaluation of the patient's baseline HRV may indicate that rather than increasing parasympathetic activity, sympathetic activity should be decreased.
  • the patient's condition may require stimulation of both parasympathetic and sympathetic activity, but wherein parasympathetic activity should be increased more than sympathetic activity.
  • the patient's ECG is first recorded before the beginning of the treatment to make a preliminary diagnosis of the patient and to determine generally what method of treatment is indicated.
  • PS power spectrum
  • the region of interest is the frequency range over 0.15 Hz, which, as stated above, is thought generally to be a reflection of parasympathetic activity.
  • the stimuli is determined based on changes in the area below the power spectrum curve in this region.
  • An increase in area reflects an increase in parasympathetic activity, and a decrease in area reflects a decrease in parasympathetic activity.
  • feedback of this information comprises modulation of light and/or sound stimuli which follows a positive feedback correlation with the measured changes in the area under the curve and in the height of the highest peak in this region.
  • the visual stimuli will comprise colors in the range between green and violet.
  • the audio stimuli will comprise sounds in the low frequency range between approximately 25 and approximately 440 Hz.
  • the background default visual stimulus is a blue-green color of 490 nm. This color is displayed on a TV or computer screen.
  • the area S under the variability curve in the region of interest is continuously monitored. If this area increases, the color of at least a part of the screen is shifted toward the violet end of the color spectrum during T A , and if this area decreases, the color is shifted toward the green end of the color spectrum during T A .
  • the SCI for color hue ( ⁇ hu) is
  • the FIV 0.8, and for every D ⁇ , the stimuli shift is determined by the product of FIV and SCI.
  • the central tendency peak or the height of the curve at its peak A in the region of interest is also continuously monitored, and as the height of the peak increases, the brightness of the color is increased, and as the height of the peak decreases, the brightness is decreased, albeit in smaller proportion than the increase in intensity.
  • the SCI for color brightness ( ⁇ br) is
  • the background auditory stimulus is a sound having a frequency of 234 Hz, which is within the low frequency range between 25 Hz and 440 Hz.
  • the frequency of the sound during T A is decreased, and as the area S decreases, the frequency of the sound is increased, albeit in smaller proportion than the decreases in pitch.
  • the loudness of the sound is decreased, and as the height of the central tendency peak above 0.15 Hz decreases, the loudness of the sound is increased, albeit in smaller proportion than the decreases in loudness. Modulation of loudness is determined in the same way as the modulation of brightness, which is described above.
  • the stimuli duration T A may be constant or may vary from approximately 0.06 seconds to approximately 0.18 seconds. The default duration is approximately 0.12 sec. As parasympathetic activity increases, T A is increased, and as parasympathetic activity decreases, T A is decreased, albeit in smaller proportion than the increases in duration. In other words, in this example there is an SCI for T A which is equal to
  • T A is not synchronized with any biological cycle, and T(off) may be any arbitrary predefined time value.
  • the stimuli is applied during a predefined phase of a physiological cycle.
  • the stimuli may be synchronized with the cardiac cycle, using the R fiducial point of the QRS complex as a reference.
  • the stimuli will be applied during the latter or diastolic phase of the cardiac cycle.
  • the delay (T W ) between the fiducial point in the QRS complex of the ECG signal and the starting point of the stimuli is calculated according to equations (1), (4), (5), and (6) explained above.
  • a subject is treated for migraines and/or poor mental concentration by assisting the autonomic nervous system in stimulating or increasing sympathetic activity. Again, alternative treatments may be possible.
  • the subject's baseline heart rate variability is initially determined by time domain and spectral analysis as described above.
  • the power spectrum region of interest is below 0.15 Hz, and the stimuli is determined based on changes in the area and in the height of the highest peak in this region.
  • the visual stimulus will comprise colors ranging from green to red, and the audio stimulus will range from approximately 440 Hz to approximately 1320 Hz.
  • the background visual stimulus is an orange color of approximately 600 nm, which is a color between green and red. This color is displayed on a TV or computer screen.
  • the area S under the variability curve in the region of interest is continuously monitored. If this area increases, the color of at least part of the screen is shifted toward the red end of the color spectrum, and if this area decreases, the color of at least part of the screen is shifted toward the green end of the color spectrum. Color shifts toward the green end of the spectrum are smaller than color shifts toward the red end.
  • the SCI for color hue ( ⁇ hu) is
  • /S, and FIV 0.7.
  • the height of the central tendency peak A in the region of interest is also continuously monitored, and as the height of the peak increases, the brightness or intensity of the color is increased, and as the height of the peak decreases, the intensity is decreased, albeit in smaller proportion than the increase in intensity.
  • SCI for color brightness ( ⁇ br) is
  • the background auditory stimulus is a sound having a frequency of 880 Hz.
  • the SCI value for sound frequency ( ⁇ sf) is
  • the stimuli duration T A may be constant or may vary from approximately 0.12 seconds to approximately 0.16 seconds. The default duration is approximately 0.14 sec.
  • T A is decreased, and as sympathetic activity decreases, T A is increased, albeit in smaller proportion than the decreases in duration.
  • the SCI for T A is
  • T(off) which represents the time period between consecutive stimuli, may be determined in at least one of the following two ways.
  • T A is not synchronized with any biological cycle, and T(off) may have any arbitrary predetermined time value.
  • the stimuli is applied during a predefined phase of a physiological cycle, such as the cardiac cycle.
  • a physiological cycle such as the cardiac cycle.
  • the stimuli will be applied during the early or systolic phase of the cardiac cycle, where the delay T W between the fiducial point in the QRS complex of the ECG signal and the starting point of the stimuli will be calculated according to equations (1) and (2) explained above.

Abstract

A method and apparatus for affecting the autonomic nervous system of a subject using stimuli based on separate analysis of the sympathetic and/or parasympathetic branches of the subject's autonomic nervous system. The present invention also relates to a method and apparatus for affecting the autonomic nervous system, wherein stimuli is applied in coordination with cyclical activities of the subjects body such as respiration or cardiac cycle.

Description

    FIELD OF INVENTION
  • The present invention relates generally to a method and apparatus for affecting a subject's health or condition by using information regarding the sympathetic and/or parasympathetic branches of the subject's autonomic nervous system to apply and/or modulate stimuli to the subject. The present invention also relates to a method and apparatus for affecting the autonomic nervous system, wherein stimuli is applied in coordination with cyclical activities of the subject's body. [0001]
  • BACKGROUND
  • The autonomic nervous system (ANS) and its role in health and pathology is a field of medicine that has been explored and written about at great length. There are also several prior art methods and devices that use the concept of applying sensory stimuli to a subject's body to affect the subject's health or condition. [0002]
  • I. Application of Stimuli to Affect the Body [0003]
  • The prior art includes methods and devices for applying stimuli to a patient's body. For example, U.S. Pat. No. 5,577,990 to Widjaja et al. discloses a device that directs light and sound toward a patient, apparently eliciting a relaxation response from the patient. It is also known in the art to combine stimuli with feedback from the body. For example, U.S. Pat. No. 5,562,719 to Lopez-Claros discloses a method and apparatus for treating disorders such as Seasonal Affective Disorder by preferentially directing light therapy to the non-dominant hemisphere of the brain. In addition, U.S. Pat. No. 4,289,121 to Kupriyanovich discloses a method and device for controlling the functional state of the central nervous system using audio and light signals applied according to the body's biorhythms that correspond to a stable state of the central nervous system. This reference describes modulating frequencies depending on the electroencephalogram (EEG), electrocardiogram (ECG), or the measured respiration rate of the patient, wherein the amplitude or rhythmic signals correspond to the volume of sound and brightness of light. Kupriyanovich also suggests that the device include a feedback system to automatically vary the illumination in response to certain changes in the patient's vital signs, including the patient's pulse rate, temperature, and respiratory rate. A patient's vital signs, however, do not provide complete information about the autonomic nervous system. The same vital sign reading for one patient may represent different levels of sympathetic and parasympathetic activity. Therefore, the complexity or the various dimensions of the autonomic nervous system activity is not fully reflected in a patient's vital sign measurements or in the patient's biorhythms. [0004]
  • In addition, in the prior art methods and devices, stimuli which generally affect the sympathetic branch differently than the parasympathetic branch, are applied without regard to the autonomic nervous system balance and activity. For example, U.S. Pat. No. 5,076,281 to Gavish discloses a biorhythm modulator, which produces music-like sound pattern signals based on a patient's biorhythmic activity. Although Gavish notes that certain activities of the body are associated with the sympathetic nervous system, the rhythm of the sound synthesized patterns are based simply on a patient's monitored respiration rate, and not on separate analyses of the sympathetic and parasympathetic branches of the autonomic nervous system. [0005]
  • The prior art also includes stimulation treatments and biofeedback treatments that involve the patient's cognitive awareness and involvement in the treatment. These treatments are also known to include adjustment of the stimuli through trial and error. In the Kupriyanovich reference described above, the patient chooses the initial light and audio signals. However, a patient's subjective feelings do not accurately reflect the complex interactions of the patient's autonomic nervous system. [0006]
  • It is known that, generally, colors ranging from green to blue or violet are calming colors, and that these colors have the effect of stimulating the parasympathetic branch. It is also known that, generally, colors ranging from red to yellow are rousing colors, and that these colors have the effect of stimulating the sympathetic branch. Hospitals tend to incorporate greens and blues in their interior color scheme in order to calm and soothe patients, whereas fast food restaurants are typically red, yellow, and orange in order to move customers in and out of the restaurants quickly. In addition, it is known that increasing the brightness or intensity of a color increases its stimulatory effect. [0007]
  • It is also known that, generally, sounds having a pitch below 500 cycles or Hz tend to have a calming effect, whereas sounds having a pitch above 500 cycles tend to have a rousing effect. In addition, it is known that the louder the sound, the greater the stimulation. [0008]
  • As is apparent from the discussion above, known methods or devices generally do not provide for treatment of a patient based on the full range of information which can be ascertained from the condition of the autonomic nervous system. For example, prior art treatment methods or devices generally do not take advantage of independent or separate analyses of the sympathetic and parasympathetic branches of the autonomic nervous system. [0009]
  • II. The Autonomic Nervous System and Heart Rate Variability [0010]
  • A. Nervous System [0011]
  • The nervous system comprises the central nervous system and the peripheral nervous system. The central nervous system comprises the brain and spinal cord, and the peripheral nervous system comprises a network of nerves that connects the brain and spinal cord to the rest of the body. [0012]
  • The brain, which is the site of cognitive awareness and the control center for the rest of the body, comprises the cerebrum, the brain stem, and the cerebellum. The brain coordinates the ability to move, touch, taste, smell, hear, and see. The cerebrum regulates a variety of voluntary activities of the body, including speech, thought, planning, and initiating communication or action. [0013]
  • A variety of critical body functions are automatically regulated by the brain stem. These functions include adjusting posture, regulating breathing, swallowing, and heart rate, controlling the rate at which the body burns food, and increasing alertness when needed. The autonomic nervous system is a part of the peripheral nervous system and comprises the nerves that communicate between the brain stem and the body's internal organs. [0014]
  • The autonomic nervous system comprises the sympathetic and parasympathetic branches or systems, and it functions below the conscious level through complex interactions between its two branches to respond quickly and continuously to perturbations that threaten the stability of the body's internal environment. [0015]
  • Responses to sympathetic and parasympathetic stimulation are frequently antagonistic. For example, they have opposing or antagonistic effects on heart rate. Whereas stimulation of the sympathetic branch increases heart rate, stimulation of the parasympathetic branch decreases heart rate. In addition, the body's response to activity in one branch depends on the level of activity in the other branch. [0016]
  • A useful, albeit simplistic, analogy for the parasympathetic and sympathetic branches is that the sympathetic branch functions as the body's gas pedal and the parasympathetic branch functions as the body's brakes. Sympathetic and parasympathetic activity make up a complex, dynamic system that is continuously adjusting to changing conditions in the body and in the outside environment. The autonomic nervous system strives to optimize activity in each branch and to balance the two branches at every passing moment, depending on both internal and external conditions. [0017]
  • B. Heart Rate [0018]
  • Normal rhythmic contractions of the heart occur because of spontaneous electrical pacemaker activity of cells in the sinoatrial (SA) node. The heart rate, i.e., the time interval between heartbeats, is determined by how long it takes the membranes of these pacemaker cells to spontaneously depolarize to the threshold level. The heart beats at a spontaneous or intrinsic rate, which is approximately 100 beats per minute, in the absence of outside influences. Outside influences are required to increase or decrease the heart rate from its intrinsic rate. [0019]
  • The two most important outside influences on heart rate come from the autonomic nervous system. Fibers from both the sympathetic branch and parasympathetic branch of the autonomic nervous system terminate on cells in the SA node and both can modify the intrinsic heart rate. Activating the cardiac sympathetic nerves increases cardiac sympathetic tone, thereby increasing heart rate. Increasing cardiac parasympathetic tone, on the other hand, slows the heart rate. Both sympathetic and parasympathetic nerves influence heart rate by altering the course of spontaneous depolarization of the resting potential in SA pacemaker cells. [0020]
  • C. Heart Rate Variability (HRV) [0021]
  • Heart rate variability is the amount of heart rate fluctuation around a mean heart rate. Such fluctuations reflect oscillations in sympathetic-parasympathetic balance associated with a variety of factors, including respiration, baroreceptor reflexes, vasomotor control, and thermoregulatory processes. The main periodic fluctuations found are respiratory sinus arrhythmia and baroreflex-related and thermoregulation-related heart rate variability. [0022]
  • Heart rate variability is demonstrated by every normal person's heart, regardless of that person's state of health and regardless of the presence of stress or disturbances. Even a sleeping person displays heart rate variability. Each person has a measurable baseline heart rate variability even in the absence of external stressors, such as traffic, screaming babies, and looming deadlines. [0023]
  • Heart rate variability can be used as a mirror of the cardiorespiratory control system, and it is a valuable tool to investigate the sympathetic and parasympathetic function of the autonomic nervous system. Heart rate variability provides information about sympathetic-parasympathetic interplay and balance, which includes other valuable information about the nervous system, including, for example, the risk for sudden cardiac death in patients after myocardial infarction. [0024]
  • Heart rate variability measurements are easy to perform, are noninvasive, and are easily and accurately reproducible. In addition, heart rate variability has been found to be largely unaffected by placebos. [0025]
  • Heart rate variability can be influenced by physiologic and maturational factors. Maturation of the autonomic nervous system results in an increase in heart rate variability with gestational age and during early post-natal life. Heart rate variability decreases with age, and this decline begins in childhood. In addition, heart rate variability is influenced by provocation and physical disorders. [0026]
  • SUMMARY OF THE INVENTION
  • The present invention relates generally to a method and apparatus for affecting a subject's health or condition by using information regarding the sympathetic and/or parasympathetic branch of the autonomic nervous system to modulate and/or apply stimuli to the patient. This invention may be used to treat patients with various pathologies as well as to treat healthy subjects in order to improve or refine the functioning of their autonomic nervous system. [0027]
  • One aspect of the present invention involves stimulation of the autonomic nervous system based on separate analysis of the sympathetic branch and/or the parasympathetic branch of the autonomic nervous system. In a preferred embodiment, this invention does not require any physical or mental efforts on the part of the patient, nor does it involve the patient's cognitive awareness or involvement. The present invention may even be used to treat a sleeping patient. In a preferred embodiment the autonomic nervous system is monitored and the stimuli is modulated by continuous assessment of heart rate variability. [0028]
  • Another aspect of the present invention is to provide a real-time feedback loop or system, wherein information about the autonomic nervous system is continuously monitored and then fed back or conveyed to the patient through afferent neural pathways, i.e., nerves that conduct impulses from the periphery of the body inward to the spinal cord. In a preferred embodiment, this feedback loop is carried out by application of sensory stimuli. [0029]
  • Another aspect of the present invention is to coordinate the stimuli with specific phases of cyclical activities of the body, such as systole or diastole in the cardiac cycle, and inspiration or expiration in the breathing cycle.[0030]
  • BRIEF DESCRIPTION OF THE FIGURES
  • Reference is next made to a brief description of the figures, which are intended to illustrate the apparatus and method according to the present invention. The figures and detailed description which follow are intended to be merely illustrative, and are not intended to limit the scope of the invention as set forth in the appended claims. [0031]
  • FIG. 1 is a block diagram, generally illustrating an embodiment of an apparatus of the present invention; [0032]
  • FIG. 2 is an electrocardiogram (ECG) of a typical series of heartbeats, each showing the QRS complex, from which a fiducial point is identified in order to measure heart rate and heart rate variability; [0033]
  • FIG. 3 is a schematic of an amplifier, which may be used in the apparatus of the present invention; [0034]
  • FIG. 4 is a block diagram, generally illustrating an embodiment of the method of the present invention; [0035]
  • FIG. 5 illustrates a typical power spectrum; [0036]
  • FIG. 6 is a flowchart depicting the function of the signal processor in one embodiment of the present invention; [0037]
  • FIG. 7 is a flowchart depicting the function of the stimuli modulation software in one embodiment of the present invention; and [0038]
  • FIG. 8 is a chart summarizing the stimuli parameters in one embodiment of the present invention as shown in FIG. 4.[0039]
  • DETAILED DESCRIPTION OF THE INVENTION
  • The present invention is generally directed toward affecting a patient's autonomic nervous system, whether to treat a condition or simply to enhance performance of the autonomic nervous system. In a preferred embodiment the invention involves applying one or more default sensory stimuli, such as a visual stimulus and/or an audio stimulus to the patient, obtaining separate information about the parasympathetic and/or sympathetic branches of the autonomic nervous system by monitoring the patient, which information is preferably used to continuously alter the default stimuli according to the information obtained. In a preferred method, information about the autonomic nervous system is obtained by transforming the patient's ECG electrical signal into a measurement of heart rate variability parameters as a function of frequency. [0040]
  • An apparatus according to one embodiment of the present invention is depicted schematically in FIG. 1. As shown in this figure, the apparatus generally comprises [0041] patient monitor 100, stimuli modulating unit 200, and stimuli applicator 300.
  • Patient or [0042] subject monitor 100 is utilized to monitor a condition of the patient or subject that may be analyzed to separately assess sympathetic and/or parasympathetic activity of the patient's autonomic nervous system. The monitored patient condition or parameter may be any patient activity, including physiological, cognitive, and behavioral activity of the patient or subject. In a preferred embodiment, patient monitor 100 comprises commercially available ECG machine. Other suitable patient monitors include digital cameras to quantify the amount of the body's agitation, and a skin resistance galvanometer.
  • Signals representative of the patient condition monitored are transmitted to [0043] stimuli modulating unit 200 by patient monitor 100. The apparatus is controlled by stimuli modulating unit 200, which includes central processing unit (CPU) 210, memory 240, in which data may be stored, signal processor 250, and stimuli modulation software 260, which contains the expert system comprising the algorithm to control the stimulus or stimuli. In a preferred embodiment, CPU 210 is coupled with monitor 230 and/or keyboard 220 to provide an interface with the operator of the apparatus. However, these interface elements are optional and may comprise alternative elements known in the art. Several of these components (i.e., CPU 210, keyboard 220, monitor 230, memory 240, signal processor 250) are well known in the art. In fact, these components are included in most existing computers.
  • In an exemplary embodiment, the apparatus includes an electrocardiograph connected to a standard Intel based personal computer (PC) or workstation with a color monitor, a soundboard, and an audio headset. The electrocardiograph is used to measure electrical currents associated with heart muscle activity, from which the subject's HRV may be analyzed and the subject's sympathetic and/or parasympathetic activity assessed. The algorithm is executed by the PC, the visual output is displayed on the color monitor, and the audio output is transmitted by the headset. [0044]
  • The functions performed by [0045] signal processor 250 are shown in the flowchart of FIG. 6. First, as indicated in step 510, signal processor 250 amplifies the signal received from patient monitor 100. Then, in step 520, the analog signal is converted to a digital signal. It should be noted that some of the steps shown in FIG. 6 may not be required. For example, if the signals received by signal processor 250 are digital signals, then converting the signal to a digital signal, i.e., step 520, is obviously not necessary. Next, in step 530, a time domain analysis is performed on the signal. This is followed by step 540 comprising a frequency domain analysis of the signal. The hardware necessary to perform these functions are well known in the art. Alternatively the functionality of signal processor 250 may be incorporated in software as part of stimuli modulation software 260 or separately in memory 240, as is known in the art. Ultimately, signal processor 250 produces output parameters which are indicative of sympathetic and/or parasympathetic activities of patient 400.
  • The [0046] stimulation modulation software 260 then determines the appropriate stimuli or the appropriate modulation of the stimuli. The algorithm by which this is accomplished is depicted in the flowchart of FIG. 7, and is described in detail herein.
  • [0047] Stimuli applicator 300 applies appropriate stimulus to the patient 400 based on information received from modulating unit 200. As schematically represented in FIG. 1, applicator 300 may include headphones for sound stimulus and a television or computer monitor for light stimulus. Light stimuli may be applied in any suitable manner, such as by LCD or CRT monitor, incandescent, fluorescent, or neon lighting. A variety of different types of stimuli also may be used, including pressure applied to the patient and thermal radiation. Also, virtual reality helmets may be used as applicator 300, including an open or closed helmet with glasses and earphones attached. The apparatus may be completely wireless, incorporating infrared technology.
  • According to further alternative embodiments of the present invention, it is not necessary that all components of the apparatus be located together or directly linked together. For example, a subject provided with ECG electrodes and an appropriate interface to his or her own home computer could access modulating [0048] unit 200 or both patient monitor 100 and modulating unit 200 at a remote location via the Internet or a direct dial-in connection. The multimedia features of current home computers make them ideal stimuli applicators 300, providing both sound and light stimuli. Alternatively, patient monitor 100 and stimuli applicator 300 may be incorporated together into a home unit which would communicate with modulating unit 200 via a standard telephone connection. In this embodiment, subjects without home computers or computer skills. would simply dial in to the care provider for treatment without having to make regular visits to a clinic or other treatment center. Regular, periodic dial-in sessions would create individual patient histories which may be automatically monitored for significant variations to provide automatic warnings to appropriate caregivers as may be required. In alternative embodiments, signal processor 250 and stimuli modulation software 260 may be downloaded from the Internet to home computers, such as personal computers, Web TVs, and Digital TVs.
  • According to the method of the present invention, the patient is monitored in order to ascertain the activity of the sympathetic and/or parasympathetic branch of the patient's autonomic nervous system. In the preferred embodiment, this activity is assessed through analysis of the patient's heart rate variability (HRV). Measuring heart rate variability, as opposed to simply heart rate, yields a wealth of information about the autonomic nervous system which leads to improved diagnosis of a variety of pathologies, including hypertension, cardiac ischemia, myocardial infarction, diabetic neuropathy and other autonomic dysfunctions. Improved diagnosis in turn leads to more specific, and subsequently more effective, therapy. [0049]
  • The measurement of heart rate variability contains information about several physical parameters of autonomic nervous system activity, including heart rate. Therefore, this measurement provides a multi-dimensional picture of the patient's condition that could not otherwise be obtained from simply measuring a patient's heart rate, which is simply the average number of heartbeats in a certain time interval. Information about the autonomic nervous system provided by heart rate variability as compared to heart rate is analogous to information about a moving body provided by acceleration and velocity measurements as compared to only velocity measurements. Because heart rate variability measurements provide a more complete picture, better diagnosis and treatment may be obtained based on analysis of heart rate variability as opposed to simply heart rate. [0050]
  • It should be noted, however, that any method of monitoring or assessing the activity of the autonomic nervous system, for example breathing, may be incorporated in the present invention as long as sympathetic and/or parasympathetic branch of the autonomic nervous system may be separately analyzed. [0051]
  • As shown in FIG. 7, in [0052] step 610, the signal processor output parameters are mathematically manipulated. According to the preferred method, autonomic nervous system activity is measured through statistical and spectral analysis of a series of RR intervals from a patient's ECG. In this analysis, a resting baseline ECG of the patient is first obtained by a standard electrocardiograph. Such an ECG is illustrated in FIG. 2. As shown in FIG. 2, each heartbeat includes what is referred to as a QRS complex. The average QRS complex spans approximately 75 msec.
  • Errors, resulting from noise effects, missing data, arrhythmic or ectopic beats will distort the HRV analysis. These errors may be eliminated or reduced through interpolation of previous or successive RR intervals in the ECG signal. Other algorithms known in the art such as linear regression or auto-correlation functions may also be used to reduce these errors. [0053]
  • To obtain an HRV measurement with minimal noise, the subject is instructed not to engage in certain activities for approximately four hours prior to taking the ECG, which may affect his heart rate and heart rate variability, such as taking stimulants, including chocolate. The ECG is taken while the subject is at rest in a semi-supine position. While the heart rate is being monitored, the subject is instructed to keep still and to avoid physical and mental activities such as chewing gum making mathematical calculations, or memorizing information. [0054]
  • The ECG signal, which is in the millivolt range, is then amplified into the volt range (see [0055] step 510 in FIG. 6). Several amplification methods are known in the art. Typically, amplifiers are manufactured to order, and special hardware may be designed for this particular application. A type of ECG amplifier schematic that is appropriate for the present invention is shown in FIG. 3. This circuit diagram employs known symbols for circuit elements. This ECG amplifier has two input bootstrapped buffers as an input stage and a differential amplifier as an output stage. The bootstrapped buffer makes the input impedance very high both at d.c. and at 50 Hz. The total gain for this amplifier is 180.
  • One of ordinary skill will be able to construct an amplifier based on this circuit diagram or an appropriate amplifier for the present invention. Alternatively, one may incorporate already existing software or any other suitable amplification method known in the art. [0056]
  • The analog signal is then sampled and transformed into a digital signal by a standard digital signal processing card preferably having 12-bit resolution (see [0057] step 520 in FIG. 6). The electrocardiograph itself may be provided with a digital output. For the present method, it is sufficient to sample the ECG signal at a rate of 1000 times per second. Obviously, however, sampling at a higher rate up until the Nyquist rate will improve the accuracy of the digital representation, whereas sampling at a lower rate reduces the signal to noise ratio. In any case, a sampling rate lower than 500 times per second is not recommended.
  • The digital signal is then processed using the pattern recognition method described herein, or any other suitable mathematical technique known in the art (see [0058] step 530 in FIG. 6). The pattern recognition method is preferred because it renders accurate readings even for erratic signal patterns, such as may be obtained, for example, from hyperactive kids. Through more sophisticated pattern recognition methods, increasingly stable and noise independent fiducial points may be obtained based on fewer data points.
  • According to standard pattern recognition methods known in the art, linear regression is used on the QR slope data and the RS slope data to obtain the most representative straight line for both the QR slope and the RS slope. The point of intersection of both lines represent the R[0059] n fiducial point for each QRS waveform n. The time lag between Rn and Rn+1 is the value RR, which is typically expressed in milliseconds, and a series of RR values is referred to as an RR series.
  • The RR values resulting from sinoatrial (SA) node depolarization are referred to as normal to normal (NN) intervals, and in the present invention, only the NN values are considered. The method of transforming an RR series into an NN series is well known in the art. [0060]
  • A series of NN values are measured over a segment or time window defined by N heartbeats, and the heart rate variability parameters of interest are calculated based on the ECG recording within the segment, where N=30k and k is an integer between 1 and 6, and 4 and 120 are the default values for k and N respectively. This process is repeated continuously until the end of the treatment over successive segments or time windows, wherein the starting points of each time window are separated by a predetermined interval. In the preferred embodiment, this predetermined interval is determined by M heartbeats. The preferred values for M are 4, 10, and 23, and 4 is the default value for M. [0061]
  • Typically, the subject's ECG is obtained both to measure the subject's baseline HRV and to monitor the subject during treatment. The type of stimuli to be used and the parameters of the stimuli to be modulated is initially determined based on evaluation of the patient's baseline HRV. This baseline HRV may be obtained from a single ECG segment of at least 1500 RR intervals, which typically will require no more than 24 minutes of recording. During treatment, on the other hand, the subject's ECG is continuously obtained until the end of the treatment, wherein several successive segments of heartbeats are identified, and wherein each segment comprises N heartbeats, as explained above. Stimuli is applied and modulated based on information about the autonomic nervous system obtained from these segments. [0062]
  • Heart rate variability may be determined in at least the following two ways: (1) by calculation of indices based on statistical operations on NN intervals (time domain analysis) or (2) by spectral (frequency domain) analysis of an array of NN intervals. [0063]
  • To assess baseline values of HRV or to calculate changes in ANS activity during treatment, the following time domain parameters of HRV may be used in the present method: [0064]
  • (1) mean NN intervals (mNN); [0065]
  • (2) mean heart rate as derived from the mean NN intervals; [0066]
  • (3) difference between the minimal and maximal NN interval; [0067]
  • (4) standard deviations of the NN intervals (SDNN), which is the root of variance; [0068]
  • (5) root mean square of successive differences in the NN intervals (rMSSD); and [0069]
  • (6) percentile of intervals in which the absolute value of the difference between them is greater than 50 msec (Pnn50). [0070]
  • In alternative embodiments of the present invention, these parameters and other time domain parameters of ANS activity can provide information to apply or modulate sensory stimuli to the patient. [0071]
  • Frequency domain analysis or power spectral density (PSD) analysis is a preferred method for measuring heart rate variability according to the present invention (see [0072] step 540 in FIG. 6). This analysis provides a measurement of power as a function of frequency. PSD analysis enables one to evaluate the contribution to variance of frequency-specific oscillations. Therefore, one can measure not only the amount of variability, but also its distribution in relation to oscillation frequency (i.e., number of heart rate or NN value fluctuations per second). The point of PSD analysis is to pass from a continuous time function, i.e., the ECG tracing or curve, to a discrete-time function, which is representative of the instantaneous heart rate.
  • There are a number of methods of performing frequency domain analysis. Two popular methods include fast Fourier transformation and autoregressive (AR) modeling. Both methods are well known in the art, and both yield similar results. Fast Fourier transformation spectra, however, are characterized by discrete peaks for the several frequency components, whereas the autoregressive method results in continuous smooth spectra of activity. Fast Fourier transformation analysis is particularly useful in determining the area under the spectrum, whereas AR modeling is particularly useful in determining central tendencies or peaks in given frequency ranges. [0073]
  • Other methods of performing frequency domain analysis includes coarse-graining spectral analysis, useful in extracting the harmonic components from a broad band noise spectrum and wavelet analysis, capable of processing data at different scales or resolutions, a method well suited for approximating data containing sharp discontinuities. [0074]
  • Fourier analysis is a simple, widely used technique that involves decomposing the series of sequential RR intervals into a sum of sinusoidal functions of different amplitudes and frequencies. The resulting power spectrum can be graphed as the magnitude of variability as a function of frequency. The power spectrum therefore reflects the amplitude of the heart rate fluctuations present at different oscillation frequencies. Fourier analysis may be performed on a short ECG recordings of 0.5 minute to several minutes to much longer recordings. [0075]
  • The autoregressive model is a linear prediction formula that attempts to predict an output y[n] of a system based on previous outputs (y[n−1], y[n−2], etc.) and inputs (x[n], x[n−1], x[n−2], etc.). Deriving the linear prediction formula involves determining the coefficients a1, a2, and b0, b1, b2, etc. in the following equation: [0076]
  • y[n](estimated)=a1*y[n−1]+a2*y[n−2] . . . +b0*x[n]+b1*x[n−1]+. . .
  • In other words, the system is determined by calculating a set of coefficients that yield an accurate prediction y[n]. The model will differ depending on the order of the model (n) chosen, higher orders corresponding to larger n values. Higher order models will produce higher resolution and more noise, whereas lower order models will produce less noise, or a smoother result, comprising less information. [0077]
  • Other than the time and frequency analyses of a series of NN intervals outlined above, geometrical methods can also provide information about HRV. The most common geometrical methods include: a) sample density distribution of NN interval duration; b) sample density distribution of differences between adjacent NN intervals; and c) Lorenz plot of NN or RR intervals (i.e., XY scattergram). These analyses yield shapes or patterns (e.g., the Lorenz plot characteristically yield linear, triangular, or elliptical shapes), and the geometrical and/or graphical properties of these resulting patterns provide information about various aspects of HRV. Again, any method that provides separate information about the sympathetic and/or parasympathetic branches of the ANS may be incorporated in the present invention. [0078]
  • Periodic variations in heart rate at different frequency ranges reflect different aspects of autonomic nervous system activity. The method of the present invention generally focuses on short-term ECG segments. For these short term segments, the following spectral bands may be identified: a) a very low frequency (VLF) band from 0.003 Hz to 0.04 Hz; b) a low frequency (LF) band from 0.04 Hz to 0.15 Hz; and c) a high frequency (HF) band from 0.15 Hz to 0.45 Hz in adults and to 0.5 Hz in children. Because the VLF band may contain many non-periodic components and may be significantly affected by certain algorithms used in the HRV analysis, it is not considered in a preferred embodiment of the present invention. [0079]
  • In a preferred embodiment of the present invention, PSD analysis is used to measure heart rate variability. Specifically, the 16th order transformation of the AR model is used. An example of a power spectrum obtained by such an analysis is shown in FIG. 5, depicting a graph of variability (i.e., power) as a function of frequency (Hz). [0080]
  • Variability measurements below 0.15 Hz on the variability curve is generally a reflection of sympathetic activity, and variability measurements above 0.15 Hz is generally a reflection of parasympathetic activity. Separate evaluation of the parasympathetic and sympathetic branches of the ANS is possible through PSD analysis. The power spectrum may also be divided into several other frequency bands or regions, in addition to or other than these two classical frequency regions. For example, the power spectrum above 0.04 Hz may be divided into three regions, wherein the lower region is below approximately 0.08 Hz, the middle region is located between approximately 0.08 Hz and 0.15 Hz, and the high region is above 0.15 Hz. [0081]
  • Characterizing each frequency region are two parameters of particular interest: the area under the curve of the power spectrum (S), which represents variance and is measured in msec[0082] 2; and the height of the curve at its peak (A), representing the maximum value of variability within the frequency region and is measured in msec2/Hz. The methods of the present invention also includes other parameters which are derived from these parameters, including: (a) absolute value in the change in S over consecutive time windows (|ΔS|); (b) the absolute value in the change in A over consecutive time windows (|ΔA|); (c) the absolute fractional change of S (|ΔS|/S); and (d) the absolute fractional change of A (|ΔA|/A).
  • According to the present invention, the stimuli applied to the patient or subject may comprise any stimuli, including the sensory stimuli of taste, smell, and touch. In the preferred method of the present method, a background stimuli is first applied to the patient. The background stimuli includes a visual background stimulus and/or an audio background stimulus. Then, as shown in FIG. 7, the stimuli is modulated according to the activity of the sympathetic and/or parasympathetic branches of the subject's autonomic nervous system (see step [0083] 620).
  • The background stimuli may have default values for certain parameters, depending on the type of stimuli. For light and sound stimuli, frequency default values may be given. Background values for color brightness and sound loudness will preferably be adjusted to a minimum, and according to each individual subject, and his level of perception. [0084]
  • In the preferred embodiment, the visible color spectrum frequency range is taken to be approximately 380 nm to approximately 770 nm. This visible frequency range is divided into a low frequency color spectrum of approximately 380 nm (red) to approximately 530 nm (green) and a high frequency color spectrum of approximately 530 nm to approximately 770 nm (blue). The visual background stimulus generally consists of one of two colors. One default color has a hue within the high frequency color spectrum (i.e., approximately 490 nm), and the other default color has a hue within the low frequency color spectrum (i.e., approximately 600 nm). [0085]
  • In the preferred embodiment, the audible frequency range is taken to be approximately 25 Hz to approximately 1320 Hz. This audible frequency range is divided into a low frequency range from approximately 25 Hz to approximately 440 Hz, and a high frequency range from approximately 440 Hz to approximately 1320 Hz. One default sound has a frequency of approximately 880 Hz, which falls within the high frequency range, and the other default sound has a frequency of approximately 234 Hz, which falls within the low frequency range. In addition, the audio background stimulus comprises any type of sound, such as pink noise or violin sounds and has a given loudness. [0086]
  • In alternative embodiments, the background values for the sound stimulus is total silence (i.e., frequencies below 16 Hz) and total darkness. [0087]
  • In the present method, the stimuli is applied during a predefined time interval, which lasts at least 60 msec, referred to as an “action time” (T[0088] A). When the stimuli consists of light, TA should not be shorter than 60 msec to avoid fundamental frequencies above 16 Hz because these may trigger epileptic seizures. Other types of stimuli may have similar limitations.
  • The action time T[0089] A repeats itself after an interval designated as time off or TOFF. during which the stimulus will return to its background values. In one embodiment, the stimuli modulation follows a cycle time equal to TA+TOFF. These time values are previously defined for each type of stimulus at the expert system in the stimuli modulation software 260.
  • When the stimuli is to be applied in synchronization with a biological cycle, there will be a time interval, referred to as “waiting time” (T[0090] W), between a fiducial point of the biological cycle and the beginning of TA. For example, when the stimuli is to be applied in synchronization with the cardiac cycle, TW is the interval between the RN fiducial point and the beginning of TA. TA may end before the occurrence of the next RN+1 fiducial point, resulting in an additional “gap time” (TG) between the end of an action time and the beginning of the next waiting time, starting at the next fiducial point. For example, if an NN interval is 783 msec, with a TW of a 200 ms; and TA of 180 msec, then TG=TNN−(TW+TA)=403 msec.
  • It should be noted that, when a single stimulus is applied and modulated in synchronization with a biological cycle, the total cycle time for that particular stimulus will be equal to T[0091] W+TA+TG, where TW and TA have been previously defined, for each type of stimulus and/or its parameters within stimuli modulation software 260.
  • During each action time T[0092] A, the stimuli are altered away from their default values (see step 620 in FIG. 7), according to two parameters referred to as the Stimulus Control Index (SCI) and the Directional Index (D). The Stimulus Control Index (SCI) determines the degree of change from the default for any stimulus parameter, and the Direction Index determines the direction of change from the default. In this example, there are four SCIs, including an SCI for color hue (θhu), an SCI for color brightness (θbr), an SCI for sound volume (θlo), and an SCI for sound frequency (θsf). The SCI is either the absolute fractional change of S or the absolute fractional change of A, and will be a number between 0 and 1. The Directional Index is either D+ or D−, depending on the increase or decrease of S or A. D(s) is D+ if S increases, D(s) is D− if S decreases, D(a) is D+ if A increases, and D(a) is D− if A decreases. These values are explained in more detail in the following example, which is explained with reference to FIG. 8.
  • In this example, the values for S, |ΔS|, A, |ΔA| are first calculated. As summarized in FIG. 8, the absolute fractional change of S is calculated, and this value determines the degree of the change from the default value of color hue and sound loudness. In other words, in this example, the SCI for color hue (θhu) and sound loudness (θlo) is |ΔS|/S. [0093]
  • The direction of this change is determined by the Directional Index for S. The Directional Index is D(s)+ if S increases from one time window to the next, and the Directional Index is D(s)− if S decreases from one time window to the next. If the light stimulus default value is between green and blue, D(s)+ indicates a change in color toward the blue region of the color spectrum, and D(s)− indicates a change in color toward the green region of the color spectrum. In addition, D(s)+ indicates an louder audio signal, and D(s)− indicates a quieter audio signal. [0094]
  • The absolute fractional change of A is also calculated, and, as summarized in FIG. 8, this value determines the value for the change from the default value of color brightness and sound frequency. In other words, in this example, the SCI for color brightness (θbr) and sound frequency (θsf) is |ΔA|/A. [0095]
  • The direction of this change is determined by the Directional Index for A. [0096]
  • The Directional Index is D(a)+ if A increases from one time window to the next, and the Directional Index is D(a)− if A decreases from one time window to the next. D(a)+ indicates an increase in brightness, and D(a)− indicates a decrease in brightness. In addition, D(a)+ indicates an increase in sound frequency, and D(a)− indicates a decrease in sound frequency. [0097]
  • Depending on these parameters, the stimuli shifts away from the default during T[0098] A. Specifically, as set forth above, one default audio stimuli is a sound having a frequency of 234 Hz. The SCI for sound frequency, or θsf, is used to determine the change in frequency required. This is determined by multiplying θsf with the appropriate sound frequency scale range. The appropriate sound frequency scale range represents the maximum possible range in which the sound frequency may shift.
  • In this example, the sound frequency is limited to the low frequency range of 25 Hz to 440 Hz. Therefore, the sound frequency may shift from 234 Hz up to a maximum of 440 Hz (a scale range of 206 Hz) or may shift down from 234 Hz to a minimum of 25 Hz (a scale range of 209 Hz). Given θsf=0.5, therefore, two results may occur. If D(a)=D+ (due to an increase in A from the previous time window), the change in frequency would be a positive shift of 0.5 (206 Hz) or 103 Hz. In other words, the sound frequency will shift up by 103 Hz to a sound frequency of 337 Hz, which is applied during T[0099] A. If, on the other hand, D(a)=D− (due to a decrease in A from the previous time window), the change in frequency would be a negative shift of 0.5 (209 Hz) or approximately 104 Hz. In other words, the sound frequency will shift down by 104 Hz to a sound frequency of 130 Hz, which is applied during TA.
  • In this example, as shown in FIG. 8, the SCI for color hue (θhu) and sound loudness (θlo) are determined by the absolute fractional change in S, and the SCI for color brightness (θbr) and sound frequency (θsf) are determined by the absolute fractional change in A. However, in alternative embodiments, the four SCIs may also be determined by any combination of the absolute fractional changes in S and A. Similarly, the Directional Index for any of the stimuli parameters may also be determined by any combination of changes in S and A. Appropriate scale ranges must be determined for each stimuli parameter to be modulated, and these scale ranges are stored within [0100] stimuli modulation software 260.
  • In alternative embodiments of the present invention, stimuli modulation is determined by further multiplying the SCI with a feedback index value (FIV) between 0 and 1 (0<FIV<1) depending on the effect intended. For example, if a positive feedback effect is intended, the stimuli modulation will be determined by SCI*FIV (e.g., FIV=0.8) for each D− that appears. In other words, stimuli shifts in the negative direction will be less than stimuli shifts in the positive direction. Similarly, if a negative feedback effect is intended, the stimuli modulation will be determined by SCI*FIV for each D+ that appears. An FIV of 1 indicates equal shifts in both directions. [0101]
  • FIV values will depend on the effect desired on sympathetic or parasympathetic activity. Initial assessment of a patient, based on baseline HRV measurement, together with accumulated experience and empirical data, will assist in determining FIV values at the expert system in [0102] stimuli modulation software 260 of the apparatus. Other alternative means of introducing a negative or positive feedback component in the applied and/or modulated stimuli may be incorporated in the present invention.
  • In an alternative method of the present invention, the stimuli is applied in coordination with cyclical activity of the body, such as with systole or diastole of the cardiac cycle and/or with inspiration or expiration of the breathing cycle. [0103]
  • It has been found that the balance between sympathetic and parasympathetic activity shifts or changes during periodic biological cycles, such as the menstrual and breathing cycles. There is also a circadian rhythm to ANS activity, wherein parasympathetic activity increases during the night and sympathetic activity increases during the day. In short biological cycles, like the heartbeat, a rise in sympathetic tone increases myocardial contractility and may prolong systolic ejection time. In addition, the QT interval in the electrocardiogram corresponding to the depolarization-repolarization cycle of the ventricles also increases if sympathetic activity increases. Sympathetic activity is related to the duration of systole and inspiration, and the diurnal period. Similarly, parasympathetic activity is related to the duration of diastole and expiration, and the nocturnal period. Due to these relationships, an alternative method of modulating the stimuli comprises establishing a veto for predefined sections or phases of a selected biological cycle. Halting or withdrawing the stimulus during predefined phases of a biological cycle optimizes the chances of selectively influencing the sympathetic and/or the parasympathetic branches of the ANS. [0104]
  • Stimuli may therefore be applied in coordination with systole to stimulate the sympathetic branch, while stimuli may be applied in coordination with diastole to stimulate the parasympathetic branch. [0105]
  • This synchronization is obtained by the following algorithm: [0106]
  • 1) The mean value of a normal cardiac cycle (mNN) is obtained in milliseconds (msec), from 30 prior successive NN intervals. [0107]
  • 2) The time period between the fiducial point, calculated by the method shown above, and the end of the mechanical systole, is the systolic interval or T[0108] S. This systolic interval can be calculated to a precision of a few milliseconds by the following equation:
  • T s=512−(123.000/mNN)  (1)
  • 3) To coordinate T[0109] A with TS, the following equations are used:
  • T W +T A =T S −Z  (2)
  • T A≧60  (3) (for light stimuli)
  • where T[0110] W is the waiting time in msec, between the fiducial point and the 10 beginning of TA. Any combinations for TW and TA values are possible.
  • 4) To coordinate T[0111] A with the cardiac interval following TS, known as diastole, the following equations are used:
  • T W =T S +Z  (4)
  • T A≦mNN−(T W+180)  (5)
  • T A≦60  (6) (for light stimuli)
  • The value of Z, given in msec, may differ for different stimuli acting on the same biological cycle. In the preferred method, the value of Z is at least 40 msec. [0112]
  • Alternatively, stimuli may be applied in coordination with inspiration to stimulate the sympathetic branch, while stimuli is applied in coordination with expiration to stimulate the parasympathetic branch. [0113]
  • Stimuli coordination with the breathing cycle may require monitoring the lung movements in order to determine the beginning and the end of the inspiration period, as well as the total duration of the respiratory cycle. On average, the breathing cycle T[0114] B, is four times longer than the cardiac cycle and generally lasts more than 3000 msec. The inspiration period generally lasts around 1000 msec. Therefore, even if the breathing monitoring device may not be highly accurate, synchronizing stimuli application and/or modulation lasting around 200 or 300 msc, for example, with relatively long inspiration or expiration periods, will not present technical difficulties. Devices for monitoring the breathing cycle include chest belts with, for example, piezo-electric sensors. Another type of monitor includes finger photopletismographic devices which are sensitive to venous pressure oscillations coupled with the mechanical wave generated by the lung's periodical volume changes during the breathing cycle. After eliminating noise by proper filtering, the first derivative of the breathing wave form curve will identify the beginning and end of the inspiration slope. Once the respiratory fiducial point at the beginning of the inspiration phase is obtained, average values for the breathing cycle (mTB) together with its two main phases, can then be calculated. Stimuli coordination with breathing can be achieved with the following algorithm, with values given in msec. To coordinate TA with the inspiration period T1 the following equations are used:
  • T W +T A =T 1 −Z  (7)
  • T A≧60  (8) (for light stimuli)
  • To coordinate T[0115] A with the expiration period the following equations are used:
  • T W =T 1 +Z  (9)
  • T A ≦mTB−(TW+Z)  (10)
  • T A≦60  (11) (for light stimuli)
  • where T[0116] W is the time interval between the respiratory fiducial point and the beginning of the action time TA., in msec and mTB is the average TB value of at least 5 breathing cycles, in msec. The value of Z may change according to the case and in the preferred method Z is at least 100 msec.
  • Simultaneous coordination with the cardiac and respiratory cycle is also possible. In order to influence sympathetic activity equations (1), (2), and (7) are used. T[0117] A will be given the same value in the above algorithms. These particular combinations could be proven to be very useful in patients suffering from asthma. Similarly, in order to influence the parasympathetic activity equations (1), (4), (5), (9), and (10) are used, giving TA the same value.
  • One aspect of the invention relates to the timing of application of the stimuli. In a further alternative method of the present invention, for example, the stimuli may be applied at certain times that are significant to sympathetic or parasympathetic activities. As stated above, it is known that the sympathetic branch of the autonomic nervous system is more active during the day, and that the parasympathetic branch of the autonomic nervous system is more active at night. Therefore, one aspect of the invention is to apply sensory stimuli either during the day or during the night depending on which branch of the autonomic nervous system to be stimulated or the particular health situation to be addressed. [0118]
  • It should also be noted that the time of day also has an effect on HRV measurement. HRV measurements for clinical evaluation may require short time ECG recordings, including at least nine hundred RR intervals. Although these recordings may be taken at any time, they should be renormalized for a predefined hour, such as 12 am. To compare HRV changes in one patient before and after treatment, more accurate results will be obtained if the ECG recordings of that patient are taken having the same length and at the same hour in the day. [0119]
  • The following examples are intended to illustrate the method according to the present invention. However, they are not intended to limit the scope of the invention as set forth in the appended claims. [0120]
  • EXAMPLE 1
  • In this example, a subject is treated for insomnia and/or stress by stimulating or increasing parasympathetic activity of the autonomic nervous system. It should be noted that an evaluation of the patient's baseline HRV may indicate that rather than increasing parasympathetic activity, sympathetic activity should be decreased. Alternatively, the patient's condition may require stimulation of both parasympathetic and sympathetic activity, but wherein parasympathetic activity should be increased more than sympathetic activity. [0121]
  • As described above, the patient's ECG is first recorded before the beginning of the treatment to make a preliminary diagnosis of the patient and to determine generally what method of treatment is indicated. During treatment, the patient's ECG is continuously obtained from which successive power spectrum (PS) variability curves are continuously derived from successive and overlapping time windows, wherein each time window spans 120 heartbeats (i.e., k=4 (default value) and N=30k=120 heartbeats) or approximately 2 minutes of the ECG recording. The beginning of each window is separated from beginning of the previous window by M=10 heartbeats or approximately 10-second intervals. [0122]
  • In this example, the region of interest is the frequency range over 0.15 Hz, which, as stated above, is thought generally to be a reflection of parasympathetic activity. The stimuli is determined based on changes in the area below the power spectrum curve in this region. An increase in area reflects an increase in parasympathetic activity, and a decrease in area reflects a decrease in parasympathetic activity. In one embodiment of the present invention, feedback of this information comprises modulation of light and/or sound stimuli which follows a positive feedback correlation with the measured changes in the area under the curve and in the height of the highest peak in this region. [0123]
  • Because a parasympathetic effect is intended, the visual stimuli will comprise colors in the range between green and violet. The audio stimuli will comprise sounds in the low frequency range between approximately 25 and approximately 440 Hz. [0124]
  • The background default visual stimulus is a blue-green color of 490 nm. This color is displayed on a TV or computer screen. The area S under the variability curve in the region of interest is continuously monitored. If this area increases, the color of at least a part of the screen is shifted toward the violet end of the color spectrum during T[0125] A, and if this area decreases, the color is shifted toward the green end of the color spectrum during TA. The SCI for color hue (θhu) is |ΔS|/S. Because a positive feedback effect is desired, color shifts toward the green end of the spectrum will be smaller than color shifts toward the violet end. In this example, the FIV=0.8, and for every D−, the stimuli shift is determined by the product of FIV and SCI.
  • The central tendency peak or the height of the curve at its peak A in the region of interest is also continuously monitored, and as the height of the peak increases, the brightness of the color is increased, and as the height of the peak decreases, the brightness is decreased, albeit in smaller proportion than the increase in intensity. The SCI for color brightness (θbr) is |ΔA|/A. Again, for every D−, the stimuli shift is determined by the product of FIV and SCI. [0126]
  • The background auditory stimulus is a sound having a frequency of 234 Hz, which is within the low frequency range between 25 Hz and 440 Hz. As the area S under the variability curve in the region of interest increases, the frequency of the sound during T[0127] A is decreased, and as the area S decreases, the frequency of the sound is increased, albeit in smaller proportion than the decreases in pitch. In addition, as the height of the central tendency peak above 0.15 Hz increases, the loudness of the sound is decreased, and as the height of the central tendency peak above 0.15 Hz decreases, the loudness of the sound is increased, albeit in smaller proportion than the decreases in loudness. Modulation of loudness is determined in the same way as the modulation of brightness, which is described above.
  • In addition, the stimuli duration T[0128] A may be constant or may vary from approximately 0.06 seconds to approximately 0.18 seconds. The default duration is approximately 0.12 sec. As parasympathetic activity increases, TA is increased, and as parasympathetic activity decreases, TA is decreased, albeit in smaller proportion than the increases in duration. In other words, in this example there is an SCI for TA which is equal to |ΔS|/S. T(off), which represents the time period between consecutive auditory stimuli, may be determined in at least one of the following two ways.
  • In one embodiment of the method, T[0129] A is not synchronized with any biological cycle, and T(off) may be any arbitrary predefined time value.
  • In an alternative embodiment, the stimuli is applied during a predefined phase of a physiological cycle. As explained previously for example, the stimuli may be synchronized with the cardiac cycle, using the R fiducial point of the QRS complex as a reference. To increase parasympathetic activity, the stimuli will be applied during the latter or diastolic phase of the cardiac cycle. The delay (T[0130] W) between the fiducial point in the QRS complex of the ECG signal and the starting point of the stimuli is calculated according to equations (1), (4), (5), and (6) explained above.
  • EXAMPLE 2
  • In this example, a subject is treated for migraines and/or poor mental concentration by assisting the autonomic nervous system in stimulating or increasing sympathetic activity. Again, alternative treatments may be possible. The subject's baseline heart rate variability is initially determined by time domain and spectral analysis as described above. [0131]
  • From the patient's ECG, successive PS variability curves are obtained from successive and overlapping time windows, wherein each time window spans 90 heartbeats (i.e., k=3 and N=30k=90 heartbeats) or approximately 1.5 minutes of the ECG recording. The beginning of each window is separated from the beginning of the previous one by M=4 heartbeats, which is approximately 4 seconds. In this example, the power spectrum region of interest is below 0.15 Hz, and the stimuli is determined based on changes in the area and in the height of the highest peak in this region. [0132]
  • Because a sympathetic effect is intended, the visual stimulus will comprise colors ranging from green to red, and the audio stimulus will range from approximately 440 Hz to approximately 1320 Hz. [0133]
  • The background visual stimulus is an orange color of approximately 600 nm, which is a color between green and red. This color is displayed on a TV or computer screen. The area S under the variability curve in the region of interest is continuously monitored. If this area increases, the color of at least part of the screen is shifted toward the red end of the color spectrum, and if this area decreases, the color of at least part of the screen is shifted toward the green end of the color spectrum. Color shifts toward the green end of the spectrum are smaller than color shifts toward the red end. In this example, the SCI for color hue (θhu) is |ΔS|/S, and FIV=0.7. [0134]
  • The height of the central tendency peak A in the region of interest is also continuously monitored, and as the height of the peak increases, the brightness or intensity of the color is increased, and as the height of the peak decreases, the intensity is decreased, albeit in smaller proportion than the increase in intensity. In this example, SCI for color brightness (θbr) is |ΔA|/A. [0135]
  • The background auditory stimulus is a sound having a frequency of 880 Hz. As sympathetic activity increases, the frequency of the sound is increased, and as sympathetic activity decreases, the frequency of the sound decreases, albeit in smaller proportion than the increases in pitch. The SCI value for sound frequency (θsf) is |ΔS|/S. [0136]
  • In addition, as the height of the central tendency peak A below 0.15 Hz increases, the loudness of the sound is increased, and as the height of the central tendency peak below 0.15 Hz decreases, the loudness of the sound is decreased, albeit in small proportion than the decreases in loudness. The SCI value for sound loudness (θlo) is |ΔA|/A. [0137]
  • In addition, the stimuli duration T[0138] A may be constant or may vary from approximately 0.12 seconds to approximately 0.16 seconds. The default duration is approximately 0.14 sec. As sympathetic activity increases, TA is decreased, and as sympathetic activity decreases, TA is increased, albeit in smaller proportion than the decreases in duration. In other words, the SCI for TA is |ΔS|/S. T(off), which represents the time period between consecutive stimuli, may be determined in at least one of the following two ways.
  • As described above in Example 1, in one embodiment of the present invention, T[0139] A is not synchronized with any biological cycle, and T(off) may have any arbitrary predetermined time value.
  • In another embodiment of the present method, the stimuli is applied during a predefined phase of a physiological cycle, such as the cardiac cycle. To increase sympathetic activity, the stimuli will be applied during the early or systolic phase of the cardiac cycle, where the delay T[0140] W between the fiducial point in the QRS complex of the ECG signal and the starting point of the stimuli will be calculated according to equations (1) and (2) explained above.
  • It will be readily apparent to those in the art that variations of this method and modifications of this apparatus are possible which fall within the scope of the appended claims. [0141]

Claims (36)

1-65. (canceled)
66. A method for non-invasive treatment of a condition by affecting an autonomic nervous system of a human subject having a heart cycle wherein the method comprises:
applying a plurality of stimuli to the human subject wherein the stimuli are of a same stimulus type and have a same stimulus quality modulated by a stimulus parameter and wherein the plurality of stimuli comprise: a first stimulus applied at a time T1 and having a stimulus parameter S1; a second stimulus applied at a time T2 and having a stimulus parameter S2; and a third stimulus applied at a time T3 and having a stimulus parameter S3;
monitoring the heart cycle of the human subject during a first time period to generate first heart data wherein the first time period ends after T1;
analyzing the first heart data to assess a first activity level (A1) of the autonomic nervous system;
calculating S2 based on A1;
monitoring the heart cycle of the human subject during a second time period to generate second heart data wherein the second time period ends after T2;
analyzing the second heart data to assess a second activity level (A2) of the autonomic nervous system; and
calculating S3 based on A2.
67. The method of claim 66, wherein the stimulus type is a sensory stimulus type.
68. The method of claim 66, wherein the stimulus type is selected from the group consisting of: light, sound, taste, smell, temperature, balance and touch.
69. The method of claim 68, wherein the stimulus quality is a quality selected from the group consisting of: of brightness, intensity, loudness, color, wavelength, frequency, flavor, aroma and pitch.
70. The method of claim 68, wherein the plurality of stimuli are synchronized to the heart cycle of the human subject.
71. The method of claim 70, wherein each of analyzing the first heart data and analyzing the second heart data comprises assessing a heart rate variability.
72. The method of claim 71, wherein assessing a heart rate variability comprises performing frequency domain analysis.
73. The method of claim 72, wherein performing frequency domain analysis is accomplished using at least one method selected form the group consisting of: Fast Fourier Transform analysis, autoregressive analysis and wavelet analysis.
74. The method of claim 73, wherein each of A1 and A2 is calculated so as to reflect the activity of a same branch of the autonomic nervous system selected from the group consisting of: the parasympathetic branch and sympathetic branch.
75. The method of claim 71, wherein the first time period overlaps the second time period and the second time period overlaps the third time period.
76. A method for affecting the functioning of an autonomic nervous system of a human subject comprising:
applying a first stimulus commencing at a first time (T1) to the human subject wherein the stimulus has a stimulus characteristic determined by a first value;
monitoring a cyclic activity of the human subject during a first time period to generate first activity data wherein the first time period ends after T1;
analyzing the first activity data to assess a first activity level of the autonomic nervous system;
calculating a second value based on the first activity level;
applying a second stimulus commencing at a second time (T2) having a stimulus characteristic determined by the second value to the human subject;
monitoring the cyclic activity of the human subject during a second time period to generate second activity data wherein the second time period ends after T2;
analyzing the second activity data to assess a second activity level of the autonomic nervous system;
calculating a third value based on the second activity level; and
applying a third stimulus commencing at a third time (T3) having a stimulus characteristic defined by the third value to the human subject.
77. The method of claim 76, wherein each of the first stimulus, second stimulus and third stimulus is a sensory stimulus of a same type.
78. The method of claim 77, wherein the type is selected from the group consisting of: light, sound, taste, smell, temperature, balance and touch.
79. The method of claim 78, wherein each stimulus characteristic is a stimulus characteristic of a same type selected from the group consisting of: of brightness, intensity, loudness, color, wavelength, frequency, and pitch.
80. The method of claim 79, wherein each stimulus is of a same type selected from the group consisting of: light and sound.
81. The method of claim 78, wherein T1, T2 and T3 are synchronized to the cyclic activity of the human subject.
82. The method of claim 81, wherein the cyclic activity is a cardiopulmonary activity.
83. The method of claim 82, wherein the cyclic activity is breathing.
84. The method of claim 83, wherein each of the first activity level and second activity level is calculated so as to reflect the activity of a branch of the autonomic nervous system selected from the group consisting of: the parasympathetic branch and sympathetic branch.
85. The method of claim 76, wherein each of the first activity level and second activity level is calculated so as to reflect the activity of any one branch of the autonomic nervous system selected from the group consisting of: the parasympathetic branch and sympathetic branch.
86. The method of claim 85, wherein the first time period overlaps the second time period and the second time period overlaps the third time period.
87. Apparatus for treating a condition by affecting the autonomic nervous system of a human subject comprising:
stimuli applicator means for applying a first stimulus having a stimulus characteristic determined by a first parameter (S1) to the human subject at a first time (T1), applying a second stimulus having a stimulus characteristic determined by a second parameter (S2) to the human subject at a second time (T2) and applying a stimulus having a stimulus characteristic determined by a third parameter (S3) to the human subject at a third time (T3);
activity monitor means for monitoring a cyclic activity of the human subject during a first period ending after T1 to generate first activity data and monitoring the cyclic activity of the human subject during a second time period ending after T2 to generate second activity data; and
stimuli modulator means for analyzing the first activity data to assess a first performance (A1) of the autonomic nervous system during the first time period and analyzing the second activity data to assess a second performance (A2) of the autonomic nervous system during the second period; and wherein the stimuli modulator means calculates S2 based on A1 and calculates S3 based on A2.
88. The method of claim 87, wherein each of the first stimulus, second stimulus and third stimulus is a sensory stimulus of a same type.
89. The method of claim 87, wherein each of the first stimulus, second stimulus and third stimulus is a stimulus of a same type selected from the group consisting of: light, sound, taste, smell, temperature, balance and touch.
90. The apparatus of claim 89, wherein each stimulus characteristic is a stimulus characteristic of a same type selected from the group consisting of: of brightness, intensity, loudness, color, wavelength, flavor, aroma, frequency, and pitch.
91. The apparatus of claim 90, wherein each stimulus is of a same type selected from the group consisting of: light and sound.
92. The apparatus of claim 89, wherein the stimuli applicator means comprises means for synchronizing T1, T2 and T3 to the cyclic activity of the human subject.
93. The apparatus of claim 92, wherein the cyclic activity is a cardiopulmonary activity.
94. The apparatus of claim 93, wherein the cyclic activity is breathing.
95. The apparatus of claim 92, wherein A1 and A2 are calculated so as to reflect the performance of any one branch of the autonomic nervous system selected from the group consisting of: the parasympathetic branch and sympathetic branch.
96. The apparatus of claim 95, wherein the stimuli modulator means comprises means for performing frequency domain analysis of the first activity data and second activity data.
97. The apparatus of claim 92, where the cyclic activity is cardiovascular activity.
98. The apparatus of claim 97, wherein the activity monitor means comprises an ECG monitor.
99. The apparatus of claim 98, wherein the stimuli modulator means comprises means for performing heart rate variability analysis of the first activity data and second activity data.
100. The apparatus of claim 99, wherein the stimuli applicator means comprises an audiovisual display device.
US10/778,306 1998-10-21 2004-02-17 Method and apparatus for affecting the autonomic nervous system Abandoned US20040230252A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US10/778,306 US20040230252A1 (en) 1998-10-21 2004-02-17 Method and apparatus for affecting the autonomic nervous system
US12/169,914 US8442632B2 (en) 1998-10-21 2008-07-09 Method and apparatus for affecting the autonomic nervous system

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US17656698A 1998-10-21 1998-10-21
US71504600A 2000-11-20 2000-11-20
US10/778,306 US20040230252A1 (en) 1998-10-21 2004-02-17 Method and apparatus for affecting the autonomic nervous system

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US71504600A Continuation 1998-10-21 2000-11-20

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/169,914 Continuation US8442632B2 (en) 1998-10-21 2008-07-09 Method and apparatus for affecting the autonomic nervous system

Publications (1)

Publication Number Publication Date
US20040230252A1 true US20040230252A1 (en) 2004-11-18

Family

ID=46300870

Family Applications (2)

Application Number Title Priority Date Filing Date
US10/778,306 Abandoned US20040230252A1 (en) 1998-10-21 2004-02-17 Method and apparatus for affecting the autonomic nervous system
US12/169,914 Expired - Fee Related US8442632B2 (en) 1998-10-21 2008-07-09 Method and apparatus for affecting the autonomic nervous system

Family Applications After (1)

Application Number Title Priority Date Filing Date
US12/169,914 Expired - Fee Related US8442632B2 (en) 1998-10-21 2008-07-09 Method and apparatus for affecting the autonomic nervous system

Country Status (1)

Country Link
US (2) US20040230252A1 (en)

Cited By (72)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020049395A1 (en) * 2000-08-24 2002-04-25 Timi 3 Systems for applying ultrasound energy to the thoracic cavity
US20020072691A1 (en) * 2000-08-24 2002-06-13 Timi 3 Systems, Inc. Systems and methods for applying ultrasonic energy to the thoracic cavity
US20020082529A1 (en) * 2000-08-24 2002-06-27 Timi 3 Systems, Inc. Systems and methods for applying pulsed ultrasonic energy
US20020091339A1 (en) * 2000-08-24 2002-07-11 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to stimulating circulatory activity in a targeted body region of an individual
US20030050560A1 (en) * 2000-08-24 2003-03-13 Timi 3 Systems, Inc. Systems and methods for delivering ultrasonic energy
US20030050576A1 (en) * 2000-08-24 2003-03-13 Timi 3 Systems, Inc. Systems and methods for delivering ultrasound energy at an output power level that remains essentially constant despite variations in transducer impedance
US20030055363A1 (en) * 2000-08-24 2003-03-20 Timi 3 Systems, Inc. Systems and methods for monitoring and enabling use of a medical instrument
US20030069526A1 (en) * 2000-08-24 2003-04-10 Timi 3 Systems, Inc. Applicators that house and support ultrasound transducers for transcutaneous delivery of ultrasound energy
US20040073115A1 (en) * 2000-08-24 2004-04-15 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to increase tissue perfusion and/or vasodilation without substantial deep heating of tissue
US20040153009A1 (en) * 2003-02-05 2004-08-05 Timi 3 Systems, Inc. Systems and methods for applying audible acoustic energy to increase tissue perfusion and/or vasodilation
US20050153268A1 (en) * 2004-01-14 2005-07-14 Junkin William H. Brain-based processing skill enhancement
US20060211955A1 (en) * 2000-08-24 2006-09-21 Horzewski Michael J Systems and methods for applying ultrasound energy to stimulating circulatory activity in a targeted body region of an individual
US20070056582A1 (en) * 2005-03-18 2007-03-15 Michael Wood Methods and devices for relieving stress
US20070073361A1 (en) * 2005-09-23 2007-03-29 Bioq, Inc. Medical device for restoration of autonomic and immune functions impaired by neuropathy
US7245967B1 (en) 2002-06-12 2007-07-17 Pacesetter, Inc. Parasympathetic nerve stimulation for termination of supraventricular arrhythmias
US20070173903A1 (en) * 2005-08-30 2007-07-26 Bioq, Inc. Medical device for restoration of neurological function impaired by peripheral neuropathy
US20070185533A1 (en) * 2006-02-03 2007-08-09 Brainstate Technologies, Inc. Neurofeedback conditioning method
US20070276202A1 (en) * 2004-07-01 2007-11-29 Heikki Raisanen Method and Apparatus for Measuring and Monitoring Vital Functions or Presence
WO2008009978A1 (en) * 2006-07-20 2008-01-24 Jon Clare A medical hypnosis device for controlling the administration of a hypnosis experience
US7403819B1 (en) 2002-06-12 2008-07-22 Pacesetter, Inc. Parasympathetic nerve stimulation for control of AV conduction
US20080208084A1 (en) * 2003-02-05 2008-08-28 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to increase tissue perfusion and/or vasodilation without substantial deep heating of tissue
WO2008147958A1 (en) * 2007-05-25 2008-12-04 Thomas Jefferson University Method of treatment to improve concentration and/or memory in a subject in need thereof
WO2009039799A2 (en) 2007-09-28 2009-04-02 Oliver Zafiris Method and system for determining a reaction signal for a selected point in a data processing system following the effect of at least one input signal
US7620611B1 (en) * 1999-07-20 2009-11-17 Trilithic, Inc. Neural networks for ingress monitoring
US20090318813A1 (en) * 2000-08-24 2009-12-24 Timi 3 Systems, Inc. Applicators that house and support ultrasound transducers for transcutaneou delivery of ultrasound energy
EP2148725A1 (en) * 2007-05-11 2010-02-03 Clarimedix Visible light modulation of mitochondrial function in hypoxia and disease
US20100036267A1 (en) * 2008-08-07 2010-02-11 Wen-Chien Liao Analysis device and analyzing method for autonomic nerve state
US7691049B2 (en) 2004-03-18 2010-04-06 Respironics, Inc. Methods and devices for relieving stress
US7711430B2 (en) 2006-02-10 2010-05-04 Electrocore Llc Methods and apparatus for treating anaphylaxis using electrical modulation
US7725188B2 (en) 2006-02-10 2010-05-25 Electrocore Llc Electrical stimulation treatment of hypotension
US7747324B2 (en) 2005-11-10 2010-06-29 Electrocore Llc Electrical stimulation treatment of bronchial constriction
US20100277164A1 (en) * 2006-09-01 2010-11-04 Commonwealth Scientific And Industrial Research Organisation Method and apparatus for signal recovery
US7837719B2 (en) 2002-05-09 2010-11-23 Daemen College Electrical stimulation unit and waterbath system
US8041428B2 (en) 2006-02-10 2011-10-18 Electrocore Llc Electrical stimulation treatment of hypotension
US8088127B2 (en) 2008-05-09 2012-01-03 Innovative Pulmonary Solutions, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8172827B2 (en) 2003-05-13 2012-05-08 Innovative Pulmonary Solutions, Inc. Apparatus for treating asthma using neurotoxin
US20120164618A1 (en) * 2010-12-22 2012-06-28 Brightstar Learning Monotonous game-like task to promote effortless automatic recognition of sight words
WO2012122574A1 (en) * 2011-03-10 2012-09-13 Dawn Nowlin Method of training a living body to not react to substances
WO2012176098A1 (en) * 2011-06-20 2012-12-27 Koninklijke Philips Electronics N.V. Adapting patient room ambient stimuli to patient healing status
WO2013087895A1 (en) 2011-12-16 2013-06-20 Chordate Medical Ag Ans stimulation
US8483831B1 (en) 2008-02-15 2013-07-09 Holaira, Inc. System and method for bronchial dilation
US8740895B2 (en) 2009-10-27 2014-06-03 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US8812112B2 (en) 2005-11-10 2014-08-19 ElectroCore, LLC Electrical treatment of bronchial constriction
US8840537B2 (en) 2005-11-10 2014-09-23 ElectroCore, LLC Non-invasive treatment of bronchial constriction
US8911439B2 (en) 2009-11-11 2014-12-16 Holaira, Inc. Non-invasive and minimally invasive denervation methods and systems for performing the same
US8932218B1 (en) * 2009-02-19 2015-01-13 Jeffrey D. Thompson Methodology, use and benefits of neuroacoustic frequencies for assessing and improving the health and well-being of living organisms
US20150031003A1 (en) * 2013-07-24 2015-01-29 Aspen Performance Technologies Neuroperformance
WO2015038852A1 (en) * 2013-09-13 2015-03-19 Xhale, Inc. Unified display of cardiorespiratory parameters
US20150086952A1 (en) * 2012-05-09 2015-03-26 Koninklijke Philips N.V. Device and method for supporting a behavior change of a person
US9011355B2 (en) 2011-10-19 2015-04-21 Sympara Medical, Inc. Methods and devices for treating hypertension
US20150165201A1 (en) * 2012-09-07 2015-06-18 Hyogo College Of Medicine Swallowing assist device
US9149328B2 (en) 2009-11-11 2015-10-06 Holaira, Inc. Systems, apparatuses, and methods for treating tissue and controlling stenosis
US9198618B2 (en) 2011-12-16 2015-12-01 Chordate Medical Ab Pressure sensing system and method
US20150359888A1 (en) * 2014-06-13 2015-12-17 Palo Alto Investors Methods and Compositions for Restoring Homeostatic Capacity of a Subject
US20160023017A1 (en) * 2013-04-04 2016-01-28 Circadian Zirclight Inc. Lighting system for protecting circadian neuroendocrine function
US9398933B2 (en) 2012-12-27 2016-07-26 Holaira, Inc. Methods for improving drug efficacy including a combination of drug administration and nerve modulation
US20160256108A1 (en) * 2015-03-05 2016-09-08 Palo Alto Investors Homeostatic Capacity Evaluation
US9474684B2 (en) 2012-03-20 2016-10-25 Chordate Medical Ab Electroactive vibration method
US9579247B2 (en) 2011-12-16 2017-02-28 Chordate Medical Ab Treatment of headache disorders
US9782320B2 (en) 2011-12-16 2017-10-10 Chordate Medical Ab Double stimulation
US9872814B2 (en) 2012-03-20 2018-01-23 Chordate Medical Ag Vibration pattern for vibration stimulation
US9895279B2 (en) 2011-12-16 2018-02-20 Chordate Medical Ab Stimulation of hypothalamus
WO2018160749A1 (en) * 2017-02-28 2018-09-07 Kokon, Inc. Multisensory technology for stress reduction and stress resilience
US20190290882A1 (en) * 2018-03-22 2019-09-26 Industrial Technology Research Institute Light health care system and light health care method
US20190354815A1 (en) * 2018-05-15 2019-11-21 Palo Alto Investors LP Homeostatic Capacity Evaluation of Artificial Intelligence Systems
US11020051B2 (en) 2015-11-30 2021-06-01 Palo Alto Investors Methods of enhancing homeostatic capacity in a subject by modulating homeostatic system synchrony, and devices for use in practicing the same
US11045140B2 (en) * 2015-03-05 2021-06-29 Palo Alto Investors Homeostatic capacity evaluation
US11382513B2 (en) 2016-11-08 2022-07-12 Palo Alto Investors Methods and compositions for treating a condition in a subject
US11399727B2 (en) * 2017-11-15 2022-08-02 Technion Research & Development Foundation Limited Non-invasive analysis of sinoatrial node and autonomic nervous input to heart function
US11478605B2 (en) * 2017-09-28 2022-10-25 Panasonic Intellectual Property Management Co., Ltd. Autonomic nerve control device, autonomic nerve control system, and autonomic nerve control method
US20230099854A1 (en) * 2013-12-12 2023-03-30 Alivecor, Inc. Methods and systems for arrhythmia tracking and scoring
US11813456B2 (en) 2015-11-30 2023-11-14 Palo Alto Investors LP Methods of enhancing homeostatic capacity in a subject by increasing homeostatic system component responsiveness, and devices for use in practicing the same

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9962098B2 (en) 2006-06-02 2018-05-08 Global Cardiac Monitors, Inc. Heart monitor electrode system
US7711415B1 (en) * 2006-11-08 2010-05-04 Pacesetter, Inc. Implantable devices, and methods for use therewith, for monitoring sympathetic and parasympathetic influences on the heart
US8382667B2 (en) * 2010-10-01 2013-02-26 Flint Hills Scientific, Llc Detecting, quantifying, and/or classifying seizures using multimodal data
US20100042012A1 (en) * 2008-08-15 2010-02-18 Karim Alhussiny Diagnostic device for remote sensing and transmitting biophysiological signals
US20150201858A1 (en) * 2008-08-15 2015-07-23 Global Cardiac Monitors, Inc. Diagnostic device for remote sensing and transmitting biophysiological signals
US9943461B1 (en) 2012-02-29 2018-04-17 Frederick Muench Systems, devices, components and methods for triggering or inducing resonance or high amplitude oscillations in a cardiovascular system of a patient
US10973736B2 (en) 2012-02-29 2021-04-13 Frederick J. Muench Systems, devices, components and methods for triggering or inducing resonance or high amplitude oscillations in a cardiovascular system of a patient
US10098810B1 (en) 2013-02-27 2018-10-16 Frederick Muench Systems, devices, components and methods for triggering or inducing resonance or high amplitude oscillations in a cardiovascular system of a patient
US10632040B2 (en) 2012-02-29 2020-04-28 Frederick Muench Systems, devices, components and methods for triggering or inducing resonance or high amplitude oscillations in a cardiovascular system of a patient
US9610443B1 (en) 2013-02-27 2017-04-04 Steven G Dean Methods to trigger high amplitude oscillations or resonance in the cardiovascular system of a patient using electrical stimulation
US10765884B1 (en) 2013-02-27 2020-09-08 Steven G Dean Methods to trigger high amplitude oscillations or resonance in the cardiovascular system of a patient using electromagnetic stimulation
US10251790B2 (en) 2013-06-28 2019-04-09 Nocira, Llc Method for external ear canal pressure regulation to alleviate disorder symptoms
US9039639B2 (en) 2013-06-28 2015-05-26 Gbs Ventures Llc External ear canal pressure regulation system
FR3037808B1 (en) * 2015-06-24 2017-07-28 Codesna METHOD AND SYSTEM FOR REGULATING THE AUTONOMOUS NERVOUS SYSTEM OF A SUBJECT
US10760566B2 (en) 2016-07-22 2020-09-01 Nocira, Llc Magnetically driven pressure generator
EP3551145A4 (en) * 2016-12-06 2020-11-18 Nocira, LLC Systems and methods for treating neurological disorders
WO2018157143A1 (en) 2017-02-27 2018-08-30 Nocira, Llc Ear pumps
CN113208610B (en) * 2020-01-17 2022-06-28 广西医科大学 Phase position point estimation method and device and closed-loop nerve stimulation system

Citations (86)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3032029A (en) * 1958-07-09 1962-05-01 Thompson Ramo Wooldridge Inc System controlling apparatus and method
US3762396A (en) * 1970-10-30 1973-10-02 United Biscuits Ltd Method and apparatus for inducing sleep by applying electrical pulses to plural portions of the head
US3773049A (en) * 1970-11-13 1973-11-20 L Rabichev Apparatus for the treatment of neuropsychic and somatic diseases with heat, light, sound and vhf electromagnetic radiation
US3837331A (en) * 1972-10-24 1974-09-24 S Ross System and method for controlling the nervous system of a living organism
US3900034A (en) * 1974-04-10 1975-08-19 Us Energy Photochemical stimulation of nerves
US3967616A (en) * 1972-10-24 1976-07-06 Ross Sidney A Multichannel system for and a multifactorial method of controlling the nervous system of a living organism
US3989051A (en) * 1974-09-30 1976-11-02 Valentin Matveevich Nozhnikov Apparatus for current pulse action upon central nervous system
US4126128A (en) * 1977-04-04 1978-11-21 Haruo Takahashi Electric device for thermally stimulating autonomic nervous system
US4195626A (en) * 1977-03-29 1980-04-01 Schweizer Helgi Jon Device for the production and application of body stimuli structures
US4203452A (en) * 1977-08-10 1980-05-20 Cohen David B Efficiency of human learning employing the electroencephalograph and a long-term learning diagnostic-remedial process
US4232678A (en) * 1977-05-16 1980-11-11 Joseph Skovajsa Device for the local treatment of a patient, and more particularly applicable in acupuncture and auriculotheraphy
US4289121A (en) * 1978-11-26 1981-09-15 Kupriyanovich Leonid I Method for controlling functional state of central nervous system and device for effecting same
US4315502A (en) * 1979-10-11 1982-02-16 Gorges Denis E Learning-relaxation device
US4436097A (en) * 1982-06-07 1984-03-13 Cunningham Patrick J Cardiovascular exercise apparatus
US4503863A (en) * 1979-06-29 1985-03-12 Katims Jefferson J Method and apparatus for transcutaneous electrical stimulation
US4508105A (en) * 1983-02-08 1985-04-02 Whitten Glen A Shadow generating apparatus
US4510944A (en) * 1982-12-30 1985-04-16 Porges Stephen W Method and apparatus for evaluating rhythmic oscillations in aperiodic physiological response systems
US4573449A (en) * 1983-03-08 1986-03-04 Warnke Egon F Method for stimulating the falling asleep and/or relaxing behavior of a person and an arrangement therefor
US4616659A (en) * 1985-05-06 1986-10-14 At&T Bell Laboratories Heart rate detection utilizing autoregressive analysis
US4762131A (en) * 1985-09-02 1988-08-09 Shigeru Okuda Driller providing light and sound stimuli for stimulating the visuosensory and audiosensory nerves, for use in the rehabilitation of postapoplectic paralytic patients
US4765337A (en) * 1982-08-09 1988-08-23 Edward B. Hunter Method of improving egg production in poultry hens
US4777960A (en) * 1986-08-18 1988-10-18 Massachusetts Institute Of Technology Method and apparatus for the assessment of autonomic response by broad-band excitation
US4794934A (en) * 1985-06-27 1989-01-03 Hiroshi Motoyama Apparatus for diagnosing the functions of the internal organs and autonomic nervous system of the living body
US4834701A (en) * 1984-08-24 1989-05-30 Ken Hayashibara Apparatus for inducing frequency reduction in brain wave
US4858609A (en) * 1987-12-04 1989-08-22 Cole Roger J Bright light mask
US4862361A (en) * 1985-06-05 1989-08-29 Massachusetts Institute Of Technology Methods and apparatus for monitoring cardiovascular regulation using heart rate power spectral analysis
US4882535A (en) * 1985-12-03 1989-11-21 Benjamin Gavish Device and method for monitoring small desplacements of a peak in a frequency spectrum
US4883067A (en) * 1987-05-15 1989-11-28 Neurosonics, Inc. Method and apparatus for translating the EEG into music to induce and control various psychological and physiological states and to control a musical instrument
US4892106A (en) * 1987-10-19 1990-01-09 Gleeson Iii William J Multiple afferent sensory stimulation device
US4895149A (en) * 1983-08-29 1990-01-23 Morez Jean Bernard Apparatus for non-invasive stimulation of acupuncture points
US4902274A (en) * 1988-10-11 1990-02-20 Gleeson Iii William J Multiple afferent sensory stimulation device
US4930504A (en) * 1987-11-13 1990-06-05 Diamantopoulos Costas A Device for biostimulation of tissue and method for treatment of tissue
US4966164A (en) * 1987-12-29 1990-10-30 Tradatlantex Ag Combined sound generating device and electrical acupuncture device and method for using the same
US5036858A (en) * 1990-03-22 1991-08-06 Carter John L Method and apparatus for changing brain wave frequency
US5046494A (en) * 1990-08-27 1991-09-10 John Searfoss Phototherapy method
US5047006A (en) * 1989-11-06 1991-09-10 Howard Brandston Personal integrating sphere system
US5076281A (en) * 1988-05-31 1991-12-31 Benjamin Gavish Device and method for effecting rhythmic body activity
US5079682A (en) * 1987-08-07 1992-01-07 Roberts Seth D Arising aid
US5108361A (en) * 1987-11-10 1992-04-28 Piet Hein Device for influencing an organism
US5135468A (en) * 1990-08-02 1992-08-04 Meissner Juergen P Method and apparatus of varying the brain state of a person by means of an audio signal
US5137018A (en) * 1989-09-15 1992-08-11 Chuprikov Anatoly P Method for treating the emotional condition of an individual
US5163426A (en) * 1987-06-26 1992-11-17 Brigham And Women's Hospital Assessment and modification of a subject's endogenous circadian cycle
US5167228A (en) * 1987-06-26 1992-12-01 Brigham And Women's Hospital Assessment and modification of endogenous circadian phase and amplitude
US5197941A (en) * 1991-07-30 1993-03-30 Barbara Whitaker Portable device for controlling Circadian Rhythm disorders
US5203326A (en) * 1991-12-18 1993-04-20 Telectronics Pacing Systems, Inc. Antiarrhythmia pacer using antiarrhythmia pacing and autonomic nerve stimulation therapy
US5219322A (en) * 1992-06-01 1993-06-15 Weathers Lawrence R Psychotherapy apparatus and method for treating undesirable emotional arousal of a patient
US5242376A (en) * 1990-10-09 1993-09-07 21St Century Holdings, Inc. Relaxation device and method
US5241967A (en) * 1988-12-23 1993-09-07 Pioneer Electronic Corporation System for evoking electroencephalogram signals
US5265598A (en) * 1990-08-27 1993-11-30 Energy Spectrum Foundation Phototherapy method
US5291400A (en) * 1992-04-09 1994-03-01 Spacelabs Medical, Inc. System for heart rate variability analysis
US5292345A (en) * 1987-08-24 1994-03-08 Gerardo Ernesto A Portable photoneuronic energizer
US5295491A (en) * 1991-09-26 1994-03-22 Sam Technology, Inc. Non-invasive human neurocognitive performance capability testing method and system
US5419338A (en) * 1994-02-22 1995-05-30 City Of Hope Autonomic nervous system testing by bi-variate spectral analysis of heart period and QT interval variability
US5423328A (en) * 1993-01-20 1995-06-13 Gavish; Benjamin Stress detecting device and method for monitoring breathing
US5423325A (en) * 1993-03-12 1995-06-13 Hewlett-Packard Corporation Methods for enhancement of HRV and late potentials measurements
US5447527A (en) * 1989-04-27 1995-09-05 Suzanne Maureen Waldman Therapeutic light method
US5447528A (en) * 1989-10-30 1995-09-05 Gerardo; Ernesto Method of treating seasonal affective disorder
US5465729A (en) * 1992-03-13 1995-11-14 Mindscope Incorporated Method and apparatus for biofeedback
US5503637A (en) * 1987-06-26 1996-04-02 Light Sciences, Inc. Apparatus for producing and delivering high-intensity light to a subject
US5522854A (en) * 1994-05-19 1996-06-04 Duke University Method and apparatus for the prevention of arrhythmia by nerve stimulation
US5545192A (en) * 1987-06-26 1996-08-13 Brigham And Women's Hospital Intermittent use of bright light to modify the circadian phase
US5562719A (en) * 1995-03-06 1996-10-08 Lopez-Claros; Marcelo E. Light therapy method and apparatus
US5577990A (en) * 1993-12-28 1996-11-26 Nusa Widjaja Trophotropic response method
US5596993A (en) * 1994-09-21 1997-01-28 Beth Israel Hospital Fetal data processing system and method
US5603331A (en) * 1996-02-12 1997-02-18 Cardiac Pacemakers, Inc. Data logging system for implantable cardiac device
US5643173A (en) * 1995-09-01 1997-07-01 Welles; William F. Method and apparatus for stress relief
US5662117A (en) * 1992-03-13 1997-09-02 Mindscope Incorporated Biofeedback methods and controls
US5682901A (en) * 1993-08-03 1997-11-04 Kamen; Peter Walter Method and apparatus for measuring autonomic activity of a patient
US5694939A (en) * 1995-10-03 1997-12-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Autogenic-feedback training exercise (AFTE) method and system
US5700282A (en) * 1995-10-13 1997-12-23 Zabara; Jacob Heart rhythm stabilization using a neurocybernetic prosthesis
US5709645A (en) * 1996-01-30 1998-01-20 Comptronic Devices Limited Independent field photic stimulator
US5718235A (en) * 1992-10-06 1998-02-17 Gw Scientific, Inc. Detection of abnormal and induction of normal heart rate variability
US5727556A (en) * 1993-02-10 1998-03-17 Weth; Gosbert Method for pain therapy and/or for influencing the vegetative nervous system
US5797840A (en) * 1994-09-14 1998-08-25 Ramot University Authority For Applied Research & Industrial Development Ltd. Apparatus and method for time dependent power spectrum analysis of physiological signals
US5800481A (en) * 1995-12-28 1998-09-01 Loos; Hendricus G. Thermal excitation of sensory resonances
US5913876A (en) * 1996-02-20 1999-06-22 Cardiothoracic Systems, Inc. Method and apparatus for using vagus nerve stimulation in surgery
US5919217A (en) * 1987-12-08 1999-07-06 Medic-Light, Inc. Portable phototherapy unit
US5974262A (en) * 1997-08-15 1999-10-26 Fuller Research Corporation System for generating output based on involuntary and voluntary user input without providing output information to induce user to alter involuntary input
US6007569A (en) * 1997-02-18 1999-12-28 Frenkel; Richard E. Quantifying stress reduction and medical treatment as a result of colored light therapies
US6017302A (en) * 1997-10-31 2000-01-25 Loos; Hendricus G. Subliminal acoustic manipulation of nervous systems
US6035233A (en) * 1995-12-11 2000-03-07 Intermedics Inc. Implantable medical device responsive to heart rate variability analysis
US6063108A (en) * 1997-01-06 2000-05-16 Salansky; Norman Method and apparatus for localized low energy photon therapy (LEPT)
US6129748A (en) * 1996-03-22 2000-10-10 Kamei; Tsutomu Apparatus for applying pulsed light to the forehead of a user
US6238424B1 (en) * 1996-06-07 2001-05-29 Biolight Patent Holding Ab Device for external treatment with pulsating light of high duty cycle
US6238425B1 (en) * 1996-06-07 2001-05-29 Biolight Patent Holding Ab Device for external medical treatment with monochromatic light
US6662032B1 (en) * 1999-07-06 2003-12-09 Intercure Ltd. Interventive-diagnostic device

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4832038A (en) * 1985-06-05 1989-05-23 The Board Of Trustees Of University Of Illinois Apparatus for monitoring cardiovascular regulation using heart rate power spectral analysis
JPH0832261B2 (en) * 1988-06-10 1996-03-29 三菱電機株式会社 Mental tension monitor
JP2000070373A (en) 1998-06-19 2000-03-07 Matsushita Electric Works Ltd Biological signal feedback control device and recording medium used therefor
TW548383B (en) 1998-08-21 2003-08-21 Matsushita Electric Ind Co Ltd Radiation energy irradiation device

Patent Citations (90)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3032029A (en) * 1958-07-09 1962-05-01 Thompson Ramo Wooldridge Inc System controlling apparatus and method
US3762396A (en) * 1970-10-30 1973-10-02 United Biscuits Ltd Method and apparatus for inducing sleep by applying electrical pulses to plural portions of the head
US3773049A (en) * 1970-11-13 1973-11-20 L Rabichev Apparatus for the treatment of neuropsychic and somatic diseases with heat, light, sound and vhf electromagnetic radiation
US3837331A (en) * 1972-10-24 1974-09-24 S Ross System and method for controlling the nervous system of a living organism
US3967616A (en) * 1972-10-24 1976-07-06 Ross Sidney A Multichannel system for and a multifactorial method of controlling the nervous system of a living organism
US3900034A (en) * 1974-04-10 1975-08-19 Us Energy Photochemical stimulation of nerves
US3989051A (en) * 1974-09-30 1976-11-02 Valentin Matveevich Nozhnikov Apparatus for current pulse action upon central nervous system
US4195626A (en) * 1977-03-29 1980-04-01 Schweizer Helgi Jon Device for the production and application of body stimuli structures
US4126128A (en) * 1977-04-04 1978-11-21 Haruo Takahashi Electric device for thermally stimulating autonomic nervous system
US4232678A (en) * 1977-05-16 1980-11-11 Joseph Skovajsa Device for the local treatment of a patient, and more particularly applicable in acupuncture and auriculotheraphy
US4203452A (en) * 1977-08-10 1980-05-20 Cohen David B Efficiency of human learning employing the electroencephalograph and a long-term learning diagnostic-remedial process
US4289121A (en) * 1978-11-26 1981-09-15 Kupriyanovich Leonid I Method for controlling functional state of central nervous system and device for effecting same
US4503863A (en) * 1979-06-29 1985-03-12 Katims Jefferson J Method and apparatus for transcutaneous electrical stimulation
US4315502A (en) * 1979-10-11 1982-02-16 Gorges Denis E Learning-relaxation device
US4436097A (en) * 1982-06-07 1984-03-13 Cunningham Patrick J Cardiovascular exercise apparatus
US4765337A (en) * 1982-08-09 1988-08-23 Edward B. Hunter Method of improving egg production in poultry hens
US4510944A (en) * 1982-12-30 1985-04-16 Porges Stephen W Method and apparatus for evaluating rhythmic oscillations in aperiodic physiological response systems
US4508105A (en) * 1983-02-08 1985-04-02 Whitten Glen A Shadow generating apparatus
US4573449A (en) * 1983-03-08 1986-03-04 Warnke Egon F Method for stimulating the falling asleep and/or relaxing behavior of a person and an arrangement therefor
US4895149A (en) * 1983-08-29 1990-01-23 Morez Jean Bernard Apparatus for non-invasive stimulation of acupuncture points
US4834701A (en) * 1984-08-24 1989-05-30 Ken Hayashibara Apparatus for inducing frequency reduction in brain wave
US4616659A (en) * 1985-05-06 1986-10-14 At&T Bell Laboratories Heart rate detection utilizing autoregressive analysis
US4862361A (en) * 1985-06-05 1989-08-29 Massachusetts Institute Of Technology Methods and apparatus for monitoring cardiovascular regulation using heart rate power spectral analysis
US4794934A (en) * 1985-06-27 1989-01-03 Hiroshi Motoyama Apparatus for diagnosing the functions of the internal organs and autonomic nervous system of the living body
US4762131A (en) * 1985-09-02 1988-08-09 Shigeru Okuda Driller providing light and sound stimuli for stimulating the visuosensory and audiosensory nerves, for use in the rehabilitation of postapoplectic paralytic patients
US4882535A (en) * 1985-12-03 1989-11-21 Benjamin Gavish Device and method for monitoring small desplacements of a peak in a frequency spectrum
US4777960A (en) * 1986-08-18 1988-10-18 Massachusetts Institute Of Technology Method and apparatus for the assessment of autonomic response by broad-band excitation
US4883067A (en) * 1987-05-15 1989-11-28 Neurosonics, Inc. Method and apparatus for translating the EEG into music to induce and control various psychological and physiological states and to control a musical instrument
US5545192A (en) * 1987-06-26 1996-08-13 Brigham And Women's Hospital Intermittent use of bright light to modify the circadian phase
US5503637A (en) * 1987-06-26 1996-04-02 Light Sciences, Inc. Apparatus for producing and delivering high-intensity light to a subject
US5176133A (en) * 1987-06-26 1993-01-05 Brigham And Women's Hospital Assessment and modification of circadian phase and amplitude
US5304212A (en) * 1987-06-26 1994-04-19 Brigham And Women's Hospital Assessment and modification of a human subject's circadian cycle
US5167228A (en) * 1987-06-26 1992-12-01 Brigham And Women's Hospital Assessment and modification of endogenous circadian phase and amplitude
US5163426A (en) * 1987-06-26 1992-11-17 Brigham And Women's Hospital Assessment and modification of a subject's endogenous circadian cycle
US5079682A (en) * 1987-08-07 1992-01-07 Roberts Seth D Arising aid
US5292345A (en) * 1987-08-24 1994-03-08 Gerardo Ernesto A Portable photoneuronic energizer
US4892106A (en) * 1987-10-19 1990-01-09 Gleeson Iii William J Multiple afferent sensory stimulation device
US5108361A (en) * 1987-11-10 1992-04-28 Piet Hein Device for influencing an organism
US4930504A (en) * 1987-11-13 1990-06-05 Diamantopoulos Costas A Device for biostimulation of tissue and method for treatment of tissue
US4858609A (en) * 1987-12-04 1989-08-22 Cole Roger J Bright light mask
US5919217A (en) * 1987-12-08 1999-07-06 Medic-Light, Inc. Portable phototherapy unit
US4966164A (en) * 1987-12-29 1990-10-30 Tradatlantex Ag Combined sound generating device and electrical acupuncture device and method for using the same
US5076281A (en) * 1988-05-31 1991-12-31 Benjamin Gavish Device and method for effecting rhythmic body activity
US4902274A (en) * 1988-10-11 1990-02-20 Gleeson Iii William J Multiple afferent sensory stimulation device
US5241967A (en) * 1988-12-23 1993-09-07 Pioneer Electronic Corporation System for evoking electroencephalogram signals
US5447527A (en) * 1989-04-27 1995-09-05 Suzanne Maureen Waldman Therapeutic light method
US5137018A (en) * 1989-09-15 1992-08-11 Chuprikov Anatoly P Method for treating the emotional condition of an individual
US5447528A (en) * 1989-10-30 1995-09-05 Gerardo; Ernesto Method of treating seasonal affective disorder
US5047006A (en) * 1989-11-06 1991-09-10 Howard Brandston Personal integrating sphere system
US5036858A (en) * 1990-03-22 1991-08-06 Carter John L Method and apparatus for changing brain wave frequency
US5135468A (en) * 1990-08-02 1992-08-04 Meissner Juergen P Method and apparatus of varying the brain state of a person by means of an audio signal
US5265598A (en) * 1990-08-27 1993-11-30 Energy Spectrum Foundation Phototherapy method
US5046494A (en) * 1990-08-27 1991-09-10 John Searfoss Phototherapy method
US5242376A (en) * 1990-10-09 1993-09-07 21St Century Holdings, Inc. Relaxation device and method
US5197941A (en) * 1991-07-30 1993-03-30 Barbara Whitaker Portable device for controlling Circadian Rhythm disorders
US5295491A (en) * 1991-09-26 1994-03-22 Sam Technology, Inc. Non-invasive human neurocognitive performance capability testing method and system
US5203326A (en) * 1991-12-18 1993-04-20 Telectronics Pacing Systems, Inc. Antiarrhythmia pacer using antiarrhythmia pacing and autonomic nerve stimulation therapy
US5662117A (en) * 1992-03-13 1997-09-02 Mindscope Incorporated Biofeedback methods and controls
US5465729A (en) * 1992-03-13 1995-11-14 Mindscope Incorporated Method and apparatus for biofeedback
US5291400A (en) * 1992-04-09 1994-03-01 Spacelabs Medical, Inc. System for heart rate variability analysis
US5219322A (en) * 1992-06-01 1993-06-15 Weathers Lawrence R Psychotherapy apparatus and method for treating undesirable emotional arousal of a patient
US5718235A (en) * 1992-10-06 1998-02-17 Gw Scientific, Inc. Detection of abnormal and induction of normal heart rate variability
US5423328A (en) * 1993-01-20 1995-06-13 Gavish; Benjamin Stress detecting device and method for monitoring breathing
US5727556A (en) * 1993-02-10 1998-03-17 Weth; Gosbert Method for pain therapy and/or for influencing the vegetative nervous system
US5423325A (en) * 1993-03-12 1995-06-13 Hewlett-Packard Corporation Methods for enhancement of HRV and late potentials measurements
US5503159A (en) * 1993-03-12 1996-04-02 Hewlett-Packard Company Method for enhancement of late potentials measurements
US5682901A (en) * 1993-08-03 1997-11-04 Kamen; Peter Walter Method and apparatus for measuring autonomic activity of a patient
US5577990A (en) * 1993-12-28 1996-11-26 Nusa Widjaja Trophotropic response method
US5599274A (en) * 1993-12-28 1997-02-04 Nusa Widjaja Trophotropic response system
US5419338A (en) * 1994-02-22 1995-05-30 City Of Hope Autonomic nervous system testing by bi-variate spectral analysis of heart period and QT interval variability
US5522854A (en) * 1994-05-19 1996-06-04 Duke University Method and apparatus for the prevention of arrhythmia by nerve stimulation
US5797840A (en) * 1994-09-14 1998-08-25 Ramot University Authority For Applied Research & Industrial Development Ltd. Apparatus and method for time dependent power spectrum analysis of physiological signals
US5596993A (en) * 1994-09-21 1997-01-28 Beth Israel Hospital Fetal data processing system and method
US5562719A (en) * 1995-03-06 1996-10-08 Lopez-Claros; Marcelo E. Light therapy method and apparatus
US5643173A (en) * 1995-09-01 1997-07-01 Welles; William F. Method and apparatus for stress relief
US5694939A (en) * 1995-10-03 1997-12-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Autogenic-feedback training exercise (AFTE) method and system
US5700282A (en) * 1995-10-13 1997-12-23 Zabara; Jacob Heart rhythm stabilization using a neurocybernetic prosthesis
US6035233A (en) * 1995-12-11 2000-03-07 Intermedics Inc. Implantable medical device responsive to heart rate variability analysis
US5800481A (en) * 1995-12-28 1998-09-01 Loos; Hendricus G. Thermal excitation of sensory resonances
US5709645A (en) * 1996-01-30 1998-01-20 Comptronic Devices Limited Independent field photic stimulator
US5603331A (en) * 1996-02-12 1997-02-18 Cardiac Pacemakers, Inc. Data logging system for implantable cardiac device
US5913876A (en) * 1996-02-20 1999-06-22 Cardiothoracic Systems, Inc. Method and apparatus for using vagus nerve stimulation in surgery
US6129748A (en) * 1996-03-22 2000-10-10 Kamei; Tsutomu Apparatus for applying pulsed light to the forehead of a user
US6238425B1 (en) * 1996-06-07 2001-05-29 Biolight Patent Holding Ab Device for external medical treatment with monochromatic light
US6238424B1 (en) * 1996-06-07 2001-05-29 Biolight Patent Holding Ab Device for external treatment with pulsating light of high duty cycle
US6063108A (en) * 1997-01-06 2000-05-16 Salansky; Norman Method and apparatus for localized low energy photon therapy (LEPT)
US6007569A (en) * 1997-02-18 1999-12-28 Frenkel; Richard E. Quantifying stress reduction and medical treatment as a result of colored light therapies
US5974262A (en) * 1997-08-15 1999-10-26 Fuller Research Corporation System for generating output based on involuntary and voluntary user input without providing output information to induce user to alter involuntary input
US6017302A (en) * 1997-10-31 2000-01-25 Loos; Hendricus G. Subliminal acoustic manipulation of nervous systems
US6662032B1 (en) * 1999-07-06 2003-12-09 Intercure Ltd. Interventive-diagnostic device

Cited By (149)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090328126A1 (en) * 1999-07-20 2009-12-31 Sinde Gary W Neural networks for ingress monitoring
US7620611B1 (en) * 1999-07-20 2009-11-17 Trilithic, Inc. Neural networks for ingress monitoring
US20030050560A1 (en) * 2000-08-24 2003-03-13 Timi 3 Systems, Inc. Systems and methods for delivering ultrasonic energy
US20080167556A1 (en) * 2000-08-24 2008-07-10 Timi 3 Systems, Inc. Systems and methods for delivering ultrasound energy
US20020091339A1 (en) * 2000-08-24 2002-07-11 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to stimulating circulatory activity in a targeted body region of an individual
US20020193708A1 (en) * 2000-08-24 2002-12-19 Thompson Todd A. Systems to promote dissolution of thrombi in the thoracic cavity
US20070265601A1 (en) * 2000-08-24 2007-11-15 Timi 3 Systems, Inc. Systems and methods for monitoring and enabling use of a medical instrument
US20030050576A1 (en) * 2000-08-24 2003-03-13 Timi 3 Systems, Inc. Systems and methods for delivering ultrasound energy at an output power level that remains essentially constant despite variations in transducer impedance
US20030055363A1 (en) * 2000-08-24 2003-03-20 Timi 3 Systems, Inc. Systems and methods for monitoring and enabling use of a medical instrument
US20030069526A1 (en) * 2000-08-24 2003-04-10 Timi 3 Systems, Inc. Applicators that house and support ultrasound transducers for transcutaneous delivery of ultrasound energy
US20040073115A1 (en) * 2000-08-24 2004-04-15 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to increase tissue perfusion and/or vasodilation without substantial deep heating of tissue
US20020049395A1 (en) * 2000-08-24 2002-04-25 Timi 3 Systems for applying ultrasound energy to the thoracic cavity
US7335169B2 (en) 2000-08-24 2008-02-26 Timi 3 Systems, Inc. Systems and methods for delivering ultrasound energy at an output power level that remains essentially constant despite variations in transducer impedance
US20020082529A1 (en) * 2000-08-24 2002-06-27 Timi 3 Systems, Inc. Systems and methods for applying pulsed ultrasonic energy
US20090318813A1 (en) * 2000-08-24 2009-12-24 Timi 3 Systems, Inc. Applicators that house and support ultrasound transducers for transcutaneou delivery of ultrasound energy
US20020072690A1 (en) * 2000-08-24 2002-06-13 Timi 3 Transportable systems for applying ultrasound energy to the thoracic cavity
US20080064991A1 (en) * 2000-08-24 2008-03-13 Timi 3 Systems, Inc. Systems and methods for delivering ultrasonic energy
US7220232B2 (en) 2000-08-24 2007-05-22 Timi 3 Systems, Inc. Method for delivering ultrasonic energy
US20020072691A1 (en) * 2000-08-24 2002-06-13 Timi 3 Systems, Inc. Systems and methods for applying ultrasonic energy to the thoracic cavity
US20080195002A1 (en) * 2000-08-24 2008-08-14 Timi 3 Systems, Inc. Systems and methods for applying ultrasonic energy to the thoracic cavity
US7241270B2 (en) 2000-08-24 2007-07-10 Timi 3 Systems Inc. Systems and methods for monitoring and enabling use of a medical instrument
US20060211955A1 (en) * 2000-08-24 2006-09-21 Horzewski Michael J Systems and methods for applying ultrasound energy to stimulating circulatory activity in a targeted body region of an individual
US7837719B2 (en) 2002-05-09 2010-11-23 Daemen College Electrical stimulation unit and waterbath system
US7403819B1 (en) 2002-06-12 2008-07-22 Pacesetter, Inc. Parasympathetic nerve stimulation for control of AV conduction
US7245967B1 (en) 2002-06-12 2007-07-17 Pacesetter, Inc. Parasympathetic nerve stimulation for termination of supraventricular arrhythmias
US20100049100A1 (en) * 2002-07-24 2010-02-25 Timi 3 Systems, Inc. Systems and methods for applying audible acoustic energy to increase tissue perfusion and/or vasodilation
US20080208084A1 (en) * 2003-02-05 2008-08-28 Timi 3 Systems, Inc. Systems and methods for applying ultrasound energy to increase tissue perfusion and/or vasodilation without substantial deep heating of tissue
US20040153009A1 (en) * 2003-02-05 2004-08-05 Timi 3 Systems, Inc. Systems and methods for applying audible acoustic energy to increase tissue perfusion and/or vasodilation
US20070244415A1 (en) * 2003-02-05 2007-10-18 Timi 3 Systems, Inc. Systems and methods for applying audible acoustic energy to increase tissue perfusion and/or vasodilation
US7229423B2 (en) 2003-02-05 2007-06-12 Timi 3 System, Inc Systems and methods for applying audible acoustic energy to increase tissue perfusion and/or vasodilation
US8172827B2 (en) 2003-05-13 2012-05-08 Innovative Pulmonary Solutions, Inc. Apparatus for treating asthma using neurotoxin
US10953170B2 (en) 2003-05-13 2021-03-23 Nuvaira, Inc. Apparatus for treating asthma using neurotoxin
US9339618B2 (en) 2003-05-13 2016-05-17 Holaira, Inc. Method and apparatus for controlling narrowing of at least one airway
US20050153268A1 (en) * 2004-01-14 2005-07-14 Junkin William H. Brain-based processing skill enhancement
US20100174200A1 (en) * 2004-03-18 2010-07-08 Respironics, Inc. Methods and devices for relieving stress
US7691049B2 (en) 2004-03-18 2010-04-06 Respironics, Inc. Methods and devices for relieving stress
US8938288B2 (en) 2004-03-18 2015-01-20 Respironics, Inc. Methods and devices for relieving stress
US20070276202A1 (en) * 2004-07-01 2007-11-29 Heikki Raisanen Method and Apparatus for Measuring and Monitoring Vital Functions or Presence
US8428702B2 (en) 2005-03-18 2013-04-23 Respironics, Inc. Methods and devices for relieving stress
US8002711B2 (en) 2005-03-18 2011-08-23 Respironics, Inc. Methods and devices for relieving stress
US20070056582A1 (en) * 2005-03-18 2007-03-15 Michael Wood Methods and devices for relieving stress
US20070173903A1 (en) * 2005-08-30 2007-07-26 Bioq, Inc. Medical device for restoration of neurological function impaired by peripheral neuropathy
US20070203533A1 (en) * 2005-08-30 2007-08-30 Bioq, Inc. Implantable medical device for restoration of neurological function impaired by peripheral neuropathy
WO2007038432A3 (en) * 2005-09-23 2007-07-05 Bioq Inc Medical device for restoration of autonomic and immune functions impaired by neuropathy
US20070073361A1 (en) * 2005-09-23 2007-03-29 Bioq, Inc. Medical device for restoration of autonomic and immune functions impaired by neuropathy
WO2007038432A2 (en) * 2005-09-23 2007-04-05 Bioq, Inc. Medical device for restoration of autonomic and immune functions impaired by neuropathy
US8812112B2 (en) 2005-11-10 2014-08-19 ElectroCore, LLC Electrical treatment of bronchial constriction
US8840537B2 (en) 2005-11-10 2014-09-23 ElectroCore, LLC Non-invasive treatment of bronchial constriction
US7747324B2 (en) 2005-11-10 2010-06-29 Electrocore Llc Electrical stimulation treatment of bronchial constriction
US9037247B2 (en) 2005-11-10 2015-05-19 ElectroCore, LLC Non-invasive treatment of bronchial constriction
US20070185533A1 (en) * 2006-02-03 2007-08-09 Brainstate Technologies, Inc. Neurofeedback conditioning method
US7725188B2 (en) 2006-02-10 2010-05-25 Electrocore Llc Electrical stimulation treatment of hypotension
US8612004B2 (en) 2006-02-10 2013-12-17 ElectroCore, LLC Electrical stimulation treatment of hypotension
US8483835B2 (en) 2006-02-10 2013-07-09 ElectroCore, LLC Methods and apparatus for treating anaphylaxis using electrical modulation
US7869879B2 (en) 2006-02-10 2011-01-11 Electrocore Llc Electrical stimulation treatment of hypotension
US8233988B2 (en) 2006-02-10 2012-07-31 Electrocore Llc Electrical stimulation treatment of hypotension
US8204598B2 (en) 2006-02-10 2012-06-19 Electrocore Llc Methods and apparatus for treating bronchial restriction using electrical modulation
US8010197B2 (en) 2006-02-10 2011-08-30 Electrocore Llc Methods and apparatus for treating anaphylaxis using electrical modulation
US8041428B2 (en) 2006-02-10 2011-10-18 Electrocore Llc Electrical stimulation treatment of hypotension
US7711430B2 (en) 2006-02-10 2010-05-04 Electrocore Llc Methods and apparatus for treating anaphylaxis using electrical modulation
US8099167B1 (en) 2006-02-10 2012-01-17 Electrocore Llc Methods and apparatus for treating anaphylaxis using electrical modulation
US20100010289A1 (en) * 2006-07-20 2010-01-14 Jon Clare Medical Hypnosis Device For Controlling The Administration Of A Hypnosis Experience
WO2008009978A1 (en) * 2006-07-20 2008-01-24 Jon Clare A medical hypnosis device for controlling the administration of a hypnosis experience
US8517912B2 (en) 2006-07-20 2013-08-27 Neurocoach Limited Medical hypnosis device for controlling the administration of a hypnosis experience
US20100277164A1 (en) * 2006-09-01 2010-11-04 Commonwealth Scientific And Industrial Research Organisation Method and apparatus for signal recovery
EP2148725A1 (en) * 2007-05-11 2010-02-03 Clarimedix Visible light modulation of mitochondrial function in hypoxia and disease
EP2148725A4 (en) * 2007-05-11 2010-08-25 Clarimedix Visible light modulation of mitochondrial function in hypoxia and disease
US8845704B2 (en) 2007-05-11 2014-09-30 Clarimedix Inc. Visible light modulation of mitochondrial function in hypoxia and disease
US20100331928A1 (en) * 2007-05-11 2010-12-30 Clarimedix Visible light modulation of mitochondrial function in hypoxia and disease
EP2653190A1 (en) * 2007-05-25 2013-10-23 Thomas Jefferson University Method of improving concentration and/or memory in a subject in need thereof
US8734317B2 (en) 2007-05-25 2014-05-27 Thomas Jefferson University Method to improve concentration and/or memory in a subject in need thereof
WO2008147958A1 (en) * 2007-05-25 2008-12-04 Thomas Jefferson University Method of treatment to improve concentration and/or memory in a subject in need thereof
US20090023977A1 (en) * 2007-05-25 2009-01-22 Thomas Jefferson University Method of Treatment to Improve Concentration and/or Memory in a Subject in Need Thereof
DE112007003329B4 (en) * 2007-09-28 2013-03-07 Oliver Zafiris Determining a response signal for a selected location in an information processing system upon the action of at least one input signal
WO2009039799A2 (en) 2007-09-28 2009-04-02 Oliver Zafiris Method and system for determining a reaction signal for a selected point in a data processing system following the effect of at least one input signal
US8600696B2 (en) 2007-09-28 2013-12-03 Oliver Zafiris Method and system for determining a reaction signal for a selected location in an information processing system following the effect of at least one input signal
US20110010125A1 (en) * 2007-09-28 2011-01-13 Oliver Zafiris Method and system for determining a reaction signal for a selected location in an information processing system following the effect of at least one input signal
US11058879B2 (en) 2008-02-15 2021-07-13 Nuvaira, Inc. System and method for bronchial dilation
US8483831B1 (en) 2008-02-15 2013-07-09 Holaira, Inc. System and method for bronchial dilation
US8489192B1 (en) 2008-02-15 2013-07-16 Holaira, Inc. System and method for bronchial dilation
US9125643B2 (en) 2008-02-15 2015-09-08 Holaira, Inc. System and method for bronchial dilation
US8731672B2 (en) 2008-02-15 2014-05-20 Holaira, Inc. System and method for bronchial dilation
US8088127B2 (en) 2008-05-09 2012-01-03 Innovative Pulmonary Solutions, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8808280B2 (en) 2008-05-09 2014-08-19 Holaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8821489B2 (en) 2008-05-09 2014-09-02 Holaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8961508B2 (en) 2008-05-09 2015-02-24 Holaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8226638B2 (en) 2008-05-09 2012-07-24 Innovative Pulmonary Solutions, Inc. Systems, assemblies, and methods for treating a bronchial tree
US9668809B2 (en) 2008-05-09 2017-06-06 Holaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US10149714B2 (en) 2008-05-09 2018-12-11 Nuvaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US8961507B2 (en) 2008-05-09 2015-02-24 Holaira, Inc. Systems, assemblies, and methods for treating a bronchial tree
US20100036267A1 (en) * 2008-08-07 2010-02-11 Wen-Chien Liao Analysis device and analyzing method for autonomic nerve state
US9480812B1 (en) * 2009-02-19 2016-11-01 Jeffrey D. Thompson Methodology, system, use, and benefits of neuroacoustic frequencies for assessing and improving the health and well-being of living organisms
US8932218B1 (en) * 2009-02-19 2015-01-13 Jeffrey D. Thompson Methodology, use and benefits of neuroacoustic frequencies for assessing and improving the health and well-being of living organisms
US9017324B2 (en) 2009-10-27 2015-04-28 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US8777943B2 (en) 2009-10-27 2014-07-15 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US9931162B2 (en) 2009-10-27 2018-04-03 Nuvaira, Inc. Delivery devices with coolable energy emitting assemblies
US9005195B2 (en) 2009-10-27 2015-04-14 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US9675412B2 (en) 2009-10-27 2017-06-13 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US9649153B2 (en) 2009-10-27 2017-05-16 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US8932289B2 (en) 2009-10-27 2015-01-13 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US8740895B2 (en) 2009-10-27 2014-06-03 Holaira, Inc. Delivery devices with coolable energy emitting assemblies
US9149328B2 (en) 2009-11-11 2015-10-06 Holaira, Inc. Systems, apparatuses, and methods for treating tissue and controlling stenosis
US10610283B2 (en) 2009-11-11 2020-04-07 Nuvaira, Inc. Non-invasive and minimally invasive denervation methods and systems for performing the same
US11389233B2 (en) 2009-11-11 2022-07-19 Nuvaira, Inc. Systems, apparatuses, and methods for treating tissue and controlling stenosis
US8911439B2 (en) 2009-11-11 2014-12-16 Holaira, Inc. Non-invasive and minimally invasive denervation methods and systems for performing the same
US9649154B2 (en) 2009-11-11 2017-05-16 Holaira, Inc. Non-invasive and minimally invasive denervation methods and systems for performing the same
US11712283B2 (en) 2009-11-11 2023-08-01 Nuvaira, Inc. Non-invasive and minimally invasive denervation methods and systems for performing the same
US20120164618A1 (en) * 2010-12-22 2012-06-28 Brightstar Learning Monotonous game-like task to promote effortless automatic recognition of sight words
US9691289B2 (en) * 2010-12-22 2017-06-27 Brightstar Learning Monotonous game-like task to promote effortless automatic recognition of sight words
WO2012122574A1 (en) * 2011-03-10 2012-09-13 Dawn Nowlin Method of training a living body to not react to substances
WO2012176098A1 (en) * 2011-06-20 2012-12-27 Koninklijke Philips Electronics N.V. Adapting patient room ambient stimuli to patient healing status
US9011355B2 (en) 2011-10-19 2015-04-21 Sympara Medical, Inc. Methods and devices for treating hypertension
US9782320B2 (en) 2011-12-16 2017-10-10 Chordate Medical Ab Double stimulation
WO2013087895A1 (en) 2011-12-16 2013-06-20 Chordate Medical Ag Ans stimulation
US9579247B2 (en) 2011-12-16 2017-02-28 Chordate Medical Ab Treatment of headache disorders
US9198618B2 (en) 2011-12-16 2015-12-01 Chordate Medical Ab Pressure sensing system and method
US9451889B2 (en) 2011-12-16 2016-09-27 Chordate Medical Ab Pressure sensing system and method
US9895279B2 (en) 2011-12-16 2018-02-20 Chordate Medical Ab Stimulation of hypothalamus
US10758446B2 (en) 2011-12-16 2020-09-01 Chordate Medical Ab Treatment of headache disorders
US9486381B2 (en) 2011-12-16 2016-11-08 Chordate Medical Ab ALS treatment
US11452666B2 (en) 2011-12-16 2022-09-27 Chordate Medical Ab Treatment of headache disorders
US9872814B2 (en) 2012-03-20 2018-01-23 Chordate Medical Ag Vibration pattern for vibration stimulation
US9474684B2 (en) 2012-03-20 2016-10-25 Chordate Medical Ab Electroactive vibration method
US20150086952A1 (en) * 2012-05-09 2015-03-26 Koninklijke Philips N.V. Device and method for supporting a behavior change of a person
US20150165201A1 (en) * 2012-09-07 2015-06-18 Hyogo College Of Medicine Swallowing assist device
US9440075B2 (en) * 2012-09-07 2016-09-13 Hyogo College Of Medicine Swallowing assist device
US9398933B2 (en) 2012-12-27 2016-07-26 Holaira, Inc. Methods for improving drug efficacy including a combination of drug administration and nerve modulation
US20160023017A1 (en) * 2013-04-04 2016-01-28 Circadian Zirclight Inc. Lighting system for protecting circadian neuroendocrine function
US11577091B2 (en) 2013-04-04 2023-02-14 Korrus, Inc. Lighting system for protecting circadian neuroendocrine function
US9827440B2 (en) * 2013-04-04 2017-11-28 Circadian Zirclight Inc. Lighting system for protecting circadian neuroendocrine function
US10786685B2 (en) 2013-04-04 2020-09-29 Circadian Zirclight Inc. Lighting system for protecting circadian neuroendocrine function
US20150031003A1 (en) * 2013-07-24 2015-01-29 Aspen Performance Technologies Neuroperformance
WO2015038852A1 (en) * 2013-09-13 2015-03-19 Xhale, Inc. Unified display of cardiorespiratory parameters
US20230099854A1 (en) * 2013-12-12 2023-03-30 Alivecor, Inc. Methods and systems for arrhythmia tracking and scoring
US10835134B2 (en) * 2014-06-13 2020-11-17 Palo Alto Investors Methods and compositions for restoring homeostatic capacity of a subject
US20150359888A1 (en) * 2014-06-13 2015-12-17 Palo Alto Investors Methods and Compositions for Restoring Homeostatic Capacity of a Subject
US20160256108A1 (en) * 2015-03-05 2016-09-08 Palo Alto Investors Homeostatic Capacity Evaluation
US11045140B2 (en) * 2015-03-05 2021-06-29 Palo Alto Investors Homeostatic capacity evaluation
US10631781B2 (en) * 2015-03-05 2020-04-28 Palo Alto Investors Homeostatic capacity evaluation
US11813456B2 (en) 2015-11-30 2023-11-14 Palo Alto Investors LP Methods of enhancing homeostatic capacity in a subject by increasing homeostatic system component responsiveness, and devices for use in practicing the same
US11020051B2 (en) 2015-11-30 2021-06-01 Palo Alto Investors Methods of enhancing homeostatic capacity in a subject by modulating homeostatic system synchrony, and devices for use in practicing the same
US11382513B2 (en) 2016-11-08 2022-07-12 Palo Alto Investors Methods and compositions for treating a condition in a subject
WO2018160749A1 (en) * 2017-02-28 2018-09-07 Kokon, Inc. Multisensory technology for stress reduction and stress resilience
US11779227B2 (en) 2017-02-28 2023-10-10 Kokon, Inc. Multisensory technology for stress reduction and stress resilience
US11478605B2 (en) * 2017-09-28 2022-10-25 Panasonic Intellectual Property Management Co., Ltd. Autonomic nerve control device, autonomic nerve control system, and autonomic nerve control method
US11399727B2 (en) * 2017-11-15 2022-08-02 Technion Research & Development Foundation Limited Non-invasive analysis of sinoatrial node and autonomic nervous input to heart function
US11197975B2 (en) * 2018-03-22 2021-12-14 Industrial Technology Research Institute Light health care system and light health care method
US20190290882A1 (en) * 2018-03-22 2019-09-26 Industrial Technology Research Institute Light health care system and light health care method
US20190354815A1 (en) * 2018-05-15 2019-11-21 Palo Alto Investors LP Homeostatic Capacity Evaluation of Artificial Intelligence Systems

Also Published As

Publication number Publication date
US20080269821A1 (en) 2008-10-30
US8442632B2 (en) 2013-05-14

Similar Documents

Publication Publication Date Title
US8442632B2 (en) Method and apparatus for affecting the autonomic nervous system
Lehrer et al. Resonant frequency biofeedback training to increase cardiac variability: Rationale and manual for training
US10531827B2 (en) Apparatus and method for beneficial modification of biorhythmic activity
Hayano et al. Effects of respiratory interval on vagal modulation of heart rate
Trinder et al. Autonomic activity during human sleep as a function of time and sleep stage
US7117032B2 (en) Systems and methods for facilitating physiological coherence using respiration training
EP0872255B1 (en) Relax guiding device and biofeedback guiding device
Kar et al. Effect of sleep deprivation on functional connectivity of EEG channels
JPH0257933B2 (en)
PT1156851E (en) Method and apparatus for facilitating physiological coherence and autonomic balance
JP2007501657A (en) Electrophysiological intuition indicator
KR100287549B1 (en) Brain development apparatus with brain wave biofeedback and method thereof
Gurel et al. Automatic detection of target engagement in transcutaneous cervical vagal nerve stimulation for traumatic stress triggers
KR20140103560A (en) Electroencephalogram detection apparatus of hair band type providing sleep electroencephalogram learning memory consolidation analysis function of brain and method thereof
US20210205574A1 (en) Systems and methods for delivering sensory stimulation to facilitate sleep onset
US20050154251A1 (en) Method and apparatus for correction of functional state of person
Blasi et al. Autonomic cardiovascular control following transient arousal from sleep: a time-varying closed-loop model
Cui et al. A wearable system for cardiopulmonary assessment and personalized respiratory training
Ichimaru et al. Effect of sleep stage on the relationship between respiration and heart rate variability
Rodrigues et al. The influence of music stimulation on heart rate variability: Preliminary results
WO2006021820A1 (en) System and methods to treat headache and migraine
McClain et al. An acute bout of a controlled breathing frequency lowers sympathetic neural outflow but not blood pressure in healthy normotensive subjects
Celka et al. Photoplethysmogram combined with the Art of Eastern Pulse Analysis: a Novel Way to Assess Body and Mind Energy Balance
JP7146196B2 (en) Score calculation device and method, and score calculation device control program
US11633152B2 (en) Method of monitoring respiratory rate in a health monitoring device

Legal Events

Date Code Title Description
AS Assignment

Owner name: AUTO BALANCE LTD., ISRAEL

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KULLOK, SAUL;KULLOK, JOSE R.;REEL/FRAME:015471/0925

Effective date: 20021219

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION