Closed loop control algorithm for an artificial pancreas

1. A method to control an insulin pump responsive to a controller which receives data from a glucose sensor, the method comprising performing the following steps automatically via the controller:

receiving a plurality of glucose level measurements from the glucose sensor;

calculating an estimated plasma insulin concentration in a blood stream in a body at a predetermined time interval based on an amount of an insulin delivered by the insulin pump, and a lag time from a time of the insulin being delivered to a time of the insulin entering the blood stream;

constructing an autoregressive model configured to describe a relationship between the estimated plasma insulin concentration at a time t and a difference between two successive glucose level measurements at the times t and t-k, wherein k is equal to the predetermined time interval;

computing a plurality of parameters of the autoregressive model by using a plurality of sets of the estimated plasma insulin concentration and the glucose level measurements at a plurality of times;

calculating an amount of insulin to be delivered at a current time using the autoregressive model and a proportional-integral-derivative (PID) controller respectively, wherein the PID controller calculates the amount of insulin to be delivered based on a difference between the glucose level measurement of the current time and a desirable glucose level and a plurality of predetermined gains;

tuning the plurality of parameters of the autoregressive model and the plurality of predetermined gains of the PlD controller respectively until calculation results of the amount of insulin to be delivered calculated by using the autoregressive model and the amount of insulin to be delivered calculated by using the PID controller are the same, wherein if the calculation results are not the same, the calculation results of the autoregressive model and the PID controller are replaced by average of the calculation results, respectively, to recompute the plurality of parameters of the autoregressive model and the plurality of predetermined gains of the PID controller; and

determining the amount of insulin to be delivered at the current time according to final calculation results from the step of tuning to create a desired future glucose level.

2. The method according to claim 1 , wherein the plurality of times comprises a time A and a time B; the autoregressive model comprises relating a blood glucose level of the time B to a plasma insulin of the time A, wherein the insulin delivered at the time A enters the blood stream from the time B.

3. The method according to claim 1 , wherein tuning the plurality of parameters of the autoregressive model and the plurality of predetermined gains of the MD controller comprises:

comparing the amount of insulin to be delivered calculated by using the autoregressive model to the amount of insulin to be delivered calculated by using the PID controller; and

repeating the steps of tuning and comparing until the difference is eliminated.

4. The method according to claim 1 , wherein automatically performing the steps in claim 1 for each of a plurality of discrete time intervals with updated sensor measurements via the controller.

5. An artificial pancreas using a dosed loop control, comprising:

a) a glucose sensor configured to continually measure respective glucose levels at discrete time intervals and provide respective glucose measurement data;

b) an insulin pump configured to deliver insulin in response to a delivery control signal; and

c) a controller configured to, for each of a plurality of the discrete time intervals:

i) receiving a plurality of glucose level measurements from the glucose sensor;

ii) calculating an estimated plasma insulin concentration in a blood stream in a body at a predetermined time interval based on an amount of an insulin delivered by the insulin pump, and a lag time from a time of the insulin being delivered to a time of the insulin entering the blood stream;

iii) constructing an autoregressive model configured to describe a relationship between the estimated plasma insulin concentration at a time t and a difference between two glucose level successive measurements at the times t and t-k, wherein k is equal to the predetermined time interval;

iv) computing a plurality of parameters of the autoregressive model by using a plurality respective sets of the estimated plasma insulin concentration and the plurality of glucose level measurements at a plurality of times;

v) calculating an amount of insulin to be delivered at a current time using the autoregressive model and a proportional-integral-derivative (PID) controller respectively, wherein the PID controller calculates the amount of insulin to be delivered based on a difference between the glucose level measurement of the current time and a desirable glucose level and a plurality of predetermined gains;

vi) tuning the plurality of parameters of the autoregressive model and the plurality of predetermined gains of the PID controller respectively until calculation results of the amount of insulin to be delivered calculated by using the autoregressive model and the amount of insulin to be delivered calculated by using the PID controller are the same, wherein if the calculation results are not the same, the calculation results of the autoregressive model and the PID controller are replaced by average of the calculation results, respectively, to recompute the plurality of parameters of the autoregressive model and the plurality of predetermined gains of the PID controller; and

vii) determining the amount of insulin to be delivered according to a final calculation result from the tuning step.

6. The artificial pancreas according to claim 5 , wherein the controller is one of a processer in the glucose sensor, a processer in the insulin pump, a processer in a handset, or a processing module of a smart device.