US20030195404A1  Method and apparatus for realtime control of physiological parameters  Google Patents
Method and apparatus for realtime control of physiological parameters Download PDFInfo
 Publication number
 US20030195404A1 US20030195404A1 US10417438 US41743803A US2003195404A1 US 20030195404 A1 US20030195404 A1 US 20030195404A1 US 10417438 US10417438 US 10417438 US 41743803 A US41743803 A US 41743803A US 2003195404 A1 US2003195404 A1 US 2003195404A1
 Authority
 US
 Grant status
 Application
 Patent type
 Prior art keywords
 time
 glucose
 estimator
 embodiment
 optimal
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Abandoned
Links
Images
Classifications

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
 A61B5/00—Detecting, measuring or recording for diagnostic purposes; Identification of persons
 A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
 A61B5/00—Detecting, measuring or recording for diagnostic purposes; Identification of persons
 A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
 A61B5/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
 A61B5/00—Detecting, measuring or recording for diagnostic purposes; Identification of persons
 A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
 A61B5/7235—Details of waveform analysis
 A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F19/00—Digital computing or data processing equipment or methods, specially adapted for specific applications
 G06F19/30—Medical informatics, i.e. computerbased analysis or dissemination of patient or disease data
 G06F19/34—Computerassisted medical diagnosis or treatment, e.g. computerised prescription or delivery of medication or diets, computerised local control of medical devices, medical expert systems or telemedicine
 G06F19/3456—Computerassisted prescription or delivery of medication, e.g. prescription filling or compliance checking

 G—PHYSICS
 G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
 G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
 G16H10/00—ICT specially adapted for the handling or processing of patientrelated medical or healthcare data
 G16H10/40—ICT specially adapted for the handling or processing of patientrelated medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis

 G—PHYSICS
 G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
 G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
 G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
 G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
 A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
 A61B2560/02—Operational features
 A61B2560/0242—Operational features adapted to measure environmental factors, e.g. temperature, pollution

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
 A61M2230/00—Measuring parameters of the user
 A61M2230/20—Blood composition characteristics
 A61M2230/201—Glucose concentration

 A—HUMAN NECESSITIES
 A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
 A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
 A61M5/00—Devices for bringing media into the body in a subcutaneous, intravascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, armrests
 A61M5/14—Infusion devices, e.g. infusing by gravity; Blood infusion; Accessories therefor
 A61M5/142—Pressure infusion, e.g. using pumps

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06F—ELECTRIC DIGITAL DATA PROCESSING
 G06F19/00—Digital computing or data processing equipment or methods, specially adapted for specific applications
 G06F19/30—Medical informatics, i.e. computerbased analysis or dissemination of patient or disease data
 G06F19/34—Computerassisted medical diagnosis or treatment, e.g. computerised prescription or delivery of medication or diets, computerised local control of medical devices, medical expert systems or telemedicine
 G06F19/3481—Computerassisted prescription or delivery of treatment by physical action, e.g. surgery or physical exercise

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
 Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
 Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
 Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change.
 Y02A90/20—Information and communication technologies [ICT] supporting adaptation to climate change. specially adapted for the handling or processing of medical or healthcare data, relating to climate change
 Y02A90/26—Information and communication technologies [ICT] supporting adaptation to climate change. specially adapted for the handling or processing of medical or healthcare data, relating to climate change for diagnosis or treatment, for medical simulation or for handling medical devices

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
 Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y10S128/00—Surgery
 Y10S128/92—Computer assisted medical diagnostics
 Y10S128/923—Computer assisted medical diagnostics by comparison of patient data to other data

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
 Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y10S128/00—Surgery
 Y10S128/92—Computer assisted medical diagnostics
 Y10S128/923—Computer assisted medical diagnostics by comparison of patient data to other data
 Y10S128/924—Computer assisted medical diagnostics by comparison of patient data to other data using artificial intelligence
