# Oscillometric measurement of systolic and diastolic blood pressures validated in a physiologic mathematical model

- Charles F Babbs
^{1}Email author

**11**:56

**DOI: **10.1186/1475-925X-11-56

© Babbs; licensee BioMed Central Ltd. 2012

**Received: **28 June 2012

**Accepted: **3 August 2012

**Published: **22 August 2012

## Abstract

### Background

The oscillometric method of measuring blood pressure with an automated cuff yields valid estimates of mean pressure but questionable estimates of systolic and diastolic pressures. Existing algorithms are sensitive to differences in pulse pressure and artery stiffness. Some are closely guarded trade secrets. Accurate extraction of systolic and diastolic pressures from the envelope of cuff pressure oscillations remains an open problem in biomedical engineering.

### Methods

A new analysis of relevant anatomy, physiology and physics reveals the mechanisms underlying the production of cuff pressure oscillations as well as a way to extract systolic and diastolic pressures from the envelope of oscillations in any individual subject. Stiffness characteristics of the compressed artery segment can be extracted from the envelope shape to create an individualized mathematical model. The model is tested with a matrix of possible systolic and diastolic pressure values, and the minimum least squares difference between observed and predicted envelope functions indicates the best fit choices of systolic and diastolic pressure within the test matrix.

### Results

The model reproduces realistic cuff pressure oscillations. The regression procedure extracts systolic and diastolic pressures accurately in the face of varying pulse pressure and arterial stiffness. The root mean squared error in extracted systolic and diastolic pressures over a range of challenging test scenarios is 0.3 mmHg.

### Conclusions

A new algorithm based on physics and physiology allows accurate extraction of systolic and diastolic pressures from cuff pressure oscillations in a way that can be validated, criticized, and updated in the public domain.

## Background

The ejection of blood from the left ventricle of the heart into the aorta produces pulsatile blood pressure in arteries. Systolic blood pressure is the maximum pulsatile pressure and diastolic pressure is the minimum pulsatile pressure in the arteries, the minimum occurring just before the next ventricular contraction. Normal systolic/diastolic values are near 120/80 mmHg. Normal mean arterial pressure is about 95 mmHg [1].

Blood pressure is measured noninvasively by occluding a major artery (typically the brachial artery in the arm) with an external pneumatic cuff. When the pressure in the cuff is higher than the blood pressure inside the artery, the artery collapses. As the pressure in the external cuff is slowly decreased by venting through a bleed valve, cuff pressure drops below systolic blood pressure, and blood will begin to spurt through the artery. These spurts cause the artery in the cuffed region to expand with each pulse and also cause the famous characteristic sounds called Korotkoff sounds. The pressure in the cuff when blood first passes through the cuffed region of the artery is an estimate of systolic pressure. The pressure in the cuff when blood first starts to flow continuously is an estimate of diastolic pressure. There are several ways to detect pulsatile blood flow as the cuff is deflated: palpation, auscultation over the artery with a stethoscope to hear the Korotkoff sounds, and recording cuff pressure oscillations. These correspond to the three main techniques for measuring blood pressure using a cuff [2].

In the palpatory method the appearance of a distal pulse indicates that cuff pressure has just fallen below systolic arterial pressure. In the auscultatory method the appearance of the Korotkoff sounds similarly denotes systolic pressure, and disappearance or muffling of the sounds denotes diastolic pressure. In the oscillometric method the cuff pressure is high pass filtered to extract the small oscillations at the cardiac frequency and the envelope of these oscillations is computed, for example as the area obtained by integrating each pulse [3]. These oscillations in cuff pressure increase in amplitude as cuff pressure falls between systolic and mean arterial pressure. The oscillations then decrease in amplitude as cuff pressure falls below mean arterial pressure. The corresponding oscillation envelope function is interpreted by computer aided analysis to extract estimates of blood pressure.

The point of maximal oscillations corresponds closely to mean arterial pressure [4–6]. Points on the envelope corresponding to systolic and diastolic pressure, however, are less well established. Frequently a version of the maximum amplitude algorithm [7] is used to estimate systolic and diastolic pressure values. The point of maximal oscillations is used to divide the envelope into rising and falling phases. Then characteristic ratios or fractions of the peak amplitude are used to find points corresponding to systolic pressure on the rising phase of the envelope and to diastolic pressure on the falling phase of the envelope.

