CA1126343A - Derivation of steady values of blood pressures - Google Patents

Derivation of steady values of blood pressures

Info

Publication number
CA1126343A
CA1126343A CA321,401A CA321401A CA1126343A CA 1126343 A CA1126343 A CA 1126343A CA 321401 A CA321401 A CA 321401A CA 1126343 A CA1126343 A CA 1126343A
Authority
CA
Canada
Prior art keywords
signal
value
input
output
difference
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.)
Expired
Application number
CA321,401A
Other languages
French (fr)
Inventor
Allan P. Sherman
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.)
HP Inc
Original Assignee
Hewlett Packard Co
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 Hewlett Packard Co filed Critical Hewlett Packard Co
Application granted granted Critical
Publication of CA1126343A publication Critical patent/CA1126343A/en
Expired legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/021Measuring pressure in heart or blood vessels
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/90Blood pressure recorder

Abstract

ABSTRACT

Apparatus for deriving an average of data points in which the contribution of each data point to the average is generally inversely related to the difference between its amplitude and the amplitudes of a few previous data points.

Description

Background of the Invention In monitoring the systolic and diastolic blood pressure of a patient by directly measuring the pressures of the maximum and minimum peaks of the blood pressure signal, the presence of equipment noise can cause fluctuations in the peak pressures that are not physiological in origin, and the presence of variations in the pressure of the thoracic cavity due to respiration can cause the peak values to slowly un-dulate. If the peak values are displayed on a strip chart recording or the like, the effects of such changes in their values can be discounted by the careful analysis of a skilled observer, but it is often desired to avoid the need for such analysis by displaying the systolic and diastolic values in numerical form. Fluctuations in the numerically displayed peak values would make it difficult for even a skilled observer to determine the useful systolic and diastolic pressures.
This can be even more difficult if pulmonary artery pressures are being monitored because the amplitude of the undulation in pressure due to respiration can be comparable to the pulmonary artery pressures being measured.
Steady values of the systolic and diastolic blood pressures can be derived by respectively applying the fluc-tuating maximum and minimum peak values to suitable low pass filters, but the systolic and diastolic values derived in this manner are as much affected by peaks of noise as by peaks representing heart action. Furthermore, if the upper limit of the filter is low enough to attenuate the slcw fluctuation in the peak values due to respiration, the response may be too slow to reveal sudden changes in the peak levels that are of physiological significance.
, ~ .

~1~6~343 Brief Discussion of the Invention Although this invention will be described and illus-trated as it would be used in deriving values of systolic and diastolic blood pressures from a blood pressure signal in such manner that the values are relatively free from the effects of artifacts and respiration, it can be used to derive a stable indicator of the level of any signal that changes more slowly than undesired perturbations. While this indicator re-sponds quickly to sustained changes in the nature of the signal, it is relatively insensitive to spurious artifacts. Further~
more, the indicator is an average of the signal in which each constituent portion of the signal affects the average in a generally inverse relationship to the amount of its variation from the immediately preceding constituents of the signal. The immediately preceding constituents can be the portion of the signal occurring within a predetermined prior period, the por-tion of the signal represented by a predetermined number of samples, or an exponential average of the entire previous signal whether it be continuous or represented by samples.
However the maximum and minimum peaks of the blood pressure signal may be detected, the systolic and diastolic values to be displayed are derived from them in accordance with the invention in such manner as to yield values that are steady and physiologically significant. In general, this is accom-plished by separately producing moving averages of the maximum and minimum peak values in such manner as to give greater weight to those peaks having amplitudes closer to the average amplitudes of their predecessors than to those having ampli-tudes that are significantly different from the amplitudes of their predecessors. Accordingly, a large positive peak that 1126~43 results from the detection of the amplitude of a single relatively large spike of noise will have very little, if any, effect on the average value that is displayed as the systolic pressure, and a large negative peak that results from the detection of the amplitude of a single relatively large spike of noise will have very little, if any, effect on the average value that is displayed as the diastolic pressure. In either case, the average values to be displayed will gradually change if there is a number of similar noise spikes occurring in succession. Thus, the displayed values can follow a step change in the signal level. On the other hand, if the amplitudes of the peaks gradually change, as during respiration cycles, the displayed values will be very close to the average of the peaks within the cycle. This is especially important in monitoring pulmonary artery pressures as the variation in the numerically displayed peaks is much reduced. In fact, it can be reduced to a value within the resolution of the display so that respiration has no effect.
According to one aspect of this invention there is provided ~paratus for processing an input signal, comprising a signal input terminal, first means for deriving at its output a first signal representative of an average of at least some of the values of a signal applied to its input, means providing a value at the input of said first means representing the value of the signal at said signaI input terminal, r ~ 3~3 second means for deriving at its output a second signal representative of the absolute value of the relative deviation between the values represented by signals applied to first and second inputs thereof, means for providing a value at the first input of said second means representing the value at said signal input terminal, and means providing a value at the second input of said second means representing the signal value at the output of said first means, third means for deriving at its output a third signal representative of the difference between the values represented by signals coupled to different inputs thereof, a signal output terminal having means for main-taining the value of a signal applied to it, means providing a value at one input of said third means representing the value at said signal input terminal, and means providing a value at the other input of said *ird means representing the value at said signal output terminal, fourth means coupled to the outputs of said second and third means for deriving a :ourth signal at its output representing a portion of the value represented by the third signal, the portion being larger when said second signal becomes smaller and smaller when said second signal becomes larger, control means coupled to said fourth means for changing the value represented by the signal at said signal output by an amount proportional to said fourth signal~ and means for selectively activating said first, 0 second, third and fourth means.
- 4a -According to another aspect of this invention there is provided apparatus ~or processing an input signalrepresenting the variation in blood pressure with time in which there are systolic and diastolic peaks and peaks due to noise as well as undulations due to respiration in such manner as to produce an output signal having a value repre-senting one of the systolic and diastolic peak pressures in which the effects of noise and respiration undulations are reduced, comprising an input to which a signal representing the variation in blood pressure with time may be applied, first means coupled to said input for deriving a first signal approximately representing the peaks of said one of the systolic and diastolic pressures of a blood pressure signal when present in such manner that the effects of the undulations and noise may be present, second means coupled to the output of said first means for deriving a second signal representing an average of some previous peaks represented by said first signal, third means responsive to said first and second signals for dexiving a control signal representative of the absolute value of the deviation between the first and second signals, an output from which an output signal may be derived, fourth means responsive to said first signal and a signal at said output for providing a fourth signal representative of the difference in value between the value represented by the first signal and the value represented - 4b -.~

by the signal at said output, and means for changing the signal at said output by a portion of the value represented by said fourth signal, the said portion increasing as said control signal gets smaller and decreasing as said control signal gets larger.
According to another aspect of this invention there is provided a method of processing an input signal repre-senting variations in systolic or diastolic blood pressures in such manner as to derive an output signal in which the effects of undulations contained in the input signal that are caused by respiration as well as the effects of false values contained in the input signal that are caused by noise are reduced, comprising sampling the amplitude of the input signal at successive instants of time, deriving an average of the values of at least some of the prior samples, deriving a first difference between the value of each sample and the value of the average amplitude of samples prior to it, determining a second difference between the value of each sample and the value of the output signal, and changing the value represented by the output signal in the direction of the input signal by a portion of the second difference that increases as the first difference decreases and decreases as the first difference increases.
According to another aspect of this invention there is provided a method of processing an input signal in such - 4c -r~

6~43 manner as to derive an output signal which is more responsive to a slow change in value of the input signal than it is to a fast change in the value of the input signal, comprising sampling the amplitude of the input signal at successive instants of time, deriving an average of the values of at least some of the prior samples, deriving a first difference between the value of each sample and the value of the average amplitude of 0 samples prior to it, determining a second difference between the value of each sample and the value of the output signal, and changing the value represented by the output signal in the direction of the input signal by a portion of the second difference that increases as the first difference decreases and decreases as the first difference increases.

- 4d -~ ,.

