`Eckerle et al.
`
`ANA01
`
`US005243992A
`[11] Patent Number:
`[45] Date of Patent:
`
`5,243,992
`Sep. 14, 1993
`
`[54] PULSE RATE SENSOR SYSTEM
`[75]
`Inventors:
`Joseph S, Eckerle, Redwood City;
`Dale W. Ploeger, San Francisco;
`Steven T. Holmes, Palo Alto; Thomas
`P, Low, Woodside; Rudolf Elbrecht,
`Los Altos; Philip R. Jeuck, III,
`Menlo Park; Ronald E, Pelrine,
`Menlo Park; Victor T. Newton, Jr.,
`Menlo Park, all of Calif.
`
`[73] Assignee: Colin Electronics Co., Ltd., Aichi,
`Japan
`[21] Appl. No.: 502,028
`[22] Filed:
`Mar, 30, 1990
`
`4,239,048 12/1980 Steuer srcssecccececcenecee 128/666
`. 128/687
`4,307,728 12/1981 Walton..
`
`4,353,372 10/1982 Ayer......
`128/640
`
`4,409,983 10/1983 Albert .......
`128/690
`6/1984 Hirano etal. .
`364/417
`4,456,959
`
`3/1987 Berger etal. .
`4,646,749
`128/678
`
`4,667,680 5/1987 Ellis ......
`128/672
`364/417
`4,712,179 12/1987 Heimer..
`4,735,213 4/1988 Shirasaki
`128/681
`4,799,491
`1/1989 Eckerle.....
`«. 128/672
`4,802,488 2/1989 Eckerle
`«128/372
`4,836,213
`6/1989 Wenzel et al. eccsccseneesenene 128/672
`
`
`
`Primary Examiner—Lee S. Cohen
`Assistant Examiner—Marianne H. Parker
`Attorney, Agent, or Firm—Oliff & Berridge
`
`[57]
`ABSTRACT
`(SU)
`Tent CUS ssssssscccs ccccasceasarecasiqveanstseansieencts A61B 5/02
`
`B52) RISaccerssnsseversnnanersnvesnyneansys 128/690; 128/672
`A pulse rate sensor system is packaged in a wristwatch
`[58] Field of Search............... 128/687, 722, 699, 689,
`sized assembly and is worn by the user to provide an
`128/677, 681, 683, 690
`accurate determination of pulse rate. A tonometer sen-
`sor is provided to detect heartbeat pressure waves pro-
`duced by a superficial artery. A microcomputer manip-
`ulates the unprocessed tonometer sensor elementsignals
`using multiple algorithms to determine an accurate
`pulse rate.
`
`[56]
`
`References Cited
`U.S. PATENT DOCUMENTS
`3,999,537 12/1976 Nooiles siu.ceccserescesseescesseeeccsees 128/687
`4,058,118 11/1977 Stupayetal.
`128/2.05 T
`4,086,916 5/1978 Freemanetal.
`+ 128/690
`
`we 128/689
`4,181,134
`1/1980 Mason et al.
`4,202,350 5/1980 Walton ......0...cee 128/690
`
`36 Claims, 9 Drawing Sheets
`
`0001
`
`Apple Inc.
`APL1047
`U.S. Patent No. 8,923,941
`FITBIT, Ex. 1047
`
`Apple Inc.
`APL1047
`U.S. Patent No. 8,923,941
`
`0001
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 1 of 9
`
`5,243,992
`
`FIG.|
`
`FITBIT, Ex. 1047
`
`0002
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 2 of 9
`
`5,243,992
`
`FLEXIBLE
`PC. BOARD
`
`
`
`
`SENSOR 2
`
`
`___AXIS #2 9a
`
`GIMBAL 9
`ane #1
`MOUNT TO SPRING {1
`
`b
`
`FIG.2
`
` \
`
`FIG. 3C
`
`FIG. 3A
`
`FIG 3B
`
`0003
`
`FITBIT, Ex. 1047
`
`0003
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 3 of 9
`
`5,243,992
`
`SOIVNY
`
`
`
`Wiigig(==4yg4yangEeaF\ssvdEE]SSdF<hi||MO[iin=fHOIHea
`
`YSLYIANOD
`
`AV1dS10
`
`92
`
`0004
`
`FITBIT, Ex. 1047
`
`0004
`
`FITBIT, Ex. 1047
`
`
`
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 4 of 9
`
`5,243,992
`
`
`
`cannaA|
`
`
`
`
`
`
`INITIALIZE
`
`CLOCK TICK
`OCCURED
`
`Y
`UPDATE TIMERS
`
`|
`
`| |
`
`|
`
`DATA
`
`
`
`AVERAGE ALL
`
`ce ee ee ee eee
`
`
`ACCUMULATE
`
`CORRELATION
`SUMMATIONS
`
`
`
`0005
`
`FITBIT, Ex. 1047
`
`0005
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 5 of 9
`
`5,243,992
`
`y
`
`DECREASE
`SLOPE
`THRESHOLD BY 1/2
`
`RESETIC
`
`SL
`
`SLOP
`AX SLOPE
`
`ve”
`
`6d
`
`aueoument
`MAX SLOPE
`
`
`
`§j "6k
`
`TIMER TOO ©
`
`,
`
`FOUND SYSTOLIC SET] 6f
`SIGNAL VALID SET
`SYSTOLIC FOUND
`
`SLOPE THRESHOLD =
`
`serene’ f
`
`SAVE SYSTOLIC
`AND DIASTOLIC
`
`|,
`
`SAVE SLOPE
`AS MAX SLOPE
`
`DOWNHILL FLAG
`
`}-6N
`
`‘
`
`SET FLAG THAT
`SLOPE THRESHOLD| &P
`eaeeree
`
`RESET ABSOLUTE
`MINIMUM AND
`)
`MAXIMUM STORAGE
`eee ee eeee5
`
`6!
`
`0006
`
`FITBIT, Ex. 1047
`
`0006
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 6 of 9
`
`5,243,992
`
`Proto rs =
`GMENT
`@)
` SERMENT
`!
`|
`|
`
`||
`
`i
`|
`
`|
`
`|
`
`||
`
`|
`I
`|
`
`(PULSE AMPLITUDE)
`
`AVAILABLE
`
`lI
`
`|
`|
`|
`
`
`
`SET ELEMENT OF
`INTEREST TO ELEMENT |
`
`
`WITH MOST MAXIMUMS
`
`
`INLAST 5 TRIGGERS
`
`
`eere es: Se ae ee
`
`|
`!
`
`||D
`
`WTMAM170 CeSEGMENT#8|
`
`ont Ba
`ETERMINE ELEMENT|=|
`PERIOD >0.7
`MUM
`WITH MAXI
`Tb
`POSITIVE PULSE
`AMPLITUDE IF
`
`
`
`| | | | | | || | | | |1 | | | | | ||| | | | | |
`
`| ||D
`
`RRENT P
`
`+ PREVIOUS PERIOD
`
`\
`
`:
`
`1
`
`\
`
`l
`
`
`
`It
`
`i
`1
`|
`I
`
`|
`
`INCREASE C TO A
`
`|
`
`B
`
`a
`
`|
`SIGNAL]
`acu oe
`VALID= Tj
`FALSE]
`
`e Sst StepsSaas a
`SAVE CURRENT
`on
`SYSTOLIC PERIOD
`=
`RAGE)
`|
`[X4-C)+(THISPULSE)| 190
`L—-~~-~-]-------- -
`Priereeececey mri 7
`
`ee
`:
`is
`|
`EWAVERAG
`
`
`
`
`
`| Naor as|VDToto averace -®) |
`
`
`!
`Fo ern ier7
`I
`|
`1
`!
`1
`!
`FIG. 5C
`|
`
`DECREASE C TOA
`
`{0c
`
`I |
`
`[MINIMUM OF0
`)
`hosiceeces taser aia aceFe|
`
`|
`
`|
`
`0007
`
`FITBIT, Ex. 1047
`
`0007
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 7 of 9
`
`5,243,992
`
`oo Fe eee 5
`|
`$ i!
`!
`wee)
`|
`1b
`Bosse cee
`ea ee ee =
`roa SS SEMEN
`|
`[CALCULATE
`~oee—|'
`[COEFFICIENT
`7
`(2c
`12b
`I
`
`|
`SIGNAL|
`0.6<
`|!
`SEGMENT
`
` <2.0? FALSE| |#14 >>*-elyaLip=) |—COEFFICIENT (8)}
`
`
`
`
`
`
`
`iL ___Moc gee cee J
`14b
`{4a
` VALID
`SETa
`|
`ca
`CALCULATE
`PULSEIN LAST
`PLAY ||!
`30SEC.? = i!|PULSERATE
`
`GMENT
`
`N
`
`{SIGNAL
`
`
`
`SN
`
`|1\
`
`—-—-—ee
`
`[SET POST SYSTOLIC
`|
`Lonm -1 DELAY TO 025
`OF PULSE TIME
`
`I3b
`
`t3c
`
`| | :
`
`RESET SYSTOLIC
`| TIMER
`
`| |
`
`UPDATE DISPLAY
`[13d
`|IF VALID PULSE FOUND
`|
`Cae eee eee 4
`
`0008
`
`FITBIT, Ex. 1047
`
`0008
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 8 of 9
`
`5,243,992
`
`(Y091)¥73409WNSNwnsddé
`02002O0Ob02
`02oseeSeesOLS-OSESépS-
`
`
`“WNOISGJONVHN3“TWNOIS
`
`ATIVILNSYSISIO~GQ3I9OVYSAV
`SIO}
`Sl-0}0Ol
`2Old
`oOI4
`(YOdI)¥T3N09=NOD
`0S£0€-éOo}|(}-33D)
`AVY=LINSWI13
`G2€G-2Olt
`
`‘TWNOIS
`
`
`
`
`
`($30)
`
`
`
`
`
`0009
`
`FITBIT, Ex. 1047
`
`0009
`
`FITBIT, Ex. 1047
`
`
`
`U.S. Patent
`
`Sep. 14, 1993
`
`Sheet 9 of 9
`
`5,243,992
`
`
`
`0010
`
`FITBIT, Ex. 1047
`
`0010
`
`FITBIT, Ex. 1047
`
`
`
`5,243,992
`
`1
`
`PULSE 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
`tonometer 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 meansof a novel
`pressure sensing array and multiple methodsfor identifi-
`cation and elimination of artifacts.
`Other methods and apparatus are known for measur-
`ing pulse rates 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 motionartifacts are disclosed in U.S. Pat. Nos.
`4,307,728, 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 enhancedsig-
`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 compensatefor
`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.
`Thepulse 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 element(s) 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. Theresult 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 twofalls 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,488. These principles are also de-
`scribed in several publications including anarticle enti-
`tled “Tonometry, Arterial,” in Volume 4 of the Ency-
`clopedia of Medical Devices and Instruments. (J. G.
`Webster, Editor, John Wiley & Sons, 1988). 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 ofthe 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 ofthe arterial wall is
`not necessary and a lower hold downforce can be used.
`This results in greater comfort for the wearer.
`SUMMARY OF THE INVENTION
`
`20
`
`40
`
`Accordingly, the present invention has been devel-
`oped to overcome the foregoing shortcomingsofexist-
`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
`50
`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 inventionis 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 elementsignals
`by subtracting a value from all these signals based on a
`spatially weighted average ofthese signals.
`An additional object of the present invention is to
`cancelout artifacts from a sensor element which exceed
`a level predetermined to be the maximum level of a
`valid blood pressure signal.
`
`60
`
`0011
`
`FITBIT, Ex. 1047
`
`0011
`
`FITBIT, Ex. 1047
`
`
`
`5,243,992
`
`3
`Still another object of the present invention is to
`accurately process inverted waveforms caused by mis-
`alignmentor 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 elementsignal corresponding to one of the sen-
`sor elements, the selected sensor elementsignal 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 atleast two of
`the pulse rates, each of the two pulse rates having the
`value within a predetermined 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 onat
`least one of the pressure sensing elements and produc-
`ing a plurality of sensor elementsignals, at least one of
`the sensor elementsignals 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 meansfor deter-
`mining a pulse rate based onat 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 meansat 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 arrangementofa 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 gimbal assembly;
`FIGS. 3A, B and C show low, medium and high
`curvature spring profiles, respectively.
`FIG.4 is a block diagram ofthe 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 embodimentofthe present inven-
`tion;
`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
`
`4
`FIGS. 8A and B are graphical representations of
`normal and inverted blood pressure waveforms, respec-
`tively.
`DESCRIPTION OF THE PREFERRED
`EMBODIMENTS
`
`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
`which in tum is connected to case 3 by a spring 4.
`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-
`tached 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 pinned 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 10 connects the tonometer
`sensor 2 to circuitry (shown in FIG. 4)in case 3.
`Tonometer sensor 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 6 individual
`sensor elements are necessary for good accuracy in
`pulse rate measurementbut 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 remainderofthe 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 donsthe pulse rate moni-
`tor 1, the tonometer sensor 2 is positioned abovea 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 8 of a strap 8 and locking the strap in posi-
`tion by meansof a buckle 7. Strap 8 is prevented from
`directly contacting the tonometer sensor 2, the gimbal
`assembly 9 and the spring 4 by a protective band 6
`whichis part of the strap 8. Protective band 6 has a box
`shaped cutout section whichallowsit 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 4 and allows the tonometersensor 2 to pivot
`so as to lie flat against the wrist while being worn. Gim-
`bal assembly 9 allows about 20 degrees ofrotation about
`each of its two axes. The position on spring 4 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 mustbe 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
`
`35
`
`45
`
`55
`
`65
`
`0012
`
`FITBIT, Ex. 1047
`
`0012
`
`FITBIT, Ex. 1047
`
`
`
`9,243,992
`
`_ A
`
`20
`
`5
`6
`each individual but ranges from about 100 grams to
`element signals, counts, flag settings and calculated
`about 500 grams. Some wearers may require a higher
`values. A ROM 72 contains the operating program
`hold down force to obtain a reliable pulse signal from
`software for CPU 68. CPU 68, through I/O 74,receives
`the tonometer sensor 2. Other wearers may be sensitive
`digital data from ADC 66 and outputs pulse rate calcu-
`to the hold down force of the tonometer sensor 2
`lations to display 76.
`against their wrists and desire the lowest possible hold
`The operation of CPU 68 is best understood byrefer-
`downforce.
`ring to the flow chart of FIGS. 5A, B, C and D.
`Ideally, since the hold down force of the tonometer
`When the pulse rate sensor system is turned on or
`sensor 2 is controlled only by the deflection of the
`reset by the wearer, the system goes through an initial-
`spring 4 and since the size and shape of the wrist can
`ization routine to clear stored data. At this point, a
`vary greatly between individuals, each spring 4 would
`clock (not shown)starts and a clock signal is provided
`have to be custom fit for each wearer. In general, three
`to the CPU 68 to set a predetermined sampling period.
`configurations for spring 4, shown in FIGS. 3A, B and
`Referring to FIG. 5A, while performing segment 1 of
`C, will suffice to cover the majority of the population.
`the program, sub-step la checks to see if a clock tick has
`Springs 4A, B and C, shown in FIGS. 3A, B and C,
`occurred. If a clock tick is not detected at sub-step la,
`provide adequate length and curvature adjustment to
`the cycle repeats. If a clock tick is detected, the CPU 68
`cover the general population. The springs 4A, 4B and
`executes sub-step 15 and samples the outputs ofall sen-
`4C differ only in their radius of curvature.
`sor elements of the tonometer sensor 2, In addition,
`The discussion below of the computation ofthe pulse
`system timers are updated. Control then passes to pro-
`gram segment 2.
`rate can best be understood by first understanding the
`important features of a normal blood pressure wave-
`During program segment 2, each data element for
`each sensor element is checked to see if it exceeds a
`form. Referring to FIG. 8A a normal blood pressure
`waveform 100 with an average blood pressure 102 is
`predetermined maximum value. If an actual value ex-
`shown. The point 104 on the waveform 100 where
`ceeds the predetermined maximumvalue,i.e. the prede-
`blood pressure is maximum is referred to as systole,
`termined maximum amplitude, the data from that sensor
`while the point 106 where blood pressure is a minimum
`element is set to zero. The program also sets a flag so
`is referred to as diastole. It is known from scientific
`that all data from that sensor elementis set to zero for
`studies that the maximum absolute value of the slope of
`the next five seconds. This eliminates unusually large
`wave form 100 occurs just prior to systole, in region
`signals which are usually noise. Program control then
`108, when the blood pressure is rising from diastole. A
`passes to program segment 3.
`dichrotic notch 110 is present in the blood pressure
`Differential enhancement ofthe signal occurs during
`waveform 100 of some subjects. Referring to FIG. 8B,
`execution of program segment3. During sub-step3a,all
`an inverted waveform 100'
`is shown having systole
`signals from sensor elements not previously set to zero
`(104), diastole (106),
`the region (108) and dichrotic
`are averaged to produce a spatially weighted average
`notch (110).
`signal. In the preferred embodiment, the weighted fac-
`Computation of pulse rate is performed by electronic
`tor is about 1.0, but can be adjusted as described below.
`circuitry shown in FIG. 4 and located in case 3,
`in
`At sub-step 34, the spatially weighted average signal is
`accordance with the flow chart shown in FIGS. 5A, B,
`subtracted from each of the actual sensor elementsig-
`nals.
`C and D.
`Referringfirst to FIG.4, circuitry to process tonome-
`Thedifferential enhancement algorithm reduces mo-
`ter signals and compute pulse rate in accordance with
`tion artifacts in the tonometersensorsignals. This algo-
`the present invention comprises: a preamplifier 58; a
`rithm can only be used when multiple sensor elements
`high passfilter 60; a low pass filter 62; an amplifier 64;
`are employed simultaneously. Since all of the sensor
`an analog-to-digital digital converter (ADC) 66; a cen-
`elementsignals are often affected equally by a motion
`tral processing unit (CPU) 68; random access memory
`artifact, the differential enhancement algorithm aids in
`(RAM)70; read only memory (ROM)72; input/output
`distinguishing between artifacts and blood pressure
`unit (I/O) 74 and display 76. As shown in FIG.4, the
`signals. For example, motionartifacts such as footsteps
`preamplifier 58, filters 60 and 62, amplifier 64 and ADC
`usually affect all sensor elements in the same way.
`66 each include inputs corresponding to the individual
`Blood pressure signals, on the other hand,affect only
`sensor elements.
`sensor elements which are directly over or very near to
`the artery.
`During operation, the outputofall sensor elements of
`tonometer sensor 2 are routed to preamplifier 58 which
`Differential enhancement addsall of the signals from
`amplifies the signals before filtering. The output of pre-
`all of the sensor elements together and formsa spatially
`amplifier 58 is the input to high pass filter 60 which
`weighted average signal. If each sensor elementis af-
`removes the DC and very low frequency components
`fected in the same way by a motion, the spatially
`of the signals. The outputof the high passfilter 60 is the
`weighted average signal will be an accurate representa-
`input to the low pass filter 62 which removes some of
`tion of the motion artifact. This spatially weighted aver-
`age signal is then subtracted from eachindividual sensor
`the high frequency components ofthe signal. (High pass
`filter 60 coupled with low pass filter 62 effectively act as
`elementsignal to form a differentially enhanced signal
`a band pass filter.) The preferred bandwidth of the
`for each sensor elementsignal. (I.e. an approximation of
`the motion artifact is subtracted from the raw sensor
`effective band pass filer is about 0.1-30 Hz. After filter-
`ing, the signals are amplified by amplifier 64 to a level
`element signals to produce corrected sensor element
`compatible with the ADC 66. ADC 66 multiplexes the
`signals.) The raw sensor element signals are expected to
`signals and converts them to digital data whichis sent to
`be either motion artifact signals or motionartifact sig-
`I/O 74. CPU 68 reads the digital data form I/O 74 and
`nals plus blood pressure signals. For example, FIG. 6
`stores the digital data in RAM 70.
`shows the effect of the differential enhancement algo-
`RAM 70 is segmented to formaplurality of data
`rithm onthe signals from a three elementsensorarray.
`Signals for only one clock tick are shown.
`buffers for storage of digital data representing sensor
`
`30
`
`35
`
`40
`
`60
`
`nantA
`
`0013
`
`FITBIT, Ex. 1047
`
`0013
`
`FITBIT, Ex. 1047
`
`
`
`7
`Differential enhancementis part of a larger class of
`algorithms where the processed output for each valueis
`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 1 = (2/3\raw value 1) +
`(—1/3)(raw value 2) + (—1/3)raw value 3).
`
`Other extensions of this procedure are possible and
`different weighting factors than those used by the basic
`differential enhancementalgorithm 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 elementsignals located near to a selected sensor
`element.
`Returning to FIG. 5A, after performing program
`segment3, 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 numberof clock “ticks” (i.e. the elapsed
`time) that the processor has spent looking for a systole.
`At the start, ICORis 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 ICORis set
`equal to 1, 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 whenthe processorstarts looking for a systole and
`the time when it actually finds a (presumed valid) sys-
`tole is stored in an array, called CORELA.
`ICORis
`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 CORELA(ICOR).
`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, CO-
`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
`waveforms for the current and previous heartbeats is
`mathematically defined as:
`
`.
`(CORELA(/)*CORE.
`COR =F (SORELCORELA—
`
`a)
`
`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 waveformsare not similar and at least one of them
`is probably distorted by a movementartifact. The pro-
`gram accepts values of COR between about0.6 and 2.0
`
`3,243,992
`
`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 ofthe
`above example) have exactly the same duration. In
`other words, Eg. (1) 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 if the durationsofthe
`two heartbeats are not exactly equal. Eq. (1) 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
`thoughit is not used in a mathematically rigorous way.
`In the preferred embodiment, the array, CORELB,of
`the aboveillustration is not used. Instead, a running
`summation of the numerator of Eg. (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-
`sured by a selected element at time ICOR,the following
`sequenceis executed:
`
`30
`
`la5
`
`a iy
`
`55
`
`60
`
`10 DSUM = DSUM + (BP(ELEMENT) x BAELEMENT))
`
`20 NSUM = NSUM + BRELEMENT) x CORELA(ICOR)
`3 CORELA(ICOR) = BRELEMENT)
`
`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
`CORELA(ICOR)still contains the pressure data from
`the previous heartbeat. Line 30 updates CORELAfor
`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 memoryre-
`quirements. However, Eg. (1) defines the correlation
`coefficient so it may be calculated by other procedures
`without departing from the teachingsofthis disclosure.
`The table shown in FIG, 7 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/DSUM= 3850/5325=0.723 L 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 =(CORELA(j/)*CORELB(/)/(AMPA)x(AMPB)
`
`where AMPA=sqrt(DSUM), AMPB=sqrt(DSUM’)),
`and DSUM'is the value of DSUM for the previous
`heartbeat. This definition makes the correlation coeffi-
`cient independent of any overall gain change between
`one heartbeat and the next.
`After executing sub-step 4a, the program checks to
`see if a systolic peak has been detected within a prede-
`
`0014
`
`FITBIT, Ex. 1047
`
`0014
`
`FITBIT, Ex. 1047
`
`
`
`9
`termined. number of sampling periods, e.g., the last 0.25
`seconds, at sub-step 44. If the answer is yes, program
`segment4 loops back to program segment1 to accumu-
`late sensor elementsignals. If the answer is no, program
`control passes to program segment5,
`A normal, non-inverted pulse waveform will have a
`fairly large negative slope as the blood pressure dr