Abstract
A realtime controller operating as an artificial pancreas uses a Kalman control algorithm to control glucose level of a patient in real time. The realtime controller receives an estimate of the patient glucose level and a reference glucose level. The estimate of the patient glucose level can be provided by an optimal estimator implemented using a linearized Kalman filter. The estimated glucose level and the reference glucose level are processed by the Kalman control algorithm to determine a control command in real time. The Kalman control algorithm has a dynamic process forced by the control command a cost function determining a relative level of control. The control command is provided to a dispenser which secretes insulin or glucagon in response to the control command to correct a relatively high glucose level or a relatively low glucose level.
Description
 [0001]This application is a continuation of U.S. application Ser. No. 09/960,855, filed on Sep. 21, 2001, which claims the benefit of U.S. Provisional Application No. 60/234,632, filed on Sep. 22, 2000, and entitled Real Time Estimation & Control Of Biological Process.
 [0002]A portion of the disclosure of this patent document contains material that is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or of the patent disclosure as it appears in the Patent and Trademark Office patent files or records, but otherwise reserves all copyright rights whatsoever.
 [0003]1. Field of the Invention
 [0004]This present invention relates generally to a method and apparatus for controlling physiological parameters and more particularly to an optimal controller for controlling glucose levels in a patient.
 [0005]2. Description of the Related Art
 [0006]Different types of sensors (e.g., optical sensors) are available for monitoring of physiological parameters (e.g., glucose concentration). Glucose monitoring is typically performed by people with diabetes mellitus which is a medical condition involving a body's inability to produce the quantity or quality of insulin needed to maintain a normal circulating blood glucose. Frequent monitoring of glucose is generally necessary to provide effective treatment and to prevent long term complications of diabetes (e.g., blindness, kidney failure, heart failure, etc.). New methods of monitoring glucose are fast, painless and convenient alternatives to the typical capillary blood glucose (CBG) measurements which involve finger pricks that are painful, inconvenient and difficult to perform for long term.
 [0007]Optical measurement of glucose is performed by focusing a beam of light onto the body. Optical sensors determine glucose concentration by analyzing optical signal changes in wavelength, polarization or intensity of light. However, many factors other than glucose concentration also contribute to the optical signal changes. For example, sensor characteristics (e.g., aging), environmental variations (e.g., changes in temperature, humidity, skin hydration, pH, etc.), and physiological variations (e.g., changes in tissue fluid due to activity, diet, medication or hormone fluctuations) affect sensor measurements.
 [0008]Various methods are used to improve the accuracy of the sensor measurements. One method (e.g., multivariate spectral analysis) utilizes calibration models developed by initially measuring known glucose concentrations to correct subsequent sensor measurements. The calibration models become inaccurate over time due to dynamic changes in physiological processes. Another method (e.g., adaptive noise canceling) utilizes signal processing to cancel portions of the sensor measurements unrelated to glucose concentration. For example, two substantially simultaneous sensor measurements at different wavelengths make up a composite signal which can be processed to cancel its unknown and erratic portions. However, many sensors do not provide substantially simultaneous measurements at two different wavelengths.
 [0009]The present invention solves these and other problems by providing a method and apparatus for making optimal estimates of a physiological parameter (e.g., glucose level), assessing reliability of the optimal estimates, and/or providing optimal control of the physiological parameter in real time using one or more sensor measurements at each measurement time epoch (or interval). The sensor measurements can be timebased (e.g., every five minutes) to provide continuous monitoring and/or regulation of the physiological parameter. The sensor measurements are a function of the physiological parameter within specified uncertainties.
 [0010]An optimal estimator provides an accurate estimate of glucose level in real time using a sensor with at least one output. In one embodiment, the optimal estimator is integrated with the sensor and an output display to be a compact glucose monitoring device which can be worn by a patient for continuous monitoring and realtime display of glucose level. In an alternate embodiment, the optimal estimator is a separate unit which can interface with different types of sensors and provide one or more outputs for display, further processing by another device, or storage on a memory device.
 [0011]In one embodiment, the optimal estimator employs a priori deterministic dynamic models developed with stochastic variables and uncertain parameters to make estimates of glucose level. For example, glucose level is defined as one of the stochastic (or random) variables. Dynamic mathematical models define process propagation (i.e., how physiological and sensor parameters change in time) and measurement relationship (i.e., how physiological and sensor parameters relate to environmental conditions). Environmental conditions (e.g., temperature, humidity, pH, patient activity, etc.) can be provided to the optimal estimator intermittently or periodically via environment sensors and/or data entries by a patient or a doctor.
 [0012]The optimal estimator uses dynamic models to propagate estimates of respective stochastic variables, error variances, and error covariances forward in time. At each measurement time epoch, the optimal estimator generates realtime estimates of the stochastic variables using one or more sensor outputs and any ancillary input related to environmental conditions. In one embodiment, the optimal estimator employs a linearized Kalman filter to perform optimal estimation of the stochastic variables (e.g., glucose level). In particular, an extended Kalman filter is used to accommodate nonlinear stochastic models.
 [0013]Before making realtime estimates, the optimal estimator is initialized by providing initial values for the stochastic variables, error variances, and error covariances. For example, a CBG measurement or another direct glucose measurement is performed at initialization to provide a starting value for the stochastic variable corresponding to glucose level.
 [0014]In one embodiment, the optimal estimator provides one or more outputs to a patient health monitor which is capable of optimized realtime decisions and displays. The patient health monitor evaluates system performance by assessing the performance of the sensor and/or optimal estimator in real time. For example, the patient health monitor applies statistical testing to determine the reliability of the realtime estimates of the stochastic variables by the optimal estimator. The statistical testing is performed in real time on residual errors of the optimal estimator to establish performance measures.
 [0015]In one embodiment, the patient health monitor acts as an input/output interface between the patient or medical staff (e.g., a doctor, nurse, or other healthcare provider) and the optimal estimator. For example, environmental conditions can be provided to the patient health monitor for forwarding to the optimal estimator. Optimal estimator outputs can be provided to the patient health monitor for display or forwarding to an external device (e.g., a computer or a data storage device).
 [0016]In one embodiment, the optimal estimator provides one or more outputs to an optimal controller which can regulate in real time the physiological parameter being monitored. For example, an optimal controller responds to realtime optimal estimator outputs and provides an output to operate an actuator. In the case of glucose control, the actuator can be a dispenser or a pump which secretes insulin to correct a relatively high glucose level and glucagon to correct a relatively low glucose level. The optimal controller takes advantage of a priori information regarding the statistical characteristics of the actuator and is able to control the output of the actuator to be within specified uncertainties.
 [0017]In one embodiment, the optimal estimator and the optimal controller form an optimal closedloop system. For example, a glucose sensor, an optimal estimator, an optimal controller, and an insulin/glucagon dispenser work together as an artificial pancreas to continuously regulate glucose level. The glucose sensor can be internal or external to a patient's body. The optimal controller provides a control feedback to the optimal estimator to account for delivery of the insulin/glucagon.
 [0018]The optimal closedloop system is effective in a variety of biomedical applications. For example, cardiovascular functions can be continuously regulated by using sensors to detect blood pressure, blood oxygen level, physical activity and the like, an optimal estimator to process the sensor measurements and make realtime estimates of heart function parameters, and an optimal controller to control operations of an artificial device (e.g., a pacemaker) in real time based on the realtime estimates from the optimal estimator to achieve a set of desired heart function parameter values. Other artificial devices (e.g., artificial limbs, bionic ears, and bionic eyes) can be part of similar optimal closedloop systems with sensors detecting nerve signals or other appropriate signals.
 [0019]The optimal closedloop system is also effective in optimal treatment of chronic illnesses (e.g., HIV). Some medications for treatment of chronic illnesses are relatively toxic to the body. Over delivery of medication generally has adverse effects on the patient. The optimal closedloop system is capable of providing effective and safe treatment for the patient. For example, an optimal estimator provides realtime estimates of key physiological parameters using one or more sensors, and an optimal controller controls a slow infusion of medication in real time based on the realtime estimates from the optimal estimator to obtain desirable values for the key physiological parameters.
 [0020]In one embodiment, the optimal estimator, patient health monitor, and optimal controller are software algorithms which can be implemented using respective microprocessors. New information regarding process propagation or measurement relationship can be easily incorporated by modifying, reconfiguring, and/or adding to the software algorithms. The optimal estimator, patient health monitor, and optimal controller can be implemented as one joint algorithm or separate respective algorithms which function together to provide an optimal closedloop system.
 [0021][0021]FIG. 1A is a block diagram of one embodiment of an estimator.
 [0022][0022]FIG. 1B is a block diagram of another embodiment of a glucosemonitoring device.
 [0023][0023]FIG. 2 illustrates one embodiment of modeling physiological processes in a linearized Kalman filter application.
 [0024][0024]FIG. 3A describes one embodiment of an estimation function which depicts a linearized Kalman filter formulation.
 [0025][0025]FIG. 3B illustrates an initialization in one embodiment of an optimal estimator.
 [0026][0026]FIG. 3C illustrates a timeupdate cycle in accordance with one embodiment of an optimal estimator.
 [0027][0027]FIG. 3D illustrates a measurementupdate cycle in accordance with one embodiment of an optimal estimator.
 [0028][0028]FIGS. 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
 [0029][0029]FIGS. 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
 [0030][0030]FIG. 6A is a block diagram of one embodiment of a patient health monitor.
 [0031][0031]FIG. 6B illustrates one embodiment of a residual test data process in the patient health monitor.
 [0032][0032]FIG. 6C illustrates one embodiment of a statistical test process in the patient health monitor.
 [0033][0033]FIG. 6D illustrates one embodiment of an input/output interface in the patient health monitor.
 [0034][0034]FIGS. 7A and 7B illustrate time history plots of residual test data with respect to measurement time.
 [0035][0035]FIG. 8 is a functional diagram of one embodiment of an artificial pancreas.
 [0036][0036]FIG. 9A illustrates one embodiment of an optimal controller for a closed loop system.
 [0037][0037]FIG. 9B illustrates one embodiment of a control model for a controller in accordance with the present invention.
 [0038][0038]FIG. 9C illustrates one embodiment of a control algorithm.
 [0039]The present invention involves application of realtime optimal estimation, optimized realtime decisions and displays (e.g., a patient health monitor), or optimal realtime control to physiological processes. In one embodiment, the realtime optimal estimation, the optimized realtime decisions and displays, and the optimal realtime control are implemented as separate modules which can be combined functionally. In an alternate embodiment, the realtime optimal estimation, the optimized realtime decisions and displays, and the optimal realtime control are implemented as one joint algorithm.
 [0040]In one embodiment, input to an optimal estimator is provided by a physical sensor (or a plurality of sensors) which measures some arbitrary, but known, function (or functions) of variables or parameters to be estimated to within specified uncertainties and whose statistical characteristics are known. In one embodiment, an output of a realtime controller is provided to a physical controllable dispenser, or actuator, whose output is some known function of parameters to be controlled within specified uncertainties and whose statistical characteristics are known. In one embodiment, a decision and display function utilizes statistical testing of estimator residual errors using internally computed, and updated, estimator variances and covariances.
 [0041]In one embodiment, the present invention is implemented as a software algorithm. The present invention uses models (e.g., dynamic process and measurement models). For best performance, the models should reflect the latest and most complete information available. As new and more complete information is developed, performance can be improved through incorporation of this information by simply modifying the software algorithm of the present invention.
 [0042]Embodiments of the present invention will be described hereinafter with reference to the drawings. FIG. 1A is a block diagram of one embodiment of an estimator. The estimator uses a linearized Kalman filter. The linearized Kalman filter accommodates nonlinear process models and/or nonlinear measurement models. In one embodiment, the linearized Kalman filter is a discrete extended Kalman filter which is linearized after each update using best estimates.
 [0043]A general formulation of a continuousdiscrete extended Kalman filter is provided in Table 1. In an alternate embodiment, the linearized Kalman filter is linearized about a nominal set for which a general formulation is given in Table 2.
 [0044]The estimator computes an estimator gain based on time updated and measurement updated error variable variances and covariances. In one embodiment, the estimator is implemented using discrete formulations. In an alternate embodiment, the estimator is implemented using continuous formulations. In the actual development of algorithms, one can choose a Covariance formulation or an Information formulation depending on initialization uncertainty considerations. Further, the use of Bierman factorization techniques (UDU{circumflex over ( )}T) in the implementation leads to numerically stable algorithms which are excellent for operation over very long time periods.
 [0045]In one embodiment, the estimator of FIG. 1A is applied to the problem of monitoring patient glucose levels. Any type of physical sensor which measures some function of glucose can be used. In one embodiment, one or more capillary blood glucose (CBG) measurements are obtained on initialization of the estimator and when estimated glucose deviations exceed computed variance levels.
 [0046]In another embodiment, additional capillary glucose values are obtained on a periodic basis to assess sensor function. As an example, CBG measurements may be obtained once or twice a day pre and one hour post prandial or when the estimator determines that the glucose values are out of range of predetermined limits.
 [0047]In one embodiment, the estimator is contained in a small portable package. The estimator can be operated by a patient. Alternatively, the estimator can be operated by medical staff in a hospital or clinic. The estimator (e.g., Kalman estimator) uses a Kalman filter to provided realtime estimates of patient glucose levels based upon glucose sensor measurements. The Kalman estimator can also use additional information to more closely predict time propagation of glucose levels, such as exercise, food intake, insulin administration, or other factors which influence glucose levels or the function of the sensor (e.g., local pH, temperature or oxygen tension near the physical sensor).
 [0048][0048]FIG. 1B is a block diagram of one embodiment of a glucosemonitoring device 100. The glucosemonitoring device 100 includes a glucose sensor 108, an ancillary sensor 110, a glucose estimator 114, and a patient health monitor 114. The glucosemonitoring device 100 is used to provide realtime glucose estimates of a patient 102. In one embodiment, the glucosemonitoring device 100 is an integrated unit which is portable by the patient 102 to provide continuous glucose monitoring and realtime displays.
 [0049]The glucose sensor 108 (e.g., a probe, patch, infrared, or laser sensor) is coupled to the patient 102 and outputs a measurement ƒ(g) which is a function of the glucose level of the patient 102. The glucose sensor 108 makes measurements within specified uncertainties and has known statistical characteristics. The glucose sensor 108 provides the measurement ƒ(g) to the glucose estimator 112. In one embodiment, the glucose sensor 108 makes measurements periodically. In an alternate embodiment, the glucose sensor 108 makes measurements intermittently or upon command.
 [0050]The ancillary sensor 110 is coupled to the patient 102, the glucose sensor 108, and/or surroundings of the patient 102 and/or glucose sensor 108 to provide information regarding environmental and/or glucose sensor conditions (e.g., temperature, humidity, local pH, etc.) which affect the measurement. One or more outputs of the ancillary sensor 110 are provided to the glucose estimator 112. In one embodiment, the ancillary sensor 110 provides outputs to the patient health monitor 114 which can process the information for display or for forwarding to the glucose estimator 112. In an alternate embodiment, the ancillary sensor 110 is not a part of the glucosemonitoring device 100.
 [0051]The patient 102 and/or a medical staff 104 (i.e., a user) can provide information on the environmental and glucose sensor conditions as well as other information affecting the measurement. In one embodiment, the patient 102 and/or medical staff 104 inputs information (e.g., exercise activity, food intake, insulin administration, etc.) using the patient health monitor 114. The patient health monitor 114 acts as an input/output interface or a means for the user to configure the glucose estimator 112. The patient health monitor 114 forwards the information to the glucose estimator 112.
 [0052]In one embodiment, the patient health monitor 114 is a display device. For example, an output of the glucose estimator 112 (e.g., an optimal realtime estimate of glucose) is provided to the patient health monitor 114 which displays the information in a comprehensible format for the patient 102 and/or the medical staff 104. The patient health monitor 114 can also contemporaneously display information provided by the glucose sensor 108, the ancillary sensor 110, the patient 102 and/or the medical staff 104 which affects the measurement used to make the optimal glucose estimate.
 [0053]In another embodiment, the patient health monitor 114 is a status indicator. For example, the glucose estimator 112 provides outputs (e.g., residuals and variances) to the patient health monitor 114 which applies statistical testing to determine the reliability of the glucose sensor 108 and/or the glucose estimator 112. The patient health monitor 114 provides a warning when poor performance is detected.
 [0054]In one embodiment, the glucose estimator 112 is a linearized Kalman filter (e.g., a discrete extended Kalman filter) to account for a nonlinear process model and/or measurement model. The glucose estimator (or Kalman estimator) 112 provides realtime estimates of the glucose level in the patient 102. An initialization measurement 106 (e.g., a CBG measurement) is obtained from the patient 102 and provided to the Kalman estimator 112 to initialize the Kalman estimator 112.
 [0055][0055]FIG. 2 is a flow chart illustrating one embodiment of a modeling process for physiological processes. In one embodiment, physiological processes are described by nonlinear stochastic models. In one embodiment, the flow chart of FIG. 2 illustrates one method of developing a dynamic model for an optimal estimator. The method includes steps for state vector development, nominal dynamic/measurement model development, linearized model development, and uncertainty modeling and nominal model verification.
 [0056][0056]FIG. 2 illustrates modeling of physiological processes in a linearized Kalman filter application. As an example in glucose estimation, in a first block of FIG. 2, glucose is an estimation variable. The time rate of change of glucose might become another estimation variable. An uncertain parameter might be the glucose sensor scale factor, and additionally, the rate of change of scale factor over time could be another.
 [0057]In one embodiment of a second block of FIG. 2, the way in which variables and parameters nominally propagate over time may change with conditions. In the example of glucose estimation, inputs can be used to identify patient related activities: eating; exercising; sleeping; insulin injection; etc. With these identified patient related activities, additional state variables can be identified and modeled. For example, if patient eating can be related to a change in glucose level over some specified time interval, dynamics can be implemented within the estimator model which will propagate (or extrapolate) a rise in the glucose level over that time interval. This rise may be modeled by appropriate functions whose variables contain uncertainties, which may be added as elements of the state vector. In an analogous way, decreases in glucose levels (e.g., due to insulin injections) can be modeled. These models may be general in nature, or they may be patient specific. Consequently, patient related activities that have a significant impact on glucose levels, or the rate of change of glucose levels, can be accounted for within the dynamic process model.
 [0058]Insofar as sensor modeling is concerned, experience shows that, for example, infrared sensor measurement bias errors vary with, among other things, temperature. If this variation with temperature can be modeled, and included in the process model, then a temperature measurement will improve estimator performance. A particular physical glucose sensor may have a scale factor which has a characteristic decay in sensitivity over time as discussed in further detail below. There may be other variables that can be measured which will affect the physical function of a sensor such as local pH, oxygen tension, etc.
 [0059]In one embodiment, the nominal modeling is comprised of three types: 1) predictable characteristics of a physical sensor function over time (e.g., a fixed rate in decline of sensor output), 2) other measurable physical variables which may affect sensor function (e.g., local temperature, pH, etc.), and 3) predictable changes in the model which occur with patient related activities (e.g., exercise, eating or insulin administration). Changes to the dynamic process model may add variables and/or uncertain parameters to the state vector and changes to measurement models. As a result, the activities indicated in blocks 1 and 2 of FIG. 2 constitute an iterative process.
 [0060]In one embodiment of block 3 of FIG. 2, the relationship between variables, parameters, and measurements determines which parts of the processes and measurements are nonlinear and are linearized. Further, the relationship between variables, parameters, and measurements determines if the variables and parameters are observable and can be estimated. In certain cases, observability can be enhanced through introduction of additional modeling information. For the glucose example, the ability to estimate sensor scale factor and/or detect sensor failure may be improved by modeling glucose propagation changes due to insulin injections or ingestion of sugar. For example, tracking known changes enhances scale factor observability through estimator generated correlations. Further, if it is known that glucose levels vary and the sensor measurement does not change accordingly within prescribed levels of uncertainty, a sensor problem is indicated.
 [0061]In one embodiment of block 4 of FIG. 2, the development of a relatively large database is used to empirically verify and/or modify the nominal nonlinear dynamic process/measurement models and derive uncertainty levels associated with the variables, parameters, and measurements. In one embodiment, the empirical data can be fitted to nonlinear functions using a nonlinear regression package contained in a commercially available software application program such as “Mathematica” from Wolfram Research, Inc. Analytical functions may be added or modified using the test database. The repeatability of the fit over nominal ranges of the patient environment determines the uncertain parameters and the variations establish uncertainty levels. The more accurate the dynamic process and measurement models, the more the uncertainties are reduced, and the greater the estimator performance. As the process evolves, the modeling becomes better defined through iterations between blocks 4 and 2. Certain portions of the models may be developed on an individual basis (e.g., for a specific patient).
 [0062]A database is used to empirically develop and verify models. An embodiment discussed below uses the database to develop two separate dynamic process models. For example, after processing a number of data sets from a physical sensor that was used to monitor the glucose level of various patients, it was observed that the sensor scale factor was equally likely to move up or down over the first fifteen to twentyfive hours. However, the scale factor tended to decay for the remaining life of the sensor after this period of time, usually three or four days. These observations are incorporated in an embodiment of the estimator discussed in further detail below.
 [0063][0063]FIG. 3A describes an estimation function which depicts a linearized Kalman filter formulation. In one embodiment, the particular form is that of a discrete extended Kalman filter which linearizes after each update using best estimates. In the formulation, a vector whose elements comprise variables and/or parameters (with uncertainties) for making estimates defines the state of a system.
 [0064]There is a distinction between variables (e.g., random variables) and parameters with uncertainty. Variables, such as glucose and rate of change of glucose, are estimated and can be controlled (if control is implemented). Parameters with uncertainty are part of the model structure not known precisely and are estimated and updated (like variables) but not controlled, e.g., a glucose sensor scale factor or insulin dispensing controller scale factor.
 [0065]In one embodiment, realtime variable and parameter estimates are used to relinearize the model following each update. Inputs to the estimator can consist of any measurement which is related to, or can be correlated with, any element in the state vector. In the case of glucose estimation, other inputs can consist of dynamic process configuration control based on patient related activities or other circumstances. Following initialization, the time update and measurement update cycles form a recursive loop. The time update period is the time interval between the receipt of measurements. This is a function of the sensor and of acceptable latency in the estimates.
 [0066][0066]FIGS. 3B, 3C and 3D illustrate, respectively, Initialization and the recursive Time Update and Measurement Update cycles according to one embodiment of the present invention. Table 3 defines, in more detail, symbols used in these figures. In one embodiment, dynamic process and measurement modeling is contained in software algorithms with parameter and structure updates in real time.
 [0067]One embodiment of an initialization of the estimation process is described in FIG. 3B. A state vector estimate, Xeo, contains the initial estimates of the process variables and model parameters with uncertainty, while the covariance matrix, Po, contains the initial variances and covariances associated with the Xeo elements.
 [0068][0068]FIG. 3B shows an example embodiment with a twoelement state vector which is based on empirical observations described above. As an example, the point in time when the scale factor begins to decay was chosen as 20 hours, a nominal value over the database. An exponential decay was chosen to model this decay rate and is consistent with the first derivative of scale factor equal to a parameter, alpha, multiplied by the variable scale factor. In this case, alpha is not modeled as an uncertain parameter; a value of 0.012 was chosen as a nominal value, over the database, for the fiveminute cycle time of this physical sensor.
 [0069]The initial estimate of glucose is 150 mg/dl (state vector element) with an initial uncertainty variance of 100 mg/dl squared (a covariance matrix element). When processing the data, the glucose element was initialized by setting it substantially equal to the first capillary blood glucose measurement. The nominal initial scale factor value, for the physical sensor, is 0.25 nanoamps/(mg/dl) with a parameter uncertainty (variance) of 0.1 nanoamps/(mg/dl) squared.
 [0070]The initial covariance between the glucose variable and the scale factor parameter is zero. Correlation between glucose and scale factor will develop as the estimator processes the sensor measurements. Each type of sensor will have its own model and characteristics. The measurement uncertainty is 5 nanoamps squared, an element of the R matrix. In one embodiment, this is a scalar measurement and the R matrix contains a single element.
 [0071]A second measurement in this example embodiment is an occasional direct measurement, such as a capillary blood glucose measurement (CBG), with a unity scale factor, and a measurement error uncertainty of 15 mg/dl squared (a single element in a second R matrix). This is probably better modeled as 15% of the measured glucose value. The growth in uncertainty of glucose and scale factor from measurement to measurement is, respectively, 20 mg/dl squared and 0.002 nanoamps squared (elements of the process noise matrix, Q).
 [0072][0072]FIG. 3C illustrates a timeupdate process of the recursive process in one embodiment of an optimal estimator. In this figure, a negative superscript indicates a time update while a positive superscript indicates a measurement update. In the brackets, a letter “i” indicates time at the ith interval and (i−1) indicates time at the previous time interval. In an update process, a state vector is first updated since these elements are used to update matrices and to bring the time epoch of the estimated measurement to be consistent with that of the next measurement to be received.
 [0073]In an example embodiment with a twoelement state vector, the dynamic process is linear. With no patient inputs, the first derivative of the measured physiologic variable is zero, corresponding to the case when the level (on average) is as likely to either go up or to go down. No additional a priori information is assumed about the time propagation of glucose.
 [0074]In one embodiment of estimating glucose, the solution to the scale factor propagation after the first 20 hours is defined by the exponential shown in the second column. Consequently, for this embodiment, the dynamic process function (f), is linear, is not a function of state vector elements, and, from linear system theory, the transition matrix (A) is the 2 by 2 identity matrix for the first 20 hours and thereafter is defined by the 2 by 2 matrix in the second column.
 [0075]In one embodiment, the measurement function (h) for the sensor measurement is nonlinear in the state vector elements. If Ge and Se are used to denote glucose and scale factor estimates, respectively, then Ye=Se*Ge. Definitions of the above terms are provided in Table 3. When linearized using best estimates, the linearized measurement matrix H=[Se Ge] and is of the same functional form both before and after 20 hours.
 [0076][0076]FIG. 3D illustrates a measurementupdate process of one embodiment of an optimal estimator. The measurement update sequence begins with the computation of the gain matrix, K(i). The difference between the actual sensor measurement and the best estimate of the measurement is computed: y(i)=Ym(i)−Se(i)*Ge(i). This difference, or residual, when multiplied by the gain matrix and added to the timeupdated estimate produces the measurementupdated estimate.
 [0077]The covariance matrix is then measurement updated, reflecting the level of uncertainty in the estimates following the processing of a measurement. In the glucose example, if a second measurement is available, such as a CBG, then the measurement sequence is again cycled through, starting with a new gain computation, and using the appropriate measurement matrix and new best estimate of the next measurement. Following the processing of all available measurements at the ith time epoch, the updated state vector and covariance matrix are then used to start the time update for the (i+1)th time epoch which begins the next cycle.
 [0078]One embodiment of an estimation algorithm illustrating the initialization, the time update process, and the measurement update process discussed above is provided in Table 4. Table 4 is an algorithm programmed in the MATLAB language (from Math Works). This printout defines a working program and has been used to process a significant number of data sets. In the glucose example, the estimation results from the process of two data files, both gathered from the same patient, and taken about a month apart, are shown in FIGS. 4A, 4B, 5A, and 5B.
 [0079]In one embodiment, sensor inputs are provided and processed every 5 minutes. Occasional CBGs are also provided. For example, two CBGs per day were processed by an estimator; and additional CBG values were used to judge estimator performance by comparing glucose estimates with actual CBG values not used by the estimator.
 [0080][0080]FIGS. 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time. FIG. 4A shows a time history of the realtime estimates of glucose (e.g., every 5 minutes) along with all available discrete CBGs.
 [0081][0081]FIG. 4B shows a time history of the glucose estimates along with CBGs that were processed. In this figure, the estimated glucose value took several rapid swings between approximately 48 hours and 58 hours. FIG. 4B indicates that only sensor measurements were processed during that interval and no CBGs were processed. However, the CBGs plotted in FIG. 4A indicate that the glucose estimates did tend to follow the excursions of the patient glucose levels.
 [0082][0082]FIGS. 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time. The time duration for these runs is about 4 days without patient inputs. Over this time period, the patient ate, slept, exercised, and took insulin injections. Dynamic models to account for these activities, much like the decaying scale factor, could be implemented and called into use upon command.
 [0083][0083]FIG. 6A is a block diagram of one embodiment of a patient health monitor. On example of an algorithm for the patient health monitor is in Table 4 which generates statistical test data based on Kalman filter residuals as well as test displays. In one embodiment, the patient health monitor generates realtime decisions and displays which are integrated with a Kalman filter. The patient health monitor allows the patient or medical staff to interact with a Kalman estimator and/or Kalman controller described herein.
 [0084]In one embodiment, the patient health monitor provides insight into how well a Kalman filter is working through the filter residual that is the difference between the estimate of the measurement at the time the measurement is received and the actual measurement. In another embodiment, other checks are used from time to time, such as the CBGs in the glucose example described above. An example of another check on the status of a sensor is through the use of patient inputs or signals indicating that something is changing in a prescribed way and then noting whether or not the sensor is observing this change within prescribed uncertainties.
 [0085]In one embodiment, the realtime displays and decisions of the patient health monitor uses some occasional outside checks but relies substantially on results of statistical testing performed on filter residuals. If estimates of the measurements are, on average, good (e.g., residuals are small and unbiased), then the estimator is generally working well, and vice versa. More specifically, elements of the covariance matrix can be used to construct statistical test applications.
 [0086]For hospital applications, realtime displays of glucose estimates along with realtime displays of estimator performance test results can be important visual inputs to the medical staff. Requests for additional CBG measurements or the sounding of an alarm in the event glucose estimates exceed critical limits may also prove useful. In one embodiment, a reduced number of outputs is provided in relatively small estimators for individual use.
 [0087][0087]FIG. 6B illustrates one embodiment of a residual test data process in the patient health monitor. For example, the residual (y), covariance matrix (P), measurement matrix (H), and measurement noise matrix (R) are available from an estimator algorithm at each time epoch (i). If the estimator is operating properly, the sequence of residuals has the property of zero mean, white noise sequence, i.e., any two residuals taken at different times are uncorrelated (E[y(i)*y(j)]=0, for all j not equal to i). This condition provides a unique means for constructing statistical tests.
 [0088]Visually, a time history plot of the residuals, Sy(i) in FIG. 6B, should appear random, zero mean, and unbaised. If they are summed over time, the deviation of the sum from zero should, on average, grow as the square root of time, as should its absolute value, ASy(i). The fact that E[y(i)*y(j)]=0, for i not equal to j, also means that the sum of the variances, SV(i) in FIG. 6B, which is easily computed, is equal to the variance of the sum of residuals for a properly performing estimator. The standard deviation of the sum of residuals, StdSV(i) also grows as the square root of time.
 [0089][0089]FIG. 6C illustrates one embodiment of a statistical test process in the patient health monitor. FIG. 6C defines statistical tests which can be constructed based on filter residual test data. One of these compares the absolute value of the sum of the residuals with the standard deviation of the sum of residuals. On average, the ASy(i) should be bounded by the StdSV(i). If not, this indicates that the deviations of Sy(i) are growing faster than that of a white noise sequence, implying degraded estimator performance.
 [0090]Realtime displays of ASy(i) and StdSV(i) histories can provide a visual picture of estimator performance. An example is provided in FIGS. 7A and 7B which illustrate time history plots of residual test data with respect to measurement time. These two plots are one form of residual tests for the two sets of glucose estimation results provided in FIGS. 4A, 4B, 5A, and 5B respectively. In both cases, the sum of residuals are well behaved and was bounded by the Standard Deviation, Std, of the sum of residuals.
 [0091]The example embodiment for the glucose application of test data, test generation, and test result display described above is implemented by the algorithm in Table 4. Other realtime quantitative tests can also be constructed using these data. Tests on individual residuals can be performed using individual variances. If a measurement is received which causes the residual to exceed a four signal level, for example, the action might be to emit a warning and request an immediate CBG measurement. Other tests are identified in FIG. 6C and will be evaluated as the system develops.
 [0092][0092]FIG. 6D illustrates one embodiment of an input/output interface in the patient health monitor. The patient health monitor provides a capability for the patient and/or staff to communicate and interact in real time with estimation and control processes using simple commands, visual displays, and audio outputs described above. Confidence levels in estimator performance can be established and a genuine interface established whereby the estimation and control processes could request additional information to check and insure confidence in estimates, physical sensors, and physical controllers. The staff can provide useful realtime inputs to augment this process.
 [0093][0093]FIG. 8 is a block diagram of one embodiment of a Kalman optimal stochastic control solution as applied to physiological processes. In one embodiment, the application includes both optimal stochastic regulator and optimal stochastic tracking control solutions. Tracking control involves a controlled variable following a reference value, constant or dynamic, as closely as possible. Control is applied to physiological processes, wherein a control gain is computed based upon optimization criteria which minimizes controlled variable errors while minimizing application of control based on cost weightings.
 [0094]There is a duality between the computation of a Kalman estimator gain and a Kalman control gain. The Kalman estimator gain minimizes estimation error variances. The Kalman control gain minimizes variances of error between the controlled variables and the reference variables while minimizing the level of control applied. As a result, the optimal control function includes specification of the controlled variables and their associated costs as well as costs associated with the amount of control to be applied. For example, smaller control variable error costs and larger control application costs will allow the controlled variable to deviate farther from the reference, but with reduced application of control.
 [0095]In one embodiment, linearization techniques described above in association with the Kalman estimator is applied to physiological nonlinear stochastic processes. For example, linearization about nominal values or about best estimates are provided by the Kalman estimator. Uncertain parameters associated with a controllable dispenser, or actuator, are included in an estimator state vector.
 [0096]The optimal stochastic controller can use linear or nonlinear formulations and discrete or continuous time formulations. In one embodiment, the optimal stochastic controller is used with an optimal estimator described herein and/or an optimized decision and display function also described herein to form a closed loop system. The closed loop system works as an artificial pancreas when applied to a glucose problem in one embodiment.
 [0097][0097]FIG. 8 is a functional diagram showing a controllable dispenser, or pump, with the capability to secrete insulin and glucagon to control high and low glucose levels, respectively. In one embodiment of a closed loop system, an estimator and a controller share a state vector wherein the estimator estimates it and the controller controls designated elements of it. Unlike the estimator, the time varying gain computations for the controller is computationally intensive and may not be used in all applications. The controller includes time varying or steady state gain formulations.
 [0098][0098]FIG. 9A illustrates one embodiment of an optimal controller for the artificial pancreas. The glucose control problem is a tracking problem since the glucose level is controlled to a desired level which may either be constant or a function of time. In one embodiment, the dynamic process model includes a quantitative description of how glucose levels propagate, in time, as a function of insulin/glucagon secreted by the controllable pump. The dynamic process model is described by a system of first order differential or difference equations. Excluding the pump, much of the modeling for the glucose control is available through the estimator modeling development.
 [0099][0099]FIG. 9B illustrates one embodiment of a control model for a controller. In an example embodiment, the dynamic process model is described by three first order differential equations that are forced by a control variable, u. The dynamics of this process include a first order time lag and a scale factor associated with the pump, and a first order time lag and a scale factor associated with the physiological process. The model for the glucose sensors is the same as in the example glucose estimator embodiment. In the definition of the elements of the transition matrix, alpha g, s, and d are the inverse of the first order time lags associated with glucose, glucose sensor scale factor, and pump, respectively. Beta g is the glucose scale factor multiplied by alpha g. Delta t is the time interval between measurement/control application epochs. The controlled error is the difference between the estimate of glucose (provided by the estimator) and the glucose control point input.
 [0100]In a controller command, Beta d is a pump scale factor multiplied by alpha d. The measurements are the same as for the example glucose estimator, except that the measurements contain an additional zero since the measurements are not functions of the variable insulin/glucagon.
 [0101]The cost function, which is minimized by the optimal stochastic control, contains costs associated with the glucose error and the application of control. Choosing a value of Cu that is much larger than Cg will result in a relative gentle application of control. Another example embodiment would utilize a higher order model in which the first derivative of glucose would be included in the state vector, and included as a control variable. If a relatively large cost is associated with the first derivative relative to the glucose control point error, then the control will be very active when rapid changes in patient glucose occur while relatively gentle otherwise. Other example embodiments would include pump scale factor as a state vector element (uncertain parameter) as well as insulin/glucagon measurements to the estimator.
 [0102][0102]FIG. 9C illustrates one embodiment of a control algorithm based on the model defined in FIG. 9B. Using either a time varying gain or a precomputed steady state gain, control is applied at each epoch based upon a difference between an estimate of a patient glucose level and a glucose control point. A control, u, is applied to a pump over each time interval, and a state vector is time updated as indicated in this Figure using the control variable, u, as a forcing function. Other time updates are performed in accordance with embodiments discussed above in association with estimator equations.
 [0103]A potential problem in the application of closed loop control to that of physiological processes is due to potentially long time delays that may be destabilizing. These delays can be in the form of transport lags. A transport lag is the time between when control is applied and when the process action begins. An optimal technique for control using delay states is discussed in “Optimal Control of Linear Stochastic Systems with Process and Observation Time Delays” by E. J. Knobbe (Academic Press, Inc., 1989) and is hereby incorporated herein in its entirety by reference thereto. A discussion of principles developed for optimal control with process and observation time delays is provided in Table 5.
 [0104]Although described above in connection with particular embodiments of the present invention, it should be understood that the descriptions of the embodiments are illustrative of the invention and are not intended to be limiting. Various modifications and applications may occur to those skilled in the art without departing from the true spirit and scope of the invention.