~1~6343 The Drawings FIGURE 1 is a graphical representation of the output signals provided in accordance with this invention from data values that are hypothetically selected for purposes of illus-tration;
- FIGURE 2 includes a graph of an actual pulmonary artery pressure signal having a significant dwell between breaths as well as graphical representations of the values of systolic, mean, and diastolic pressures produced by apparatus constructed in accordance with the invention;
FIGURE 3 includes a graph of an actual pulmonary ar-tery pressure signal of a patient whose breathing is aided by a -4~-:~2~ 43 ventilator as well as graphical representations of the values of systolic, mean and diastolic pressures produced by apparatus constructed in accordance with the invention;
FIGURE 4 includes a graph of an actual pulmonary artery pressure signal in which the respiratory pattern is such that neither systolic nor diastolic peaks have consistent ampli-tudes, as well as a graphical representation of the values of systolic, mean and diastolic pressures produced by apparatus constructed in accordance with this invention;
FIGURE 5 includes a graph of an actual pulmonary artery pressure signal having spaced bursts of noise, as well as graphical representations of the systolic, mean and diastolic pressures produced by apparatus constructed in accordance with this invention;
FIGURE 6 includes a graph of an actual pulmonary artery pressure signal having closely spaced bursts of noise, as well as graphical representations of the systolic, mean and diastolic pressures produced by apparatus constructed in accordance with this invention;
FIGURE 7 is a diagram of the invented circuit and FIGURE 7A includes graphs used in explaining FIGURE 7;
FIGURE 8 is a block diagram of the functional parts of a digital circuit for processing data in accordance with the invention; and 25 . FIGURES 9A, 9B, 9C and 9D are flow charts for the program included in the specification that sets forth instruc-tional steps by which a Hewlett-Packard 21MX series computer can be made to process data in accordance with the invention.
Before considering either analog or digital apparatus for processing a blood pressure signal in accordance with this invention, reference is made to the graphs of FIGURE 1 for a brief discussion of the functions to be performed and an illustration of the advantage to be derived. In the drawing, the X's represent either maximum or minimum peaks of a pulmonary artery pressure and the solid line SSAVG represents a weighted average derived by apparatus constructed in accordance with this invention. Assume that the peaks have been at a constant ampli-tude long enough for the weighted average SSAVG to equal the constant amplitude of the first three peaks, and that the fourth peak has a significantly greater amplitude that is maintained through the peak number 13. The wave SSAVG increases in steps until it reaches the new amplitude. The height of successive steps increases until the one occurring at the time of the peak number 9 and then decreases. The increase is due to the fact that the differences between successive peaks 5 through 9 and an exponentially weighted average of the previous ~eaks are progressively less, so that the peaks are given greater weight in computing the height of the next step. The height of the steps at the peaks lO through 13 gradually decreases because of limits that are arbitrarily selected to limit the amount of correction due to any change. Note in particular the peak num-ber 14 that is at a much greater pressure than the peak number 13 and the fact that SSAVG changes very little in response to it.
Thus, if this peak represented a large spike of noise, it would have very little effect on the value o~ SSAVG that is to be displayed. A similar result would be obtained if the single peak 14 were much less than the peak 13.
Starting with the peak number l9, it has been assumed for purpose of illustration that the peaks vary in a nearly sinusoidal fashion between the pressures of 20.8 and 22.6 mm Hg.

.63~3 Even larger variations would normally result from respiration.
If the resolution of numerical display for the peaks is l mm Hg, it will vary from a reading of 20 to 22 mm Hg. Note, however, that the peak-to-peak amplitude of the wave SSAVG that is S produced during the sinusoidally varying peaks is about one-third of the peak-to-peak amplitude of the peaks, so that the numerical indicator would remain at 22 mm ~g and not fluctuate.
This represents the kind of reaction the invented apparatus has to undulating variations of low frequency that occur, for example, as a result of respiration. Actually, the constant variation in pressure due to respiration is larger than that shown, so that an indicator may exhibit changes in numerical value, but in any event they will be much smaller with the invention than without it.
FIGURES 2 through 6 show actual results obtained by the invention under various conditions. In all of these figures, the blood pressuxe signal is indicated by PDATA; the top horizontal line segments generally indicated by SP are the systolic pressures attained; the intermediate horizontal line
2~ segments generally indicated by MP are the mean pressures;
and the lowest horizontal line segments generally indicated by DP are the diastolic pressures. It will be noted that any changes in level of the line segments occur at regular inter-vals, rather than after each peak as in FIGURE l. This result is obtained by sampling the averages at spaced intervals.
FIGURE 2 illustrates a blood pressure signal PDATA
in which there is a significant dwell between breaths during which both maximum and minimum peaks are sufficiently con-sistent to cause them to be entered into their respective averages with a strong weighting factor. The sporadic peaks 3a~3 on either side of this dwell depart from the average so as to have very little effect. As a result, the values of SP, MP
and DP that are to be indicated by a numerical display vary only slightly. A similar result is attained in FIGURE 3 that illustrates a blood pressure signal obtained from a patient on a ventilator.
When, as in FIGURE 4, the blood pressure signal P~ATA exhibits a respiration pattern such that neither the systolic nor the diastolic peaks have consistent amplitudes, the pressure values of SP, MP and DP tend to be the averages of all the respective peaks.
Intermittent bursts of noise, such as illustrated in FIGURE 5, have little effect on the values of SP, ~P and DP, but a rapid succession of bursts of noise, such as illus-trated in FIGURE 6, does have some effect on these values.
FIGURE 7 is an analog circuit for processing the signals in accordance with this invention. A blood pressure signal PDATA, such as illustrated in FIGURE 7A, is derived from a patient P by coupling the blood pressure at an appro-priate point in his circulatory system to a transducer 2 with a catheter 4. The output signal of the transducer 2 corre-sponds to the signal PDATA and is coupled via an amplifier 6 to means within the dotted rectangle 7 for deriving signals SMAX, each of which corresponds to the amplitude of the most recent maximum peak PS of the signal PDATA, as well as signals DMIN, each of which corresponds to the amplitude of the most recent minimum peak PD. Various means may be used for this purpose, but the particular one illustrated i~ a simplified form of the system disclosed in my U.S. Patent No. 4,223,682, issued September 23, 1980, and entitled "Beat-to Beat Systolic and Diastolic Indicator". For reasons .
that will be explained, the means 7 also provides pulses respec-tively indicating when SMAX or DMIN have attained a value to be processed. These pulses could be provided by other means, such as a QRS detector operating upon the ~CG signal.
The circuits within the dotted rectangle 7 may be described as follows. The signal PDATA at the output of the amplifier 6 is applied to the non-inverting input of a com-parator 8 and to a low pass filter for deriving the approximate mean BMAVG of the signal PDATA. The filter is comprised of a resistor 10 and a capacitor 12 connected in series to ground.
The signal BMAVG appears at their junction 14 and is coupled to the inverting input of the comparator 8. Whenever PDATA is greater than BMAVG, the output of the comparator 8 is high, and when PDATA is less than BMAVG, the output of the comparator 8 is low, as indicated by the waveform C of FIGURE 7A. The out-put of the comparator 8 is connected to one end of a relay coil 15, the other end being connected to a low voltage -V.
When the output of the comparator 8 is high, the coil 15 is energized so as to pull a switch sO from the normal position shown, where it contacts a terminal 15' that is connected to the input of a minimum peak detector 16, to a position where it contacts a terminal 15" that is connected to the input of a maximum peak detector 17. The switch sO is connected to the output of the amplifier 6 so that it applies the signal PDATA
to one peak detector or the other.
The output of the comparator 8 is connected to a positive edge triggered single-shot multivibrator 18 that pro-duces pulses P18 whenever PDATA crosses BMAVG in an upward direction. These pulses are applied to set a sample-and-hold circuit 19 and via a delay 20 to clear the minimum peak 1~ t;3~3 detector 16. Thus, the signals DMIN of FIGURE 2A appear at the output of the sample-and-hold circuit 19.
The output of the comparator 8 is also connected to a negative edge triggered one-shot multivibrator 21 that produces pulses P21 whenever PDATA crosses BMAVG in a downward direction. These pulses are applied to set a sample-and-hold circuit 22 and via a delay 24 to clear the maximum peak detector 17. Thus, the signals SMAX of FIGURE 2A appear at the output of the sample-and-hold circuit 22.
Processing of a maximum peak signal SMAX is accom-plished by the circuits within a dotted line 25 in the following way. The signal SMAX at the output of the sample-and-hold circuit 22 is coupled to an input terminal 26 that, in turn, is connected to the A input of a subtractor 28. The terminal 26 is also connected to means for deriving a voltage FSAVG, fast systolic average, equal to an average of the most recent values of SMAX. In this particular circuit, the means is com-prised of a normally open switch sl connected between the input terminal 26 and a low pass filter comprised of a series combination of a resistor 29 and a capacitor 30. The pulses P21 at the output of the multivibrator 21 are applied to a relay coil 32 that operates the switch 51 When it closes, it effectively samples SMAX and applies the sample to the filter. The values of the resistor 29 and capacitor 30 in combination with the duration of the pulse P21 determine the speed with which the voltage FSAVG at the junction 34 of the capacitor 30 and the resistor 29 follows the changes in the amplitudes of the signal SMAX. Values chosen so as to make the change in voltage equal half the voltage across the resistor 29 at the beginning of each pulse have heen found ~L2~