The characteristic ratios (also known as oscillation ratios or systolic and diastolic detection ratios [8]) have been obtained experimentally by measuring cuff oscillation amplitudes at independently determined systolic or diastolic points, divided by the maximum cuff oscillation amplitude. The systolic point is found at about 50% of the peak height on the rising phase of the envelope. The diastolic point is found at about 70 percent of the peak height on the falling phase of the envelope [7]. These empirical ratios are sensitive however to changes in physiological conditions, including most importantly the pulse pressure (systolic minus diastolic blood pressure) and the degree of arterial stiffness [9, 10]. Moreover, a rational physical explanation for any particular ratio has been lacking. Since cuff pressure oscillations continue when cuff pressure falls beneath diastolic blood pressure, the endpoint for diastolic pressure is indistinct. Most practical algorithms used in commercially available devices are closely guarded trade secrets that are not subject to independent critique and validation. Hence the best way to determine systolic and diastolic arterial pressures from cuff pressure oscillations remains an open scientific problem.

The present study addresses this problem with a new approach based upon the underlying physics, anatomy, and physiology. This task requires modeling the cuff and arm and the dynamics of a partially occluded artery within the arm during cuff deflation. A second phase of the problem is the development of a regression procedure for analysis of recorded cuff pressure oscillations to extract model parameters and predict the unique systolic and diastolic pressure levels that would produce the observed cuff pressure oscillations.

### Methods Part 1: Modeling cuff pressure oscillations

#### Model of the cuff and arm

Let P_{0} be the maximal inflation pressure of the cuff at the beginning of a run. The pressure is bled down slowly at rate, r mmHg/sec (about 3 mmHg/sec [11]). During the brief period of one heartbeat the amount of air inside the cuff is roughly constant. In addition to smooth cuff deflation, small cuff pressure oscillations are caused by pulsatile expansion of the artery and the corresponding compression of the air in the cuff. One can model the cuff as a pressure vessel having nearly fixed volume, V_{0} − ΔV_{a}, where V_{0} is cuff volume between heartbeats and ΔV_{a} is the small incremental volume of blood in the artery beneath the cuff as it expands with the arterial pulse.

_{a}, in the occluded artery segment it is necessary to know the compliance of the cuff, C = ΔV/ΔP, which is obtainable from Boyle’s law as follows. Boyle’s law is PV = nRT, where P is the absolute pressure (760 mmHg plus cuff pressure with respect to atmospheric), V is the volume of air within the cuff, n is the number of moles of gas, R is the universal gas constant, and T is the absolute temperature. During the time of one heartbeat, n, R, and T are constants and n is roughly constant owing to the slow rate of cuff deflation. Hence to relate the change in cuff pressure, ΔP to the small change in cuff volume, ΔV, from artery expansion we may write $PV\approx \left(P+\Delta P\right)\cdot \left(V+\Delta V\right)\approx PV+P\Delta V+V\Delta P$, for absolute cuff pressure P. So

_{cuff}, or more precisely the time-varying and pressure-varying dynamic compliance of the sealed air inside the cuff, is

In this problem as cuff pressure is slowly released, even as cuff volume remains nearly constant, the dynamic compliance of the cuff increases significantly and its stiffness decreases. Hence a suitably exact statement of the physics requires a differential equation (1a), rather than the constant compliance approximation $P={P}_{0}-rt+{V}_{a}\left(t\right)/{C}_{\mathit{cuff}}$. However, equation (1a) may be integrated numerically to obtain a sufficiently exact representation of cuff pressure changes with superimposed cardiogenic oscillations.

#### Model of the artery segment

_{a}/dt, one can regard the artery as an elastic tube with a dynamic compliance, C

_{a}, which varies with volume and with internal minus external pressure. The dynamic compliance C

_{a}= dV

_{a}/d(P

_{a}– P

_{o}), where P

_{t}= P

_{a}– P

_{o}is the transmural pressure or the difference between pressure inside the artery and outside the artery. Then

where the artery “feels” the prevailing difference between internal blood pressure and external cuff pressure, neglecting the small cuff pressure oscillations. The time derivative of arterial pressure can be determined from a characteristic blood pressure waveform and the known rate, r, of cuff deflation. Hence, the crucial variable to be specified next is the dynamic arterial compliance, C_{a}.