TABLE 1 SUMMARY OF CONTINUOUSDISCRETE EXTENDED KALMAN FILTER System Model $\stackrel{.}{\underset{\_}{x}}\ue8a0\left(t\right)=\underset{\_}{f}\ue8a0\left(\underset{\_}{x}\ue8a0\left(t\right),t\right)+\underset{\_}{w}\ue8a0\left(t\right);\underset{\_}{w}\ue8a0\left(t\right)~N\ue8a0\left(\underset{\_}{0},Q\ue8a0\left(t\right)\right)$ Measurement Model ${\underset{\_}{z}}_{k}={\underset{\_}{h}}_{k}\ue8a0\left(\underset{\_}{x}\ue8a0\left({t}_{\underset{\_}{k}}\right)\right)+{\underset{\_}{v}}_{k};\text{\hspace{1em}}\ue89ek=1,2,\text{\hspace{1em}}\ue89e\dots \ue89e\text{\hspace{1em}};\text{\hspace{1em}}\ue89e{\underset{\_}{v}}_{k}~N\ue8a0\left(\underset{\_}{0},{R}_{k}\right)$ Initial Conditions $\underset{\_}{x}\ue8a0\left(0\right)~N\ue8a0\left({\hat{\underset{\_}{x}}}_{o},{P}_{o}\right)$ Other Assumptions $E\ue8a0\left[\underset{\_}{w}\ue8a0\left(t\right)\ue89e{v}_{k}^{T}\right]=0\ue89e\text{\hspace{1em}}\ue89e\mathrm{for}\ue89e\text{\hspace{1em}}\ue89e\mathrm{all}\ue89e\text{\hspace{1em}}\ue89ek\ue89e\text{\hspace{1em}}\ue89e\mathrm{and}\ue89e\text{\hspace{1em}}\ue89e\mathrm{all}\ue89e\text{\hspace{1em}}\ue89et$ State Estimate Propagation $\stackrel{.}{\hat{\underset{\_}{x}}}\ue8a0\left(t\right)=\underset{\_}{f}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)$ Error Covariance Propagation $\stackrel{.}{P}\ue8a0\left(t\right)=\hat{F}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)\ue89eP\ue8a0\left(t\right)+P\ue8a0\left(t\right)\ue89e{F}^{T}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)+Q\ue8a0\left(t\right)$ State Estimate Update ${\hat{\underset{\_}{x}}}_{k}\ue8a0(+)={\hat{\underset{\_}{x}}}_{k}\ue8a0()+{K}_{k}\ue8a0\left[{\underset{\_}{z}}_{k}{\underset{\_}{h}}_{k}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\right]$ Error Covariance Update ${P}_{k}\ue8a0(+)=\left[I{K}_{k}\ue89e{H}_{k}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\right]\ue89e{P}_{k}\ue8a0()$ Gain Matrix ${K}_{k}=\text{}\ue89e{P}_{k}\ue8a0()\ue89e{{H}_{k}^{T}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\ue8a0\left[{H}_{k}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\ue89e\text{\hspace{1em}}\ue89e{P}_{k}\ue8a0()\ue89e{H}_{k}^{T}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)+{R}_{k}\right]}^{1}$ Definitions $\begin{array}{c}{F\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)=\frac{\partial \underset{\_}{f}\ue8a0\left(\underset{\_}{x}\ue8a0\left(t\right),t\right)}{\partial \underset{\_}{x}\ue8a0\left(t\right)}\uf604}_{\underset{\_}{x}\ue8a0\left(t\right)=\hat{\underset{\_}{x}}\ue8a0\left(t\right)}\\ {{H}_{k}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)=\frac{\partial {\underset{\_}{h}}_{k}\ue8a0\left(\underset{\_}{x}\ue8a0\left({t}_{k}\right)\right)}{\partial \underset{\_}{x}\ue8a0\left({t}_{k}\right)}\uf604}_{\underset{\_}{x}\ue8a0\left({t}_{k}\right)={\hat{\underset{\_}{x}}}_{k}\ue8a0()}\end{array}\hspace{1em}$  [0105][0105]
TABLE 2 SUMMARY OF CONTINUOUSDISCRETE LINEARIZED KALMAN FILTER System Model $\stackrel{.}{\underset{\_}{x}}\ue8a0\left(t\right)=\underset{\_}{f}\ue8a0\left(\underset{\_}{x}\ue8a0\left(t\right),t\right)+\underset{\_}{w}\ue8a0\left(t\right);\underset{\_}{w}\ue8a0\left(t\right)~N\ue8a0\left(\underset{\_}{0},Q\ue8a0\left(t\right)\right)$ Measurement Model ${\underset{\_}{z}}_{k}={\underset{\_}{h}}_{k}\ue8a0\left(\underset{\_}{x}\ue8a0\left({t}_{\underset{\_}{k}}\right)\right)+{\underset{\_}{v}}_{k};\text{\hspace{1em}}\ue89ek=1,2,\text{\hspace{1em}}\ue89e\dots \ue89e\text{\hspace{1em}};\text{\hspace{1em}}\ue89e{\underset{\_}{v}}_{k}~N\ue8a0\left(\underset{\_}{0},{R}_{k}\right)$ Initial Conditions $\underset{\_}{x}\ue8a0\left(0\right)~N\ue8a0\left({\hat{\underset{\_}{x}}}_{o},{P}_{o}\right)$ Other Assumptions $\begin{array}{c}E\ue8a0\left[\underset{\_}{w}\ue8a0\left(t\right)\ue89e{v}_{k}^{T}\right]=0\ue89e\text{\hspace{1em}}\ue89e\mathrm{for}\ue89e\text{\hspace{1em}}\ue89e\mathrm{all}\ue89e\text{\hspace{1em}}\ue89ek\ue89e\text{\hspace{1em}}\ue89e\mathrm{and}\ue89e\text{\hspace{1em}}\ue89e\mathrm{all}\ue89e\text{\hspace{1em}}\ue89et\\ \mathrm{Nominal}\ue89e\text{\hspace{1em}}\ue89e\mathrm{trajectory}\ue89e\text{\hspace{1em}}\ue89e\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right)\ue89e\text{\hspace{1em}}\ue89e\mathrm{is}\ue89e\text{\hspace{1em}}\ue89e\mathrm{available}\end{array}\hspace{1em}$ State Estimate Propagation $\stackrel{.}{\hat{\underset{\_}{x}}}\ue8a0\left(t\right)=\underset{\_}{f}\ue8a0\left(\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right),t\right)+P\ue8a0\left(\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right),t\right)\ue8a0\left[\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right)\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right)\right]$ Error Covariance Propagation $\stackrel{.}{P}\ue8a0\left(t\right)=\hat{F}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)\ue89eP\ue8a0\left(t\right)+P\ue8a0\left(t\right)\ue89e{F}^{T}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)+Q\ue8a0\left(t\right)$ State Estimate Update ${\hat{\underset{\_}{x}}}_{k}\ue8a0(+)={\hat{\underset{\_}{x}}}_{k}\ue8a0()+{K}_{k}\ue8a0\left[{\underset{\_}{z}}_{k}{\underset{\_}{b}}_{k}\ue8a0\left(\stackrel{\_}{\underset{\_}{x}}\ue8a0\left({t}_{k}\right)\right){H}_{k}\ue8a0\left(\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right)\right)\ue8a0\left[{\hat{\underset{\_}{x}}}_{k}\ue8a0()\stackrel{\_}{\underset{\_}{x}}\ue8a0\left({t}_{k}\right)\right]\right]$ Error Covariance Update ${P}_{k}\ue8a0(+)=\left[I{K}_{k}\ue89e{H}_{k}\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left({t}_{k}\right)\right)\right]\ue89e{P}_{k}\ue8a0()$ Gain Matrix ${K}_{k}=\text{}\ue89e{P}_{k}\ue8a0()\ue89e{{H}_{k}^{T}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\ue8a0\left[{H}_{k}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)\ue89e\text{\hspace{1em}}\ue89e{P}_{k}\ue8a0()\ue89e{H}_{k}^{T}\ue8a0\left({\hat{\underset{\_}{x}}}_{k}\ue8a0()\right)+{R}_{k}\right]}^{1}$ Definitions $\begin{array}{c}{F\ue8a0\left(\hat{\underset{\_}{x}}\ue8a0\left(t\right),t\right)=\frac{\partial \underset{\_}{f}\ue8a0\left(\underset{\_}{x}\ue8a0\left(t\right),t\right)}{\partial \underset{\_}{x}\ue8a0\left(t\right)}\uf604}_{\underset{\_}{x}\ue8a0\left(t\right)=\stackrel{\_}{\underset{\_}{x}}\ue8a0\left(t\right)}\\ {{H}_{k}\ue8a0\left({\stackrel{\_}{\underset{\_}{x}}}_{k}\ue8a0\left({t}_{k}\right)\right)=\frac{\partial {\underset{\_}{h}}_{k}\ue8a0\left(\underset{\_}{x}\ue8a0\left({t}_{k}\right)\right)}{\partial \underset{\_}{x}\ue8a0\left({t}_{k}\right)}\uf604}_{\underset{\_}{x}\ue8a0\left({t}_{k}\right)={\stackrel{\_}{\underset{\_}{x}}}_{k}\ue8a0\left({t}_{k}\right)}\end{array}\hspace{1em}$  [0106][0106]
TABLE 3 State Vector X, Xe, Xe(0): These symbols denote, respectively, the state vector, the best estimate (or optimal estimate) of the state vector, and the initial best estimate of the state vector. They contain both variables and parameters as defined above and, with regard to estimation, there are no distinctions. Dynamic Process Model dX = f(X,t), dXe = f(Xe,t): These symbols denote the time derivative of the state vector and estimated state vector, respectively. They represent a system of first order differential equations (or difference equations) which describe the manner in which the state elements propagate in time. The function, f, may be linear or nonlinear. Measurement Model Ym = h(X) + e, Ye = h(Xe): These symbols denote the actual measurement and the best estimate of the measurement, respectively. The function, h, defines the arbitrary, but known, way in which the state vector elements are related to the measurement, and e represents the sensor measurement error. Covariance Matrix P, P(0): Denotes, respectively, the matrix of variances and covariances associated with the errors of each of the state variable estimates (within the estimated state vector) and their initial values, Transition Matrix A, A(0): Denotes the transition matrix, which is used to propagate the covariance matrix forward it time (time update), and it's initial value, respectively; in the linear case it is also used to time update the state vector estimate. Process Noise Matrix Q: Denotes the matrix of variances and covariances associated with error growth uncertainty accumulated in the state variable estimates since the last measurement update, Measurement Matrix H, H(0): Denotes the measurement matrix and it's initial value, respectively. This matrix defines the linear functional relationship between the measurement and the state vector elements. If, the measurement model is linear, then H = h; if h is nonlinear, then H is defined by linearizing h using partial derivatives or perturbation techniques. Measurement Noise Matrix R: Denotes the matrix of variances and covariances associated with measurement error uncertainties. Kalman Gain Matrix K: Denotes the Kalman gain matrix which, when multiplied by the difference between the actual measurement, and the best estimate of the measurement, yields the estimated state correction. This estimated state correction, when added to the old best estimate, becomes the new best estimate. Estimated Measurement Error (Residual) y = Ym − Ye: Denotes the difference between the actual measurement and the best estimate of the measurement. This difference multiplied by the Kalman gain yields the correction to the previous state vector estimate.  [0107][0107]
TABLE 4 % Optimal Glucose Estimator. Extended Kalman Filter; Computes Statistical Test Data for Patient Health Monitor %Ed_JYP_Estimator; Twostep scalefactor model, Exp = .999 for t>20; r=15{circumflex over ( )}2 for y<10 ; x=[150;.25]; hh=[1 0];p=[100{circumflex over ( )}2 0 ;0 .1{circumflex over ( )}2];a=[1 0 ; 0 1]; q=[20{circumflex over ( )}2 0 ; 0 .002{circumflex over ( )}2]; rr=15{circumflex over ( )}2; r= 5{circumflex over ( )}2; I=[1 0;0 1]; load EdJY528.prn, [m,n]=size(EdJY528), t=EdJY528 1:, 1)/12; y=EdJY528(:,2); g=EdJY528(:,3); gm=EdJYS28(:,4); x=[gmf(1); .25]; sz=0.0; SVZ = p(1,1)*K(2){circumflex over ( )}2 + p(2,2)*x(1){circumflex over ( )}2; for i = 1:m; . . . r = 5{circumflex over ( )}2; . . . zz = 0; . . . if y(i) < 10 r = 15{circumflex over ( )}2; . . . end; if t(i) > 20 a =(1 0 ; 0 .999]; . . . end; p=a*p*a′+q; . . . x=a*x;. . . if gm(i) > 0 k=p*hh*/(hh*p*hh*+rr); . . . zz=(gm(i)−x(1)); . . . x=x+k*zz; . . . p=(I−k*hh) *p; . . . end; h=[x(2) x(1)]; . . . k=p*h′/ (h*p*hh+r); . . . x=x+k*z; . . . p=(I−k*h) *p; . . . az=abs (z); . . . sz=sz+z;. . . asz = abs (sz); . . . VZ = p(1,1) *x(2){circumflex over ( )}2 + p(2,2)*x(1){circumflex over ( )}2 + 2*p(1,2)*x(1)*x(2)+ r ; . . . sigZ = sqrt(VZ); . . . SVZ = SVZ +VZ; . . . siqSZ = sqrt(SVZ); . . . VZZ = p(1.1) + rr; . . . sigZZ = sqrt(VZZ); . . . azz = abs(zz); . . . xhistory(i,1) = x(1); . . . xhistory(i,2) = x(2); . . . gmhistory =gm; . . . ghistory =g; . . . zhistory(i) =z; . . . sigsz_history(i) = sigSZ; . . . sigz_history(i) = 2*sigZ; . . . aszhistory(i) = asz; . . . azhistory(i) = az; . . . zzhistory(i) = zz; . . . sigzz_history(i) = 2*sigZZ; . . . azzhistory(i) = azz; . . . sig_history(i,1)=sqrt(p(1,1)); . . . sig_history(i,2)=sqrt(p(2,2)); . . . yhistory=y;. . . end; figure (1) plot(t,xhistory(:,2)); xlabel (′ Measurement Time (hrs) ′); ylabel(′ Scale Factor Estimate; (ISIG Units/(mg/dl))′); title(′Scale Factor Estimate vs Time′); figure (2) plot (t, gmhistory (:), t,xhistory (:,1)); xlabel(′Time (hrs)′); Ylabel(′Glucose Meas & Estimated Glucose; (mg/dl)′); title(′All Glucose Measurements & Estimated Glucose vs Time′); figure (3) plot (t,yhistory (:), t, zhistory(:)); xlabel(′Time (hrs)′); Ylabel(′ISIG Meas & z Residual; (ISIG Units) ′); title(′ISIG Measurement & z residual vs Time′); figure (4) plot(t,xhistory(:,1),t,ghistory(:,1)); xlabel(′ Time (hrs)′); ylabel(′Glucose & Estimated Glucose; (mg/dl)′); title(′Glucose Estimate & All Glucose CBGs vs Time′); figure (5) plot(t,sigzz_history(:),t,azzhistory(:)); xlabel (′Time (hrs) ′); ylabel(′Abs(zz) Residual & two sigma(zz); (mg/dl) ′); title(′Residual Test: 2 Siqma(zz) & Abs(zz) Residual vs Time′); figure (6) plot(t,azhistory(:),t,sigz_history(:)); xlabel(′Measurement Time (hrs)′); ylabel(′z Residual & Standard Deviation; (ISIG Units)′); title(′Residual Test Abs(z) & 2 Sigma(z) vs Time′); figure (7) plot(t,aszhistory(:),t,sigsz_history(:)); xlabel(′ Measurement Time (hrs)′); ylabel(′ sz Residual & Sigma; (ISIG Units) ′); title(′Residual Test : Abs(sz) & Sigma(sz) vs Time′);  [0108][0108]
TABLE 5 PROBLEM DEFINITION The general, discretetime system representation with explicit process and observation time delays is given by $x\ue8a0\left(k+1\right)=\sum _{i=1}^{q}\ue89e{A}_{i}\ue8a0\left(k\right)\ue89eX\ue8a0\left(k+1i\right)+B\ue8a0\left(k\right)\ue89eU\ue8a0\left(k\right)+{W}_{1}\ue8a0\left(k\right),$ (1) $Y\ue8a0\left(k\right)=\sum _{j=1}^{p}\ue89e{C}_{j}\ue8a0\left(k\right)\ue89eX\ue8a0\left(k+1j\right)+{W}_{2}\ue8a0\left(k\right),$ (2) where p and q are integers > 1; W_{1 }and W_{2 }are zeromean whitenoise sequences such that $\begin{array}{c}E\ue89e\left\{{W}_{1}\ue8a0\left(k\right)\ue89e{W}_{1}^{T}\ue8a0\left(k\right)\right\}={V}_{1}\ue8a0\left(k\right)\\ E\ue89e\left\{{W}_{2}\ue8a0\left(k\right)\ue89e{W}_{2}^{T}\ue8a0\left(k\right)\right\}={V}_{2}\ue8a0\left(k\right)\\ E\ue89e\left\{{W}_{1}\ue8a0\left(i\right)\ue89e{W}_{1}^{T}\ue8a0\left(j\right)\right\}=0\end{array}\}\ue89e\text{\hspace{1em}}\ue89e{\forall}_{i,j,k};\hspace{1em}$ Y is an m × 1 observation vector; X is an n × 1 random state vector whose initial uncertainty is uncorrelated with W_{1 }and W_{2 }and with initial covariance Q_{0}; and U is the control input vector. The objective is to find the control function (functional) U(k) for k = 1, 2, . . . that minimizes an expected quadratic cost function for the linear stochastic regulator defined by Eqs. (1) and (2). Because linear stochastic tracking problems can be formulated as linear stochastic regulator problems by combining the reference and plant models in an augmented system [3], this system representation applies equally to the tracking problem. PROBLEM FORMULATION AND SOLUTION The system representation defined by Eqs. (1) and (2) can be cast in stochastic regulator form by augmenting the state vector with the timedelayed states, that is, $\stackrel{\_}{X}\ue8a0\left(k+1\right)=\stackrel{\_}{A}\ue8a0\left(k\right)\ue89e\stackrel{\_}{X}\ue8a0\left(k\right)+{\stackrel{\_}{B}\ue8a0\left(k\right)}^{*}\ue89eU\ue8a0\left(k\right)+{\stackrel{\_}{W}}_{1}\ue8a0\left(k\right),$ (3) and $Y\ue8a0\left(k\right)=\stackrel{\_}{C}\ue8a0\left(k\right)\ue89e\stackrel{\_}{X}\ue8a0\left(k\right)+{W}_{2}\ue8a0\left(k\right)$ (4) $\begin{array}{c}\stackrel{\_}{A}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{ccccc}0& I& 0& \cdots & 0\\ \text{\hspace{1em}}& 0& \text{\hspace{1em}}& \ddots & \text{\hspace{1em}}\\ \vdots & \text{\hspace{1em}}& \ddots & \text{\hspace{1em}}& \vdots \\ \text{\hspace{1em}}& \ddots & \text{\hspace{1em}}& \text{\hspace{1em}}& \text{\hspace{1em}}\\ 0& \text{\hspace{1em}}& \cdots & 0& I\\ {A}_{k}\ue8a0\left(k\right)& \text{\hspace{1em}}& \cdots & \text{\hspace{1em}}& {A}_{1}\ue8a0\left(k\right)\end{array}\right],\stackrel{\_}{X}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{c}X\ue8a0\left(kh+1\right)\\ X\ue8a0\left(kh+2\right)\\ \vdots \\ X\ue8a0\left(k1\right)\\ X\ue8a0\left(k\right)\end{array}\right],\\ \stackrel{\_}{B}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{c}0\\ \vdots \\ 0\\ B\ue8a0\left(k\right)\end{array}\right],\stackrel{\_}{W}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{c}0\\ \vdots \\ 0\\ {W}_{1}\ue8a0\left(k\right)\end{array}\right],{\stackrel{\_}{C}\ue8a0\left(k\right)}^{T}\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{c}{C}_{h}\ue8a0\left(k\right)\\ {C}_{h1}\ue8a0\left(k\right)\\ \vdots \\ {C}_{2}\ue8a0\left(k\right)\\ {C}_{1}\ue8a0\left(k\right)\end{array}\right],\end{array}\hspace{1em}$ and the dimensionality of the system is defined by h, where h = max(p, q) In Eq. (3), the control remains the same as in Eq. (1) because; physically, past states cannot be changed or controlled. In Eq. (4) the observation Y and the observation noise W_{2 }also remain unchanged. In Eqs. (1) and (2), p and q are not, in general. equal. For example. if the observation is not a linear function of all the delayed states contained in {overscore (X)}. then h = q > p. and the appropriate (q − p) submatrices in {overscore (C)}. of Eq. (4). are set equal to zero. If the process evolution is not a linear function of all the delayed states contained in {overscore (X)}, then h = p > q, and the appropriate submatrices in {overscore (A)}. of Eq. (3). are set equal to zero. The dimensions of {overscore (B)} and {overscore (W)}_{1},in Eq. (3), must also be consistent with the integer h. Therefore, without loss of generality, the system defined by Eqs. (1) and (2) can be represented by the augmented system model described by Eqs. (3) and (4). Given that the control is to minimize the expected value of a quadratic cost function of the form, $E\ue89e\left\{\sum _{k={k}_{0}}^{{k}_{1}1}\ue89e\left[{\stackrel{\_}{X}\ue8a0\left(k\right)}^{T}\ue89e{\stackrel{\_}{R}}_{1}\ue8a0\left(k\right)\ue89e\stackrel{\_}{X}\ue8a0\left(k\right)+{U}^{T}\ue8a0\left(k\right)\ue89e{R}_{2}\ue8a0\left(k\right)\ue89eU\ue8a0\left(k\right)\right]+\stackrel{\_}{X}\ue8a0\left({k}_{1}\right)\ue89e\stackrel{\_}{P}\ue8a0\left({k}_{1}\right)\ue89e\stackrel{\_}{X}\ue8a0\left({k}_{1}\right)\right\},$ (5) then the separation principle applies [3]. Also, the optimal linear stochastic control of the augmented system is given by a deterministic, optimal linear controller with state input {overscore (X)} (or estimated state feedback), which is provided by an optimal onestep predictor using the augmented model [4]. That is, $U\ue8a0\left(k\right)=\stackrel{\_}{F}\ue8a0\left(k\right)\ue89e\hat{\stackrel{\_}{X}}\ue8a0\left(k\right);\text{\hspace{1em}}\ue89ek={k}_{0},{k}_{0}+1,\dots ,{k}_{1},$ (6) where the control gain {overscore (F)} satisfies $\begin{array}{c}\stackrel{\_}{F}\ue8a0\left(k\right)=\ue89e\{{R}_{2}\ue8a0\left(k\right)+{\stackrel{\_}{B}}^{T}\ue8a0\left(k\right)[{\stackrel{\_}{R}}_{1}\ue8a0\left(k+1\right)+\stackrel{\_}{P}(k+\\ {\ue89e1]\ue89eB\ue8a0\left(k\right)\}}^{1}\ue89e{B}^{T}\ue8a0\left(k\right)\ue8a0\left[\stackrel{\_}{R}\ue8a0\left(k+1\right)+P\ue8a0\left(k+1\right)\right]\ue89e\stackrel{\_}{A}\ue8a0\left(k\right);\end{array}\hspace{1em}$ (7) the matrix {overscore (P)} satisfies the recursive matrix Riccati equation, $\stackrel{\_}{P}\ue8a0\left(k\right)={\stackrel{\_}{A}}^{1}\ue8a0\left(k\right)[{\stackrel{\_}{R}}_{1}\ue8a0\left(k+1\right)+\stackrel{\_}{P}\left(k+1\right]\ue8a0\left[\stackrel{\_}{A}\ue8a0\left(k\right)\stackrel{\_}{B}\ue8a0\left(k\right)\ue89e\stackrel{\_}{F}\ue8a0\left(k\right)\right];$ (8) $\mathrm{the}\ue89e\text{\hspace{1em}}\ue89e\mathrm{one}\ue89e\text{}\ue89e\mathrm{step}\ue89e\text{\hspace{1em}}\ue89e\mathrm{predictor}\ue89e\text{\hspace{1em}}\ue89e\mathrm{output}\ue89e\text{\hspace{1em}}\ue89e\hat{\stackrel{\_}{X}}\ue89e\text{\hspace{1em}}\ue89e\mathrm{is}\ue89e\text{\hspace{1em}}\ue89e\mathrm{defined}\ue89e\text{\hspace{1em}}\ue89e\mathrm{by}$ $\hat{\stackrel{\_}{X}}\ue8a0\left(k+1\right)=\stackrel{\_}{A}\ue8a0\left(k\right)\ue89e\hat{\stackrel{\_}{X}}\ue8a0\left(k\right)+\stackrel{\_}{B}\ue8a0\left(k\right)\ue89eU\ue8a0\left(k\right)+\stackrel{\_}{K}\ue8a0\left(k\right)\ue8a0\left[Y\ue8a0\left(k\right)\stackrel{\_}{C}\ue8a0\left(k\right)\ue89e\hat{\stackrel{\_}{X}}\ue8a0\left(k\right)\right];$ (9) the estimator gain {overscore (K)} satisfies $\stackrel{\_}{K}\ue8a0\left(k\right)=\stackrel{\_}{A}\ue8a0\left(k\right)\ue89e\stackrel{\_}{Q}\ue8a0\left(k\right)\ue89e{{\stackrel{\_}{C}}^{T}\ue8a0\left(k\right)\ue8a0\left[\stackrel{\_}{C}\ue8a0\left(k\right)\ue89e\stackrel{\_}{Q}\ue8a0\left(k\right)\ue89e{\stackrel{\_}{C}}^{T}\ue8a0\left(k\right)+{V}_{2}\ue8a0\left(k\right)\right]}^{1};$ (10) and the state estimation error covariance matrix {overscore (Q)} satisfies the recursive matrix Riccati equation, $\stackrel{\_}{Q}\ue8a0\left(k+1\right)=\left[\stackrel{\_}{A}\ue8a0\left(k\right)\stackrel{\_}{K}\ue8a0\left(k\right)\ue89e\stackrel{\_}{C}\ue8a0\left(k\right)\right]\ue89e\stackrel{\_}{Q}\ue8a0\left(k\right)\ue89e{\stackrel{\_}{A}}^{T}\ue8a0\left(k\right)+{\stackrel{\_}{V}}_{1}\ue8a0\left(k\right).$ (11) The final value of {overscore (P)} used to “initialize” Eq. (8) (which is solved backward in time) is the final value defined in the quadratic cost function of Eq. (5), that is, {overscore (P)}(k_{1}) = {overscore (P)}_{1}. The initial value of {overscore (Q)} used to initialize Eq. (11) is the error covariance of the initial estimate of X, that is, {overscore (Q)}(0) = {overscore (Q)}_{0}. If the system statistics defined in Eqs. (1) and (2) are gaussian, then the above solution is the optimal solution without qualification; if not, then it is the optimal linear control solution. The expected system performance is determined by analyzing the augmented system as a linear, stochastic regulator problem PRACTICAL APPLICATIONS From control gain Eqs. (7) and (8) and estimator gain Eqs. (10) and (11), it can be shown that the required dimensions of the controller and estimator are not, in general, equal. The dimensions of the control matrix Riccati equation (8) is determined by the number of delayed states in the process evolution, that is, {overscore (P)} has the dimensions of n · q × n · q. Consequently, for the special case of measurement delays only (q = 1), the controller implementation is unaffected. However, the dimensions of the estimator matrix Riccati equation (11) are determined by the maximum of p and q. Since h = max(p, q), then {overscore (Q)} has the dimensions n · h × n · h. For the general case, the optimal control is given by $U\ue8a0\left(k\right)=\sum _{i=1}^{q}\ue89e{F}_{i}\ue8a0\left(k\right)\ue89e\hat{X}\ue8a0\left(k+1i\right),$ (12) where $\stackrel{\_}{F}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{ccccc}{F}_{q}\ue8a0\left(k\right)& {F}_{q1}\ue8a0\left(k\right)& \dots & {F}_{2}\ue8a0\left(k\right)& {F}_{1}\ue8a0\left(k\right)\end{array}\right],$ (13) and ${\hat{\stackrel{\_}{X}}}^{T}\ue8a0\left(k\right)\ue89e\stackrel{\Delta}{=}\ue89e\left[\begin{array}{cccccc}{\hat{X}}^{T}\ue8a0\left(k+1h\right)& \dots & {\hat{X}}^{T}\ue8a0\left(k+1q\right)& \dots & {\hat{X}}^{T}\ue8a0\left(k1\right)& {\hat{X}}^{T}\ue8a0\left(k\right)\end{array}\right].$ (14) We note that, in Eq. (14), {circumflex over (X)}(k) is the onestep predicted value of the original system state vector; {circumflex over (X)}(k − 1) is the filtered value; {circumflex over (X)}(k − 2) is the onestep smoothed value; and, finally, {circumflex over (X)}(k + 1 − h) is the (h − 2)th smoothed value. Typically, the realtime computational requirements associated with the implementation of timevarying optimal control are always stressing due to the “backwardintime” recursion that is required to obtain solutions {overscore (P)} of Eq. (8). However, in many highaccuracy applications, the {overscore (A)}, {overscore (B)}, {overscore (R)}_{1}, and R_{2} matrices of the augmented system and cost function can be treated as time invariant over the time intervals of interest. Further, if the augmented system satisfies the relatively minor requirements of stabilizability and detectability, then the control gain {overscore (F)} will converge to a unique value such that the steady state optimal control law is time invariant, asymptotically stable, and mini mizes the quadratic cost function of Eq (5) as k_{1 }→ ∞. For this case, the steadystate gain matrix {overscore (F)}_{ss }can be computed offline and stored for real time use; the realtime computations required to implement this steadystate control law are negligible. Usually, for tracking and regulator problems, the steadystate control gains and timevarying gains are such that the initial value of the timevarying gain is equal to the steadystate gain but is less than or equal to the steadystate gain as time progresses, that is, ${\stackrel{\_}{F}}_{\mathrm{ss}}\ge \stackrel{\_}{F}\ue8a0\left(k\right)\ue89e\text{\hspace{1em}}\ue89e\mathrm{for}\ue89e\text{\hspace{1em}}\ue89e{k}_{o}\le k\le {k}_{1}.$ As a consequence, the steadystate gain tends to maintain the controlled state closer to the estimated state, but at the cost of control energy. If accuracy is the significant criterion, then the steadystate gain is not only easy to implement but also provides essentially equivalent or better accuracy. This, however, is not the case for the estimator gain, even though in most cases where linear steadystate optimal stochastic control is implemented, both the steadystate control and the estimator gains are used. For the es timator, and the same general conditions as before, the initial timevarying estimator gain is usually significantly larger than the steadystate gain, to account for initial uncertainties in the knowledge of the system state. As time progresses, the timevarying gain converges to the steadystate gain, that is, ${\stackrel{\_}{K}}_{\mathrm{ss}}\le \stackrel{\_}{K}\ue8a0\left(k\right)\ue89e\text{\hspace{1em}}\ue89e\mathrm{for}\ue89e\text{\hspace{1em}}\ue89e{k}_{o}\le k\le {k}_{1}.$ Consequently, initial system performance is significantly degraded if K_{ss }is used. In fact, if applied to a “linearized” system, the estimator can actually diverge given the initial small steadystate filter gains. Hence a good compromise between performance and computational complexity is to choose the steadystate controller with the timevarying estimator. Further, because the estimator matrix Riccati equation is solved forward in time, the computations associated with the timevarying filter gain are orders of magnitude less than with the timevarying control gain and can usually be implemented in real time. Although the computations associated with the augmented system es timator are significantly increased because of the increased dimensions, some simplifications can be made. The estimator Riccati equation can be separated into a time update and a measurement update, where the time update for the augmented system becomes primarily one of data transfer. If {overscore (Q)}, from Eq. (11), is defined as $Q\ue8a0\left(k\right)=\left[\begin{array}{cccc}{Q}_{1,1}\ue8a0\left(k\right)& {Q}_{1,2}\ue8a0\left(k\right)& \cdots & {Q}_{1,h}\ue8a0\left(k\right)\\ {Q}_{2,1}\ue8a0\left(k\right)& {Q}_{2,2}\ue8a0\left(k\right)& \cdots & \vdots \\ \vdots & \text{\hspace{1em}}& \ddots & \text{\hspace{1em}}\\ {Q}_{h,1}\ue8a0\left(k\right)& \cdots & \text{\hspace{1em}}& {Q}_{h,h}\ue8a0\left(k\right)\end{array}\right],\hspace{1em}$ (15) then, for the special case where q = 1 and p ≧ 2, the time update becomes $\begin{array}{c}\begin{array}{c}{Q}_{i,j}\ue8a0\left(k+1\right)={Q}_{i+1,j+1}\ue8a0\left(k\right)\\ {Q}_{i,P}\ue8a0\left(k+1\right)={Q}_{i+1,P}\ue8a0\left(k\right)\ue89e{A}_{1}^{T}\ue8a0\left(k\right)\end{array}\}\ue89e\text{\hspace{1em}}\ue89ei,j=1,\dots \ue89e\text{\hspace{1em}},p1\\ {Q}_{p,p}\ue8a0\left(k+1\right)={A}_{1}\ue8a0\left(k\right)\ue89e{Q}_{p,p}\ue8a0\left(k\right)\ue89e{A}_{1}^{T}\ue8a0\left(k\right)+{V}_{1}\ue8a0\left(k\right)\end{array}\hspace{1em}\}.$ (16) For the worst case, where q > p > I, the time update is given by $\begin{array}{c}{Q}_{i,j}\ue8a0\left(k+1\right)={Q}_{i+1,j+1}\ue8a0\left(k\right),\text{\hspace{1em}}\ue89ei,j=1,\dots \ue89e\text{\hspace{1em}},q1,\\ {Q}_{i,q}\ue8a0\left(k+1\right)=\sum _{j=1}^{q}\ue89e{Q}_{i,j}\ue8a0\left(k\right)\ue89e{A}_{\left(q+1\right)j}^{T}\ue8a0\left(k\right),\text{\hspace{1em}}\ue89ei=1,\dots \ue89e\text{\hspace{1em}},q1,\\ {Q}_{q,q}\ue8a0\left(k+1\right)=\sum _{i=1}^{q}\ue89e\sum _{j=1}^{q}\ue89e{A}_{\left(q+1\right)i}\ue8a0\left(k\right)\ue89e{Q}_{i,j}\ue8a0\left(k\right)\ue89e{A}_{\left(q+1\right)j}^{T}\ue8a0\left(k\right)+{V}_{1}\ue8a0\left(k\right)\end{array}\}.$ (17) In practice, another computational simplification results because rarely, if ever, is the actual system measurement a function of every timedelayed state element, and likewise for the actual system process model. Hence the submatrices A_{i }and C_{i }(for i > 1) in {overscore (A)} and {overscore (C)}, respectively, are usually of significantly reduced dimension. This, in turn, significantly reduces the dimensions of the augmented system model. We note that the iteration interval Δt (implied in the discretetime system representation) or some integer number of Δt's should be set equal to the time delay. For variable time delays, it may be advantageous to use (a) variable iteration intervals for the estimator, (b) fixed iteration intervals for the controller, and (c) variabletime updates of the state estimate to time synch the estimator output with the controller input. Finally, for the analogous continuoustime problem, the filter equation becomes a partial differential equation with a boundary condition, and the covanance equation becomes a partial differential matrix equation with three boundary conditions (see [2]). Consequently, the most practical control solution is obtained by discretizing the continuoustime system representation and then applying the approach of Section III. One technique for discretizing the continuoustime representation is by using the Z transform method; see, for example, [5]. Using a scalar differential equation with one delayed slate as an example, we have $\stackrel{.}{X}\ue8a0\left(t\right)=\sum _{i=0}^{1}\ue89e{a}_{i}\ue89eX\ue8a0\left(ti\ue89e\text{\hspace{1em}}\ue89e\mathrm{\Delta t}\right)+U\ue8a0\left(t\right).$ (18) If we assume that the output solution X(t) is sampled in discrete time, and if the system process can be reasonably approximated by a continuous system in which {dot over (X)} is driven by the output of a high frequency sampling of the right side of Eq. (18), then the discretetime process model is given by $X\ue8a0\left(k+1\right)=\left(1+{a}_{0}\ue89e\mathrm{\Delta t}\right)\ue89eX\ue8a0\left(k\right)+\left({a}_{1}\ue89e\mathrm{\Delta t}\right)\ue89eX\ue8a0\left(k\mathrm{\Delta t}\right)+\left(\mathrm{\Delta t}\right)\ue89eU\ue8a0\left(k\right).$ (19) The approximate discretetime solution to Eq. (18), which is given by Eq. (19), is now in a form consistent with Eq. (1) of Section II, and the ap proach of Section III can be applied.
Claims (4)
 1. A method for closeloop control of a physiological parameter comprising the acts of:obtaining a measurement of the physiological parameter from a patient;providing the measurement to an optimal estimator in real time, wherein the optimal estimator outputs a best estimate of the physiological parameter in real time based on the measurement;providing the best estimate of the physiological parameter to an optimal controller in real time, wherein the optimal controller outputs a control command in real time based on the best estimate of the physiological parameter and a control reference; andproviding the control command to an actuator, wherein the actuator provides an output to adjust the physiological parameter.
 2. The method of