to work well, but other values can be used. In any case, the value of FSAVG is the average of very recent values of SMAX. The voltage FSAVG at the junction 34 is applied to the B input of the subtractor 28. The effect of each sample of SMAX on the voltage across the capacitor 30 decreases exponentially with the number of succeeding samples, so that theoretically the voltage on the capacitor 30 is the result of all samples that have occurred since the apparatus was turned on, but the most recent sample has the greatest effect, and the effect of each successive previous sample is less.
The output of the subtractor 28 is A-B and therefore equal to the value of SMAX minus the value of FSAVG, which is the difference between the amplitude of the current value of SMAX and the average of the amplitudes of a few recent values of SMAX. Whether the output of the subtractor 28 is positive or negative, it is made to be a positive voltage of the same numerical value by an absolute value circuit 35. The output of the absolute value circuit 35 is a divisor control voltage and is applied via a diode dl to a divisor input Y of a divider 37. A resistor 36 is connected between the Y input and ground. In order to provide a minimum value of positive voltage at the divisor input Y, a diode d2 is connected between a tap 38 on a potentiometer resistor 39 and the divisor input Y. The resistor 39 is connected between a point of positive voltage and ground.

~J63~3 The manner in which the voltage applied to the dividend input X of the divider 37 is derived will be ex-plained at a later point. The output of the divider 37, X/Y, is applied via a voltage controlled current source 40 to an output terminal 42 so as to charge or discharge an output capacitor 44 that is connected between the terminal 42 and ground by an amount proportional to X/Y. The voltage of the terminal 42, SSAVG, is the weighted moving average of the values of SMAX that is to be used for display purposes.
The voltage at the dividend input X of the divider 37 is derived as follows. A connection is made between the output terminal 42 and an input B' of a subtractor 46. A
normally open switch s2 that is operated by the relay coil 32 is connected between the input terminal 26 and the other input A' of the subtractor 46. A resistor 50 is connected between the inputs A' and B'. In between pulses P~, when the relay coil 32 is not energized and the switch s2 is open, the voltage SSAVG at the output terminal 42 is applied to both inputs A' and B' of the subtractor 46, so that its output is zero. During the pulses P~,the switch s2 is closed so as to apply the cur-rent value of SMAX at the input terminal 26 to the A' input of the subtractor ~6. Its output A'-B' during the pulses P~ is therefore equal to the difference between the current maxlmum SMAX and the old value of the moving average SSAVG. This out-put is applied to the dividend input X of the divider 37.
Thus, if a maximum detected peak of PDATA abruptly increases to a much larger value so that a much larger value of SMAX appears at the input terminal 26, the difference be-tween this value of SMAX and the moving average SSAVG is divided by an amount proportional to the difference between the value SMAX and the average amplitude of a small number of previous peaks determined by the time constant of the resistor 29 and the capacitor 30 and t~,e duration of the pulses P2].
In certain cases, this is necessarily modified by the appli-cation of the minimum positive voltage to the input Y of the divider 37 in order to prevent the divisor from being zero.
C Referring again to FIGURE 1, assume that the ~1~ in-dicate values of data peaks having corresponding values of SMAX.
As SSAVG gradually approaches the level of peaks 4 through 9, A-B, or SMAX-F~AVG, gets smaller so that the amount of charge placed on the output capacitor 44 during successive pulses P21 increases. From peak 10 through peak 13, the increase becomes less because the dividend A'-B', or SMAX-SSAVG, gets smaller and the divisor is limited to the value set by potentiometer 39.
It can be seen that the moving average SSAVG does not change in response to the much larger peak number 4 because, even though SM~X-SSAVG is large, the value of A-B, or SMAX-FSAVG, that determines the magnitude of the divisor is also large.
Even though the value of the moving average SSAVG
at the output terminal 42 can change at the time of the pulses P21, the changes to be applied to an indicator 51 can be made to occur at a slower rate by coupling it via a sample-and-hold circuit 52 to the output terminal 42 and triggering it at suitable times with pulses from a timer 54. This accounts for the fact that changes in the amplitudes of the horizontal line segments SP, MP and DP of FIGURES 2 through 6 occur less fre~uently than the peaks in the signal PDATA.
Initialization In order to initialize the circuit, power is applied by means not shown, and a switch s3 is momentarily closed so 6~43 as to energize a relay coil 56 and close normally open switches s4, s5, s6 and S7. The closure of the switch S4 applies the PDATA signal from the amplifier 6 to the capacitor 30; the closure of the switch s5 applies the PDATA signal to the out-put capacitor 44; and the closure of the switches s6 and S7 applies the PDATA signal to the same points in a circuit not shown but contained within a rectangle 58 for deriving the current value of the diastolic pressure in the same manner as the systolic voltage is derived. The pulses Pl8 from the multivibrator 18 and the voltage DMIN at the output of the sample-and-hold circuit 19 are also applied to the circuit 58.
The diastolic circuit contained in the rectangle 58 is assumed to be the same as those portions of the systolic circuit con-tained within the dotted enclosure line 25. Its output may be applied to the indicator 51 via a sample-and-hold circuit 60 that is triggered by pulses from the timer 54.
Operation From the description above, it can be seen that a new current value of SSAVG is developed across the output capacitor 44 in response to each of the pulses P21. The amount by which it is changed depends on the output of the divider 37. This output is equal to the value of the voltage applied to its dividend input X, which is equal to the difference between the value of the current systolic peak SMAX and the old value of SSAVG divided by the greater of the voltages at the tap 38 or the absolute value of SMAX-FSAVG. Thus, (1) SSAVG (new) = SSAVG (old) + SMA_ - SSAVG (old) The potentiometer 39 and its tap 38 insure that F will never be below some given value.
3~3 The value of F is proportional to the absolute value of the difference between the amplitude of the current systolic peak SMAX and the former systolic average FSAVG that is de-rived by the subtractor 28. Therefore, we can say (2) F = ¦SMAX--FSAVG
where K is determined by the gains of the subtractor 28 and the absolute value circuit 35. The value of F is held to a minimum by the potentiometer 39 so as to limit the change in value of SSAVG for each beat. This also avoids dividing by zero. Substituting, (3) SSAVG (new) = SSAVG (old) + K ISMAX-SSAVG (old) where the minimum value of the denominator SMAX-FSAVG is established by the setting of the tap 39 on the potentiometer resistor 38.
Expressing these values in millimeters of mercury, appropriate values have been determined to be:
- ¦SMAX-FSAVG¦minimum Radial Artery Pressure 1/2 _ _ .
Pulmonary Artery Pressure 1/10 1/5 During the time that 52 is closed, the output of the subtractor 46 can be seen to represent an "error signal" or the disagreement between the signals at the input terminal 26 and the output terminal ~2. When s2 is opened, the signal SSAVG is applied by the resistor 50 to the A and B inputs of f~ ~ r- r ~ r ~ the subtractor 46 so as to cause ~ signal to be zero. The ¦SMAX_FSAVG¦ is seen to control the amount o adjustment that the "error signal" is allowed to make to the signal at the output terminal 42, i.e., to the voltage across the capacitor 44.
Stated in a general way, it can be said that the 3~3 subtractor 46 is a means for providing a difference signal representing the difference between the current ~alue of SMAX and the old value of the moving average SSAVG that may be used in producing a numerical display, and that the sub-tractor 28 is a means for providing a signal indicative of the portion of the difference signal that is to be used in changing the old value of the moving average SSAVG to a new value.
The means within the dotted rectangle 7 provides signals SMAX and DMIN as well as pulses P18 and P21 that re-spectively indicate when new va~ues of SMAX and DMIN have been respectively determined, i.e., a pulse P18 occurs just after the search for the value DMIN of a diastslic peak PD has been completed and a pulse P21 occurs just after the search for the value SMAX of a systolic peak PS has been completed. The values of DMIN and SMAX at the outputs of the sample-and-hold circuits 19 and 22 correspond to the amplitudes of these peaks. By effectively sampling SMAX and DMIN with the pulses P21 and P18l the voltage FSAVG on the capacitor 30 is un-affected by the spacing between peaks. If the spacing is constant and known, the pulses could be dispensed with and the low pass filter comprised of the resistor 29 and the capacitor 30 could be connected directly to the terminal 26 and could be made to have a time constant such as to provide the same value of FSAVG~ The time constant could be varied as the spacing varied. Whereas peaks due to systolic and diastolic pressures may have a fairly consistent spacing, it can be seen from FIGURES 5 and 6 that noise peaks do not, so that the use of a filter alone would provide values of FSAVG that are depen-dent on spacing rather than on amplitude alone. For this ~L2~;343 reason, it is preferable that some means be provided for sampling new values of SMAX and DMIN. The diastolic pea~s are processed in an identical fashion by means not shown within the rectangle 58.
Comment In the circuit of FIGURE 7, the pulses P21 are applied to the relay coil 32 so as to close the switch sl during each pulse, thereby obtaining a sample of each new value of the analog stepped signal SMAX. During the same time, the switch s2 is closed so that the subtractor 46 can derive the difference between the sample and the former value of SSAVG at the output terminal 42. Operation in this manner relates the average SSAVG to prior samples of SMAX regard-less of how far apart they may be. If it is desired to relate SSAVG to the time average of SMAX, the switches sl and s2 could remain permanently closed and the relay 32 could be dispensed with, along with the pulses applied thereto.
In either case, one can consider that the low pass filter comprised of the resistor 29 and the capacitor 30, the subtractor 28 and the absolute value circuit 35 is a differ-entiating means for determining the absolute value of the slope of a signal; that the divider 37, the VCCS circuit 40 and the capacitor 44 effectively constitute a low pass filter having a variabl~ time constant determined by the differenti-ating means; and that the subtractor 46 is a means for deter-mining the difference between the current value of the signal and the output of the low pass filter ~ust described. Thus, when the absolute value of the slope of the signal is large, as in the case of sample numbers 3 and 14 of FIGURE 1, the effective time constant of the low pass filter is made large so that the output of the filter is changed very little and, i3~3 conversely, when the slope of the signal is low, as in the case of samples 4 through 13 of FIGURE 1, the effective time constant of the low pass filter is small so that the output of the filter is changed by a greater amount.
In the case considered, the signal that is processed is the variation of the maximum peaks or SMAX rather than the entire blood pressure signal PDATA, but the circuit could be used to continuously process a signal that is expected to have slopes within a given range so as to reduce the effect of por-tions of the signal having slopes that are outside that range.
Description of a Digital Signal Processor Direct processing of the blood pressure signal PDATA
at the output of the amplifier 6 of FIGURE 7 so as to derive maximum and minimum peaks and the production of stable values thereof in accordance with this invention can be carried out by digital logic circuits as well as by the analog circuits of FIGURE 7. Circuits operating in this manner are described in connection with the ~ewlett-Packard 21MX series computer, but it will be apparent to those skilled in the art that other com-puters could be used. It is contemplated that the computer will be an integral part of the monitoring equipment of which this invention forms a part.
FIGURE 8 is a block diagram of functional digital units required for carrying out the invention. Control of the flow of data as well as the operations performed on the data are determined by instructions stored at different locations in a ROM or instruction memory 62. Instructions are fetched from the memory 62 by a memory controller 64 in response to commands from a program counter or P register 66 and applied to a controller/sequencer 68.