_{t}<0, is easy to imagine. For the negative transmural pressure domain artery volume ${V}_{a}={V}_{a0}{e}^{a{P}_{t}}$ for positive constant, a, and zero pressure volume, V

_{a0}, of the artery. Here the dynamic compliance is clearly

_{0}). This means that for positive constant, b, and typically b < a,

These two exponential functions can be used to characterize the dynamic compliance of the artery model in terms of easily obtained data, including the collapse pressure, P_{c} < 0, defined as the pressure when the artery volume is reduced to 0.1V_{a0} (for example, P_{c} = −20 mmHg) and the normal pressure arterial compliance, C_{n}, measured at normal mid-level arterial pressure, P_{mid} , halfway between systolic and diastolic pressure.

The zero pressure volume, V_{a0}, can be known from anatomy if necessary, but as shown later is not needed if one is interested only in the relative amplitude of cuff pressure oscillations.

_{t}< 0

_{t}≥ 0

#### Forcing function—the time domain blood pressure waveform

_{a}, as a function of time, t, with ω being the angular frequency of the heartbeat, that is ω = 2πf for cardiac frequency, f, in Hz. Here SBP is systolic blood pressure, DBP is diastolic blood pressure, and PP is pulse pressure (SBP − DBP). In turn, the derivative of the arterial pressure waveform is

_{0}, r, V

_{a}, C

_{cuff}, and time. We must work with the time derivative of cuff pressure, rather than absolute cuff pressure, because the compliance of the cuff and also the form of the artery volume vs. pressure function vary with time and pressure during a run. Cuff pressure can then be computed numerically by integrating equation (1a),

with artery pressure, P_{a}, and its time derivative given by equations (5). Combining equations (1a) and (6) gives a precise model for cuff pressure oscillations.

#### Initial conditions

*Artery dimensions:*As a standard normal model consider a brachial artery with internal radius of 0.1 cm under zero distending pressure. The resting artery volume is V

_{a0}= πr

^{2}L or

*Stiffness constant a:*For collapse to 10 percent at −20 mmHg transmural pressure we have

*Stiffness constant b:*It is easy to estimate the normal pressure compliance of the brachial artery in humans, C

_{n}, from experiments using ultrasound. For example, using the data of Mai and Insana [14], the brachial artery strain (Δr/r) during a normal pulse is 4 percent for a blood pressure of 130/70 mmHg with pulse pressure 60 mmHg. In turn the volume of expansion during a pulse is 2πrΔrL, where r is the radius and L is the length of the compressed artery segment. Hence for a normal pressure radius of 0.2 cm the change in volume would be

The normal pressure compliance for the artery segment is the volume change divided by pulse pressure or

C_{n} = 0.10 ml / 60 mmHg = 0.0016 ml/mmHg.

_{mid}, would be 100 mmHg, so

Jeon et al. [13] working with a similar model used a = 0.09 mmHg^{-1.}, b = 0.03 mmHg.

#### Numerical methods

In this model equations (1), (5), and (6) govern the evolution of cuff pressure as a function of time during cuff deflation. Equation (1) can be integrated numerically using techniques such as the simple Euler method coded in Microsoft Visual Basic, Matlab, or “C”. In the results that follow cuff deflation is started from a maximal level of 150 mmHg and continues over a period of 40 sec. Pressures are plotted every 1/20^{th} second. To extract the small oscillations from the larger cuff pressure signal, as would be done in an automatic instrument by an analog high pass filter, cuff pressure at time, t, is subtracted from the average of pressures recorded between times t − Δt/2 and t + Δt/2 , where Δt is the period of the pulse. For simplicity, filtered oscillations are not computed for time points that are Δt/2 seconds from the beginning or from the end of the time domain sample.

### Methods Part 2: Interpreting cuff pressure oscillations

Given this model and the associated insight into the physics of cuff pressure oscillations, one can also devise a scheme for estimating true systolic and diastolic blood pressures from an observed time domain record of cuff pressure and filtered cuff pressure oscillations. The method is based upon the ability, just described, to predict the amplitude of pulse pressure oscillations for a given diastolic pressure and pulse pressure and the ability to deduce exponential constants, a and b, from the rising and falling regions of the oscillation amplitude envelope. Details are as follows.

#### Artery motion during cuff deflation