claim 1 , wherein the measurement is obtained using a sensor.  3. The method of
claim 1 , wherein the optimal estimator is implemented using a linearized Kalman algorithm.  4. The method of
claim 1 , wherein the control reference is provided by a patient health monitor.
Priority Applications (3)
Application Number  Priority Date  Filing Date  Title 

US23463200 true  20000922  20000922  
US09960855 US6572545B2 (en)  20000922  20010921  Method and apparatus for realtime control of physiological parameters 
US10417438 US20030195404A1 (en)  20000922  20030416  Method and apparatus for realtime control of physiological parameters 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US10417438 US20030195404A1 (en)  20000922  20030416  Method and apparatus for realtime control of physiological parameters 
Related Parent Applications (1)
Application Number  Title  Priority Date  Filing Date  

US09960855 Continuation US6572545B2 (en)  20000922  20010921  Method and apparatus for realtime control of physiological parameters 
Publications (1)
Publication Number  Publication Date 

US20030195404A1 true true US20030195404A1 (en)  20031016 
Family
ID=22882159
Family Applications (3)
Application Number  Title  Priority Date  Filing Date 

US09960855 Expired  Fee Related US6572545B2 (en)  20000922  20010921  Method and apparatus for realtime control of physiological parameters 
US09960846 Expired  Fee Related US6575905B2 (en)  20000922  20010921  Method and apparatus for realtime estimation of physiological parameters 
US10417438 Abandoned US20030195404A1 (en)  20000922  20030416  Method and apparatus for realtime control of physiological parameters 
Family Applications Before (2)
Application Number  Title  Priority Date  Filing Date 

