`
`[19]
`
`11
`
`||||I|||l|||l|||||||||||||||||ll|||||l|||||l||||||||l||||||||||||||||||||||l
`
`US00524-3992A
`Patent Number:
`
`9
`3
`5 243 992
`
`Eckerle et al.
`
`[451 Date of Patent:
`
`Sep. 14, 1993
`
`123/690
`364/-H‘!
`123/673
`
`364-I41?
`
`4,239,048 ‘I2/1980 Steuer ............................... .. 128/666
`4,3C|'.-‘J23 12/1931 W l
`4.353.372 10/1982 A;:troll.....
`4.409.933 10/1933 Albert
`4,456,959
`6/1984 Hirano et al.
`4.646.749
`3/1937 Berger 151 81-
`4,66?,680
`.......... ..
`4,7l2,lT9 12/1937 He_J.r|1er :..
`‘
`figgfig
`4'ao2'-133 2/1939 Eek:-I: IL...
`4.836.2I3
`6/1989 Wenulet ai.
`.
`.
`Primary Exammer—Lee S. Cohen
`Assistant Exam:'ner—Ma.|'ianne H. Parker
`Attorney, Agent. or FJrmA0I1fl' dc Berndge
`[57]
`ABSTRACT
`A pulse rate Sensor system is packaged in a wristwatch
`sized assembly and is worn by the user to provide an
`accurate determination of pulse rate. A tonometer sen~
`sor is provided to detect heartbeat pressure waves pro-
`duced by a superficial artery. A microcomputer ruanip-
`ulates the unprocessed tonorneter sensor element signals
`using multiple algorm-‘ms to determine an accurate
`pulse me
`‘
`
`[54] PULSE RATE SENSOR SYSTEM
`[75]
`Inventors:
`Joseph 5- Eckerle. Redwooql City;
`Dale W. Ploezer. San Francisco;
`Steven T. Holmes, Palo Alto; 'I'|:oma.5
`P. Low, Woodside; Rudolf Elbrecht,
`L05
`R_ Jenck, In,
`Mcnio park: Ronald E_ Peking,
`Menlu Park; Victor T. Newton, Jr.,
`Men“-’ Park‘ “'1 °f Calif‘
`com. Eyecu-onics C9" Ltd” Agchi,
`Japan
`
`[73] Assignee,
`
`[21] Appl. No.: 502,028
`[22] Filed:
`M". 30' 1990
`[511
`Int. Cl.5 ................................................ A613 5/02
`
`[521 U.S. Cl. ............... ..
`123/690; l23/672
`[58] Field of Search ............... 128/687, 722, 699. 689,
`123/677' 631' 633* 690
`References Cited
`U_S' PATENT DOCUMENTS
`
`[55]
`
`3,999,531‘ 12/1976 Noiles ................................ .. 123/687
`4.053.113 11/1971 Stupay etal.
`. I28/2.{J5T
`4.086.916
`5/19'.-'8 Freeman ei al.
`..
`.... .. 128/690
`4.l8l,l34 1/1980 Mason et al. . . . . . .
`. . . . .. 128/689
`4.202.350 5/ 1980 Walton .............................. .. 128.3690
`
`
`
`
`
`123/3:
`.................... .. I28/672
`
`36 Claims, 9 Drawing Sheets
`
`°°°‘
`
`U.S. Patent No. 89239
`
`Apple I
`APL10
`
`Apple Inc.
`APL1047
`U.S. Patent No. 8,923,941
`
`0001
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 1 of 9
`
`5,243,992
`
`FIG.I
`
`0002
`
`
`
`US. Patent
`
`Sep. 14, 1993
`
`Sheet 2 of 9
`
`5,243,992
`
`
`
`AxIs?1
`
`MOUNT TO smms 11
`
`FlG.2
`
`MOUNTING
`
`PLANE
`
`FIG. 3A
`
`FIG 3B
`
`0003
`
`0003
`
`
`
`US. Patent
`
`Sep. 14, 1993
`
`Sheet 3 of 9
`
`5,243,992
`
`
`
`PROCESSING68
`
`...|
`
`Q2I
`
`-ZI
`
`L!
`L)
`
`DIGITALCONVERTER
`
`0004
`
`0004
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 4 of 9
`
`5,243,992
`
`INITIALIZE
`
`
`’§I3:?;IA'EN—TaT=1" ' "'
`
`l— — — - — — — — — — -— - —
`ISEGIIENT
`I#2
`?v‘I~E+EL'3LII§Es'%'Ts
`1F
`SIGNAL TD ZERO
`
`ISEGMENT
`|#3
`
`
`
`"1
`
`II I I I I I
`
`I
`'}
`
`I JI I III I I
`
`‘1
`
`
`
`SUBTRACT
`AVERAGE FRQM
`
`- LL ELEMENTS
`
` ACCUMULATE
`
`
`CORRELATION
`SUMMATIONS
`
`
`MARKING
`TIME AFTER
`SYSJOLIC
`
`0005
`
`0005
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 5 of 9
`
`5,243,992
`
`SLOPE
`
`THRESHOLD at 1/2
`
` ‘'
`
`
`RESET
`svsrouc
`man m o
`
`- l
`
`
`
`POUND svsrouc sn 5;
`SIGNAL wtun str
`SYSTDLIC FOUWJ
`
`SLOPE n-mssnoun 5
`551; onmxsaoee
`
`°
`
`save svsrouc
`AND IIASTOLIC
`
`RESETABSGJJTE
`uluuuu AND
`MIXIMUU STORAGE
`
`5,,
`
`6!
`
`0006
`
`
`
`
`
`
`
`
`
`6m
`
`
`ABSOLUTE
`VALUE OF CURREN'l'
`SLOPE EJCEEDS
`
`Y
`5595 smpg
`AS max SLOPE
`
`5"
`
`'"""3"5§ M5 5.
`nowuuau me
`
`5151' rue rm-r
`SIDFE mass» - GP
`“C5595”
`
`0
`
`0006
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 6 of 9
`
`5,243,992
`
`IZETERMINE ELEMENT
`um-H mxmw
`{PULSE AMPLITLDEJ
`
`— — —- ———————— — -
`sE<sIIIENT#e“
`
`"
`
`H MST STRIGGERS
`
`
`
`0007
`
`
`
`AVERAGE PULSE
`RATE ' (OLD AVERKSB
`
`CURRENT PEFHOD
`ICURRENT PERIOD
`+ PREVIOUS PERIOD
`
`r -------- - - ‘I
`SEGIENT
`I
`I
`9 #7
`I
`I
`i
`I
`I
`I
`I
`I
`I
`I
`I
`I
`'
`{
`'
`Tb .
`:
`.
`I
`I
`=
`l
`I
`I
`I
`I
`
`DETERMINE ELEMENT
`WITH HAXIMUM
`POSITIVE PLI.SE
`AMPLITUE lF
`AVAILABLE
`
`SET ELEHENTO-F
`
`
`INTEREST T0 ELEMENT '3:
`WITH MOSTMAXJMJMS
`
`
`
`I I
`
`'
`|
`I
`
`0007
`
`
`
`US. Patent
`
`Sep. 14, 1993
`
`Sheet 7 of 9
`
`5,243,992
`
`
`
`I""—_ — _ - _ _ — --
`sasuem "'
`I CALCULATE
`I
`
`: I
`
`-129
`
`,_-___
`
`__________ __ I
`
`F
`
`I N
`
`“Z I
`coRI'II-:LATIo;I
`I
`.11 COEF ICE
`12c
`',.
`12b
`
`'I
`I 0
`SLGLTQE
`N
`I‘
`COEFII-‘l6C‘IENT
`,
`FALSE
`I:
`«2.o9
`IL_ _ _ __ _1___,_____; 2 _ __ _I
`'4”
`Ir * -H‘
`Sm“?!
`SET
`DISPLAY II CM-CUUITE
`#13 I
`
`
`
`
`
`0
`
`SEGMENT
`
`#"
`
`14a
`PULSE IN LAST
`
`N
`
`III I II II I I I II
`
`I I I I I I I I
`
`I
`1'
`
`I I I
`
`13°
`
`.3,
`
`we
`
`Ed
`
` vnuo
`
`305Ec_ ? FT!)BLANK .: PULSERATE
`
`I I I
`
`: SET POST SYSTDLIC
`I. ____________ — — -JI DELAY TO .2125
`, or PULSE TIME
`I
`
`' RESET svsrouc
`I TIMER
`0 I
`I UPDATE DISPLAY
`I
`IFVALIO PULSE FUJND
`
`0008
`
`0008
`
`
`
`S”U
`
`B
`
`W
`
`5
`
`293,
`
`
`4..._<zo_mBuzaxzu._$._o_m._<_,_o_m
`
`
`
`
`Dom»._._q:z.”___.rmEaEo<fi._>q3%x.Ezm__.m_._m
`
`
`
`
`
` ,1....WmoI228n9:O.n-Nmo22_
`
`
`
`
`
`
`
`
` 288.mum»E8m9n-canonem-8-N.mONcomo8on0._M:.32;.2mE8:<.E_8:32:58.5Eoo_E.n._mou200..
`
`mN0_g
`
`
`
`0009
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 9 of 9
`
`5,243,992
`
`0010
`
`0010
`
`
`
`5, 243,992
`
`5
`
`I0
`
`25
`
`30
`
`35
`
`45
`
`I
`
`PUISE RATE SENSOR SYSTEM
`
`FIELD OF THE INVENTION
`
`The present invention relates to pulse monitors with
`visual readouts of pulse rate and more particularly to a
`tonorneter sensor pulse rate monitor which employs
`multiple noise and motion artifact rejection methods to
`determine an accurate pulse rate.
`BACKGROUND OF THE INVENTION
`
`The present invention relates generally to a method
`and apparatus for measuring and displaying pulse rate
`with increased accuracy. More specifically, the present
`invention provides a method for increasing the accu-
`racy of a pulse rate sensing system by means of a novel
`pressure sensing array and multiple methods for identifi-
`cation and elimination of artifacts.
`Other methods and apparatus are known for measur-
`ing pulse rales and for rejecting pulse artifacts. For
`example, U.S. Pat. No. 4,409,983 shows a pulse measur-
`ing device which employs multiple transducers con-
`nected to averaging circuits and differential amplifiers.
`This invention helps separate signals corresponding to
`motion artifacts from the signal corresponding to a
`heartbeat pulse. Other apparatus and methods for re-
`moving motion artifacts are disclosed in U.S. Pat. Nos.
`4,307,723, 4,202,350, 4,667,680, 4,239,048, 4,181,134 and
`4,456,959. Methods used to reduce signal errors include
`the use of windowing and averaging techniques and
`auto correlation algorithms.
`The pulse rate sensor systems described above are
`subject to several sources of inaccuracies. First, it is
`difficult to reject motion and noise artifacts in many of
`these systems. This is especially true for systems em-
`ploying a single sensor element. (See U.S. Pat. Nos.
`4,202,350 and 4,239,048.) These systems have no physi-
`cal means for receiving both a pulse-plus-artifact signal
`and a separate artifact signal. Other means are required
`to compensate for, or eliminate, the error caused by
`artifacts such as motion artifacts. Signal processing
`techniques such as filtering and windowing are often
`used.
`Even those systems or methods which employ multi-
`ple sensor elements inaccurately measure pulse rate
`because only a single method is used for enhanced sig-
`nal processing. For example, different types of motion
`artifacts can occur simultaneously, and with other per-
`tubations, on the pulse sensor. It is also not unusual for
`signal errors to be interpreted as pulses or for actual
`pulses to be missed by the pulse sensor. Methods of
`pulse rate determination which do not compensate for
`these errors are inherently inaccurate under real-world
`conditions ere artifacts are present.
`For example. if a pulse rate system detects a “pulse“
`caused by noise. several adverse results may be seen.
`The pulse rate system could use the noise as the basis for
`windowing the signal. The pulse rate system could
`simply use this "pulse" as part of the overall pulse rate
`calculation. In addition, the pulse rate system could
`recognize the noise as noise and subtract out the noise.
`in some cases subtracting out a valid signal as well.
`Another source of inaccuracy that occurs using pulse
`measuring devices that measure pressure variations
`caused by a subject’s pulse (see U.S. Pat. No. 4,409.983
`for example) is inverted pulse waveforms. An inverted
`waveform can occur when the housing that holds the
`pressure sensitive element(s) is located on the artery,
`
`2
`but the pressure sensitive elernentfs) itself is located off
`the artery. In this case the subject's pulse can push up on
`the housing and lessen the pressure on the pressure
`sensitive element. The result is that a pulse waveform is
`still received, but it is inverted and shows a negative
`relative pressure. Pulse measuring devices which rely
`on pressure measurements but can correctly interpret
`only positive pressure waveforms must be placed and
`held accurately on the artery, creating additional de-
`mands on attachment of the device and/or lowering
`comfort to the user.
`Additionally, pronounced dichrotic notches can be
`found in the pulse of many people. When dichrotic
`notches are present there are two rises and two falls in
`blood pressure during a single heartbeat. These can be
`mistakenly interpreted as two heartbeats. leading to a
`major inaccuracy in pulse rate measurement.
`The present invention overcomes the problems en-
`countered with other pulse rate sensors by applying the
`principles of arterial tonometry for signal acquisition
`for a pulse rate sensor. In the invention. multiple algo-
`rithms are used in signal processing and pulse rate cal-
`culation to compensate for multiple signal errors which
`could occur during pulse rate measurement.
`The principles of arterial tonometry are described in
`several U.S. Patents including: U.S. Pat. Nos. 3,219,035;
`4,799,491 and 4,802,438. These principles are also de-
`scribed in several publications including an article enti-
`tled "Tonometry, Arterial,“ in Volume 4 of the Ency-
`clopedia of Medical Devices and Instruments. (J. G.
`Webster, Editor, John Wiley & Sons, I933). All of these
`references discuss arterial tonometry as used for the
`measurement of blood pressure.
`For blood pressure measurement, it is desirable to
`flatten a section of the arterial wall as described in these
`
`references. Flattening is produced by exerting an appro-
`priate hold down force on the tonometer sensor. For
`pulse sensing, significant flattening of the arterial wall is
`not necessary and a lower hold down force can be used.
`This results in greater comfort for the wearer.
`SUMMARY OF THE INVENTION
`
`Accordingly, the present invention has been devel-
`oped to overcome the foregoing shortcomings of exist-
`ing pulse rate sensor systems.
`It is therefore an object of the present invention to
`_ provide a method and an apparatus for measuring pulse
`rates using arterial tonometry techniques including a
`sensor array with multiple sensing elements disposed in
`an array, in order to provide increased accuracy in the
`determination of pulse rate.
`Another object of the present invention is to increase
`the accuracy of the displayed pulse rate by calculating
`the displayed pulse rate using only pulse rates deter-
`mined to be valid.
`A further object of the present invention is to deter-
`mine whether pulses detected are valid, based on the
`correlation between the present pulse and the previous
`pulse.
`A still further object of the present invention is to
`remove motion artifacts from the sensor element signals
`by subtracting a value from all these signals based on a
`spatially weighted average of these signals.
`An additional object of the present invention is to
`cancel out artifacts from a sensor element which exceed
`a level predetermined to be the maximum level of a
`valid blood pressure signal.
`
`50
`
`55
`
`65
`
`0011
`
`0011
`
`
`
`3
`Still another object of the present invention is to
`accurately process inverted waveforms caused by mis-
`alignment or shifting of the sensor elements relative to
`an underlying artery.
`These and other objects and advantages are achieved
`in accordance with the present invention by the steps
`of: sensing at least one blood pressure waveform signal
`at a predetermined sampling period using a tonometer
`sensor having a plurality of sensor elements disposed in
`an array; producing a plurality of sensor element sig-
`nals, at least one of the sensor element signals corre-
`sponding to the at least one blood pressure signal; cor-
`recting the sensor element signals using a correction
`factor based on one characteristic of the sensor element
`signals; calculating a plurality of slopes based on the
`corrected sensor element signals; selecting a corrected
`sensor element signal corresponding to one of the sen-
`sor elements, the selected sensor element signal having
`slopes greater than a predetermined slope threshold;
`determining a plurality of pulse rates based on the se-
`lected sensor element signal; computing a value corre-
`sponding to the autocorrelation of the corrected sensor
`element signal over a predetermined time period; and
`calculating a display pulse rate based on at least two of
`the pulse rates, each of the two pulse rates having the
`value within a predeterrnined range.
`These and other objects and advantages are achieved
`in accordance with the preferred embodiment of the
`present invention comprising: a tonometer sensor means
`having a plurality of pressure sensing elements disposed
`in an array, for sensing a blood pressure waveform on at
`least one of the pressure sensing elements and produc-
`ing a plurality of sensor element signals. at least one of
`the sensor element signals being indicative of the blood
`pressure acting on at least one of the pressure sensing
`elements; means for pivoting the tonometer sensor
`means about a pair of axes; means for pressing the to-
`nometer sensor means against a radial artery of a sub-
`ject; means for anchoring the pressing means on a dorsal
`side of the subject; central processing means for deter-
`mining a pulse rate based on at least one of the sensor
`element signals received from the tonometer sensor
`means; means for displaying the pulse rate; and means
`for holding the anchoring means on the subject, the
`holding means at no time contacting the pressing means.
`BRIEF DESCRIPTION OF THE DRAWINGS
`
`The preferred embodiments are described with refer-
`ence to the drawings in which:
`FIG. 1 is a general arrangement of a pulse rate sensor
`system illustrating a pulse rate sensor connected to a
`case assembly;
`FIG. 2 is a perspective view of the tonometer sensor
`supported in a girnbal assembly;
`FIGS. 3A. B and C show low, medium and high
`curvature spring profiles, respectively.
`FIG. 4 is a block diagram of the pulse rate processing
`circuitry in accordance with the preferred embodiment
`of the present invention;
`FIGS. 5A, B, C and D are a flowchart describing a
`pulse rate measurement and compensation method in
`accordance with one embodiment of the present inven-
`tton;
`FIG. 6 is a table illustrating the effect of differential
`enhancement on raw sensor element signal levels;
`FIG. 7 is a table illustrating the calculations of the
`correlation algorithm; and
`
`0012
`
`5,243,992
`
`4
`FIGS. 8A and B are graphical representations of
`normal and inverted blood pressure waveforms, respec-
`tively.
`DESCRIPTION OF THE PREFERRED
`EMBODIMENTS
`
`10
`
`15
`
`20
`
`25
`
`30
`
`35
`
`45
`
`55
`
`65
`
`Referring to FIG. 1, a pulse rate monitor 1 in accor-
`dance with the preferred embodiment of the present
`invention is shown having a tonometer sensor 2 and a
`case 3 housing processing and display circuitry. The
`tonometer sensor 2 is mounted in a gimbal assembly 9
`whichinturnisconnectedtocaae3byaspring4.
`Spring 4 includes a sensor position adjustment 5 section
`at the connection point between gimbal assembly 9 and
`spring 4.
`Further details of this arrangement are shown in FIG.
`2, which uses the same component designations found in
`FIG. 1, where possible. The tonometer sensor 2 is at-
`tsched to a sensor adapter 12 which; in turn is pinned to
`the gimbal assembly 9 at axis number 2. A spring mount-
`ing pad 11 is pirmed along axis number 1 of gimbal
`assembly 9. mounting pad 11 being the point at which
`sensor position adjustment 5 connects to gimbal assem-
`bly 9. Flexible printed circuit ll] connects the tonometer
`sensor 2 to circuitry (shown in FIG. 4) in case 3.
`Tonorneter senmr 2 is an array of pressure or force
`sensitive elements fabricated into a single structure.
`Standard photolithographic manufacturing techniques
`can be used to construct the tonometer sensor. Experi-
`mental testing indicates that about 3 to I5 individual
`sensor elements are necessary for good accuracy in
`pulse rate measurement but even a single sensor element
`adapted for use with the method and apparatus of the
`present invention will produce more accurate pulse rate
`determinations.
`
`Referring again to FIG. 1., the remainder of the pulse
`rate monitor 1 will be described in terms of wearing the
`pulse rate monitor 1. In the preferred embodiment, the
`pulse rate monitor is worn on the operator's wrist like a
`wrist watch. When the wearer dons the pulse rate moni-
`tor 1, the tonometer sensor 2 is positioned above a radial
`artery and the case 3 is positioned on the opposite side
`of the wrist. The case 3, which anchors one end of
`spring 4, is held in place by strap 8 by cinching an flexi-
`ble portion 8a of a strap 8 and locking the strap in posi-
`tion by means of a buckle T. Strap 8 is prevented from
`directly contacting the tonometer sensor 2. the gimbal
`assembly 9 and the spring 4 by a protective band 6
`which is part of the strap 8. Protective band 6 has a box
`shaped cutout section which allows it to fit around the
`tonometer sensor 2,
`the gimbal assembly 9 and the
`spring 4. without contacting any of these elements.
`The tonometer sensor 2 is held against the artery with
`a thin cantilever spring 4 which is attached to the case
`3. The spring reaches around the wrist to position the
`tonometer sensor 2. A low profile gimbal assembly 9
`(See FIGS. 1 and 2) connects the tonometer sensor 2 to
`the spring It and allows the tonometer sensor 2 to pivot
`so as to lie flat against the wrist while being worn. Gim-
`bal assembly 9 allows about 20 degrees of rotation about
`each of its two axes. The position on spring -I of the
`tonometer sensor-gimbal assembly 2. 9 is adjustable by
`means of the sensor position adjustment 5 section of
`spring 4.
`Sensor hold down force must be controlled to pro-
`vide enough pressure to partially flatten the radial ar-
`tery but not enough pressure to cause discomfort to the
`wearer. The optimum hold down force is unique to
`
`0012
`
`
`
`5,243,992
`
`6
`element signals, counts. flag settings and calculated
`values. A ROM 72 contains the operating program
`software for CPU 68. CPU 68, through 1/0 74. receives
`digital data from ADC 66 and outputs pulse rate calcu-
`lations to display 76.
`The operation of CPU 68 is best understood by refer-
`ring to the flow chart of FIGS. 5A, B, C and D.
`When the pulse rate sensor system is turned on or
`reset by the wearer, the system goes through an initial-
`ization routine to clear stored data. At this point. a
`clock (not shown) starts and a clock signal is provided
`to the CPU 68 to set a predetermined sampling period.
`Referring to FIG. 5A, while performing segment 1 of
`the program, sub-step la checks to see if a cloclt ticl-' has
`occurred. If a clock tick is not detected at sub-step la.
`the cycle repeats. If a clock tick is detected, the CPU 68
`executes sub-step lb and samples the outputs of all sen-
`sor elements of the tonometer sensor 2. In addition,
`system timers are updated. Control then passes to pro-
`gram segment 2.
`During program segment 2, each data element for
`each sensor element is checked to see if it exceeds a
`predetermined maximum value. If an actual value ex-
`ceeds the predetermined maximum value, i.e. the prede-
`termined rnaxirnurn amplitude. the data from that sensor
`element is set to zero. The program also sets a flag so
`that all data from that sensor element is set to zero for
`the next five seconds. This eliminates unusually large
`signals which are usually noise. Program control then
`passes to program segment 3.
`Differential enhancement of the signal occurs during
`execution of program segment 3. During sub-step 3a. all
`signals from sensor elements not previously set to zero
`are averaged to produce a spatially weighted average
`signal. In the preferred embodiment. the weighted fac-
`tor is about 1.0, but can be adjusted as described below.
`At sub-step 3b. the spatially weighted average signal is
`subtracted from each of the actual sensor element sig-
`nals.
`
`The differential enhancement algorithm reduces mo-
`tion artifacts in the tonometer sensor signals. This algo-
`rithm can only be used when multiple sensor elements
`are employed simultaneously. Since all of the sensor
`element signals are often affected equally by a motion
`artifact, the differential enhancement algorithm aids in
`distinguishing between artifacts and blood pressure
`signals. For example, motion artifacts such as footsteps
`usually affect all sensor elements in the same way.
`Blood pressure signals, on the other hand, affect only
`sensor elements which are directly over or very near to
`the artery.
`Differential enhancement adds all of the signals from
`all of the sensor elements together and forms a spatially
`weighted average signal. If each sensor element is af-
`fected in the same way by a motion,
`the spatially
`weighted average signal will be an accurate representa-
`tion of the motion artifact. This spatially weighted aver-
`age signal is then subtracted from each individual sensor
`element signal to form a differentially enhanced signal
`for each sensor element signal. (I.e. an approximation of
`the motion artifact is subtracted from the raw sensor
`element signals to produce corrected sensor element
`signals.) The raw sensor element signals are expected to
`be either motion artifact signals or motion artifact sig-
`nals plus blood pressure signals. For example. FIG. 6
`shows the effect of the differential enhancement algo-
`rithm on the signals from a three element sensor array.
`Signals for only one clock tick are shown.
`
`5
`each individual but ranges from about 100 grams to
`about 500 grams. Some wearers may require a higher
`hold down force to obtain a reliable pulse signal from
`the tonorneter sensor 2. Other wearers may be sensitive
`to the hold down force of the tonometer sensor 2
`against their wrists and desire the lowest possible hold
`down force.
`
`Ideally, since the hold down force of the tononteter
`sensor 2 is controlled only by the deflection of the
`spring 4 and since the size and shape of the wrist can
`vary greatly between individuals. each spring 4 would
`have to be custom fit for each wearer. In general. three
`configurations for spring 4, shown in FIGS. 3A, B and
`C. will suffice to cover the majority of the population.
`Springs «IA, B and C. shown in FIGS. 3A, B and C,
`provide adequate length and curvature adjustment to
`cover the general population. The springs 4A, 4B and
`4C differ only in their radius of curvature.
`The discussion below of the computation of the pulse
`rate can best be understood by first understanding the
`important features of a normal blood pressure wave-
`form. Referring to FIG. 8A a normal blood pressure
`waveform ‘[00 with an average blood pressure 102 is
`shown. The point 104 on the waveform 100 where
`blood pressure is maximum is referred to as systole.
`while the point 106 where blood pressure is a minimum
`is referred to as diastole. It is known from scientific
`
`studies that the maximum absolute value of the slope of
`wave form 10!} occurs just prior to systole, in region
`108, when the blood pressure is rising from diastole. A
`dichrotic notch 110 is present in the blood pressure
`waveform 100 of some subjects. Referring to FIG. 8B,
`an inverted waveform 100’
`is shown having systole
`(104), diastole (106),
`the region (108) and dichrotic
`notch (110).
`Computation of pulse rate is performed by electronic
`circuitry shown in FIG. 4 and located in case 3,
`in
`accordance with the flow chart shown in FIGS. 5A, B,
`C and D.
`Referring first to FIG. 4. circuitry to process tcnome-
`ter signals and compute pulse rate in accordance with
`the present invention comprises: a preamplifier 58; a
`high pass filter 60; a low pass filter 62; an amplifier 64;
`an analog-to~digital digital converter (ADC) 66; a cen-
`tral processing unit (CPU) 68; random access memory
`(RAM) 70; read only memory (ROM) 72; input/output
`unit (I/O) 74 and display 76. As shown in FIG. 4, the
`preamplifier 58. filters 60 and 62. amplifier 64 and ADC
`66 each include inputs corresponding to the individual
`sensor elements.
`
`During operation. the output of all sensor elements of
`tonometer sensor 2 are routed to preamplifier 58 which
`amplifies the signals before filtering. The output of‘ pre-
`amplifier 58 is the input to high pass filter 60 which
`removes the DC and very low frequency components
`of the signals. The output of the high pass filter 60 is the
`input to the low pass filter 62 which removes sortie of
`the high frequency components of the signal. (High pass
`filter 60 coupled with low pass filter 62 effectively act as
`a band pass filter.) The preferred bandwidth of the
`effective band pass filer is about 0.l—30 Hz. After filter-
`ing, the signals are amplified by amplifier 64 to a level
`compatible with the ADC 66. ADC 66 multiplexes the
`signals and converts them to digital data which is sent to
`I/O 1'4. CPU 68 reads the digital data form I/O '74 and
`stores the digital data in RAM 70.
`RAM 70 is segmented to form a plurality of data
`buffers for storage of digital data representing sensor
`
`10
`
`20
`
`25
`
`35
`
`45
`
`50
`
`55
`
`65
`
`0013
`
`0013
`
`
`
`5,243,992
`
`5
`
`10
`
`I5
`
`20
`
`25
`
`30
`
`35
`
`45
`
`55
`
`3
`as valid pulses, but rejects or ignores pulses that have
`correlation coefficients outside this predetermined
`range. Of course, different acceptance ranges for COR
`may be used without departing from the teachings of
`this disclosure.
`The traditional mathematical definition of the corre-
`lation coefficient applies only when the two waveforms
`being correlated (CORELA and CORELB of the
`above example) have exactly the same duration. In
`other words, Eq. (I) is mathematically rigorous only
`when both heartbeats’ durations are equal. The present
`invention departs from this mathematical constraint and
`allows COR to be computed even iftlie durations of the
`two heartbeats are not exactly equal. Eq. (I) can be used
`to advantage in the case where the two waveforms do
`not have equal duration since Eq. (1) can be used very
`effectively to recognize movement artifacts even
`though it is not used in a mathematically rigorous way.
`In the preferred embodiment, the array. CORELB. of
`the above illustration is not used. Instead, a running
`summation of the numerator of Eq. (1), called NSUM.
`and of the denominator of Eq. (1), called DSUM, are
`updated after each sample of sensor data is obtained.
`Specifically. if BP(ELEMENT) is the pressure mea-
`surcd by a selected element at time ICOR, the following
`sequence is executed:
`
`10 DSUM = osvu + (aP(ELEMt-:N1') >( aJ=rEi.Em-:N'rn
`
`20 NSUM = MUM + EFIELEMENT) x C‘ORELA(.lCOR)
`
`so com-:r..4ucoR) = BHELEMENT}
`
`The correlation coefficient is then simply calculated
`after systole is found as COR.=NSUM/DSUM. This
`procedure works because in line 20 of the above code
`CORELAUCOR) still contains the pressure data from
`the previous heartbeat. Line 30 updates CORELA for
`the calculation on the next (future) heartbeat. Use of
`this running summation (in place of the two arrays,
`CORELA and CORELB, of Eq. (1)) is advantageous
`because it reduces computer speed and memory re-
`quirements. However. Eq. (1) defines the correlation
`coefficient so it may be calculated by other procedures
`without departing from the teachings of this disclosure.
`The table shown in FIG. '1' gives an example of what
`the relevant variables hold for a few samples of hypo-
`thetical data. If systole occurred at ICOR =3 (in reality
`the waveform normally extends for many more samples
`before a systole is found) the correlation coefficient
`would be calculated as
`
`COR =NSUM'x‘DSUM= 385t}.!'5325==O.T23 I. Feed
`
`and this heartbeat would be accepted as valid by the
`acceptance criterion described above.
`_
`
`Finally. the correlation coefficient in other applica-
`tions is sometimes also defined as:
`
`COR =(COR£LA U)‘ CORELK,-)/(AMPA)a(AMPB)
`
`7
`Differential enhancement is part of a larger class of
`algorithms where the processed output for each value is
`a weighted sum of the unprocessed values. For exam-
`ple. differential enhancement of three raw values pro-
`duces a processed output for value 1 as follows:
`
`Processed Value I = (2.r‘J)(raw value 1) +
`{-1/3)(nw value 1) + (-— U3-}(rsw value 3).
`
`Other extensions of this procedure are possible and
`different weighting factors than those used by the basic
`differential enhancement algorithm are possible without
`departing from the teachings of this disclosure. For
`example. it may be advantageous to assign large nega-
`tive weights to sensor element signals far from a se-
`lected sensor element and large positive weights to
`sensor element signals located near to a selected sensor
`element.
`Returning to FIG. 5A, after performing program
`segment 3. program. control passes to program segment
`4 where the differentially enhanced signal data is pro-
`cessed by a correlation algorithm.
`The correlation algorithm computes a quantitative
`measure of the similarity between sensor waveforms
`over a predetermined time period of about two consec-
`utive heartbeats. In principle, the shape of a person's
`blood pressure waveform will be fairly constant from
`one heartbeat to the next. The correlation algorithm
`compares the blood pressure waveform for the current
`heartbeat with a previously recorded waveform for the
`preceding heartbeat.
`The processor executes the correlation algorithm as
`follows. A variable. ICOR. is used as a counter to keep
`track of the number of clock "ticks" (i.e. the elapsed
`time) that the processor has spent looking for a systole.
`At the start. ICOR is set equal to 1. ICOR is then incre-
`mented by one for each subsequent clock tick (i.e. each
`time data is sampled from the sensor) When ICOR is set
`equal to I. the variables DSUM and NSUM. described
`below. are both set equal to zero. The data from the
`selected element of the tonometer sensor between the
`time when the processor starts looking for a systole and
`the time when it actually finds a (presumed valid) sys-
`tole is stored in an array. called CORELA. ICOR is
`used as a pointer for the array CORELA. i.e. the value
`of the data from the selected element after ICOR clock
`ticks have occurred (after the start of looking for a
`systole) is stored in CORELAOCOR).
`For example, suppose a similar procedure had been
`used for the previous heartbeat. and the values from the
`selected element had been stored in another array, C0-
`RELB. After the systole is found for the present heart-
`beat, the correlation coefficient, COR. may be calcu-
`lated. The correlation coefficient between the pressure
`wavefonns for the current and previous heartbeats is
`mathematically defined as:
`
`00351.4 ‘CORE
`_.
`COR ‘-’ (E‘5*i£t..4g}-Eon£LA(;))
`
`U}
`
`is over all values of elapsed
`.
`where the summation,
`time, j. from the start of looking for the systole. If the
`current heartbeat is exactly identical to the previous
`heartbeat, the array. CORELB, will be equal to the
`array. CORELA, and the correlation coefficient, COR.
`will be equal to 1. If COR differs greatly from 1, the
`two waveforms are not similar and at least one of them
`is probably distorted by a movement artifact. The pro-
`gram accepts values of COR between about 0.6 and 2.0
`
`where AMPA=sqrt(DSUM). AMPB=si:p-t(DSUM'),
`and DSUM'
`is the value of DSUM for the previous
`heartbeat. This defmition makes the correlation coeffi-
`cient independent of any overall gain change between
`one heartbeat and the next.
`
`65
`
`After executing sub-step «la. the program checks to
`see if a systolic peak has been detected within a. prede-
`
`0014
`
`0014
`
`
`
`5,243,992
`
`9
`tennined. number of sampling periods. e.g., the last 0.25
`seconds. at sub-step Cb. If the answer is yes, program
`segment 4 loops back to program segment 1 to accumu-
`late sensor element signals. If the answer is no, program
`control passes to program segment 5.
`A normal, non-inverted pulse waveform will have a
`fairly large negative slope as the blood pressure drops
`from its peak value at systole. There can also be a pro-
`nounced dichrotic notch which has a local minimum
`pressure. See FIG. BA. These features can be mistak-
`enly interpreted by the program as an inverted wave-
`form (i.e. a waveform containing a large negative slope
`followed by a local minimum). Similar problems can
`occur with a true inverted waveform, with the dich-
`rotic notch just after systole being mistakenly inter-
`preted as a new, non—inverted pulse. See FIG. 813. Since
`the program must be