The pulsatile component of transmural pressure causes the artery to change in volume with each heartbeat. The magnitude of the change in transmural pressure is always equal to the pulse pressure (say, 40 mmHg) which is assumed to be constant during cuff deflation. As cuff pressure gradually decreases from well above systolic to well below diastolic pressure, the range of transmural pressure, P_{t}, experienced by the artery changes. At (a) cuff pressure is well above systolic and net distending pressure is always negative. There is a small change in arterial volume because the artery becomes less collapsed as each arterial pulse makes the transmural pressure less negative. As cuff pressure approaches systolic the relative unloading of negative pressure becomes more profound. Because of the exponential shape of the arterial pressure-volume curve, the amount of volume change accelerates. At (b) cuff pressure is close to systolic. After this point the volume change continues to increase but at a decelerating rate, because of the shape of the pressure-volume curve. Hence (b) is the inflection point for systolic pressure. At (c) cuff pressure is near mean arterial pressure and the volume change is maximal. At (d) cuff pressure is just below diastolic. After this point, as shown in (e), the volume change becomes less and less with each pulse as the increasingly distended artery becomes stiffer. Hence (d) is the inflection point for diastolic pressure. Thus the nonlinear compliance of arteries and the shape of the arterial pressure-volume curve govern the amplitude of cuff pressure oscillations.

The particular volume change of the artery from the nadir of diastolic pressure to the subsequent peak of systolic pressure can be specified analytically from Equations (4a) and (4b) as follows. Consider P_{t} as the transmural pressure at the diastolic nadir of the arterial blood pressure wave and let PP be the pulse pressure. One can imagine three domains of transmural pressure. In Domain (1) P_{t} + PP < 0. In Domain (2) P_{t} < 0 and P_{t} + PP ≥ 0. In Domain (3) P_{t} > 0. The largest artery volume oscillations occur in Domain (2) when transmural pressure oscillates between positive and negative values. Doman (1) represents the head of the oscillation envelope in time, and Domain (3) represents the tail.

Using equations (4), the artery volume changes during the rising phase of the arterial pulse in each of the three domains are

where for cuff pressure, P, systolic blood pressure SBP, and diastolic blood pressure DBP, the transmural pressure P_{t} = DBP − P, and the pulse pressure PP = SBP − DBP.

It is easy to show by differentiating expressions (7) for Domains (1), (2), and (3) that the systolic and diastolic pressure points correspond exactly to the maximal and minimal slopes d(ΔV_{a})/dP_{t}. Therefore a simple analysis for finding systolic and diastolic pressures points would involve taking local slopes of the oscillation envelope vs. pressure function. Slope taking, however, is vulnerable to noise in practical applications. An alternative approach that does not involve slope taking creates a model of each individual subject’s arm in terms of exponential constants a and b and then numerically finds the unique combination of systolic and diastolic arterial pressures that best reproduces the observed oscillation envelope.

#### Regression analysis for exponential constants

for constant, k_{1}, during a cuff deflation scan in which cuff pressure, P, varies and the other variables are constant. (Note that here C_{cuff} is very nearly constant because the rising phase of the pulse happens in a very short time, roughly 0.1 sec.) Hence, $ln\left(\Delta P\right)=ln\left({k}_{1}\right)-a\phantom{\rule{0.1em}{0ex}}P$, and a regression plot of the natural logarithm of the amplitude of pulse oscillations in the leading region of the envelope versus the instantaneous cuff pressure, P, yields a plot with slope − a. Thus we can obtain by linear regression an estimate of stiffness constant, a, as $\widehat{a}=slop{e}_{1}$. The range of the rising phase of the oscillation envelope from the beginning of the envelope to the first inflection point (maximal slope) can be used for the first semi-log regression. More simply, the range of the rising phase of the oscillation envelope from its beginning to one third of the peak height provides reasonable estimates of slope_{1}.

hence, $ln\left(\Delta P\right)=ln\left({k}_{3}\right)+b\phantom{\rule{0.1em}{0ex}}P$, and a regression plot of the natural logarithm of the amplitude of pulse oscillations in the envelope tail versus cuff pressure at the time of each pulse yields a plot with slope b. In turn, we can obtain by linear regression an estimate of stiffness constant, b, as $\widehat{b}=slop{e}_{3}$. The range of the falling phase of the oscillation envelope from the second inflection point (maximal negative slope) of the oscillation envelope to the end of the envelope can be used to define the range of the second semi-log regression. More simply, the range of the falling phase of the oscillation envelope from two thirds of the peak height to the end of the envelope provides reasonable estimates of slope_{3}. The slope estimates from the head and tail regions of the amplitude envelope include multiple points and so are relatively noise resistant. Other variables involved in the lumped constants, k_{1} and k_{3}, are not relevant to the estimation of exponential constants a and b.