US09960855 Expired  Fee Related US6572545B2 (en)  20000922  20010921  Method and apparatus for realtime control of physiological parameters 
US09960846 Expired  Fee Related US6575905B2 (en)  20000922  20010921  Method and apparatus for realtime estimation of physiological parameters 
Country Status (2)
Country  Link 

US (3)  US6572545B2 (en) 
WO (1)  WO2002024065A1 (en) 
Cited By (10)
Publication number  Priority date  Publication date  Assignee  Title 

WO2005081171A2 (en) *  20040219  20050901  Edward Henry Mathews  Insulin bolus calculator for mobile communication device 
WO2008057384A2 (en) *  20061102  20080515  The University Of North Carolina At Chapel Hill  Methods and systems for determining an intravenous insulin infusion rate 
US20080214919A1 (en) *  20061226  20080904  Lifescan, Inc.  System and method for implementation of glycemic control protocols 
US20100249561A1 (en) *  20070621  20100930  University Of Virginia Patent Foundation  LQG Artificial Pancreas Control System and Related Method 
US20110238350A1 (en) *  20100329  20110929  Commissariat A L'energie Atomique Et Aux Energies Alternatives  Calibrating method and system, recording medium for this method 
US20110264378A1 (en) *  20081126  20111027  University Of Virginia Patent Foundation  Method, System, and Computer Program Product For Tracking of Blood Glucose Variability in Diabetes 
US20110282321A1 (en) *  20070131  20111117  Medtronic Minimed, Inc.  Model predictive method and system for controlling and supervising insulin infusion 
US8417311B2 (en)  20080912  20130409  Optiscan Biomedical Corporation  Fluid component analysis system and method for glucose monitoring and control 
US8449524B2 (en)  20071010  20130528  Optiscan Biomedical Corporation  Fluid component analysis systems and methods for glucose monitoring and control 
WO2017136218A1 (en) *  20160201  20170810  Dexcom, Inc.  System and method for decision support using lifestyle factors 
Families Citing this family (233)
Publication number  Priority date  Publication date  Assignee  Title 

