Classifications

• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
  • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
  • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
  • A61B5/02405—Determining heart rate variability
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/48—Other medical applications
  • A61B5/4884—Other medical applications inducing physiological or psychological stress, e.g. applications for stress testing
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
  • A61B5/7235—Details of waveform analysis
  • A61B5/7253—Details of waveform analysis characterised by using transforms
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
  • A61B5/7235—Details of waveform analysis
  • A61B5/7253—Details of waveform analysis characterised by using transforms
  • A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms

Definitions

• the present invention relates in general to methods and apparatus for monitoring cardiovascular regulation and in particular to methods and apparatus for heart rate spectral analysis.
• Changes in cardiovascular regulation associated with congestive heart failure include attenuation of activity in the parasympathetic division of the autonomic nervous system, enhancement of activity in the sympathetic division of the autonomic nervous system, cardiac catecholamine depletion, down regulation of the beta-receptor system, increased renin-angiotensin system activity, and alteration of baroreceptor function. All of these regulatory changes require either specific clinical manipulations, such as a stress test, a Valsalva maneuver, or the like, and/or invasive maneuvers, such as cardiac biopsy, plasma catecholamine measurement, or the like, in order to determine the extent of regulatory dysfunction and its impact upon the clinical state of the patient and upon prognoses for the patient. These procedures are time consuming, and generally do not permit the formation of a clinical judgment and subsequent action within the timeframe of the course of treatment for critically ill patients in an Intensive Care Unit.
• Fluctuations from heartbeat to heartbeat in measured properties of the circulatory system reflect both the presence of a variety of naturally occurring physiological disturbances of the circulatory system homeostasis, and the dynamic response of cardiovascular control systems to these disturbances.
• the cyclic variation in intrathoracic pressure which accompanies breathing mechanically affects the return of venous blood to the heart and also affects blood pressure in pulmonary vessels and in the aorta.
• the variation in intrathoracic pressure is also coupled to a cyclic variation in heart rate through a neural mechanism mediated by the central nervous system.
• the baroreceptor reflex which is mediated by the autonomic nervous system.
• Disturbances in cardiovascular homeostasis also occur with fluctuations in the resistance of peripheral blood vessels as vascular beds regulate local blood flow to match supply with demand. These fluctuations in peripheral resistance may perturb central blood pressure and, through the baroreceptor reflex, may also lead to a compensatory variation in heart rate.
• the instantaneous rate-meter is perhaps the earliest such instrument. This meter measures each RR interval through analog or digital circuitry and displays the instantaneous heart rate.
• An improvement in the rate-meter is achieved by performing first order statistical evaluation on the RR-intervals.
• histogram displays of RR-interval differences may be generated along with their mean and standard deviations.
• Another technique for heart rate variability analysis involves the study of spectral content of the instantaneous heart rate time series.
• the computations are done on a computer.
• Such systems analyze data recorded on magnetic or punched tape. However, not only do these systems introduce additional errors during the recording process, they do not perform in real time. Furthermore, these systems are not multichannel in nature.
• a Sparse Discrete Fourier Transform algorithm which may be implemented on a personal computer (CBM 2016) and which may perform on-line monitoring of heart rate variability, based on a low pass filtered cardiac event series is disclosed in Rompelman, et al., IEEE Trans. Biomed. Engineering, BME-29, 503-510 (1982).
• a specialized hardware device also exists for low pass filtering the cardiac event series by a stepwise. convolution to create the low pass filtered cardiac event series. Coenen, et al., Medical and Biological Engineering and Computing, 15, 423-430 (1977). Nevertheless, these instruments posses a limited band width and a limited frequency resolution capability.
• an instrument which provides multi-channel spectral analysis of an instantaneous heart rate and of a respiratory activity time series.
• an instrument wherein such calculations are performed in real time at the bedside.
• An apparatus corrects artifacts in a series of heartbeats.
• Means for collecting a series of heartbeat samples are coupled to means for determining a mean interval between heartbeats.
• Means for identifying a mean variance among the intervals between heartbeats samples are coupled to means for establishing an acceptable of slewing rates as a function of the mean variance.
• Means for particularizing the absolute value of the slewing rate of a heartbeat sample relative to the mean interval are coupled to the means to determining and means for substituting the mean interval between heartbeats for all heartbeat interval samples having an absolute outside the range of acceptable slewing rates are coupled to the means for particularizing.
• a method corrects artifacts in a series of heartbeats.
• a series of heartbeat interval samples is collected and an appropriate interval between heartbeats is determined. Variances in the intervals between heartbeats are identified and an acceptable range of slewing rates is established as a function of a mean variance. An absolute value of the slewing rate of a heartbeat sample relative to the mean interval is particularized. An appropriate interval is substituted for all heartbeat interval samples having an absolute value outside the range of acceptable slewing rates.
• Apparatus calibrates a heart rate power spectrum monitor.
• Means for supplying a signal simulating a heart rate, means for generating a signal simulating a respiratory frequency fluctuation in heart rate and means for providing a signal simulating a low frequency fluctuation in heart rate are coupled to means for applying signals from these means to a heart rate power spectrum analyzer.
• Apparatus performs heart rate fluctuation power spectral analysis.
• Means for providing an electrocardiogram signal and means for supplying electroplethysmogram signal are coupled to means for obtaining a heart rate fluctuation power spectrum from an electrocardiogram signal and from an electroplethysmogram signal.
• Real time means for displaying a heart rate fluctuation power spectrum are coupled to the means for obtaining.
• Apparatus according to the present invention trends heart rate fluctuation power spectral data.
• Means for providing an electrocardiogram signal and the means for supplying an electroplethysmogram signal are coupled to means for obtaining a heart rate fluctuation power spectrum from an electrocardiogram signal and from an electroplethysmogram signal.
• Means for storing heart rate fluctuation power spectral data are coupled to means for obtaining.
• Addressable means for transmitting stored heart rate fluctuation power spectral data are coupled to the means for storing and means for converting heart rate fluctuation power spectral data into graphic form are coupled to the addressable means for transmitting.
• Real time means for displaying heart rate fluctuation power spectra are coupled to the means for converting.
• a method according to the present invention treats conditions related to malfunctions of the cardiovascular control system. A power spectrum of heart rate fluctuations in the patient are monitored.
• a level below about 0.1 (beats/min.) 2 in the power spectrum of heart rate fluctuations is identified at a frequency between about 0.04 and about 0.10 Hz as indicative of cardiovascular instability. Procedures are applied to treat the condition and thereby to increase the level of heart rate fluctuations at a frequency between about 0.04 and about 0.10 Hz.
• a method according to the present invention treats conditions related to malfunctions of the cardiovascular control system in a patient.
• a power spectrum of heart rate fluctuations is monitored in the patient.
• Heart rate in heart rate fluctuations at a frequency between about 0.04 to about 0.10 Hz is identified as indicative of cardiovascular stress. Procedures are applied to treat the condition and thereby to decrease the level of heart rate fluctuations between about 0.04 and about 0.10 Hz.
• Yet another method treats conditions related to malfunctions of the cardiovascular control system in a patient.
• a power spectrum of heart rate fluctuations in the patient is monitored.
• a ratio of the area under a heart rate power spectrum peak at a frequency between about 0.04 and 0.10 Hz to the area under a peak in the respiratory power spectrum centered at the mean respiratory rate about 0.1 Hz is identified as having an absolute value less than 2.0 for longer than or equal to about one hour as indicating of cardiac instability. Procedures are applied to treat the condition and thereby to increase the ratio.
• Still another method treats conditions related to malfunctions of the cardiovascular control system in a patient.
• a power spectrum of heart rate fluctuations in the patient is monitored.
• a ratio of the area under a heart rate power spectrum peak at a frequency between about 0.04 and 0.10 Hz to the area under a peak in the respiratory power spectrum centered at the mean respiratory rate about 0.1 Hz is identified as having an absolute value greater than or about 50 as indicating of cardiac instability.
• Procedures are applied to treat the condition and thereby to increase the ratio.
• Fig. 1 illustrates low frequency, mid-frequency and high frequency in the power spectrum of heart rate fluctuations in a dog according to the prior art
• Fig. 2 illustrates aspects of the cardiovascular control system according to the prior art
• Fig. 3 is a block diagram of apparatus for heart rate fluctuation power spectral analysis according to the present invention
• Fig. 4 illustrates address buffers and address decoding in a data acquisition device according to the present invention
• Fig. 5 illustrates components according to the present invention for interfacing an ECG apparatus with a personal computer according to the present invention
• Fig. 6 illustrates a digital to analog converter according to the present invention
• Fig. 7 illustrates a ECG trigger according to the present invention
• Fig. 8 illustrates a portable calibrator according to the present invention
• Figs. 9A and B are halves of a flow chart for software applicable to an embodiment of the present invention on a IBM personal computer;
• Fig. 10 illustrates a trend for a stable patient according to the present invention
• Fig. 11 illustrates a trend display for an unstable patient according to the present invention
• Fig. 12 is an illustration of an instantaneous heart rate according to the present invention.
• Fig. 13 is an illustration of an instantaneous heart rate fluctuation spectrum of the sort obtainable from apparatus according to the present invention.
• Fig. 14 is a stable patient's heart rate fluctuation power spectrum according to the present invention
• Fig. 15 is an unstable patient's heart rate fluctuation power spectrum according to the present invention
• Fig. 16 depicts distributions in LFP data obtained according to the present invention for stable and for unstable patients;
• Fig. 17 graphically depicts distributions of
• Fig. 18 graphically depicts data for LFP/RFP ratios according to the present invention for stable and for unstable patients.
• Power spectral methods may be used to analyze the frequency content of fluctuations in heart rate and other hemodynamic parameters.
• Short term (i.e., on a time scale of seconds to minutes) fluctuations in these parameters are concentrated in three principal spectral peaks as illustrated for a canine model in Fig. 1.
• the second identifiable spectral peak, the mid-frequency peak occurs typically between 0.1 and 0.15 Hz.
• the oscillations associated with this second peak occur at 6-9 cycles per minute, a considerably lower frequency than the respiratory frequency, and are related to the frequency response of the baroreceptor reflex.
• the third peak of the spectrum typically occurs in the frequency band of 0.04 to 0.10 Hz. This low frequency peak is related to thermoregulatory fluctuations in vasomotor tone.
• properties of the heart rate fluctuations in the conscious dog may be related to the activity of three cardiovascular control systems - the parasympathetic nervous system, the sympathetic nervous system and the renin-angiotensin system.
• Fig. 2 A simple model of the short term cardiovascular control system is illustrated in Fig. 2.
• heart rate is directly modulated by the sympathetic and parasympathetic nervous systems.
• the parasympathetic system may respond over a wide frequency range while the sympathetic system may only respond at relatively low frequencies below roughly 0.1 Hz.
• renin-angiotensin hormonal system senses blood pressure fluctuations and, through the elaboration of a substance called angiotensin II, plays the role of the guardian of the overall peripheral vascular resistance. Blockade of the renin-angiotensin system by a converting enzyme inhibitor, may remove this damping influence and may permit increased fluctuations in blood pressure and increased compensatory fluctuations in heart rate in the low frequency regime.
• the critically ill infant or child prior to, during, and after cardiac surgery at times exhibits marked changes in heart rate, blood pressure, and peripheral perfusion. These changes may be of no clinical consequence or they may indicate the existence of a major unrecognized pathology whose first outward manifestation may be sudden cardiac arrest.
• cardiovascular regulatory reserve permits objective assessment of a patient's cardiovascular stability as well as their response to medical and surgical interventions intended to improve cardiovascular function.
• a transitional microprocessor-based monitoring instrument which utilized a Z-80 microprocessor and a S-100 bus, was constructed along with a data acquisition system which interfaced the microprocessor with a Hewlett-Packard 78341 patient monitor.
• the regular respiratory signal of a patient on a ventilator was employed to obtain a respiratory spectrum and was similarly obtained through a VCO
• the respiratory frequency had to be manually entered in order to establish a fixed window for computing the power in the heart rate power spectrum in the respiratory peak.
• a smoothed heart rate "tachometer wave form" was created as follows: (1) the instantaneous heart rate time series was computed from the stored RR intervals; (2) a 1024 point time series of the instantaneous heart rate was computed from the stored instantaneous heart rate time series by sampling the latter at 4 Hz; (3) the mean heart rate computed from the 1024-point time series of instantaneous heart rate was subtracted from the smoothed series resulting in a "tachometer waveform".
• the heart rate power spectrum was computed from the heart rate "tachometer waveform" as follows: (1) a 1024-Point Fast Fourier Transform was computed using 1024 points of the tachometer cardiac tachometer waveform; and (2) the heart rate power spectrum was computed by squaring the absolute value of the previously calculated transform.
• the 256 second data segments were rejected if, (1) the patient was not in sinus rhythm; (2) transients and/or artifact were present on the cardiac "tachometer wave form"; and (3) the LFP/RFP ratios were greater than 2 standard deviations from the mean for the study period.
• the practical problems associated with this prototype monitoring instrument included the extremely tedious calculations required for use of the prototype with free-breathing patients and the large amount of data (as much as 50%, in some instances) which had to be discarded due to the presence of motion artifacts. These artifacts resulted from virtually any disturbance of the patient, even a disturbance so slight as holding the patient's hand.
• the prototype system had no capacity to identify or reject artifacts or to examine the data for dropped beats and premature triggers.
• a low value for LFP/RFP ( ⁇ 2) which is sustained for greater than one hour or a value greater than or about 50 is associated with a clinical course characterized by cardiac arrest and/or profound hypotension. At times this ratio may be the only clinical indicator of cardiovascular instability.
• the LFP/RFP ratio provides a sensitive and specific index of cardiovascular instability and may provide a clinically important, continuous, non-invasive probe of cardiovascular stability.
• a multipurpose microcomputer-based system including data basing, instantaneous heart rate and respiratory activity spectral monitor, was developed using a Hewlett Packard Series 200 Computer and Multiprogrammer as available from Hewlett-Packard.
• Advantages over the original design include: (1) error correcting routines which correct automatically for motion artifact and missed triggerings of the EKG, thus permitting a substantial increase (>30%) in available data; (2) automated trending of spectral densities along with the instanteous heart rate and respiratory activity time series; and (3) a data basing program which permits accurate temporal correlation of spectral densities with virtually every clinical intervention, routine ventilatory changes, hemodynamic, fluid monitoring and laboratory results.
• Software incorporating these advantages is included herein as Appendix A.
• FIG. 3 a block diagram of apparatus according to the present invention is illustrated.
• a source of an ECG signal 2 and a source of an electroplythsmogram signal 3 are contained within a patient monitor 4.
• a patient monitor for use with the present invention may be the System 2 Infant Monitor available from ARVEE, Incorporated, Battle Creek, Michigan.
• Source 2 is connected to an ECG trigger 5 which is in turn connected to a personal computer 7.
• Source 3 is connected to an analog to digital interface 6.
• Interface 6 is connected to analog converter 8 which is connected in turn to a personal computer 7.
• Personal computer 7 receives input from and provides output to interface 6.
• Personal computer 7 is connected to a display 9.
• Source 2 receives input from pregelled electrodes adhered to the chest wall and thigh of the patient.
• Source senses respiratory activity through a pair of electrodes by the impedence method.
• Personal computer 7 and display 9 are available as an IBM PC and a compatible display available from IBM, Incorporated, Armonk, New York. Elements 5, 6 and 8 are described below.
• address buffers and address decoding receive input from a PC bus 10.
• Nodes 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25 and 26 are respectively connected to address lines A0, A1, A2, A3, A4, A5, A6, A7, A8, A9, A10, A11,, A12, A13, A14 and A15 in PC bus 10.
• a first address buffer 100 has address inputs A0, A1, A2, A3, A4, A5, A6 and A7 which are respectively connected to nodes 11-18.
• Buffer 100 also has two gate inputs, 1G and' 2G, which are connected to ground along with a ground output GND of buffer 100.
• a power supply input V CC of buffer 100 is connected to a node 102 at a potential of +5 volts.
• a second address buffer 110 has address inputs A8, A9, A10, A11, A12, A13, A14 and A15 which are respectively connected to nodes 19-26.
• Buffer 110 also has two gate inputs, 1G and 2G, which are connected by way of a node 111 to ground.
• a ground GND output of buffer 110 is also connected to a common potential.
• Buffer 110 has a power supply input V CC which is connected to a node 112 at a potential of +5 volts.
• a status buffer 120 has address inputs A16, A17, A18 and A19 which are respectively connected to nodes 27, 28, 29 and 30. Nodes 27-30 are respectively connected to an address enable line AEN, a reset line RES, an input/output read line IOR and an input/output write line IOW in PC bus 10.
• Buffer 120 has two gate inputs, 1G and 2G, which are connected by way of a node 121 to ground. A ground output GND of buffer 120 is also connected to ground by way of node 121.
• a power supply input V CC of buffer 120 is connected to a node 122 at a potential of +5 volts.
• a data acquisition system board which is both reliable and compatible with a personal computer (PC) bus, preferably adheres to the timing requirements and the loading requirements supplied by the PC bus. This means that all connections to the PC bus should be buffered so that the load provided at any input or output of the bus is equivalent to 1 LS TTL load and high speed CMOS integrated circuits are provided for this purpose.
• PC personal computer
• address buffers are provided. This is done by buffers 100 and 110. Parts used for buffers 100, 110 and 120 are normally gated, but the gate enables, 1G and 2G, are tied to ground so that the gates are always enabled. Some of the status lines on the PC bus are buffered by a chip 120, in particular: the reset line RES; the read and write lines IOR and IOW respectively, for the input/output (10) channels; and the address enable AEN.
• An address decoder includes a three to eight line decoder 130.
• Decoder 130 has three line inputs A, B and C which are respectively connected to outputs B2, B3 and B4 of buffer 100.
• Decoder 130 has gate inputs G2A and G2B which are respectively connected to outputs B5 and B6 of buffer 100.
• a power supply VCC input of decoder 130 is connected to a node 131 at a potential of +5 volts while a ground GND output of decoder 130 is connected to a common potential.
• Outputs Y0, Y1, Y2, Y3, Y4, Y5, Y6 and Y7 are connected to inputs, of a NAND gate 140.
• a NAND gate 151 has an input connected to each of outputs B8, B9 and B10 of buffer 110.
• An output B100 of buffer 110 is connected to an input of an inverter 152 which has an output connected to an input of NAND gate 151.
• outputs B12, B13, B14 and B15 of buffer 110 are respectively connected to an input of each of inverters 153, 154, 155 and 156, each of which has an output connected to an input of NAND gate 151.
• NAND gate 151 has an output connected to an input of an inverter 157.
• a NAND gate 158 has an input connected to anoutput of inverter 157 and has an output connected to an input of an inverter 159.
• An inverter 160 has an input connected to an output B7 of buffer 100 and has an output connected to an input of NAND gate 158.
• an inverter 161 has an input connected to an output B16 of buffer 120 and has an output connected to an input of NAND 158.
• An output of inverter 159 is connected to a gate input Gl of decoder 130.
• address decoding is provided.
• a fixed address location location hex 700 to 71F (a total of 32 channels).
• the decoding of the fixed upper bytes in the address is provided by a combination of nine inverting gates, 152, 153, 154, 155, 156, 157, 159, 160 and 161, and NAND gates 151 and 158.
• These elements in combination with decoder 140, provide chip enable signals which can be used to select one or another of the functional chips on our board.
• Each of the eight chip enable signals correspond to a block of four channels .
• a chip select #0 from output to of decoder 130 corresponds to channels hex 700, 701, 702 and 703.
• a logic network for driving a data buffer includes a NAND gate 171, an inverter 172 and a NAND gate 173.
• An output of inverter 172 is connected to a first input of NAND gate 173 while an output of NAND gate 140 is connected by way of a node 174 to a second input of NAND 173 and to a first input of a NAND gate 175.
• a second input of NAND gate 175 is connected to an output of NAND gate 171.
• a node 181 is connected to an output BO of buffer 100.
• a node 182 is connected to an output Bl of buffer 100.
• Nodes 183 and 184 are respectively connected to output Y0 and output Y7 of decoder 130.
• Nodes 185 and 186 are respectively connected to an output of NAND gate 175 and an output of NAND gate
• a node 187 is connected to an output B17 of buffer 120.
• a node 188 is connected an output B18 of buffer 120, to a first input of NAND gate 171 and to an input of inverter 172.
• a node 189 is connected to a second input of NAND 171 and to an output B19 of buffer 120. Additional chips are used to provide logic which drives a data buffer connected to a data bus.
• the data bus is bidirectional in order to both transmit data to and from devices on the board. In order that this be accomplished, one must determine at any time whether or not data is either being read from or written to the board.
• NAND gate 171 NAND gate 173, AND gate 175 and inverter 172 which translates the read and write signals for the input/output (IO) channel into an output enable and a transmit enable for a data buffer.
• IO input/output
• the apparatus of Fig. 4 may be used to properly interface a device to the PC bus 10.
• components according to the present invention for interfacing an ECG apparatus with a personal computer include a port expander 200.
• Port expander 200 has four sets of 8 nodes each, the four sets correspond to four ports A, B, C and D.
• the outputs for port A are A0, A1, A2, A3, A4, A5, A6 and A7.
• the inputs corresponding to port B are B0, B1, B2, B3 , B4, B5, B6 and B7.
• Outputs corresponding to port C are C0, C1, C2, C3, C4, C5, C6 and C7.
• a set of outputs corresponding to port D includes D0, D1, D2, D3, D4, D5, D6 and D7.
• Expander 200 has a chip select input CS connected to node 184. Expander 200 also has a read input RD and a write input WR respectively connected to nodes 188 and 189.
• Expander 200 has two address inputs, AD0 and AD1 which are respectively connected to nodes 181 and 182.
• a reset RES input of expander 200 is connected to node 187.
• Inputs A0, Al, A2, A3, A4, A5 , A6, A7 are respectively connected to nodes 291, 292, 293, 294, 295, 296, 297 and 298.
• Outputs D0-D7 are respectively connected to nodes 208, 207, 206, 205, 204, 203, 202 and
• a power supply input V CC of expander 200 is connected to a node 209 at a potential of +5 volts.
• a ground GND output of expander 200 is connected to a common potential.
• Port expander 200 is used to overcome the low speed of the data bus on both A/D converter 260 and a digital analog converter. This permits slowing down the read and write signals inasmuch as they may be provided artifically on port C of expander 200 or as chip select signals from address decoder 130.
• Port C of expander 200 is a bit addressable register which allows one to individually select or deselect bits without affecting any of the other bits. This is accomplished by sending a one byte command to expander 200. Because expander 200 is given the control function, the address of expander 200 is the highest address in the set of channels. In other words, expander 200 occupies 10 channels hex 71C to hex 71F.
• the ports A, B and C on expander 200 are addresses 71C, 71D and 71E, respectively, and the control register internal for expander is at input/output I/O channel 71F.
• a timer 220 has two address inputs, AD0 and AD1 respectively connected to nodes 181 and 182.
• Timer 220 also has a read input RD connected to node 188, a write input WR connected to node 189 and a chip select input CS connected to node 184.
• a first gate input GO is connected to the CO of expander 200 while a second gate input G1 and a third gate input G2 are both connected by way of a node 223 to output C1 of expander 200.
• Timer 220 has three clock inputs CLK0, CLK1 and CLK2, of which CLK1 is connected by way of node 222 to an output OUT0 of timer 220 and input CLK2 is connected to an output OUT1 of timer 221 by way of a node 31.
• An interrupt request line IRQ4 within PC bus 10 is also connected to node 31.
• An output OUT2 is connected to a non-inverting input of an operational amplifier 224, an inverting input and a output of which are connected to a node 400.
• a power supply input V CC of timer 220 is connected to a node 221 which at a potential of +5 volts.
• Timer 221 has seven outputs D0, D1, D2, D3,
• Timer 220 includes three 16 bit timers which are addressed at hex locations 704, 705, 706, and 707. In other words, they are provided by chip select 1.
• the three clocks on timer 220 are connected in series which effectively converts it into a 48 bit counter. However, in the operation of the program, some of the bits in this counter are thrown away because the reset values are less than 65,536.
• the three clock registers are used in the following way.
• Counter 0, corresponding to input CLK 0, counts an onboard time base to be discussed later and provides an output which gives the minimum resolution of the heart rate counting. In other words. it provides the counter time base for measuring the heart rate.
• Counter #1 corresponding to input CLK 1, counts the heart rate counter time base and provides as an output an interrupt at IRQ4. This signal drives the sampling of the respiratory signal at a constant frequency, and is also used to measure interbeat intervals. In the standard data collecting mode, where one is interested in measuring the respiratory signal at 4 hertz intervals, this means that the counter 0 is set to generate output pulses at 11 microsec. intervals and that these pulses are in turn counted by counter 1 to generate 4 hertz pulses which are used to drive data acquisition from the respiratory signal.
• the last counter register, counter #2, corresponding to input CLK2, is used to count the number of respiratory sampling pulses which have been supplied. This functions as an overflow counter and always has the reset value of 65,536.
• the counter measuring interbeat intervals effectively overflows only every 65,536 respiratory sampling times, which is far in excess of what would be required to recover dropped beats which occur because the heart rate is not adequately detected.
• a counter 240 has an input 1A connected to a clock line PC CLK in PC bus 10 by way of a node 32.
• Counter 240 has a first output IQA connected to the CLKO input of timer 220.
• Counter 240 has a secnd output 1QB and has a third output 1QC.
• a clear input CLR1 of counter 240 and a ground output GND of timer 240 are connected to a common potential by way of a node 242.
• a data output buffer 280 has an output enable input OE connected to node 185 and has a tranfer enable input TE connected to a node 186.
• Eight data inputs, A0, A1, A2, A3, A4, A5, A6 and A7, of buffer 280 are respectively connected to nodes 208, 207, 206, 205, 204, 203, 202 and 201.
• a power supply V CC input of buffer 280 is connected to a source of potential at +5 volts.
• a ground GND output of buffer 280 is connected to a common potential.
• Outputs B0, B1, B2, B3, B4, B5, B6 and B7 of buffer 280 are respectively connected to data lines in PC bus 10 by way of nodes 33, 34, 35, 36, 37, 38, 39 and 40.
• the time base for this clock system is provided by counter 240.
• Timer 220 counts only at a rate of 2.6 MHz megahertz which is exceeded by the IBM PC bus clock of 4.77 megahertz.
• the IBM PC bus clock is divided by 2 using counter 240 and the result used to provide a time base at 2.38 megahertz for timer 220.
• the 4.77 megahertz clock is also divided by 8 to provide a 596 kilohertz clock which is used to drive an analog to digital (A/D) converter.
• A/D converter 260 uses this clock signal in order to properly execute the successive approximation scheme to convert analog inputs into digital outputs.
• A/D converter 260 has an output enable input OE connected to output C4 of expander 200.
• A/D converter 260 also has three inputs A, B and C which are respectively connected to outputs C5, C6 and C7 of expander 200.
• a clock input CLK of A/D converter 260 is connected to the 1QC output of counter 240.
• An address latch enable ALE and a start input STR of A/D converter 260 are connected to a node 261.
• a power supply V CC input and a reference voltage +V REF input of A/D converter 260 are connected to a node 262 at a potential of +5 volts.
• A/D converter 260 are connected to a common potential by way of a node 263.
• A/D converter 260 has seven outputs D0, D1, D2, D3, D4, D5, D6 and D7 which are respectively connected to inputs B0, B1, B2, B3, B4, B5, B6 and B7 of expander 200.
• A/D converter 260 has an end of count EOC output connected to a first input of the NAND gate 264, an output of which is connected to an input of an inverter 265.
• a second input of NAND gate 264 is connected to an output of an inverter 266 which has an input connected to node 187.
• An output of inverter 265 is connected to node 261.
• A/D converter 260 has a signal input IN connected to a node 267.
• An output of an operational amplifier 268 is connected to node 267 and to a first lead of a resistor 269.
• a second lead of resistor 269 is connected to a first lead of resistor 270, a second lead of which is connected to a source of potential at -5 volts.
• the first end of resistor 270 is also connected to an inverting input of amplifier 268 and to a first end of a resistor 271.
• a non-inverting input of amplifier 268 is tied to ground.
• a second end of resistor 271 is connected to a node 272 which provides an analog signal input ANA IN for the apparatus according to the present invention.
• A/D converter 260 is connected to port B of port expander 200. This A/D has built into it its own 8 channel analog multiplexer which allows the selection of one of eight analog signals to be converted. The channel select corresponding to inputs A. B and C of converter 260 is connected to port C on bytes 5, 6 and 7.
• analog input at input IN should be in the range of 0 to 5 volts or an input buffer should be supplied to alter this input range.
• input IN should always be provided with an analog buffer to provide isolation for both the computer and the instrument being monitored.
• the input buffer is provided by operational amplifier 268. This amplifier converts a bipolar analog input of plus or minus 5 volts to a single unipolar input of 0 to 5 volts at input IN. This analog input is used to monitor the respiration.
• A/D converter 260 is set up in a free running mode such that it continuously does conversions on the analog signal.
• the end-of-conversion pulse at output EOC is used to generate a start pulse for the A/D so that as soon as an end of conversion occurs it a new conversion is started. This is the reason for the two gates connected between end of conversion output EOC and the start input STR.
• the reset line at node 187 is also used to generate a start pulse. This means that the device will always function even after being powered up.
• the address latch enable ALE which is used to latch in the address value for the channel to be monitored, is re-latched at every start pulse.
• digital analog (D/A) converter 300 has inputs D0, D1, D2, D3, D4, D5, D6 and D7 which are respectively connected to nodes 298, 297, 296, 295, 294, 293, 292 and 291 as illustrated in Fig. 5.
• Converter 300 has a write WR input connected to node 183 and has a feedback input RFB.
• Converter 300 also has a power supply V CC input, a reference voltage V REF input and an input latch enable input ILE all of which are connected to a source of potential at +5 volts by way of a node 301.
• Converter 300 has an analog ground AGND and a digital ground output DGND, both of which are connected by way of a node 302 to a common potential.
• Converter 300 has a first output OUT1 and a second output OUT2 which are respectively connected to an inverting and a non-inverting input of an operational amplifier 303.
• the non-inverting input of amplifier 303 is also connected to a common potential by way of a node 305.
• .Amplifier 303 has an input connected to a node 306 at a potential of +12 volts and an input connected to a node 307 at a potential of -12 volts.
• An output of amplifier 303 is connected to a node 308 which is connected to the RFB input of converter 300 and to a first end of a variable resistor 309.
• a second lead of variable resistor 309 is connected to a first lead of a variable resistor 310, a second lead of which is connected to a node 311 at a potential of +5 volts.
• the second lead of resistor 309 is also connected to an inverting input of operational amplifier 312 and to a first lead of a resistor 313.
• a non-inverting input of amplifier 312 is connected to ground.
• a second lead of resistor 313 is connected to an output of amplifier 312 and to a node 391 which serves as an analog output for the apparatus according to the present invention.
• Port A of expander which is at location 71C, is attached to a D/A converter data bus which, includes nodes 291-298.
• the write latch signal for the D/A converter is provided by chip select #0. In other words, any dummy byte written to any of the addresses
• D/A converter 300 chip is in the form of differential currents generated at outputs OUT 1 and OUT 2.
• a system having two operational amplifiers is employed to convert these currents to a voltage.
• Amplifier 303 is a differential current to voltage converter which provides a signal from 0 to 5 volts.
• Amplifier 312 converts the signal to a bipolar plus or minus 5 volt signal.
• D/A converter 300 Feedback control for the current to voltage converter is provided in D/A converter 300 through input RFB so that in actuality three connections are made from the D/A chip to the first operational amplifier. Because the D/A converter is an 8 bit device, this provides 256 voltage levels which are linearly distributed between plus and minus 5 volts. This D/A output may be used to generate calibrating signals or other control signals. As illustrated in Fig. 7, a source of an ECG signal is connected by way of a node 400 to a non-inverting input of an operational amplifier 401 in an ECG trigger 60. An input of amplifier 401 is connected to a node 402 at a potential of plus 12 volts.
• An inverting input of amplifier 401 is connected to an output of amplifier 401 and to a non-inverting input of an operational amplifier 406.
• a first lead of each of resistors 403a, 403b, 403c, 403d, 403e, 403f, 403g, 403h and 403i is connected to the output of amplifier 401 while the second lead of resistor 403i is permanently connected and a second lead of one other of resistors 403a through h is connected to a node 410 by a jumper.
• a first lead of capacitor 404 is connected to node 410 while a second lead of capacitor 404 is connected to a node 405 at a potential of minus 12 volts.
• An inverting input of amplifier 406 is connected to a cathode of a diode 407, an anode of which is connected to an output of amplifier 406.
• the cathode of diode 407 is also connected to a first lead of capacitor 408 and a first lead of each of resistors 410a, 410b,
• a non-inverting input of an operational amplifier 412 is also connected to the cathode of diode 407 while an inverting input of amplifier 412 is connected to the output of amplifier 406.
• An input of amplifier 412 is connected to a node 413 at a potential of minus 12 volts.
• a first lead of resistor 414 is connected to the output of amplifier 412 while a second lead of resistor 414 is connected to a cathode of a diode 415 an anode of which is connected to ground.
• the cathode of diode 415 is also connected to an input of a Schmitt trigger 416 an output of which is connected to a line designated IRQ 3 in PC bus 10 by way of a node 491.
• ECG trigger 60 has an input buffer consisting of a non-inverting buffer of an amplifier 401 which isolates the ECG signal from the rest of the board. As illustrated in Fig. 5, the EKG trigger functions in the following manner.
• the R wave which is larger than any other signal in the ECG, causes capacitor 408 to charge up to a certain value corresponding to the peak of the R wave. Any values beneath the peak of the R wave will be rejected by amplifier 403 so that no output occurs.
• the voltage on capacitor 405 decays slowly with a rate given by the RC time constant of capacitor 405 and the resistance across elements 410a- f.
• the voltage on the capacitor is sent to the inverting input on amplifier 403 and is used as a threshold for the R wave of the EKG.
• the electrocardiogram is being passed to the non-inverting input of amplifier 406, the only time that the operational amplifier has a positive output is when the EKG signal is larger than the voltage on capacitor 405. Whenever this occurs, capacitor 408 is immediately charged up to the value at the EKG input. In other words, the voltage on capacitor 408 is a sort of envelope on the top of the electrocardiogram, although its decay rate is limited by the RC time constant. Diode 407 insures that the envelope function which is provided by capacitor 408 is the upper envelope and not the lower envelope. The lower envelope is provided by reversing the polarity of diode 407.
• the RC network of capacitor 405 and resistors 403a-i provides a low pass filtered ECG.
• the voltage on capacitor 405 is the baseline for the ECG, which may vary.
• the array of jumper selected resistors 410a-e allows variation of the time constant of the RC network containing resistors 406a-e and capacitor 408.
• this latter network which monitors the ECG envelope is referenced to the ECG baseline present on capacitor 404 permitting accurate tracking of the envelope and therefore better R wave detection.
• the jumpers may be replaced with analog switches controlled by the personal computer in order to give the computer control of RC time constant selection.
• An output from ECG trigger 60 is generated by connecting amplifier 412 in parallel with peak detector amplifier 406 so that the inputs are reversed.
• ECG trigger 60 may be modified to allow selection of various decay rates for the envelope and also to provide a floating threshold for the 0 point of the EKG.
• the ECG triggers if the R wave passes above 0 volts. However, it can be imagined that sometimes the baseline will drift far enough below 0 volts that the R wave does not cross 0 volts and in such a case this trigger would never detect the R wave. This is corrected by connecting the second leads of the charging capacitor 408 and On the selected discharging resistor of 406a-f may be connected to a low pass filter consisting of a capacitor 405 and a selected one of resistors 403a-f (to choose various discharge rates) which low pass filters the electrocardiograms and essentially selects out the baseline. This means that instead of measuring the R wave with respect to 0 volts, the R wave may be measured with respect to the floating baseline of the electrocardiogram.
• the jumper selected resistor selects an RC time constant much greater than the RR interval. So long as the baseline does not drift faster than one R wave in approximately 10 heart beats, this means that this trigger will successfully detect all R waves. Selecting one of resistors 410a-f allows variation of the RC time constant of elements 408 and 410a-f.
• an operational amplifier 500 has a non-inverting input connected to a first lead of each of resistors 501, 502 and 503.
• a second lead of resistor 501 is connected by way of a node 503a to a positive voltage source while a second lead of resistor 502 is connected by way of a node 504 to a negative voltage source.
• An inverting input of amplifier 500 is connected to a first lead of a capacitor 505, a second lead of which is connected by way of a node 506 to a negative voltage source.
• the inverting input of amplifier 505 is also connected to a first lead of a variable resistor 507 and to a first lead of a resistor 508 a second lead of which is connected to an output of amplifier 500.
• the output of amplifier 500 is also connected to a second lead of resistor 503.
• Amplifier 500 has an input connected by way of a node 509 to a positive voltage source and by way of a node 510 to a negative voltage source.
• a second lead of resistor 507 is connected to a non-inverting input of an amplifier 511, an inverting input of which is connected to an output of amplifier 511 by way of a node 591 which provides an output port for a simulated respiratory frequency.
• a first lead of a resistor 512 is connected to node 591 while a second lead of resistor 512 is connected to a first lead of a resistor 513 and to a first lead of a capacitor 514, a second lead of which is connected by way of a node 515 to a negative voltage source.
• a second lead of resistor 513 is connected to an output of an operational amplifier 514 and to an inverting input of amplifier 515 is connected to a first lead of a resistor 516, to a first lead of a capacitor 517 and to an inverting input of an operational amplifier 518.
• the second lead of capacitor 517 is connected by way of a node 519 to a negative voltage source.
• a non-inverting input of amplifier 518 is connected to a first lead of each of resistors 520, 521 and 522.
• a second lead of resistor 520 is connected by way of a node 523 to a positive voltage source while a second lead of resistor 521 is connected by way of a node 524 to a negative voltage source.
• a second lead of resistor 522 is connected to an output of amplifier 518 and to a second lead of resistor 516.
• An inverting input of an operational amplifier 525 is connected to the first lead of resistor 513 and to a first lead of a variable resistor 526.
• a non-inverting input of amplifier 525 is connected to a first lead of each of resistors 527, 528 and 529.
• a second lead of resistor 527 is connected to a node 530 at a positive potential while a second lead of resistor 528 is connected by way of a node 531 to a negative voltage source.
• a second lead of resistor 529 is connected to a second lead of resistor 526 and to an output of amplifier 525 at a node 592 which provides a square wave output simulating a modulated heart rate pulse.
• a first lead of a capacitor 532 is connected to node 592 while a second lead of capacitor 532 is connected by way of a node 593 to a first lead of a resistor 533, a second lead of which is connected to ground.
• Node 593 provides an output port for a spike output simulating the R wave of an EKG.
• the source of positive potential for the portable calibrator 70 may be at a voltage between about plus 5 and about plus 18 volts.
• the negative voltage source for portable calibrator 70 may be at a potential of about minus 18 volts to about minus 5 volts.
• Portable calibrator 70 provides test signal for the heart rate spectral analysis hardware which, although not of a truly calibrated nature, does allow one to evaluate whether or not the software and hardware is functional.
• Each of the output signals provided is a triangle wave which represents the respiration and a frequency modulated pulse train representing the heart rate.
• the modulation of the heart rate is provided at two frequencies which simulate a respiratory modulation and also a low frequency modulation.
• the basic circuit of calibrator 70 for providing each pulse train consists of an oscillator having one operational amplifier as typified by the respiratory frequency modulator.
• a charging capacitor 505 and a variable resistor 507 provide an RC circuit which is charged by the output of the amplifier 500. It is also discharged by the amplifier 500 when the output of the amplifier 500 is low.
• Progressive cycles of the oscillator consist of charging and discharging the capacitor at the rate prescribed by the RC circuit.
• the reference level which determines whether or not one is discharging or charging is provided at the non-inverting input of the amplifier 500.
• capacitor 505 begins as being completely discharged, then the voltage at the inverting input for the operational amplifier 500 is low.
• the output of the operational amplifier 500 is therefore high and this means that the input at the non-inverting input is 2/3 the voltage between the negative voltage source V and the positive voltage source V+.
• the capacitor 505 begins to charge.
• the threshold voltage at the non-inverting input is changed and now.becomes only 1/3 the way from the negative voltage source to the positive voltage source.
• the voltage on the charging capacitor 505 varies between 1/3 and 2/3 the difference between the negative and the positive voltage source. This determines the range of output on capacitor 505.
• the voltage at capacitor 505 is buffered by a non-inverting buffer 511 and this provides the respiratory signal at node 591.
• An identical oscillator is used to provide low frequency modulation.
• the difference in the two frequencies is obtained by adjusting the respective variable resistors, 505 and 517,which set the RC time constants.
• the outputs of these two modulators are fed by resistors 512 and 513 into the charging capacitor 514 for the heart rate.
• the heart rate oscillator is similar in design and consists of variable resistor 526 and capacitor 532 which charges and discharges in cycles with the range of voltages on the capacitor ranging between 1/3 the distance from the negative voltage source to the positive voltage source to 2/3 the voltage between the negative voltage source and the positive voltage source.
• Resistors 512 and 513 which connect the outputs of the low frequency and respiratory frequency modulators to the heart rate modulator, allow a small amount of current to flow into charging capacitor 514 of the heart rate modulator. This alters the charging rate of capacitor 514 and thereby affects the rate at which the heart rate oscillator oscillates. For example, on a positive cycle of the respiratory frequency modulator, the heart rate capacitor is charging more rapidly towards the plus side because more current is being supplied on the plus side of the cycle.
• the output of the heart rate modulator is sent through an RC filter comprising capacitor 532 and resistor 533 which converts the square wave output of the heart rate modulator into a spike output which may be sent to an R wave detector.
• the spike output includes both positive and negative spikes so that an R detector which depends on a high frequency filtering function may be discharging at twice the heart rate, inasmuch as it may trigger on both positive and negative spikes.
• a block diagram may be constructed for the main program (designated SYNCTS19) and for sub-routine modules (SYNC7s, GWINDOW3, and FGRAPH8 ) .
• This block diagram may be used in order to better interpret a complete program for heart rate fluctuation spectral analysis useful on an IBM personal computer, as illustrated in Appendix B.
• programs are provided for a Hewlett-Packard and an IBM computer herein, the software and other aspects of the present invention may be readily modified for use with other mini- and micro-computers.
• Appendix B is a routine for removing artifacts from a detected heart rate provided for by an electrocardiograph machine.
• This program computes histograms from the heart rate data in order to generate a tachometer waveform. The most common rate on the histogram is selected as the correct rate and other rates are interpreted in light of it. Specifically, in order to correct for a spurious extra trigger, where a first and a second beat are close together while a third beat is spaced at an abnormally long interval, the second beat is discarded if the first beat to second beat interval is less than a predetermined value. The resulting interval between the first and the third beats is divided by an integer in order to provide a more normal intrabeat interval.
• the intrabeat interval is divided by that multiple, most commonly two, in order to provide a more correct interval length. If the slewing rate of the heartbeat samples is outside of an acceptable range of slewing rates determined as a function of a mean variance, and the problem cannot be identified as a missed trigger or as a spurious extra trigger, or if the three previous intervals have been corrected, a determined mean interval, against which all other intervals are judged, is substituted for the inappropriate interval.
• the slew rate is calculated on a moving average of the heart rate waveform and corrects for triggers that fall within the parameters of 0.05 Hz (3 beats/min.) per beat and five times the maximum slew. This artifact-correcting routine never slews more than 10 percent of the heart rate waveform.
• a graphic routine for trending heart rate fluctuation spectral data The parameters of LFP, RFP, LFP/RFP ratio and heart rate are plotted on a graph over time to show trends in the four parameters. These trends may then be studied in order to examine the effects of various clinical interventions. Values for the parameters heart rate, LFP/RFP ratio, LFP and RFP are stored and may be called up at any point in time through a graphin routine in order to provide a graphic depiction of the course of a patient's condition. This sort of graphic depiction is illustrated for a stable patient in Fig. 10 and for an unstable patient in Fig. 11.
• a routine is provided for the segmentation of data and subsequent reanalysis.
• data from the analog to digital converter 260 is collected continuously into a buffer and is dumped to a disk in blocks of 1,024 numbers (2,048 bytes equals 1,024 words and each block is referred to as a record or EPOCH).
• the time of heartbeat occurrence as measured by the signal provided by outputs OUT1 and OUT2 of timer 220 are collected continuously into two buffers (hb buffer 1 and hb buffer 2). These times are dumped to the disk in blocks of 1,024 pairs of numbers (1,024 from each buffer which equals 2,048 bytes or 1,024 words each).
• the heart rate is less than the sample rate of A/D converter 260 as required by signal processing, there are fewer heartbeat disk dumps.
• the A/D and heartbeat data must correspond to the same time interval for the purpose of doing correlations.
• the correspondence may be determined from (1) the record number in a A/D file and (2) the absolute of the times stored in the heartbeat file (time differences used for intrabeat intervals).
• the instantaneous heart rate signal is generated backwards in time from the heartbeat corresponding to the last A/D sample in the record of interest. This means that if the heart rate signal is analyzed on a frequency scale not corresponding to the respiration data (e.g.
• the heart rate waveform extends backwards in time beyond the beginning of the present A/D record. This means that the heart rate waveform overlaps the heart rate waveform corresponding to the previous A/D records. Overlapping permits lower frequency analysis than would be possible if only data corresponding to the present record were used (as in the prototype apparatus). Also, overlapping leads to the smoothing of parameters and to the subsequent reduction of fluctuating artifacts. In addition, it becomes less critical at what point analysis begins.
• a calibration program providing a software driven calibrator which may provide more realistic spectral data than the portable calibrator of Fig. 8, is contained within the program of Appendix A for a Hewlett-Packard micro-computer.
• Appendix C is a program which, although not tested, is believed to provide the same sort of software-driven calibration for an IBM personal computer through the data acquisition system of Figs. 4 through 7.
• outputs OUT0 and OUT1 of timer 220 in Fig. 5 generate a time base used via interrupt request line IRQ4 to clock data from a buffer.to D/A converter 300.
• This buffer contains a respiratory waveform which may be a sign wave or any selected waveform as obtained by changing the contents of the buffer.
• Output OUT2 of timer 220 generates a heartbeat pulse as its output. In order to work properly, this pulse must be returned to the ECG trigger through node 400 or directly to interrupt request line IRQ3. If the latter course is chosen, however, node 491 must be disconnected from the output of Schmitt trigger 416. By returning the pulse to the ECG trigger, the computer is informed that the timer is through counting the present RR interval and needs a new interval to be loaded into a timer register of timer 220.
• a display of instantaneous heart rate as provided by an electrocardiograph machine, and as illustrated in Fig. 12, may be converted into an instantaneous heart rate fluctuation spectrum as illustrated in Fig. 14.
• a typical spectrum for a stable patient is illustrated in Fig. 14 while a typical spectrum for an unstable patient is illustrated in Fig. 15.
• Example I and Example II relate respectively to diagnosis and to treatment employing the present invention.
• Parts suitable for use in construction of the apparatus as illustrated in Figs. 4 through 9 may include those as listed in Tables I, II, III and IV.
• Heart rate spectral analysis was applied to the study of congestive heart failure in infants and children.
• Congestive heart failure is characterized by a marked alteration in cardiovascular regulation.
• many cardiovascular functions which are normally monitored in cardiac intensive care units such as: mean heart rate; arterial blood pressure; arterial blood gases; left arterial pressure and right arterial pressure; right atrial, left atrial and pulmonary artery oxygen saturations; the peripheral pulses; peripheral perfusion; and cardiac output
• the usually-monitored cardiovascular function parameters may be within a normal range immediately before a major cardiovascular crisis, such as hypotension or cardiac arrest, inasmuch as the cardiovascular regulatory system maintains these parameters within a normal range up to the point of system failure.
• Real time heart rate spectral analysis was performed on a dedicated personal computer using a 6809E Motorola Microprocessor-Based System.
• a data acquisition system interfaced the computer with a patient monitor, available from Hewlett-Packard, Palo Alto, California, as Model No. 78341.
• the heart rate power spectrum was calculated in continuous 256 second data epochs.
• a QRS synchronization pulse from the patient monitor was used to determine an RR interval sequence.
• An instantaneous heart rate signal was computed from RR interval sequence and the magnitude of the signal was set to the reciprocal of the current interbeat interval.
• the instantaneous heart rate signal was sampled at 4 Hz and the mean heart rate was substracted from the resulting one thousand twenty-four point time series.
• a power spectrum was computed by squaring the absolute value of a Fast Fourier Transform of the one thousand twenty-four point time series. Values for low frequency power (LFP) were computed by integrating the spectrum of between 0.04 and 0.1 Hz.
• LFP low frequency power
• Respiratory frequency power was computed by integrating the heart rate power spectrum over a 0.2 Hz-wide band centered at the mean respiratory frequency. Hard copies of the heart rate time series and power spectrum were printed for each 256 second epochs. Trend graphics for the LFP, the RFP, LFP/RFP ratio, mean heart rate and respiratory rate (hereinafter referred to as the study parameters) were constructed by manually entering data in data files and analyzing the entered data by means of a computer.
• the study parameters for a given patient the LFP, the RFP and the LFP/RFP ratio (hereinafter referred to as the spectral parameters) remain fairly stable.
• Group I included seventeen stable patients whose median age was one month.
• the patients in Group I were without major post-operative complications and did not need prolonged inotropic support.
• the eight patients in Group II suffered cardiac arrest and died.
• the median age for the members of Group II was one month.
• Group III there was a total of four patients each of whom was critically ill at the time of the study but later recovered. Median age of the members of Group III was one month. Of the four members of Group III, one required re-operation, one had intermittent hypotensive episodes, and two had cardiac arrests from which they were successfully resuscitated.
• Group A included data for twenty patients and Group B included data for twelve patients.
• Typical heartrate fluctuation power spectra for Group A and B are respectively illustrated in Figs. 16 and 17.
• Patient profiles for Groups I, II and III are respectively provided in Tables V, VI and VII. These profiles include age, diagnosis and operation.
• TGA Transposition of the Great Arteries
• IVS Ventricular Septal Defect
• PS Pulmonic Stenosis
• HLHS Hypoplastic Left Heart Syndrome
• SV Single Ventricle
• SEV. severe
• COAO Coarctation of the Aorta
• MULT multiple
• VSD Ventricular Septal Defect
• Supra-V. Supravalulvar
• DCRV Double Chamber Right Ventricle
• TOF Tetralogy of Fallot
• AR Aortic Regurgitation
• MR Mitral
• W/IAA is with Interrupted Aortic Arch
• DORV Double Outlet Right Ventricle
• TAPVC Total Anomalous Pulmonary Venous Connections
• CCAVC Complete Common Atrial Ventricular Canal
• S/P is Status Post
• L Left
• BTS Blailock Taussig Shunt
• PA Pulmonary Artery
• ANOM is Anomalous
• B is muscle Bundle
• PULM Pulmonary
• SYS Systemic.
• the discriminate value for the LFP/RFP ratio was two.
• the range of LFP/RFP ratios was 3 to 22 (arithmetic mean 8.77).
• the range of RFPs was 0.01 to 3.13 (arithmetic mean 0.28) and the range of LFPs was 0.09 to 13.88 (arithmetic mean 1.77).
• Group B the range of LFP/RFP ratios was 0.17 to 1.9 (arithmetic mean 0.83), the ratio of RFPs was 0.02 to 0.32 (arithmetic mean 0.1), and the range of LFPs was 0.01 to 0.1 (arithmetic mean 0.5)
• the mean value of the LFP/RFP ratio was greater than two for Group I, the ratio for the stable patients fell below two for brief periods. That which distinguishes the stable from the critical patients is the sustained value for greater than or about one hour of the LFP/RFP ratio for the critical group.
• each heavy dot A represents a geometric mean
• each light line B indicates the standard error of the geometric mean
• each heavy line C represents the standard deviation of the geometric mean.
• each heavy dot A represents an arithmetic mean
• each set of slashes B1 and B2 represents the standard error of the arithmetic mean
• each set of slashes C1 and C2 represents the standard deviation of the arithemetic mean.
• the mean LFP for Group B [0.05 (Beats per minute) 2 ] was less than the mean LFP for Group A [1.77 beats per minute) 2 ], p ⁇ 0.0001. There was no significant difference between the mean RFP between the groups.
• the initial LFP/RFP ratios for the patients with isolated coarctation of the aorta ranged up to 10,000.
• the LFP/RFP ratios observed for this group immediately after an operation to correct the condition were within the range for Group A patients.
• Two patients had LFP/RFP ratios greater than 100 before discharge from the intensive care unit. These ratios were correlated with mild to moderate congestive heart failure. One of these patients died suddenly at approximately 2-1/2 months after the operation. The other two patients remained alive and well.
• LFP/RFP ratio provided the sharpest discrimination between stable and critical patients in these studies, the LFP alone discriminated between Groups A and B, p ⁇ 0.0001. Neither respiratory frequency peak power nor mean heart rate distinguished between Groups A and B. On the other hand, LFP/RFP ratios and LFP levels low levels sustained for greater than or about one hour correlate with the course of the conditions of patients who experienced cardiac arrest or severe hypotensive episodes but later recovered.
• the increase in LFP and in the LFP/RFP ratio for patients with isolated coarctation of the aorta and moderate heart failure may result from an increased activity from the sympathetic mechanism and a decreased activity of the parasympathetic mechanism.
• the decreased level of LFP and of the LFP/RFP ratio found in critical patients may be due to non-responsiveness of the sympathetic mechanism.
• Sympathetic non-responsiveness may be.due to myocardial catecholamaine depletion alone or in combination with the observed down regulation of beta receptors from cardiac tissue in the end stage of heart failure.
• EXAMPLE 2 In patients undergoing operations, shifts in body fluid disposition during surgery may lead to changes in intervascular volume (i.e. a shift of fluid out of a circulatory tree of blood vessels). Accordingly, the availability of the method of diagnosing cardiovascular stress as described in Example 1 may be used to choose among various protocols for treatment or to justify a radical change in medical or surgical treatment.
• an anesthesiologist may non-invasively monitor intravascular volume status.
• the anesthesiologist may increase the amount of fluids administered by way of intravenous injection or may take steps to reverse effects of a particular anesthetic.
• heart rate fluctuation spectral analysis may be done in real time. This capability permits correlation of treatment administered with changes in LFP or LFP/RFP ratios.
• the present invention may also be used for monitoring cardiovascular instability in the following patients in which adjustments in cardiovascular regulation may provide a central key to understanding the efficiency and efficacy of treatment.
• Ambulatory patients with known heart disease in which sudden cardiac death is a common association one example of which would be a patient with a congestive cardiomyopathy who is being treated with vasodialator drugs and for whom the LFP/RFP ratio has changed from a normal baseline level to decreased levels may then 'subsequently be either admitted to the hospital for adjustment of medications and/or observed and monitored in the physician's office while his vasodialator drug dose is increased.
• a patient with renal disease e.g.
• a patient with moderate to severe pulmonary disease resulting in hypoxemia and/or hypercarbia who requires bronchodialator and/or supplementary oxygen and/or mechanical ventilation may be treated by adjustments in bronchodialator drugs, diuretics, and/or ventilator adjustments.
• a premature infant of very low birth weight known to be at risk for intraventricular hemorrhage may, for example, develop a slow intracranial bleed associated with an abrupt increase in LFP, which may alert physicians prior to a brisk bleed thus allowing institution of appropriate changes in medical management to limit substantially known risk factors that may predispose to such an event, or may permit recognition of the presence of unsuspected circumstances that contribute to the bleed.
• neurologic disease such as one in which a patient has sustained a major intracerebral event (e.g.
• a patient may, for example, exhibit a markedly attenuated LFP/RFP ratio, secondary to massively increased parasympathetic activity which would markedly increase RFP, at the expense of LFP, but which may or may not be associated with signs of increased intracranial pressure, and which may be treated by, for example, hyperventillation, rapid diuresis, or burr hole placement.
• a patient with severe systemic infection may exhibit shock secondary to the infection process may, for example, exhibit an elevated LFP/RFP ratio which may then be subsequently used by the physician in managing the shock state by means of pressor agents and infusion of significant volumes of fluid, thus providing the physician an indication of how effectively he is treating the shocked state above and beyond the traditional measurements such as systemic blood pressure and cardiac output.
• a patient with hematologic disease associated with anemia such as Sickle Cell Anemia, exhibits an oscillation in capillary blood flow when severly anemic at the frequency associated with LFP and may exhibit large values for LFP, and for the LFP/RFP ratio may, for example, be treated by blood transfusion which may lead to an expected decrease in LFP, LFP/RFP ratio, and thus enable the physician to monitor by means of heart rate spectral analysis appropriate timing for transfusion therapy.
• a fetus prior to delivery may for example, exhibit a marked attenuation in LFP associated with severe fetal distress, and may thus alert the physician to perform an emergency Caesarean section.
• calibrators according to the present invention may be adjusted to simulate disease states as well as normal conditions. It is also understood that the present invention is not limited to use with patients whose primary disease is of the heart but that modifications may be made for use with such patients.
• Image_wtl IMAGE "Name: ",K,XXXX,”Hosp num:
• Image_wt2 IMAGE "Age: ",K,XXX,”Wt(kg): “,K,XXXX,”Ht(cm): “ ,K,XXXX, "Diag: ",K,XXXX,”Op: ",K 1320 Next_msg: !
• Block_time Pacing_rate*1.024/3600.
• Topi,Top2,Top3,Top4,Botl,Bot2,Bot3,Bot4 500 COM /Editor/ Edit_msg$[80] 510 COM /Subject/ Sub_name$[25],Hos_num$[15],Id_ age$[10],Id_wt$[10],Id_ht$[10 ],Diag$[30], Opera$[45],Halt_pg 520 COM /Io_chart/ Io_time$(8)[10],Iv_intake(8),Fluid_ in(8),In_tot(8),Urine(8 ),Chest(8),Out_ tot(8),Net(8),Io_ptr 530 COM /Lab_chart/ Lab_time$(8)[10],Na(8),K1(8),
• Pco2(8),Bgo3(8),Be(8),Vent_ptr 550 COM /Pres_chart/ Pres_time$(20)[15],Ao_s(20),Ao_ d(20),Ao_m(20),Pa_s(20),Pa_d(20),Pa_m(20), La_m(20),Ra_m(20),Pres_ptr,Pres_in 560 COM /Heart_index/ Heart_ time$(15)[15],Ci(15),Pvri(15),Svri(15),Heart_ptr 570 COM /Drugs/ Drug_time$(40)[20],Drug_ name$(40)[40],Drug_dos$(40)[20],Drug_ptr 590 DIM Io$(5,15)[30],Io_msg$(5,15)[80] 600 DIM Msg_pad$(10)[80]
• the HP 9826/9836 flexible disk has the ! following structure 1070 ! 2 sides, 33 tracks/side, 16 sectors/track, 256 ! bytes/sector
• Fluid_in(Io_ptr) FNLval(Io_msg$(Chart_num,3))

Landscapes

• Health & Medical Sciences (AREA)
• Life Sciences & Earth Sciences (AREA)
• Engineering & Computer Science (AREA)
• Cardiology (AREA)
• Molecular Biology (AREA)
• General Health & Medical Sciences (AREA)
• Veterinary Medicine (AREA)
• Physics & Mathematics (AREA)
• Physiology (AREA)
• Biophysics (AREA)
• Pathology (AREA)
• Biomedical Technology (AREA)
• Heart & Thoracic Surgery (AREA)
• Medical Informatics (AREA)
• Public Health (AREA)
• Surgery (AREA)
• Animal Behavior & Ethology (AREA)
• Psychiatry (AREA)
• Psychology (AREA)
• Signal Processing (AREA)
• Child & Adolescent Psychology (AREA)
• Developmental Disabilities (AREA)
• Hospice & Palliative Care (AREA)
• Computer Vision & Pattern Recognition (AREA)
• Social Psychology (AREA)
• Artificial Intelligence (AREA)
• Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
• Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
• Apparatus For Radiation Diagnosis (AREA)
• Ultra Sonic Daignosis Equipment (AREA)

Applications Claiming Priority (2)

Publications (2)

Family

ID=24983432

Family Applications (1)

Country Status (15)

Families Citing this family (12)

Citations (1)

Family Cites Families (6)

• 1986
  • 1986-05-27 IL IL78931A patent/IL78931A0/xx unknown
  • 1986-05-30 AU AU59654/86A patent/AU5965486A/en not_active Abandoned
  • 1986-05-30 JP JP61503324A patent/JPS63500153A/ja active Pending
  • 1986-05-30 BR BR8606714A patent/BR8606714A/pt unknown
  • 1986-05-30 WO PCT/US1986/001193 patent/WO1986007248A1/en not_active Application Discontinuation
  • 1986-05-30 EP EP19860903945 patent/EP0223846A4/en not_active Withdrawn
  • 1986-06-03 ZA ZA864119A patent/ZA864119B/xx unknown
  • 1986-06-03 ES ES555655A patent/ES8707853A1/es not_active Expired
  • 1986-06-04 GR GR861453A patent/GR861453B/el unknown
  • 1986-06-04 CA CA000510803A patent/CA1257395A/en not_active Expired
  • 1986-06-04 CN CN198686104798A patent/CN86104798A/zh active Pending
• 1987
  • 1987-02-03 DK DK055787A patent/DK55787A/da not_active Application Discontinuation
  • 1987-02-04 FI FI870470A patent/FI870470A0/fi not_active Application Discontinuation
  • 1987-02-04 NO NO870448A patent/NO870448L/no unknown
  • 1987-02-05 KR KR870700101A patent/KR870700317A/ko not_active Application Discontinuation

Patent Citations (1)

Non-Patent Citations (3)

Also Published As

Similar Documents

Legal Events