#### Least squares analysis

Having estimated elastic constants a and b for a particular envelope of oscillations from a particular patient at a particular time, it is straightforward in a computer program to find SBP and DBP values that reproduce the observed envelope function most faithfully. Let y(P) be the observed envelope amplitude as a function of cuff pressure, P, and let y_{max}(P_{max}) be the observed peak amplitude of oscillations at cuff pressure P_{max}. Let $\widehat{y}$(P, SBP, DBP) be the simulated envelope amplitude as a function of cuff pressure, P, for a particular pulse and a particular test set of systolic and diastolic pressure levels. The values of $\widehat{y}$ are obtained from equations (7) and the prevailing cuff compliance as follows

_{max}(P

_{max}, SBP, DBP) be the predicted peak of the oscillation envelope at cuff pressure P

_{max}. A figure of merit for goodness of fit between modeled and observed oscillations for particular test values of SBP and DBP is the sum of squares over all measured pulses

The values of SBP and DBP that minimize this sum of squares are the taken as the best estimates of systolic and diastolic pressure by the oscillometric method.

Here cuff pressure, P, is the cuff pressure at the time of each oscillation. Use of the amplitude normalized ratios y/y_{max} and $\widehat{y}$/$\widehat{y}$
_{max}, means that it is not necessary to know the zero pressure volume of the artery, V_{a0} , or cuff volume V_{0}, which depend on anatomy and geometry of a particular arm and cuff and are constants. It is the shape of the amplitude envelope in the pressure domain that contains the relevant information. The least squares function, SS, includes information from all of the measured oscillations and so is relatively noise resistant.

A variety of numerical methods may be used to find the unique values of SBP and DBP corresponding to the minimum sum of squares. Here, to demonstrate proof of concept, we evaluate the sum of squares, SS, over a two-dimensional matrix of candidate systolic and diastolic pressures at 1 mmHg intervals and identify the minimum sum of squares by plotting. The values of SBP and DBP corresponding to this minimum sum of squares are the best fit estimates for a particular oscillometric pressure run. The best fit model takes into account the prevailing artery stiffness and also the prevailing pulse pressure.

## Results and discussion

### Normal model

**Standard parameters for the oscillometric blood pressure model**

Parameter | Definition | Value | Units |
---|---|---|---|

P | Cuff pressure at onset of deflation | 150 | mmHg |

r | Cuff pressure decay rate | 3 | mmHg/sec |

PP | Arterial pulse pressure | 40 | mmHg |

f | Cardiac frequency (heart rate) | 80 | beats/min |

V | Artery segment volume at zero pressure | 0.3 | ml |

a | Exponential constant | 0.03 | 1/mmHg |

b | Exponential constant | 0.11 | 1/mmHg |

C | Artery segment compliance at 100 mmHg pressure | 0.0016 | ml/mmHg |

### Varying arterial compliance

### Varying pulse pressure

### Regression analysis for systolic and diastolic pressures

**Validation of algorithm for estimation of systolic and diastolic pressures**

Scenario and parameter | Actual value | Algorithm value |
---|---|---|

Normal | ||