US6036924A (en)  19971204  20000314  HewlettPackard Company  Cassette of lancet cartridges for sampling blood 
US6391005B1 (en)  19980330  20020521  Agilent Technologies, Inc.  Apparatus and method for penetration with shaft having a sensor for sensing penetration depth 
US6175752B1 (en)  19980430  20010116  Therasense, Inc.  Analyte monitoring device and methods of use 
US8480580B2 (en)  19980430  20130709  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US9066695B2 (en)  19980430  20150630  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US8688188B2 (en)  19980430  20140401  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US8465425B2 (en)  19980430  20130618  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US8346337B2 (en)  19980430  20130101  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US8974386B2 (en)  19980430  20150310  Abbott Diabetes Care Inc.  Analyte monitoring device and methods of use 
US6430525B1 (en) *  20000605  20020806  Masimo Corporation  Variable mode averager 
WO2002024065A1 (en) *  20000922  20020328  Knobbe, Lim & Buckingham  Method and apparatus for realtime estimation and control of pysiological parameters 
US8641644B2 (en)  20001121  20140204  SanofiAventis Deutschland Gmbh  Blood testing apparatus having a rotatable cartridge with multiple lancing elements and testing means 
US6560471B1 (en)  20010102  20030506  Therasense, Inc.  Analyte monitoring device and methods of use 
US7892183B2 (en)  20020419  20110222  Pelikan Technologies, Inc.  Method and apparatus for body fluid sampling and analyte sensing 
WO2002078512A8 (en)  20010402  20041202  Therasense Inc  Blood glucose tracking apparatus and methods 
US7674232B2 (en)  20020419  20100309  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US9795334B2 (en)  20020419  20171024  SanofiAventis Deutschland Gmbh  Method and apparatus for penetrating tissue 
US7316700B2 (en)  20010612  20080108  Pelikan Technologies, Inc.  Self optimizing lancing device with adaptation means to temporal variations in cutaneous properties 
US7491178B2 (en)  20020419  20090217  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7547287B2 (en)  20020419  20090616  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7226461B2 (en)  20020419  20070605  Pelikan Technologies, Inc.  Method and apparatus for a multiuse body fluid sampling device with sterility barrier release 
US7331931B2 (en)  20020419  20080219  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7291117B2 (en)  20020419  20071106  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7232451B2 (en)  20020419  20070619  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7682318B2 (en)  20010612  20100323  Pelikan Technologies, Inc.  Blood sampling apparatus and method 
US8221334B2 (en)  20020419  20120717  SanofiAventis Deutschland Gmbh  Method and apparatus for penetrating tissue 
US8267870B2 (en)  20020419  20120918  SanofiAventis Deutschland Gmbh  Method and apparatus for body fluid sampling with hybrid actuation 
US9248267B2 (en)  20020419  20160202  SanofiAventis Deustchland Gmbh  Tissue penetration device 
US7699791B2 (en)  20010612  20100420  Pelikan Technologies, Inc.  Method and apparatus for improving success rate of blood yield from a fingerstick 
US7717863B2 (en)  20020419  20100518  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7229458B2 (en)  20020419  20070612  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US8784335B2 (en)  20020419  20140722  SanofiAventis Deutschland Gmbh  Body fluid sampling device with a capacitive sensor 
US7648468B2 (en)  20020419  20100119  Pelikon Technologies, Inc.  Method and apparatus for penetrating tissue 
US7371247B2 (en)  20020419  20080513  Pelikan Technologies, Inc  Method and apparatus for penetrating tissue 
US9314194B2 (en)  20020419  20160419  SanofiAventis Deutschland Gmbh  Tissue penetration device 
US7901362B2 (en)  20020419  20110308  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US7981056B2 (en)  20020419  20110719  Pelikan Technologies, Inc.  Methods and apparatus for lancet actuation 
WO2002100460A3 (en)  20010612  20030508  Don Alden  Electric lancet actuator 
US7041068B2 (en)  20010612  20060509  Pelikan Technologies, Inc.  Sampling module device and method 
US7909778B2 (en)  20020419  20110322  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US9226699B2 (en)  20020419  20160105  SanofiAventis Deutschland Gmbh  Body fluid sampling module with a continuous compression tissue interface surface 
US8579831B2 (en)  20020419  20131112  SanofiAventis Deutschland Gmbh  Method and apparatus for penetrating tissue 
US7297122B2 (en)  20020419  20071120  Pelikan Technologies, Inc.  Method and apparatus for penetrating tissue 
US8337419B2 (en)  20020419  20121225  SanofiAventis Deutschland Gmbh  Tissue penetration device 
US9427532B2 (en)  20010612  20160830  SanofiAventis Deutschland Gmbh  Tissue penetration device 
EP1404235A4 (en)  20010612  20080820  Pelikan Technologies Inc  Method and apparatus for lancet launching device integrated onto a bloodsampling cartridge 
US7976476B2 (en)  20020419  20110712  Pelikan Technologies, Inc.  Device and method for variable speed lancet 
US7344507B2 (en)  20020419  20080318  Pelikan Technologies, Inc.  Method and apparatus for lancet actuation 
US7451126B2 (en) *  20011018  20081111  Omron Corporation  State space navigation system, user system and business methods for machine to machine business 
US9282925B2 (en)  20020212  20160315  Dexcom, Inc.  Systems and methods for replacing signal artifacts in a glucose sensor data stream 
US8702624B2 (en)  20060929  20140422  SanofiAventis Deutschland Gmbh  Analyte measurement device with a single shot actuator 
CA2484265C (en) *  20020506  20120807  Becton, Dickinson And Company  Method and device for controlling drug pharmacokinetics 
US8423113B2 (en)  20030725  20130416  Dexcom, Inc.  Systems and methods for processing sensor data 
US8260393B2 (en)  20030725  20120904  Dexcom, Inc.  Systems and methods for replacing signal data artifacts in a glucose sensor data stream 
US8574895B2 (en)  20021230  20131105  SanofiAventis Deutschland Gmbh  Method and apparatus using optical techniques to measure analyte levels 
US7016715B2 (en) *  20030113  20060321  Nellcorpuritan Bennett Incorporated  Selection of preset filter parameters based on signal quality 
US6949816B2 (en)  20030421  20050927  Motorola, Inc.  Semiconductor component having first surface area for electrically coupling to a semiconductor chip and second surface area for electrically coupling to a substrate, and method of manufacturing same 
US20040253736A1 (en) *  20030606  20041216  Phil Stout  Analytical device with prediction module and related methods 
US7850621B2 (en)  20030606  20101214  Pelikan Technologies, Inc.  Method and apparatus for body fluid sampling and analyte sensing 
US8066639B2 (en)  20030610  20111129  Abbott Diabetes Care Inc.  Glucose measuring device for use in personal area network 
US8886273B2 (en)  20030801  20141111  Dexcom, Inc.  Analyte sensor 
US8275437B2 (en)  20030801  20120925  Dexcom, Inc.  Transcutaneous analyte sensor 
US8369919B2 (en)  20030801  20130205  Dexcom, Inc.  Systems and methods for processing sensor data 
US8160669B2 (en)  20030801  20120417  Dexcom, Inc.  Transcutaneous analyte sensor 
US7933639B2 (en)  20030801  20110426  Dexcom, Inc.  System and methods for processing analyte sensor data 
US7774145B2 (en)  20030801  20100810  Dexcom, Inc.  Transcutaneous analyte sensor 
US8761856B2 (en)  20030801  20140624  Dexcom, Inc.  System and methods for processing analyte sensor data 
US8886272B2 (en)  20040713  20141111  Dexcom, Inc.  Analyte sensor 
US8845536B2 (en)  20030801  20140930  Dexcom, Inc.  Transcutaneous analyte sensor 
US20100168542A1 (en)  20030801  20100701  Dexcom, Inc.  System and methods for processing analyte sensor data 
US7778680B2 (en)  20030801  20100817  Dexcom, Inc.  System and methods for processing analyte sensor data 
US9247901B2 (en)  20030822  20160202  Dexcom, Inc.  Systems and methods for replacing signal artifacts in a glucose sensor data stream 
US8010174B2 (en)  20030822  20110830  Dexcom, Inc.  Systems and methods for replacing signal artifacts in a glucose sensor data stream 
US8233959B2 (en)  20030822  20120731  Dexcom, Inc.  Systems and methods for processing analyte sensor data 
DE10343863A1 (en) *  20030923  20050414  Roche Diagnostics Gmbh  Method and apparatus for continuously monitoring the concentration of an analyte 
WO2005033659A3 (en)  20030929  20070118  Pelikan Technologies Inc  Method and apparatus for an improved sample capture device 
WO2005037095A1 (en)  20031014  20050428  Pelikan Technologies, Inc.  Method and apparatus for a variable user interface 
US20050119866A1 (en) *  20031114  20050602  Zaleski John R.  Medical parameter processing system 
US7519408B2 (en)  20031119  20090414  Dexcom, Inc.  Integrated receiver for continuous analyte sensor 
US8287453B2 (en)  20031205  20121016  Dexcom, Inc.  Analyte sensor 
EP2239567B1 (en)  20031205  20150902  DexCom, Inc.  Calibration techniques for a continuous analyte sensor 
EP1711791B1 (en)  20031209  20141015  DexCom, Inc.  Signal processing for continuous analyte sensor 
WO2005065414A3 (en)  20031231  20051229  Pelikan Technologies Inc  Method and apparatus for improving fluidic flow and sample capture 
US20050182500A1 (en) *  20040217  20050818  Continuous Control Solutions, Inc.  Time delay definition 
US7591801B2 (en)  20040226  20090922  Dexcom, Inc.  Integrated delivery device for continuous glucose sensor 
US8808228B2 (en)  20040226  20140819  Dexcom, Inc.  Integrated medicament delivery device for use with continuous analyte sensor 
US8792955B2 (en)  20040503  20140729  Dexcom, Inc.  Transcutaneous analyte sensor 
EP1751546A2 (en)  20040520  20070214  Albatros Technologies GmbH & Co. KG  Printable hydrogel for biosensors 
EP1765194A4 (en)  20040603  20100929  Pelikan Technologies Inc  Method and apparatus for a fluid sampling device 
WO2006001797A1 (en)  20040614  20060105  Pelikan Technologies, Inc.  Low pain penetrating 
US7755503B2 (en) *  20040616  20100713  Dell Products L.P.  System and method for error messaging of an incorrect information handling system graphics cable connection 
US7783333B2 (en)  20040713  20100824  Dexcom, Inc.  Transcutaneous medical device with variable stiffness 
US7497827B2 (en)  20040713  20090303  Dexcom, Inc.  Transcutaneous analyte sensor 
US8565848B2 (en)  20040713  20131022  Dexcom, Inc.  Transcutaneous analyte sensor 
EP2327984B8 (en) *  20040713  20150422  DexCom, Inc.  Transcutaneous analyte sensor 
US8452368B2 (en)  20040713  20130528  Dexcom, Inc.  Transcutaneous analyte sensor 
US9247900B2 (en)  20040713  20160202  Dexcom, Inc.  Analyte sensor 
US20060016700A1 (en)  20040713  20060126  Dexcom, Inc.  Transcutaneous analyte sensor 
US20080027679A1 (en) *  20040721  20080131  Dror Shklarski  Wearable Device, System and Method for Measuring Physiological and/or Environmental Parameters 
US8396670B2 (en) *  20040816  20130312  Venture Milling, Inc.  Process, system and method for improving the determination of digestive effects upon an ingestable substance 
FR2874807B1 (en) *  20040903  20061110  Memscap Sa  An analysis of the physical and chemical properties of the skin 
WO2006024672A1 (en) *  20040903  20060309  Novo Nordisk A/S  System and method for estimating the glucose concentration in blood 
US7608042B2 (en) *  20040929  20091027  Intellidx, Inc.  Blood monitoring system 
US20070191716A1 (en) *  20040929  20070816  Daniel Goldberger  Blood monitoring system 
US20060229531A1 (en) *  20050201  20061012  Daniel Goldberger  Blood monitoring system 
US20060069436A1 (en) *  20040930  20060330  Depuy Spine, Inc.  Trial disk implant 
US8652831B2 (en)  20041230  20140218  SanofiAventis Deutschland Gmbh  Method and apparatus for analyte measurement test time 
US7822454B1 (en)  20050103  20101026  Pelikan Technologies, Inc.  Fluid sampling device with improved analyte detecting member configuration 
JP2008526443A (en)  20050113  20080724  ウェルチ・アリン・インコーポレーテッド  Vital signs monitor 
US7545272B2 (en)  20050208  20090609  Therasense, Inc.  RF tag on test strips, test strip vials and boxes 
US7920906B2 (en)  20050310  20110405  Dexcom, Inc.  System and methods for processing analyte sensor data for sensor calibration 
US7509156B2 (en) *  20050518  20090324  Clarian Health Partners, Inc.  System for managing glucose levels in patients with diabetes or hyperglycemia 
US8880138B2 (en)  20050930  20141104  Abbott Diabetes Care Inc.  Device for channeling fluid and methods of use 
US20070088271A1 (en) *  20051018  20070419  Richards Cynthia C  Medication device 
US7766829B2 (en)  20051104  20100803  Abbott Diabetes Care Inc.  Method and system for providing basal profile modification in analyte monitoring and management systems 
US20070123801A1 (en) *  20051128  20070531  Daniel Goldberger  Wearable, programmable automated blood testing system 
US20080200838A1 (en) *  20051128  20080821  Daniel Goldberger  Wearable, programmable automated blood testing system 
EP1954184B1 (en) *  20051129  20171011  PhysIQ Inc.  Residualbased monitoring of human health 
US8133178B2 (en)  20060222  20120313  Dexcom, Inc.  Analyte sensor 
US7826879B2 (en)  20060228  20101102  Abbott Diabetes Care Inc.  Analyte sensors and methods of use 
US7620438B2 (en)  20060331  20091117  Abbott Diabetes Care Inc.  Method and system for powering an electronic device 
US7801582B2 (en)  20060331  20100921  Abbott Diabetes Care Inc.  Analyte monitoring and management system and methods therefor 
US8226891B2 (en)  20060331  20120724  Abbott Diabetes Care Inc.  Analyte monitoring devices and methods therefor 
US20070243567A1 (en) *  20060413  20071018  Syhhong Chang  Beta cell mimicking control algorithm for artificial pancreas 
US8092385B2 (en) *  20060523  20120110  Intellidx, Inc.  Fluid access interface 
US20090209938A1 (en) *  20060529  20090820  Wristop Technologies Oy  Apparatus and method for dosing drug and wireless remote control of a drug pump 
US20080071157A1 (en)  20060607  20080320  Abbott Diabetes Care, Inc.  Analyte monitoring system and method 
CN101500628B (en) *  20060613  20130703  康尔福盛303有限公司  System and method for optimizing control of PCA and PCEA system 
EP1875933A1 (en) *  20060706  20080109  Debiotech S.A.  Medical device for delivery of a solution 
US7653425B2 (en)  20060809  20100126  Abbott Diabetes Care Inc.  Method and system for providing calibration of an analyte sensor in an analyte monitoring system 
US7618369B2 (en) *  20061002  20091117  Abbott Diabetes Care Inc.  Method and system for dynamically updating calibration parameters for an analyte sensor 
US8364231B2 (en)  20061004  20130129  Dexcom, Inc.  Analyte sensor 
US8388534B2 (en) *  20061011  20130305  Samsung Electronics Co., Ltd.  Apparatus for providing skin care information by measuring skin moisture content and method and medium for the same 
US20080306353A1 (en) *  20061103  20081211  Douglas Joel S  Calculation device for metabolic control of critically ill and/or diabetic patients 
US8732188B2 (en)  20070218  20140520  Abbott Diabetes Care Inc.  Method and system for providing contextual based medication dosage determination 
US8930203B2 (en)  20070218  20150106  Abbott Diabetes Care Inc.  Multifunction analyte test device and methods therefor 
US9636450B2 (en)  20070219  20170502  Udo Hoss  Pump system modular components for delivering medication and analyte sensing at seperate insertion sites 
US8123686B2 (en)  20070301  20120228  Abbott Diabetes Care Inc.  Method and apparatus for providing rolling data in communication systems 
US9204827B2 (en)  20070414  20151208  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in medical communication system 
EP2137637A4 (en) *  20070414  20120620  Abbott Diabetes Care Inc  Method and apparatus for providing data processing and control in medical communication system 
CA2683959C (en)  20070414  20170829  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in medical communication system 
US9008743B2 (en)  20070414  20150414  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in medical communication system 
US7768387B2 (en)  20070414  20100803  Abbott Diabetes Care Inc.  Method and apparatus for providing dynamic multistage signal amplification in a medical device 
US7928850B2 (en)  20070508  20110419  Abbott Diabetes Care Inc.  Analyte monitoring system and methods 
US8461985B2 (en)  20070508  20130611  Abbott Diabetes Care Inc.  Analyte monitoring system and methods 
US8665091B2 (en)  20070508  20140304  Abbott Diabetes Care Inc.  Method and device for determining elapsed sensor life 
US8456301B2 (en)  20070508  20130604  Abbott Diabetes Care Inc.  Analyte monitoring system and methods 
US8260558B2 (en)  20070514  20120904  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US8600681B2 (en)  20070514  20131203  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US9125548B2 (en)  20070514  20150908  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US8239166B2 (en)  20070514  20120807  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US8103471B2 (en)  20070514  20120124  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US8140312B2 (en)  20070514  20120320  Abbott Diabetes Care Inc.  Method and system for determining analyte levels 
US8560038B2 (en)  20070514  20131015  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
US8444560B2 (en)  20070514  20130521  Abbott Diabetes Care Inc.  Method and apparatus for providing data processing and control in a medical communication system 
EP2152350A4 (en)  20070608  20130327  Dexcom Inc  Integrated medicament delivery device for use with continuous analyte sensor 
EP2006786B1 (en) *  20070618  20180509  Roche Diabetes Care GmbH  Method and glucose monitoring system for monitoring individual metabolic response and for generating nutritional feedback 
US8834366B2 (en)  20070731  20140916  Abbott Diabetes Care Inc.  Method and apparatus for providing analyte sensor calibration 
US20090143725A1 (en) *  20070831  20090604  Abbott Diabetes Care, Inc.  Method of Optimizing Efficacy of Therapeutic Agent 
DE102007046864A1 (en) *  20070928  20090409  Endress + Hauser Conducta Gesellschaft für Mess und Regeltechnik mbH + Co. KG  Measured value plausibility and quality evaluating method for e.g. control circuit, involves forming and evaluating derivative of measured value according to time, where measured value is compared with predefined warning limit value 
US9452258B2 (en) *  20071009  20160927  Dexcom, Inc.  Integrated insulin delivery system with continuous glucose sensor 
GB0719969D0 (en) *  20071012  20071121  Cambridge Entpr Ltd  Substance monitoring and control in human or animal bodies 
US8374668B1 (en)  20071023  20130212  Abbott Diabetes Care Inc.  Analyte sensor with lag compensation 
US8409093B2 (en)  20071023  20130402  Abbott Diabetes Care Inc.  Assessing measures of glycemic variability 
US8377031B2 (en) *  20071023  20130219  Abbott Diabetes Care Inc.  Closed loop control system with safety parameters and methods 
US8417312B2 (en)  20071025  20130409  Dexcom, Inc.  Systems and methods for processing sensor data 
US9135402B2 (en)  20071217  20150915  Dexcom, Inc.  Systems and methods for processing sensor data 
US9839395B2 (en)  20071217  20171212  Dexcom, Inc.  Systems and methods for processing sensor data 
US8473022B2 (en)  20080131  20130625  Abbott Diabetes Care Inc.  Analyte sensor with time lag compensation 
CA2715628A1 (en)  20080221  20090827  Dexcom, Inc.  Systems and methods for processing, transmitting and displaying sensor data 
CA2659142C (en)  20080321  20171003  Lifescan Scotland Limited  Analyte testing method and system 
US8396528B2 (en)  20080325  20130312  Dexcom, Inc.  Analyte sensor 
US8346335B2 (en)  20080328  20130101  Abbott Diabetes Care Inc.  Analyte sensor calibration management 
US20110028803A1 (en) *  20080331  20110203  Stig Ollmar  Method and device for noninvasive determination of the concentration of a substance in a body fluid 
US9220456B2 (en)  20080404  20151229  Hygieia, Inc.  Systems, methods and devices for achieving glycemic balance 
US8457901B2 (en) *  20080404  20130604  Hygieia, Inc.  System for optimizing a patient's insulin dosage regimen 
WO2009126900A1 (en)  20080411  20091015  Pelikan Technologies, Inc.  Method and apparatus for analyte detecting device 
WO2009139846A1 (en) *  20080512  20091119  Department Of Veterans Affairs  Automated system and method for diabetes control 
KR101148747B1 (en) *  20080716  20120522  파나소닉 주식회사  Battery pack 
WO2010019919A1 (en) *  20080814  20100218  University Of Toledo  Multifunctional neural network system and uses thereof for glycemic forecasting 
US8734422B2 (en)  20080831  20140527  Abbott Diabetes Care Inc.  Closed loop control with improved alarm functions 
US20100057040A1 (en)  20080831  20100304  Abbott Diabetes Care, Inc.  Robust Closed Loop Control And Methods 
US8622988B2 (en)  20080831  20140107  Abbott Diabetes Care Inc.  Variable rate closed loop control and methods 
US9943644B2 (en)  20080831  20180417  Abbott Diabetes Care Inc.  Closed loop control with reference measurement and methods thereof 
US9392969B2 (en)  20080831  20160719  Abbott Diabetes Care Inc.  Closed loop control and signal attenuation detection 
US8986208B2 (en)  20080930  20150324  Abbott Diabetes Care Inc.  Analyte sensor sensitivity attenuation mitigation 
US20100174228A1 (en) *  20081024  20100708  Bruce Buckingham  Hypoglycemia prediction and control 
US8992464B2 (en)  20081111  20150331  Hygieia, Inc.  Apparatus and system for diabetes management 
US8103456B2 (en)  20090129  20120124  Abbott Diabetes Care Inc.  Method and device for early signal attenuation detection using blood glucose measurements 
US9375169B2 (en)  20090130  20160628  SanofiAventis Deutschland Gmbh  Cam drive for managing disposable penetrating member actions with a single motor and motor and control system 
FR2942387A1 (en) *  20090223  20100827  Ecole Superieure Electricite  regulation system of blood glucose of a living subject and diabetic product corresponding computer program. 
US20120123234A1 (en) *  20090226  20120517  Eran Atlas  Method and system for automatic monitoring of diabetes related treatments 
WO2010111660A1 (en)  20090327  20100930  Dexcom, Inc.  Methods and systems for promoting glucose management 
US8753290B2 (en) *  20090327  20140617  Intellectual Inspiration, Llc  Fluid transfer system and method 
WO2010121084A1 (en)  20090415  20101021  Abbott Diabetes Care Inc.  Analyte monitoring system having an alert 
DE102009018649A1 (en)  20090423  20101028  Fresenius Medical Care Deutschland Gmbh  Dialysis machine, method and apparatus for a dialysis machine 
WO2010127050A1 (en)  20090428  20101104  Abbott Diabetes Care Inc.  Error detection in critical repeating data in a wireless sensor system 
EP2425209A4 (en)  20090429  20130109  Abbott Diabetes Care Inc  Method and system for providing real time analyte sensor calibration with retrospective backfill 
US20100279418A1 (en) *  20090504  20101104  Loren Robert Larson  Glucose meter adaptable for use with handheld devices, and associated communication network 
GB0909948D0 (en) *  20090610  20090722  Dna Electronics Ltd  A glucagon pump controller 
CN104799866A (en)  20090723  20150729  雅培糖尿病护理公司  The analyte monitoring device 
EP2456351B1 (en)  20090723  20161012  Abbott Diabetes Care, Inc.  Real time management of data relating to physiological control of glucose levels 
US8478557B2 (en)  20090731  20130702  Abbott Diabetes Care Inc.  Method and apparatus for providing analyte monitoring system calibration accuracy 
EP2473422A4 (en)  20090831  20140917  Abbott Diabetes Care Inc  Displays for a medical device 
EP2473098A4 (en) *  20090831  20140409  Abbott Diabetes Care Inc  Analyte signal processing device and methods 
EP2482720A4 (en)  20090929  20140423  Abbott Diabetes Care Inc  Method and apparatus for providing notification function in analyte monitoring systems 
CA2776007A1 (en)  20090930  20110407  Mor Research Applications Ltd.  Monitoring device for mangement of insulin delivery 
US8538740B2 (en) *  20091014  20130917  International Business Machines Corporation  Realtime performance modeling of software systems with multiclass workload 
US8771251B2 (en)  20091217  20140708  Hospira, Inc.  Systems and methods for managing and delivering patient therapy through electronic drug delivery systems 
US20110184265A1 (en) *  20100122  20110728  Abbott Diabetes Care Inc.  Method and Apparatus for Providing Notification in Analyte Monitoring Systems 
US8843321B2 (en) *  20100126  20140923  Roche Diagnostics Operations, Inc.  Methods and systems for processing glucose data measured from a person having diabetes 
US9089292B2 (en) *  20100326  20150728  Medtronic Minimed, Inc.  Calibration of glucose monitoring sensor and/or insulin delivery system 
US9398869B2 (en) *  20100326  20160726  University Of Virginia Patent Foundation  Method, system, and computer program product for improving the accuracy of glucose sensors using insulin delivery observation in diabetes 
US8965476B2 (en)  20100416  20150224  SanofiAventis Deutschland Gmbh  Tissue penetration device 
US9795747B2 (en)  20100602  20171024  SanofiAventis Deutschland Gmbh  Methods and apparatus for lancet actuation 
US8636711B2 (en) *  20100614  20140128  Legacy Emanuel Hospital & Health Center  Stabilized glucagon solutions and uses therefor 
CN106326651A (en) *  20110826  20170111  弗吉尼亚大学专利基金会  Method and system for adaptive advisory control of diabetes 
US8710993B2 (en)  20111123  20140429  Abbott Diabetes Care Inc.  Mitigating single point failure of devices in an analyte monitoring system and methods thereof 
US9317656B2 (en)  20111123  20160419  Abbott Diabetes Care Inc.  Compatibility mechanisms for devices in a continuous analyte monitoring system and methods thereof 
US9339217B2 (en)  20111125  20160517  Abbott Diabetes Care Inc.  Analyte monitoring system and methods of use 
US9510747B2 (en) *  20111130  20161206  Amo Development, Llc  System and method for ophthalmic surface measurements based on objective quality estimation 
WO2013082466A1 (en) *  20111130  20130606  Amo Development, Llc.  System and method for ophthalmic surface measurements based on sequential estimates 
US20130231543A1 (en) *  20120305  20130905  Andrea Facchinetti  Method to improve safety monitoring in type1 diabetic patients by detecting in realtime failures of the glucose 
US9119529B2 (en)  20121030  20150901  Dexcom, Inc.  Systems and methods for dynamically and intelligently monitoring a host's glycemic condition after an alert is triggered 
WO2014068743A1 (en) *  20121101  20140508  テルモ株式会社  Sensing device and sensing method 
WO2014109898A1 (en) *  20130114  20140717  The Regents Of University Of California  Modelbased personalization scheme of an artificial pancreas for type i diabetes applications 
US9227014B2 (en) *  20130207  20160105  The Board Of Trustees Of The Laland Stanford Junior University  Kalman filter based onoff switch for insulin pump 
WO2014194089A1 (en) *  20130529  20141204  Hospira, Inc.  Infusion system which utilizes one or more sensors and additional information to make an air determination regarding the infusion system 
US9898585B2 (en)  20140131  20180220  Aseko, Inc.  Method and system for insulin management 
US9486580B2 (en)  20140131  20161108  Aseko, Inc.  Insulin management 
US9892234B2 (en)  20141027  20180213  Aseko, Inc.  Subcutaneous outpatient management 
US9886556B2 (en)  20150820  20180206  Aseko, Inc.  Diabetes management therapy advisor 
CN107924715A (en) *  20160125  20180417  普雷迪科特医疗有限公司  A plurality of noninvasive estimation substance concentration sensor 
Citations (25)
Publication number  Priority date  Publication date  Assignee  Title 