3':~

The analog blood pressure signal PDATA at the output of the amplifier 6 of FIGURE 7 is converted into its digital form by an analog-to-digital converter 70 and applied to an input of an input/output processor 72. The converter 70 samples the analog wave PDATA at intervals determined by a clocking device 84, which is distinct from the computer clock, not shown, contained in the controller/sequencer 68. The sequencex con-trols the application of the digitized PDATA to an A register accumulator 74 and a B register accumulator 76 that are part of an arithmetic processor including an arithmetic logic processor 78, and it also may cause a memory controller 80 to take data from one of the accumulators 74 or 76 or the logic processor 78 f ~ IQ <~ ~
C and ~h~e it into a read/write data storage memory 82 or vice-versa. The data storage memory 82 may share the same physical memory with the control instruction memory 62 and both memories may be controlled by the same memory controller, but separate memories and controllers are shown so as to more clearly identify the roles played by each.
All operations involving data manipulation usually take place in the accumulators 74 or 76. A program instruction may, for example, cause the controller/sequencer 68 to request the arithmetic logic processor 78 to perform some operation involving one or both of the accumulators. Whether the opera-tion includes adding, dividing, negating, multiplying, shifting, rotating, Boolean manipulation (AND, OR), or testing, the re-sults stay in an accumulator. If the value in the accumulator is to be stored in the data memory 82, a separate instruction will be given to the memory controller 80 by the controller/
sequencer 68.
Instructions may be taken from the control instruction memory 62 by the memory controller 64 with a sequence other than normal by resetting the program counter 66 with the arithmetic logic processor 78 or with the controller/sequencer 68, e.g., a test compare instruction could result in the con-tents of the program counter or P register 66 being modified so as to cause it to select an alternate branch in the stream of program instructions derived from the control instruction memory 62.
An external clock 84 provides asynchronous timing pulses to the input/output processor 72 causing it to inform the controller/sequencer 68. This results in the program counter 66 being shifted to a timing section of the program that might, for example, request the input/output processor 72 to sample the analog blood pressure signal PDATA with the A-D converter 70 and place it in the data memory 82 via one of the accumulators 74 or 76. After the timing section of the program is completed, the program counter 66 is set back to allow the computer to resume its normal processing.
When all manipulations of data have been completed so as to derive digital signals respectively representative of stable values of maximum or minimum peaks of the blood pressure signal PDATA, the controller/sequencer 68 issues a <. s ~
command signal euasl~`the input/output processor 72 to apply the digital signal to a digital-to-analog converter 86 that converts it to an analog signal corresponding to the peak in question.
Operation of the Digital Processor in Accordance with the Invention . _ The manner in which a Hewlett-Packard 21MX series computer processes signals representing samples of data so as ~63~3 to provide a moving average in which the contribution of each sample is weighted in a generally inverse relationship with respect to the deviation of its amplitude from an average o~
the amplitudes of a few previous samples will now be explained by reference to the flow charts of FI~URES 9A, 9B, 9C and 9D.
A step in the program set forth at the end of the specifica-tion that is pertinent to each block in the flow charts is set forth in nearby parentheses. Steps in the program indicated by an asterisk (*) are not actual computer steps, but are in-cluded for the purpose of explanation. Reference will be made where appropriate to circuit components of FIGURE 7 and to the block diagram of the computer shown in FIGURE 8, but identifi-cation of the figures referred to is not thought to be required in each case. It should be kept in mind that FIGURE 7 only shows the details of the portion of the analog circuit that derives the systolic moving average SSAVG, whereas the program includes steps for deriving both systolic, diastolic, and mean moving averages SSAVG and SDAVG. Statements to the effect that certain functions are done by or in a block means that they are done by the relevant portion of the program and not actually by the block per se.
Initialization The digital form of the blood pressure signal PDATA
at the output of the analog-to-digital converter 70 of FIGURE
8 is transferred into the accumulator 74 as indicated in the block 88 of FIGURE 9A. As only twelve significant bits are handled by the converter 70, the least significant four bits are set to ~ero because their values are not known. Negative blood pressures are not encountered, so that if PDATA is nega-tive, it must be due to an artifact. Therefore, if PDATA i5 ~6343 negative, its value is arbitrarily set to zero.
If the equipment has just been turned on, as is done in the circuit by closing the switch S3, this fact is detected by a test block Dl. In this event, the learning flag is cleared in block 90, which corresponds to opening the switch s3. The equipment is then initialized by carrying out the functions indicated in block 92. The outputs where the average SSAVG of the maximum peaks and the average SDAVG of the minimum peaks appear are set to the value of PDATA, a function performed in the circuit of FIGURE 7 by the closure of the switches S5 and 56 In the program, the mean pressure S~AVG is also set to the value of PDATA, but the derivation of the mean is not shown in the circuit. The peak detectors are reset, a function performed in the circuit by the pulses C from the multivibrators ~ and 24, and finally, the values of FSAVG and FDAVG are set at the value P0 of PDATA at the time to. In the circuit, this latter function is done by closure of the switches S4 and S7. BMAVG is set to PDATA, but this is not shown in the circuit of FIGURE 7 because it is not necessary.
A~ter initialization, or if the answer by the decision block D1 is to the effect that initialization is not re~uired, a decision block D2 determines whether PDATA is above or below the reference, which in this case is BMAVG. If the preferred system of my above-identified patent application is used, the reference will be BMAVG plus or minus offsets that are used to reduce the false peaks caused by certain undesired pertur-bations which may be present in the signal PDATA. If the answer of the decision block D2 is the PDATA is higher than BMAVG, go to point ~ in FIGURE 9B; but if it is lower, go to point in FIGURE 9C.

~1~6343 Systolic Processiny and Search for Minimum Value of DMIN
If D2 indicates that PDATA is below the reference BMAVG, a decision block D3 in FIGURE 9B tests whether this is the first sample below BMAVG. This is done by noting whether the value of DMIN is the unique value, +32767, that is in-serted in the memor~ 82 during initialization by the block 92 or by a bloc~ 100' of FIGURE 9C during normal operation. In the circuit, the first value of PDATA below BMAVG is identified by a pulse P21 provided by the negative edge triggered one-shot multivibrator 21. When D3 indicates that a first sample of PDATA below BMAVG has been received, the following functions are performed.
As indicated in block 94, the unique value of DMIN
in the data storage memory 82 is changed to PDATA. Succeeding samples of PDATA will therefore cause the block D3 to give a negative answer rather than a positive answer. A negative answer causes a decision block D4 and a block 96 to search for a minimum value of DMIN, just as the detector 16 in the circuit.
D4 tests each sample of PDATA to determine if it represents a pressure less than the last value of DMIN. If it is not less, DMIN remains the same, but if it is less, the block 96 updates DMIN in the memory 82. In either case, proceed to point (4) in the flow charts or GMEAN in the program, at which a new value of BMAVG is calculated, as will be explained in con-nection with FIGURE 9D.
An answer in D3 indicating that a first sample of PDATA has been received that is below the threshold also causes the functions indicated in blocks 98 and 100 to be performed so as to derive the weighted averages SSAVG and SMAVG
that are-desired. In the program, a more slowly changing mean ~263~3 of PDATA, termed SMAVG, i5 derived, although this is not represented in the circuit.
In block 98, the smoothness factor F of PDATA is calculated by subtracting the old value of FSAVG from the current value of SMAX and dividing by a constant K. This is done in the circuit by a subtractor 28 and the absolute value circuit 35. The division by a constant K, which is shown to be 16 in the block 98, is accomplished by the internal gains of the subtractor 28 and the absolute value circuit 35. The minimum value of F is established by the voltage supplied to the Y input of the divider 37 via the diode d2.
Derivation of a new value of FSAVG is done by adding the current value of SMAX to the old value of FSAVG and dividing by a constant indicated as being 2. This function is performed by the filter 29, 30 and the closure of sl.
The desired weighted average of the values of SMAX
is calculated by subtracting the old value of SSAVG from the present value of SMAX, dividing by F, and adding the result to the old value of SSAVG. In the circuit, the difference SMAX-SSAVG is derived by the subtractor 46 and the division by F is done in the divider 37. The addition of the results of this division to the old value of SSAVG is done by charging or discharging the capacitor 44. As explained, the value of SMAVG calculated in the block 98 i5 not performed in the circuit, but represents an averaging of the signal BMAVG with a time constant one-half that used for SSAVG. BMAVG is also updated by a similar process in the diastolic circuit 58, or as described in FIGURE 9C. Since it is updated twice as fre-quently as either SSAVG or SDAVG, it will have a net time constant of update which is comparable to both.

3~3 Diastolic Process~g and Search for Maximum Value of SMAX
If the test in decision block D2 of FIGURE 9A indi-cates that PDATA is above the threshold, the functions indi-cated in FIGURE 9C are performed, but because they are the same, except for the values used, the blocks in FIGURE 9C correspond-ing to the blocks in FIGURE 9B have the same designations, primed. The values used are indicated within the bIocks. The step numbers in the program are different. These functions are carried out by diastolic circuits that are assumed to be within the rectangle 58, but not shown.
Updating BMAVG
In FIGURE 9D, a new value of BMAVG is continuously calculated by dividing the difference between PDATA and the old value of BMAVG by a constant, which in the particular ~ example is 16. The result ~ then added to the old value of BMAVG. In the circuit, this function is performed by the low pass filter comprised of the resistor 10 and the capacitor 12.
Example Calculation That the hardware of the computer is controlled in a manner that is fully equivalent to the operation of the hardware represented in the schematic diagram of FIGURE 7 will now be illustrated by an example. Reference will be made to the block diagram of the computer illustrated in FIGU~E 8, when appropriate.
Let us consider the resistor 10 and capacitor 12 that filter the input value PDATA to derive its approximate mean BMAVG. The equivalent function is performed by the computer program steps 809 through 819 that occur at the label GMEAN in the program and on the flow chart in FIGURE 9D. The math~matical basis of this computation is described in lines 800 to 807 of ;34;~
the program. It may be shown that if this computation i~
performed at a regular interval, with a value of P~ATA repre-senting the input to the RC filter ~0, 12, that the value BMAVG derived by the computer will behave identically as a function of time to the voltage BMAVG on the capacitor 12.
This computation is implemented in the computer as follows. The instruction on line 809 causes the memory con-troller 80 to retrieve the numeric value representing BMAVG
from the data storage memory 82 and enter this value into the A register 74. The next instruction on line 810 causes the contents of the A register 74 to be negated. The instruction on line 811 causes the value PDATA to be retrieved from the data memory 82 and added to the contents of the A register 74, the result, PDATA-BMAVG, remaining in that register. The instructions on lines 812 and 813 cause the contents of the A
register to be shifted to the right four places. This results in the binary value, PDATA-BMAVG, being divided by 16. This operation is performed by the arithmetic and logic processor 78 of the computer and the results remain in the A register.
The next two instructions are provided to prevent deadband when the divisor has a value lying between 0 and 15, and is therefore not represented in the RC equivalent. In line 818, the instruction causes the value BMAVG to be fetched from memory 82 by the memory controller 80 and added to the A register~
accumulator. The last instruction causes the contents of the A register to be stored back into memory 82, replacing and updating the original value of BMAVG.
This process takes place for each new sample of PDATA
which is acquired by the analog-to-digital converter at an arbitrary rate. A rate of 50 samples per second has been found to be suitable.

1~L;;~6343 Program To ~e Used with a Hewlett-Packard 21MX Series ComPuter o~S5 HEO pRocess NEXT OAT~ ~ALUE

o~7 ~ .
015h ~ GET NExT O~T~ PoIuT FROn PRCSSURE ~UYFeR ..

U161 HEXT LDX PPRC~
~t62 O17~ L~A O~I GET DATA
OS75 lh8ERT OV~RRUN TEST HERE t IY NEEOED ) 0175 ~ . ~S
0178 NOP SUHP RCTRIEVAL COU~T . -.
OL79 ~ . -~lao ~
0~0~ AND ~IJ7760 STRIP 4 LSS~S
0182 ssA
01PJ CLA DONIT ALLo~ NEtATIVE
oI84 St~ POATA

o ~
058~ ~ CHECY FOR LE~RNIN6 ~T 9ET
0~89 ~ IF 9ET USE CURRENT DATA FOR ST~TUS OIjpL~y 0169 STARTSNG VALUES, 0~90 ~9~ ~
0192 LOA ~EARN
al9~ CLf 0 194 ANO PCTPT~S PC~
U~95 S2A~RSS

0~97 0l9a .Cn~ ! ' ~99 ~NO PCTPT, S CLEA~ LE~RNING BIT
u2a0 sr~ PCTPT~I
020I ST~ O

0203 LDA PDAT~
0204 ST~ SOAVG
0205 STA S~AVC
~200 STA SSAVC
02~5 LDA HIN4 02JS sr~ S~AXSET TO ENsuRe 02~7 LOA H~X~TH~T ~ ~EAT
OZ~8 STA O~IN ~I~L BE FùUND UEXT TIHE
~2~9 0240 LD~ PD~TA CURRENT ~ALuE
0241 aTA ~AVG SET StATUs VALUES
0242 STA fD~VG
024~ ST~ fS~VO
~244 0243 NDLRN CLb . SY9TO~IC OR OIASTO~IC SEARCH~
0249 SXD~ LD~ PDATA
0250 AD~ aH~vG
u2a~ JSS CHPAU A cnNTAINs PD~T~
0252 JHP srsTs ABOVE THResHoLo oz5~ JHP a~ss 6ELO~ TH ~G
~2~4 .TH IS HN~VG
0236 SPC 2 SE~RC~ FOR NEXT SYSToLIe VALUE
02~ NE~ OAT~ V~LUE IS IN ~ REC ~PDATA~ -u2Sg SYSTS LD~ SHAX PREVIOU9 HAX
0260 CP~ n~
u2 ~ ~nr FTSYS FIRST TIne IN~
0263 JSB cHp~a NE~ HAX ~ '~
oZ6~ ~ , 026~ ~T~ s~X N~
0260 JHP 6hEAN

11.~.3~3~3 ~267 HEO PROCES~INt OF hE~ D;ASTOLIC vil.UES
D268 ~ UST TIne ON s~tsro~C SIOE nF FENDE
0269 CA~L 7~IS TNE START 0~ ~ NE~ 8EAT
02~ -~2~ l ~
~272 FTSYS EOU ~ FlR9T TI~E S~S~OLIC
~2~
027- ~ SET SY9TOL~C TO CUR~ENT V~LUE
Q27~ 3TA 58AX
~27 02d~ ~
02~Z ~ Ev~LuATe 9~00TH~E59 oF CURRENT ~EAT CO~PARED TO tAST AVG.
02~J ~ , . .-~204 ~OB FDAVG DI~S AYG
02~3 tHH~lN~
020~ ~D~ DHIN D~AST~FAST AVG
028~ ~
02~S SS~ ~a~ V~LUe 02~9 ~ CHO~NG
6290 ~
0291 ~ SC~LE SHooTHNEs9 tNOT ~.ESS THAN 2~ -0293 ~R 4 O~V~OE 8Y 16 D294 CLA~INA A-l 029~ JS~ CHPAO
02~0 LDP ~D2 R~S SH~LLER
0297 ST~ SHOTH

0312 ~ UPDATE FAST AVERACE
~313 O~IS ~ TNIS RDuTINE UPDATES A RUNN~U~ ~vERAGE 9Y THE ~ETHDD
0~16 0~RAVG tl,O ~ fACTOR~ ~ R~VC ~ FACTOR ~ hE~VAL
031~ ~ -0319 ~tRAVO ~ RAVC ~ fACTOR ~tNENvAL - h~VG~
~20 0321 ~THIS HAS AN EXPIlNENTIAL TInE RESPONSE SIHILAR TO N R-C FSLTeR
a322 ~323 ~ .
03Z4 ~USE fACTOR ~ 1~2 CIVES T~HE CDNsTANT Of ~ ~EAT8 032~ ~
0326 ~ANO SIHPEIFY Tn RAVC~ tRAVG~Ne~VAL~ ~ 2 ~3zaLD8 fDAVC
0329 ADU OnSN D~AST
0350 CLE~fRH OIV oY 2 Assu~e DATA ~s PossTsve ~352 03~ ~
033~ ~ UPDAte SLo~ AVERAGE
0~5~ ~
03~7 ~ T~tlE CONSTANT IS FASTER FOR SHOOThER O~TA .
0338 ~ ~ TH~S C~VES FASTER RESPONSE TO A STEP CHANGL ,, ~5 ~341 0343 ~ THIS ROUTINE UPDATES A RUNNSNG AYER~Ge 8Y THE ~ETHOD
034~ ~
0345 RAVC tl~o ~ FACTOR~ ~ RAVG ~ ~ACTOR ~ NE~VAL
~46 034~ ~ tRAVC ! RAVG ~ FACTOR ~tNE~YAL RAVC) 03~a' ~
0349 ~ THI9 HAS AN CXPONENTSAL TIHE RESPONSE St~a~AR TO iN H~C FI~TER
0:~S~
Oi~
0352 ~ UOE ~ACTOR ~ ~/S~OTH~ GIVES T~HE CONSTANT 0~ 2 SEATS
a353 ~ ~OR SHOOTHe9T DATA
03S~ ~98X sTP REapoNsE tN 6 ~EAT~
0355 ~ SYoTH~lB CtVEa 98X 9TEP ReSPONSE 2N 62 8EATS
0~56 ~357 LDA SDAVC
~3$9 CHA~SNA . -0359 ADA D~N DSAST~SLO~ ~VC
~6~ C~
03~1 SS~
O~b2 cca EXTeND NE~ SIGN
03~3 DtY SHOTH DtAST~AVGXS~oT~
036~ ADA SDAVC TRUNC~TED, NoT ROUNDED
03~5 ST~ 5DAVC
~6 ~1~tj3~3 ~3 03~B
i309 ~SJ~tLA~LY, UPOATE ~EAN SLO~ AVERAGE
~70 ~3~2 0~73 LDA sHAvc 0304 C~A~NA
~3~5 ADA B~AVC ~EAN RUNN~NG AVC ~ ~EAN SLO~ AVO
0376 CL~
0~77 SS~
037~ tCb EXTEND NEÇ SICN
fl379 .. DXV S~OTM- NE~ ~ AVC75~0T~
03~ ~
03~ ARS DIV ~Y 2 03a2 ~
0~ ADA 9~AVG TRUNCATEC~ NOT ROUNDED -~B4 STA S~AVG
~3BS
~3~6 05Se ~ SET rl4 FOh NEXT tJ~E

0520 LOA ~AX
0521 8TA D~JN
fl522 J~P C~EAN
052~ hEU NEXT P2ASTOLIC VA~UE SEARCH
~52 - 0526 D~A55 LD~ D~IN
0527 CPU n~x~ FIRST TI~E7 ~52B J~p FTOIA YES

DS30 JSO C~FAd NE~ ~tN7 0532 JMP G~E~ NP
~53~ ~
0534 5T~ DMJN ~ES
Q53S J~P 6~EAN
~5J6 05~9 ~ FS~ST TS~e SNtO DSA5TOL~C REGION
0540 ~ TAKE PRE~IOUS ~AX AS SYSTO~SC FOR THIs d~AT
OS41 ~ RESET oFF9eT
0542 ~ RESET FOR hEXT T~E
~54 ~545 0~46 OS47 ETD~A EOU ~ F2RgT TtME ntASTOLtC
05~ ~
OS4~ ~ SET DI'ASTOL~C TO CURREhT VALUe 0550 STA O~I~

05~6 ~ . _ OS57 ~ EVALUATE 8~00THNESS OF CURRENT ~EAT Co~pARED TO FA5T AVG, ~ . -055~ ~ ~
055g LO~ FSAVC SYST AVC
0560 C~SN~
05bl AD~ SHAX SYST~FAST AVC
OS62 ~ ' 05~J 550 ~08 VALUE
OS~4 C~

0566 ~ SCALE S~OnTHNESS tNOT LESS tHAN 2 ~S67 OSO~ LSR A DJVIDE 8Y 16 OS69 CLA~INA A~
~57D JSd C~PAd OS71 LD~ n2 ~AS S~ALLER
0572 9Tb S~laTH
05?3 o5~7 ~ UPOATE FAST ~VERAGE
~5 ~509 0590 ~ ThlS ROUTINE UPDATES A RUNNING ~VER~GE ~Y ThE ~ET~OD
059~
05~2 ~ RAVC ~ tl,fl ~ fACTOR~ ~ RAVG ~ FACTOh ~ NE~V~L

059a ~ ~hAVC RAVG ~ FAtToR ~tNE~VAL ~ RAYG~
~95 0596 ~ Th~S ~AS AN XPONENT1A~ TI~E RESPONSE SI~ILA~ YO A~ R~C F~TER
~59 ~S~

~1~6~43 0599 a USE FACTOR 1t2 GIVES TIHE CONSTANT OF 2 OE~TS
~6~0 ~ .
06al ~ ~ND SI~PLIFY TO RAVG^ ~RAVG~hE~VAL) 7 2 ~000~ ADD SnAX SY9T .
06HS CLE~ERO OIV Dr 2 ASSUME OATA 2S POSTTIVE
06~o~ STE~ ESAVG
~!16~7 06~8 SPC 2 ~6~9 1 ~ UPDATE sLO~ AVENAC~

0612. a .-TtHE CONSTANT 25 ~sTER FOR SMOOTHER O~T~
0613 ~ TH~S tlVEs F~STER RESPONSE TO A SrEP CH~hCE
D6~4 à

~6~

06~5 ~ THIS ROUTINE UPD~TES A RUNNINC AVER~CE S~ THE ~ETHOO
06~9 ~ ~
0620 ~ RAVC ~ ~l.U r EACTOR) ~ h~VC ~ F~CTOR ~ NE~V~L
~62t 0622 ~ tRAVG ~ RAVG ~ fACTOR tNE~V~L ~ RAVG~

0624 ~ lHIS HAS AN E~POhENTl~L TIME RESPDNSE S~Hl~AR TO ~N R-C FI~TER

0626 ~ -D627 ~ USE ~ACTO~ ~ 1tSHOTH, GIVES S~E CONSTANT Of 2 SEATS
0625 ~ FOR SM002HES7 DiTA
~06Z9 ~ t9~X STEP REspoNsE IN 6 5EAT9~
053~- ~ S~DTn~l6.Gl~S 9ax SIEP HESPnNsE 2N b2-CEATj 063~ ~
0632 LnA SSAVG
0633 CMA~2NA
0634 ADA SnAx SYST~SLo~ AVG
06~S . CL~
~6J~ SS~

D63S D~V SMOT~ sYsT-AVGtsnnTH
0639 AD~ 55AVG TRUNC~TED~ NOT RDuhDED
004~ 5T~ SSAVC
~6~

0044 ~S2M2LARLy~ UPDATE nEAN SLO~ ~VERACE
a64s 06~b 064~ ~
064S LDA snAvG
~649 CHA,INA
oobO AD~ ~MAYG nE~N RUNNING AVC MEAN SLow AVC
065l CLb . 0652 SS~
oas3 CC6 EXTENO NEG 51tN
~65~ DIV SnOTH hE~ ~ ~VC75nu~H

06b6 ARS D2V ~Y 2 TO SLO~ AVG
065~ ~
06ba ADA SHAVC TRUNC~TED~ NoT ROUNOEO

~660 ~792 SET MAX ED~ hEXT T~ME
D794 LOA 111Nq 0795 . STA SMAX
~7~16 ~ .

~30- -~ZG343 a798 0~99 HED . BEAT TO SEAT HEA~ C~LCULATlOh u8D~ ~ THIS ROUTSNE UPDATES A NU~NINC AVER~GE 6~ THE HETHoD
DBHI 11 , -oB02 ~ ffAVG ~ t~O u FACIOR~ ~ RAVC FACTOR ~ Ne~VA~

oaD4 ~ ~RAVt ~ RAVG O FACTOR ~HE~VA~ r RAYG~

0806 ~ USE f~CTOR ~ i716 GIVES TIHE COhSTAHT OF 2n 5A~PLE9 OB~ ~ 20 ~ 2~ HSEC 0,4 SEt ,0AO9 C~EAH LDA BHAVG OLD AY~
O~IU CHA~lhA
~Blt ~DA PDATA N~ VALUE
Da 12 ~RS~ARS
~3 ARS~ARS DIVIDE BY 16 D a ~
OB15 SSA~RSS DITHER FOR VALUE5 BE~EEH O AuD 15 Dffl6 INA .
0~
oal6 ADA 80AVG

Y~82V) ~ .
oa2~ JMP HEXT NEXT D~TA ~ALuE
fl~22 SPC 4 ~a23 ~i~124 ~
Oa2g HIN~ DEC ~32768 B026 nAX~ DEC 3?.767 OB27 l~oUT DEC ~250 SET FoR S SECOHDS AT SO H2 ~lA2 ~029 ~ ~

~83 ~9~2 oa33 ~ED C~P~ aOUT~NE
oa34 ~oa3s OH3b 0833 ~ TH~S ROUTINE REtURH9 TU P~1 IF A ~ OR > 8 0630 ~ THi9 RouTlNE RETURhS TO P~2 IF A ~ a Ofl40 0041 LHPA~ NOP
oa 42 C~IN0 NoTEl L043 AD3 A E AND O RECISTER9 ~05T1 0044 CHB~95~,1NU~RSS
0045 I9Z CHPA~
~a46 ~D~ ~.
0047 J~P cr.P~
os4a

Claims (10)

1. Apparatus for processing an input signal, comprising a signal input terminal, first means for deriving at its output a first signal representative of an average of at least some of the values of a signal applied to its input, means providing a value at the input of said first means representing the value of the signal at said signal input terminal, second means for deriving at its output a second signal representative of the absolute value of the relative deviation between the values represented by signals applied to first and second inputs thereof, means for providing a value at the first input of said second means representing the value at said signal input terminal, and means providing a value at the second input of said second means representing the signal value at the output of said first means, third means for deriving at its output a third signal representative of the difference between the values represented by signals coupled to different inputs thereof, a signal output terminal having means for main-taining the value of a signal applied to it, means providing a value at one input of said third means representing the value at said signal input terminal, and means providing a value at the other input of said third means representing the value at said signal output terminal, fourth means coupled to the outputs of said second and third means for deriving a fourth signal at its output representing a portion of the value represented by the third signal, the portion being larger when said second signal becomes smaller and smaller when said second signal becomes larger, control means coupled to said fourth means for changing the value represented by the signal at said signal output by an amount proportional to said fourth signal, and means for selectively activating said first, second, third and fourth means.
2. Apparatus as set forth in Claim 1 having an information signal input and a peak detecting means coupled between said information signal input and said signal input.
3. Apparatus as set forth in Claim 2 wherein means are provided for producing a signal indicative of the occurrence of a new peak value in the signal applied to said signal input, and means applying said latter signal to said activating means whereby the value represented by the signal at said signal output may be changed in response to each new peak value applied to said signal input independent of time.
4. Apparatus as set forth in Claim 1 wherein means are provided for initializing the apparatus by setting the initial value represented by the signal at said signal output and the initial value represented by said first signal at selected values.
5. Apparatus as set forth in Claim 2 wherein means are provided for applying a signal representative of blood pressure variations with time to said information signal input and wherein said peak detecting means provides a signal representative of the systolic peaks of blood pressure in such manner as to include the effects of respiratory undulations and false peaks due to artifacts and noise, when present.
6. Apparatus as set forth in Claim 2 wherein means are provided for applying a signal representative of blood pressure variations with time to said information signal input and wherein said peak detecting means provide a signal representative of the diastolic peaks of blood pressure in such manner as to include the effects of respiratory undulations and false peaks due to artifacts and noise, when present.
7. An apparatus as set forth in Claim 2 wherein said first means derives a first signal that represents an exponential average of the peaks of the signal at said signal input.
8. Apparatus for processing an input signal representing the variation in blood pressure with time in which there are systolic and diastolic peaks and peaks due to noise as well as undulations due to respiration in such manner as to produce an output signal having a value repre-senting one of the systolic and diastolic peak pressures in which the effects of noise and respiration undulations are reduced, comprising an input to which a signal representing the variation in blood pressure with time may be applied, first means coupled to said input for deriving a first signal approximately representing the peaks of said one of the systolic and diastolic pressures of a blood pressure signal when present in such manner that the effects of the undulations and noise may be present, second means coupled to the output of said first means for deriving a second signal representing an average of some previous peaks represented by said first signal, third means responsive to said first and second signals for deriving a control signal representative of the absolute value of the deviation between the first and second signals, an output from which an output signal may be derived, fourth means responsive to said first signal and a signal at said output for providing a fourth signal representative of the difference in value between the value represented by the first signal and the value represented by the signal at said output, and means for changing the signal at said output by a portion of the value represented by said fourth signal, the said portion increasing as said control signal gets smaller and decreasing as said control signal gets larger.
17. A method of processing an input signal repre-senting variations in systolic or diastolic blood pressures in such manner as to derive an output signal in which the effects of undulations contained in the input signal that are caused by respiration as well as the effects of false values contained in the input signal that are caused by noise are reduced, comprising sampling the amplitude of the input signal at successive instants of time, deriving an average of the values of at least some of the prior samples, deriving a first difference between the value of each sample and the value of the average amplitude of samples prior to it, determining a second difference between the value of each sample and the value of the output signal, and changing the value represented by the output signal in the direction of the input signal by a portion of the second difference that increases as the first difference decreases and decreases as the first difference increases.
10. A method of processing an input signal in such manner as to derive an output signal which is more responsive to a slow change in value of the input signal than it is to a fast change in the value of the input signal, comprising sampling the amplitude of the input signal at successive instants of time, deriving an average of the values of at least some of the prior samples, deriving a first difference between the value of each sample and the value of the average amplitude of samples prior to it, determining a second difference between the value of each sample and the value of the output signal, and changing the value represented by the output signal in the direction of the input signal by a portion of the second difference that increases as the first difference decreases and decreases as the first difference increases.
CA321,401A 1978-04-10 1979-02-13 Derivation of steady values of blood pressures Expired CA1126343A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US894,679 1978-04-10
US05/894,679 US4190886A (en) 1978-04-10 1978-04-10 Derivation of steady values of blood pressures

Publications (1)

Publication Number Publication Date
CA1126343A true CA1126343A (en) 1982-06-22

Family

ID=25403393

Family Applications (1)

Application Number Title Priority Date Filing Date
CA321,401A Expired CA1126343A (en) 1978-04-10 1979-02-13 Derivation of steady values of blood pressures

Country Status (4)

Country Link
US (1) US4190886A (en)
JP (1) JPS5510981A (en)
CA (1) CA1126343A (en)
FR (1) FR2423009A1 (en)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4273136A (en) * 1977-11-24 1981-06-16 Sharp Kabushiki Kaisha Electronic sphygmomanometer
US4338949A (en) * 1979-08-06 1982-07-13 Croslin Michael E Method for performing non-invasive blood pressure and pulse rate measurements
US4331154A (en) * 1979-10-15 1982-05-25 Tech Engineering & Design Blood pressure and heart rate measuring watch
US4667680A (en) * 1983-11-14 1987-05-26 Hewlett-Packard Company Apparatus and method for reduction in respiration artifact in pulmonary artery pressure measurement
US4885711A (en) * 1984-03-14 1989-12-05 Neff Dennis B Data processing
US4671290A (en) * 1985-01-15 1987-06-09 Richards Medical Company Automatic tourniquet
CA1263186A (en) * 1985-07-05 1989-11-21 Maynard Ramsey, Iii Automated mean arterial blood pressure monitor with data enhancement
JPH06110B2 (en) * 1985-08-27 1994-01-05 オムロン株式会社 Electronic blood pressure monitor
US4799492A (en) * 1986-06-17 1989-01-24 Spacelabs, Inc. Method and apparatus for indirect blood pressure measurement
JPS6458240A (en) * 1987-08-31 1989-03-06 Canon Kk Blood flowmeter
DE4239215A1 (en) * 1992-11-21 1994-05-26 Wandel & Goltermann Averaging and data compression methods for measuring instruments
US5568815A (en) * 1994-11-21 1996-10-29 Becton Dickinson And Company Self-powered interface circuit for use with a transducer sensor
US5752919A (en) * 1996-12-17 1998-05-19 Johnson & Johnson Medical, Inc. Mitigation of respiratory artifact in blood pressure signal using line segment smoothing
US7909768B1 (en) * 2004-07-19 2011-03-22 Pacesetter, Inc. Reducing data acquisition, power and processing for hemodynamic signal sampling
US7324848B1 (en) 2004-07-19 2008-01-29 Pacesetter, Inc. Reducing data acquisition, power and processing for photoplethysmography and other applications
US7922664B2 (en) * 2007-04-05 2011-04-12 Coherence Llc Method and system for improving physiologic status and health via assessment of the dynamic respiratory arterial pressure wave using plethysmographic technique
US7920913B1 (en) 2007-06-29 2011-04-05 Pacesetter, Inc. Systems and methods for increasing implantable sensor accuracy
US9968304B2 (en) * 2011-10-07 2018-05-15 Edwards Lifesciences Corporation Detecting a vasoactive agent in the bloodstream
JP2016055093A (en) * 2014-09-12 2016-04-21 セイコーエプソン株式会社 Blood pressure measurement apparatus and blood pressure measurement method
KR102312667B1 (en) 2019-04-15 2021-10-15 삼성전자주식회사 Refrigerator
WO2020213861A1 (en) 2019-04-15 2020-10-22 Samsung Electronics Co., Ltd. Refrigerator
EP3904797B1 (en) 2019-08-20 2022-10-19 Samsung Electronics Co., Ltd. Refrigerator
KR20210125379A (en) 2020-04-08 2021-10-18 삼성전자주식회사 Refrigerator

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3506813A (en) * 1966-06-13 1970-04-14 Hewlett Packard Co Signal-to-noise ratio enhancement methods and means
US3484591A (en) * 1966-07-18 1969-12-16 Hewlett Packard Co Extended bandwidth signal-to-noise ratio enhancement methods and means
US3648688A (en) * 1970-07-20 1972-03-14 Human Factors Research Inc Pulse jitter measurement
US3771167A (en) * 1971-02-02 1973-11-06 Leeds & Northrup Co Method for calculating the average value of a noise corrupted signal
US3809070A (en) * 1971-07-01 1974-05-07 Doll Research Non-invasive electromagnetic bloodflow measuring system with rejection of noises
JPS526039B2 (en) * 1971-12-17 1977-02-18
JPS5342499B2 (en) * 1973-10-09 1978-11-11
US3996926A (en) * 1974-07-09 1976-12-14 Medtronic, Inc. Blood pressure monitoring system
US4009709A (en) * 1975-05-15 1977-03-01 American Optical Corporation Apparatus and process for determining systolic pressure

Also Published As

Publication number Publication date
JPS5510981A (en) 1980-01-25
JPS6247536B2 (en) 1987-10-08
FR2423009B1 (en) 1985-03-15
US4190886A (en) 1980-02-26
FR2423009A1 (en) 1979-11-09

Similar Documents

Publication Publication Date Title
CA1126343A (en) Derivation of steady values of blood pressures
US4408614A (en) Blood pressure measurement with Korotkov sound artifact information detection and rejection
US5253648A (en) Method and apparatus for excluding artifacts from automatic blood pressure measurements
US4203451A (en) Cardiovascular analysis, method and apparatus
CN108024740B (en) Blood pressure measuring method, blood pressure measuring device and terminal
EP2942007B1 (en) Detection of parameters in cardiac output related waveforms
US6083171A (en) Blood pressure monitoring apparatus
US4616659A (en) Heart rate detection utilizing autoregressive analysis
Lin et al. Reduction of interference in oscillometric arterial blood pressure measurement using fuzzy logic
EP2393423B1 (en) Calculating cardiovascular parameters
EP0310349A2 (en) Fetal monitoring during labour
EP3571984B1 (en) Systems for creating a model for use in detecting a peripheral arterial pressure decoupling in a subject and for detecting peripheral arterial pressure decoupling using arterial pressure waveform data and said model
JPS6242612B2 (en)
CN103156583A (en) Apparatus for evaluating vascular endothelial function
CN101138493A (en) Method and system utilizing spo
JPS5937107B2 (en) electrocardiograph computer
US4456959A (en) Period measurement system
JP3415796B2 (en) Periodic biological information measurement device
EP2613691B1 (en) Method, apparatus and program for the automatic processing of blood pressure signals
EP3785617A1 (en) Method and device for calibrating blood pressure measurements based on measured pulse arrival time
Bjurstdt et al. Dynamics of arterial oxygen tension in response to sinusoidal work load in man
US20220202298A1 (en) Unsupervised real-time classification for arterial blood pressure signals
EP0150176B1 (en) Blood pressure measurement with korotkov sound artifact information detection and rejection
JP3491076B2 (en) Blood pressure measurement device
GB2032623A (en) Derivation of Steady Values of Blood Pressures

Legal Events

Date Code Title Description
MKEX Expiry