Constant a (mmHg | 0.11 | 0.105 |

Constant b (mmHg | 0.03 | 0.0311 |

SBP (mmHg) | 120 | 119 |

DBP (mmHg) | 80 | 80 |

Twice normal stiffness | ||

Constant a (mmHg | 0.076 | 0.074 |

Constant b (mmHg | 0.021 | 0.0214 |

SBP (mmHg) | 120 | 119 |

DBP (mmHg) | 80 | 80 |

Half normal stiffness | ||

Constant a (mmHg | 0.158 | 0.155 |

Constant b (mmHg | 0.0432 | 0.043 |

SBP (mmHg) | 120 | 119 |

DBP (mmHg) | 80 | 80 |

Half normal pulse pressure | ||

Constant a (mmHg | 0.110 | 0.108 |

Constant b (mmHg | 0.030 | 0.030 |

SBP (mmHg) | 110 | 110 |

DBP (mmHg) | 90 | 89 |

Twice normal pulse pressure | ||

Constant a (mmHg | 0.11 | 0.105 |

Constant b (mmHg | 0.03 | 0.030 |

SBP (mmHg) | 140 | 138 |

DBP (mmHg) | 60 | 60 |

_{max}, shown as filled circles. A high value of test pressure (130/90) leads to equally discordant reconstructions in the opposite direction (heavy dashed curve to right). For both low and high test values the sum of squared differences is obviously large. The reconstructed model for the actual pressure (120/80) is shown as the solid curve. This illustration demonstrates the sensitivity of the least squares approach.

### Validation of the regression procedure

The cuff-arm-artery model of an oscillometric pressure measurement described in Methods Part One can be used to validate the regression and analysis procedure of Methods Part Two. An unlimited number and wide variety of test scenarios can be simulated in the model as unknowns for testing by the regression scheme, including a wide range of arterial stiffnesses and a wide range of pulse pressures, heart rates, blood pressure waveforms, cuff sizes, arm sizes, cuff lengths, artery diameters, etc. Importantly, the regression analysis assumes no prior knowledge of these model parameters or of the blood pressure used to generate the simulated oscillations. Cuff pressure oscillations and absolute cuff pressure are the only inputs to the algorithm for obtaining systolic and diastolic pressures.

The data summarized in Table 2 show the effectiveness of the proposed regression procedure in small sample of various possible test scenarios, including varying artery stiffness and varying pulse pressure. This small, systematic sample includes challenging cases for the algorithm. The accuracy is quite satisfactory, with reconstructed pressures within 0, 1, or 2 mmHg of the actual pressures in the face of varying artery stiffness and varying pulse pressure. The root mean squared error is $\sqrt{8}/10$ = 0.28 mmHg.

## Discussion

The challenge of creating a satisfactory theoretical treatment of the genesis and interpretation of cuff pressure oscillations has attracted a diverse community of thinkers [4, 5, 7–10, 16]. Nevertheless, specifying a valid method for extracting systolic and diastolic pressures from the envelope of cuff pressure oscillations remains an open problem. Here is presented a mathematical model incorporating anatomy, physiology, and biomechanics of arteries that predicts cuff pressure oscillations produced during noninvasive measurements of blood pressure using the oscillometric method. Understanding of the underlying mechanisms leads to a model-based algorithm for deducing systolic and diastolic pressures accurately from cuff pressure oscillations in the presence of varying arterial stiffness or varying pulse pressure.

The shape of the oscillation amplitude envelope dictates the stiffness parameters for the artery during both compression and distension. Semilog regression procedures give good estimates of the artery stiffness parameters that characterize each individual cuff deflation sequence. Using these parameters one can create and exercise an individualized cuff-arm-artery model for a wide range of possible systolic and diastolic pressures. The pair of systolic and diastolic pressures that best reproduces the observed oscillation envelope according to a least squares criterion constitutes the output of the algorithm.

When applied to amplitude normalized oscillation data the algorithm is insensitive to variations among subjects in zero pressure artery volume, V_{a0} , or initial cuff volume V_{0}, since these terms are constants that are eliminated by the normalization procedure. Compression of the entire length of artery underlying the cuff is not necessary. Incomplete coupling of cuff pressure to the artery near the ends of the cuff merely decreases the ratio V_{a0} / V_{0} without effecting the extracted systolic and diastolic pressures.

The cuff-arm-artery model can be used as well to test the validity of the algorithm for over a wide range of possible conditions by generating trial cuff pressure data for known arterial pressure waveforms. A stress test for the algorithm can be done by comparing systolic and diastolic pressure levels extracted from synthesized cuff pressure oscillations with the arterial pressure that generated the synthesized oscillations over a wide range of test conditions. These conditions may include extreme cases that are hard to reproduce experimentally, contamination with excessive noise, any conceivable blood pressure waveforms, cardiac arrhythmias such as atrial fibrillation, etc. Such computational experiments, in addition to future animal and clinical studies, can boost confidence in the reliability of the oscillometric method and can suggest further refinements.

Here for convenience we have used the bi-exponential model to generate cuff pressure oscillations for algorithm testing. However, the regression algorithm does not “know” where the sample data came from. It tries to extract constants a and b from the head and tail portions of the semi-log plot of oscillation amplitude versus cuff pressure. The resulting best fit values of a and b will still work for non-ideal or noise contaminated data to produce a model envelope that can be matched to the actual data. An extremely stiff artery with a linear pressure volume curve is easily accommodated by this process, since e^{x} ≈ 1 + x for small values of x. In this limiting case the exponential pressure-volume curve becomes linear. An exceptionally flabby artery, rather like dialysis tubing, is well described by larger values of a and b and a larger ratio a/b. Thus the family of bi-exponential models is very inclusive of a wide range of arterial mechanical properties, as suggested in Figure 3.

Classically the oscillometric method has been relatively well validated as a measure of mean arterial pressure, which is indicated by the peak of the oscillation amplitude envelope [4]. Automated oscillometric pressure monitors have found use in hospitals for critical care monitoring in which the goal is to detect any worrisome trend in blood pressure more so than the exact absolute value. Out of hospital use of the oscillometric method in screening for high blood pressure is more problematic, because heretofore the accuracy of systolic and diastolic end points has been questioned and doubted. For example Stork and Jilek [17] studied two published algorithms, differing in detail and based on cuff oscillation ratios of either 50% for systolic and 80% for diastolic or 40% for systolic and 55% for diastolic. Compared to a reference pressure of 122/78 mmHg the algorithmic methods applied to oscillometric data gave pressures of 135/88 and 144/81 mmHg, respectively. An advisory statement from the Council for High Blood Pressure Research, American Heart Association [18] stressed the need for caution in the selection of all instruments used for blood pressure determination and the need for continuing studies to validate their the safety and reliability.

Accurate measurements of blood pressure in routine clinic and office settings are important because systemic arterial hypertension is a major cause of serious complications, including accelerated atherosclerosis, heart attacks, strokes, kidney disease, and death. These serious complications increase smoothly with every point above the nominal 120/80 mmHg, hence even small increases in blood pressure are important to detect. In screening for hypertension systematic bias or inaccuracy in blood pressure readings of a few mmHg can be significant, since the difference between high normal (85 diastolic) and abnormal (90 diastolic) is only a few mmHg. A recent 1 million-patient meta-analysis suggests that a 3–4 mmHg increase in systolic blood pressure would translate into 20% higher stroke mortality and a 12% higher mortality from ischemic heart disease [19].

False negative readings would be problematic because untreated high blood pressure can lead to strokes, blindness, kidney failure, and lethal heart attacks. False positive readings would be undesirable because the usual drugs for hypertension must be taken every day for life and can be expensive. They also have side effects. Hence accurate readings are essential. Given a reliable algorithm for extracting systolic and diastolic pressures, an automatic oscillometric device could provide screening for high blood pressure that is performed in the same way each time without inter-observer variation. The present research could lead to a wider role for oscillometric blood pressure monitors in physicians’ offices and clinics.

## Conclusions

The analytical approach and algorithm presented here represent a solution to an open problem in biomedical engineering: how to determine systolic and diastolic blood pressures using the oscillometric method. Current algorithms for oscillometric blood pressure implemented in commercial devices may be quite valid but are closely held trade secrets and cannot be independently validated. The present paper provides a physically and physiologically reasonable approach in the public domain that can be independently criticized, tested, and refined. Future demonstration of real-world accuracy will require data comparing oscillometric and intra-arterial pressures in human beings over a range of test conditions including variable cuff size, arm diameter, cuff tightness, cuff deflation rate, etc. Further development and incorporation of this algorithm into commercial devices may lead to greater confidence in the accuracy of systolic and diastolic pressure readings obtained by the oscillometric method and, in turn, an expanded role for these devices.

## Author’s contributions

CB is the only author and is responsible for all aspects of the research and the intellectual and technical content of the manuscript.

## Declarations

## Authors’ Affiliations

## References

- Boron WF, Boulpaep EL:
*Medical Physiology*. 2nd edition. Philadelphia: Elsevier; 2005.Google Scholar - Geddes LA, Baker LE:
*Prinicples of Applied Biomedical Instrumentation*. Second edition. 1975.Google Scholar - Soueidan K, Chen S, Dajani HR, Bolic M, Groza V:
**Augmented blood pressure measurement through the noninvasive estimation of physiological arterial pressure variability.***Physiol Meas*2012,**33:**881–899. 10.1088/0967-3334/33/6/881View ArticleGoogle Scholar - Mauck GW, Smith CR, Geddes LA, Bourland JD:
**The meaning of the point of maximum oscillations in cuff pressure in the indirect measurement of blood pressure–part ii.***J Biomech Eng*1980,**102:**28–33. 10.1115/1.3138195View ArticleGoogle Scholar - Geddes LA, Voelz M, Combs C, Reiner D, Babbs CF:
**Characterization of the oscillometric method for measuring indirect blood pressure.***Ann Biomed Eng*1982,**10:**271–280. 10.1007/BF02367308View ArticleGoogle Scholar - Ursino M, Cristalli C:
**A mathematical study of some biomechanical factors affecting the oscillometric blood pressure measurement.***IEEE Trans Biomed Eng*1996,**43:**761–778. 10.1109/10.508540View ArticleGoogle Scholar - Baker PD, Westenskow DR, Kuck K:
**Theoretical analysis of non-invasive oscillometric maximum amplitude algorithm for estimating mean blood pressure.***Med Biol Eng Comput*1997,**35:**271–278. 10.1007/BF02530049View ArticleGoogle Scholar - Drzewiecki G, Hood R, Apple H:
**Theory of the oscillometric maximum and the systolic and diastolic detection ratios.***Ann Biomed Eng*1994,**22:**88–96. 10.1007/BF02368225View ArticleGoogle Scholar - Forster FK, Turney D:
**Oscillometric determination of diastolic, mean and systolic blood pressure–a numerical model.***J Biomech Eng*1986,**108:**359–364. 10.1115/1.3138629View ArticleGoogle Scholar - van Popele N, Bos W, de Beer N, van der Kuip D, Hoffman A, Diederick E, CM W:
**Arterial stiffness as underlying mechanism of disagreement between oscillometric blood pressure monitor and sphygmomanometer.***Hypertension*2000,**36:**484–488. 10.1161/01.HYP.36.4.484View ArticleGoogle Scholar - Pickering TG, Hall JE, Appel LJ, Falkner BE, Graves J, Hill MN, Jones DW, Kurtz T, Sheps SG, Roccella EJ:
**Recommendations for blood pressure measurement in humans and experimental animals: part 1: blood pressure measurement in humans: a statement for professionals from the Subcommittee of Professional and Public Education of the American Heart Association Council on High Blood Pressure Research.***Circulation*2005,**111:**697–716. 10.1161/01.CIR.0000154900.76284.F6View ArticleGoogle Scholar - Posey J, Geddes L:
**Measurement of the modulus of elasticity of the arterial wall.***Cardiovasc Res Cent Bull*1973,**11:**83–103.Google Scholar - Jeon G, Jung J, Kim I, Jeon A, Yoon S, Son J, Kim J, Ye S, Ro J, Kim D, Kim C:
**A simulation for estimation of the blood pressure using arterial pressure-volume model.***World Acad Sci Eng Technol*2007,**30:**366–371.Google Scholar - Mai JJ, Insana MF:
**Strain imaging of internal deformation.***Ultrasound Med Biol*2002,**28:**1475–1484. 10.1016/S0301-5629(02)00645-2View ArticleGoogle Scholar - Ursino M, Cristalli C:
**Mathematical modeling of noninvasive blood pressure estimation techniques–Part I: pressure transmission across the arm tissue.***J Biomech Eng*1995,**117:**107–116. 10.1115/1.2792258View ArticleGoogle Scholar - Ursino M, Cristalli C:
**Mathematical modeling of noninvasive blood pressure estimation techniques–Part II: brachial hemodynamics.***J Biomech Eng*1995,**117:**117–126. 10.1115/1.2792259View ArticleGoogle Scholar - Stork M, Jilek J:
**Cuff pressure waveforms: thier current and prospective application in biomedical instrumentation.**In*Biomedical Engineering Trends in Electronics, Communications and Software*. Edited by: Laskovski AN. InTech; 2011:193–210.Google Scholar - Jones DW, Frohlich ED, Grim CM, Grim CE, Taubert KA:
**Mercury sphygmomanometers should not be abandoned: An advisory statement from the Council for High Blood Pressure Research, American Heart Association.***Hypertension*2001,**37:**185–186. 10.1161/01.HYP.37.2.185View ArticleGoogle Scholar - Lewington S, Clarke R, Qizilbash N, Peto R, Collins R:
**Age-specific relevance of usual blood pressure to vascular mortality: a meta-analysis of individual data for one million adults in 61 prospective studies.***Lancet*2002,**360:**1903–1913.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.