US4655225A (en) *  19850418  19870407  Kurabo Industries Ltd.  Spectrophotometric method and apparatus for the noninvasive 
US5070874A (en) *  19900130  19911210  Biocontrol Technology, Inc.  Noninvasive determination of glucose concentration in body of patients 
US5086229A (en) *  19890119  19920204  Futrex, Inc.  Noninvasive measurement of blood glucose 
US5218553A (en) *  19890823  19930608  The PerkinElmer Corporation  Method for determining concentrations by means of atomic emission spectroscopy 
US5298022A (en) *  19890529  19940329  Amplifon Spa  Wearable artificial pancreas 
US5308982A (en) *  19911004  19940503  PerkinElmer Corporation  Method and apparatus for comparing spectra 
US5494032A (en) *  19910712  19960227  Sandia Corporation  Oximeter for reliable clinical determination of blood oxygen saturation in a fetus 
US5497772A (en) *  19931119  19960312  Alfred E. Mann Foundation For Scientific Research  Glucose monitoring system 
US5685299A (en) *  19910307  19971111  Masimo Corporation  Signal processing apparatus 
US5792047A (en) *  19970115  19980811  Coggins; George  Physiological parameter monitoring and biofeedback apparatus 
US5822715A (en) *  19970110  19981013  Health Hero Network  Diabetes management system and method for controlling blood glucose 
US5840020A (en) *  19960212  19981124  Nokia Mobile Phones, Ltd.  Monitoring method and a monitoring equipment 
US5995860A (en) *  19950706  19991130  Thomas Jefferson University  Implantable sensor and system for measurement and control of blood constituent levels 
US6023009A (en) *  19960223  20000208  Circe Biomedical, Inc.  Artificial pancreas 
US6036642A (en) *  19910307  20000314  Masimo Corporation  Signal processing apparatus and method 
US6061582A (en) *  19940214  20000509  University Of Iowa Research Foundation  Method and apparatus for noninvasive determination of physiological chemicals, particularly glucose 
US6081735A (en) *  19910307  20000627  Masimo Corporation  Signal processing apparatus 
US6083172A (en) *  19950807  20000704  Nellcor Puritan Bennett Incorporated  Method and apparatus for estimating physiological parameters using modelbased adaptive filtering 
US6152876A (en) *  19970418  20001128  Rio Grande Medical Technologies, Inc.  Method for noninvasive blood analyte measurement with improved optical interface 
US6272480B1 (en) *  19971017  20010807  Siemens Aktiengesellschaft  Method and arrangement for the neural modelling of a dynamic system with nonlinear stochastic behavior 
US6368272B1 (en) *  19980410  20020409  Proactive Metabolics Company  Equipment and method for contemporaneous decision supporting metabolic control 
US6413223B1 (en) *  19990601  20020702  Massachussetts Institute Of Technology  Cuffless continuous blood pressure monitor 
US6544212B2 (en) *  20010731  20030408  Roche Diagnostics Corporation  Diabetes management system 
US6558351B1 (en) *  19990603  20030506  Medtronic Minimed, Inc.  Closed loop system for controlling insulin infusion 
US6572545B2 (en) *  20000922  20030603  Knobbe, Lim & Buckingham  Method and apparatus for realtime control of physiological parameters 
Family Cites Families (6)
Publication number  Priority date  Publication date  Assignee  Title 

US5387327A (en) *  19921019  19950207  Duquesne University Of The Holy Ghost  Implantable nonenzymatic electrochemical glucose sensor 
US5687733A (en) *  19951026  19971118  Baxter International Inc.  System and method for estimating cardiac output 
DE69725621D1 (en)  19970321  20031120  Nellcor Puritan Bennett Inc  Method and apparatus for adaptively averaging data signals 
US6216690B1 (en) *  19971015  20010417  DatexOhmeda, Inc.  Method and apparatus for rapid control of set inspired gas concentration in anesthesia delivery systems 
US6175752B1 (en) *  19980430  20010116  Therasense, Inc.  Analyte monitoring device and methods of use 
EP1112021A1 (en) *  19980910  20010704  Spectrx, Inc.  Attribute compensation for analyte detection and/or continuous monitoring 
Patent Citations (27)
Publication number  Priority date  Publication date  Assignee  Title 

US4655225A (en) *  19850418  19870407  Kurabo Industries Ltd.  Spectrophotometric method and apparatus for the noninvasive 
US5086229A (en) *  19890119  19920204  Futrex, Inc.  Noninvasive measurement of blood glucose 
US5298022A (en) *  19890529  19940329  Amplifon Spa  Wearable artificial pancreas 
US5218553A (en) *  19890823  19930608  The PerkinElmer Corporation  Method for determining concentrations by means of atomic emission spectroscopy 
US5070874A (en) *  19900130  19911210  Biocontrol Technology, Inc.  Noninvasive determination of glucose concentration in body of patients 
US6036642A (en) *  19910307  20000314  Masimo Corporation  Signal processing apparatus and method 
US6263222B1 (en) *  19910307  20010717  Masimo Corporation  Signal processing apparatus 
US6236872B1 (en) *  19910307  20010522  Masimo Corporation  Signal processing apparatus 
US5685299A (en) *  19910307  19971111  Masimo Corporation  Signal processing apparatus 
US6081735A (en) *  19910307  20000627  Masimo Corporation  Signal processing apparatus 
US5494032A (en) *  19910712  19960227  Sandia Corporation  Oximeter for reliable clinical determination of blood oxygen saturation in a fetus 
US5308982A (en) *  19911004  19940503  PerkinElmer Corporation  Method and apparatus for comparing spectra 
US5497772A (en) *  19931119  19960312  Alfred E. Mann Foundation For Scientific Research  Glucose monitoring system 
US6061582A (en) *  19940214  20000509  University Of Iowa Research Foundation  Method and apparatus for noninvasive determination of physiological chemicals, particularly glucose 
US5995860A (en) *  19950706  19991130  Thomas Jefferson University  Implantable sensor and system for measurement and control of blood constituent levels 
US6083172A (en) *  19950807  20000704  Nellcor Puritan Bennett Incorporated  Method and apparatus for estimating physiological parameters using modelbased adaptive filtering 
US5840020A (en) *  19960212  19981124  Nokia Mobile Phones, Ltd.  Monitoring method and a monitoring equipment 
US6023009A (en) *  19960223  20000208  Circe Biomedical, Inc.  Artificial pancreas 
US5822715A (en) *  19970110  19981013  Health Hero Network  Diabetes management system and method for controlling blood glucose 
US5792047A (en) *  19970115  19980811  Coggins; George  Physiological parameter monitoring and biofeedback apparatus 
US6152876A (en) *  19970418  20001128  Rio Grande Medical Technologies, Inc.  Method for noninvasive blood analyte measurement with improved optical interface 
US6272480B1 (en) *  19971017  20010807  Siemens Aktiengesellschaft  Method and arrangement for the neural modelling of a dynamic system with nonlinear stochastic behavior 
US6368272B1 (en) *  19980410  20020409  Proactive Metabolics Company  Equipment and method for contemporaneous decision supporting metabolic control 
US6413223B1 (en) *  19990601  20020702  Massachussetts Institute Of Technology  Cuffless continuous blood pressure monitor 
US6558351B1 (en) *  19990603  20030506  Medtronic Minimed, Inc.  Closed loop system for controlling insulin infusion 
US6572545B2 (en) *  20000922  20030603  Knobbe, Lim & Buckingham  Method and apparatus for realtime control of physiological parameters 
US6544212B2 (en) *  20010731  20030408  Roche Diagnostics Corporation  Diabetes management system 
Cited By (18)
Publication number  Priority date  Publication date  Assignee  Title 

WO2005081171A2 (en) *  20040219  20050901  Edward Henry Mathews  Insulin bolus calculator for mobile communication device 
WO2005081171A3 (en) *  20040219  20060914  Edward Henry Mathews  Insulin bolus calculator for mobile communication device 
WO2008057384A2 (en) *  20061102  20080515  The University Of North Carolina At Chapel Hill  Methods and systems for determining an intravenous insulin infusion rate 
WO2008057384A3 (en) *  20061102  20080904  Susan S Braithwaite  Methods and systems for determining an intravenous insulin infusion rate 
US20100137788A1 (en) *  20061102  20100603  Braithwaite Susan S  Methods and systems for determining an intravenous insulin infusion rate to correct hyperglycemia of a patient, to maintain euglycemia of a patient, and to prevent hypoglycemia of a patient 
US8721584B2 (en) *  20061102  20140513  The University Of North Carolina At Chapel Hill  Methods and systems for determining an intravenous insulin infusion rate to correct hyperglycemia of a patient, to maintain euglycemia of a patient, and to prevent hypoglycemia of a patient 
US20080214919A1 (en) *  20061226  20080904  Lifescan, Inc.  System and method for implementation of glycemic control protocols 
US20110282321A1 (en) *  20070131  20111117  Medtronic Minimed, Inc.  Model predictive method and system for controlling and supervising insulin infusion 
US20100249561A1 (en) *  20070621  20100930  University Of Virginia Patent Foundation  LQG Artificial Pancreas Control System and Related Method 
US9414782B2 (en)  20071010  20160816  Optiscan Biomedical Corporation  Fluid component analysis systems and methods for glucose monitoring and control 
US8449524B2 (en)  20071010  20130528  Optiscan Biomedical Corporation  Fluid component analysis systems and methods for glucose monitoring and control 
US9302045B2 (en)  20080912  20160405  Optiscan Biomedical Corporation  Fluid component analysis system and method for glucose monitoring and control 
US8417311B2 (en)  20080912  20130409  Optiscan Biomedical Corporation  Fluid component analysis system and method for glucose monitoring and control 
US9317657B2 (en) *  20081126  20160419  University Of Virginia Patent Foundation  Method, system, and computer program product for tracking of blood glucose variability in diabetes 
US20110264378A1 (en) *  20081126  20111027  University Of Virginia Patent Foundation  Method, System, and Computer Program Product For Tracking of Blood Glucose Variability in Diabetes 
US8577635B2 (en) *  20100329  20131105  Commissariat A L'energie Atomique Et Aux Energies Alternatives  Calibrating method and system, recording medium for this method 
US20110238350A1 (en) *  20100329  20110929  Commissariat A L'energie Atomique Et Aux Energies Alternatives  Calibrating method and system, recording medium for this method 
WO2017136218A1 (en) *  20160201  20170810  Dexcom, Inc.  System and method for decision support using lifestyle factors 
Also Published As
Publication number  Publication date  Type 

US20020111547A1 (en)  20020815  application 
US20020099282A1 (en)  20020725  application 
WO2002024065A1 (en)  20020328  application 
US6572545B2 (en)  20030603  grant 
US6575905B2 (en)  20030610  grant 
Similar Documents
Publication  Publication Date  Title 

Heise et al.  Noninvasive blood glucose sensors based on near‐infrared spectroscopy  
R Brown et al.  A Bayesian semiparametric joint hierarchical model for longitudinal and survival data  
Van Buuren et al.  Multiple imputation of missing blood pressure covariates in survival analysis  
US6918874B1 (en)  Attribute compensation for analyte detection and/or continuous monitoring  
US5355880A (en)  Reliable noninvasive measurement of blood gases  
US8998809B2 (en)  Systems and methods for calibrating minimally invasive and noninvasive physiological sensor devices  
US20050197553A1 (en)  Patient management of diabetes treatment  
US7392167B2 (en)  Disease simulation system and method  
Magni et al.  Model predictive control of glucose concentration in type I diabetic patients: An in silico trial  
US6045512A (en)  System and method for continuous estimation and display of cardiac ejection fraction and end diastolic volume  
US6678548B1 (en)  Unified probabilistic framework for predicting and detecting seizure onsets in the brain and multitherapeutic device  
HolsteinRathlou et al.  Renal blood flow regulation and arterial pressure fluctuations: a case study in nonlinear dynamics  
US20070282186A1 (en)  Blood glucose monitor with an integrated data management system  
Gross et al.  Performance evaluation of the MiniMed® continuous glucose monitoring system during patient home use  
Hovorka et al.  Closing the loop: the adicol experience  
US6522912B1 (en)  Ear type thermometer  
Klonoff  A review of continuous glucose monitoring technology  
US5398688A (en)  Method, system and instrument for monitoring food intake  
US20120173200A1 (en)  Method, System and Computer Program Product for RealTime Detection of Sensitivity Decline in Analyte Sensors  
US20090112478A1 (en)  Modified Sensor Calibration Algorithm  
Mastrototaro  The MiniMed continuous glucose monitoring system  
US6882940B2 (en)  Methods and devices for prediction of hypoglycemic events  
US6879971B1 (en)  Automated method for building a model  
US6243696B1 (en)  Automated method for building a model  
US20070173761A1 (en)  Apparatus and method for controlling insulin infusion with state variable feedback 