Control of delayed nonlinear model of type 1 diabetes using an improved sliding model strategy

Document Type : Original Article

Authors

1 West Virginia University, Department of Electrical and Computer, Morgantown, USA

2 Department of Electrical and Computer, Science and Research Branch, Islamic Azad University, Tehran, Iran

10.22034/jbr.2021.295685.1042

Abstract

Diabetes has been identified as the sixth leading cause of death in the world. For this reason, diabetes is one of the new and important topics in control science. The goal of diabetes control and treatment is to keep blood glucose at a normal level in the body. Insulin must be injected into the body to control blood glucose in type 1 diabetes. It is so obvious that this injection must be completely controlled. Because an excessive increase in insulin injection can cause human death. In fact, it is not possible to generalize many common controllers for latency systems to delayed systems, and researchers are studying ways to design these controllers for delayed systems. In general, the number of differential equations is reduced in expressing the system behaviors in a delayed nonlinear form, while the method of controlling them will be much more complex. In this paper, the improved sliding mode control technique is used. The proposed improved sliding mode controller is suitable for pumping and delivering insulin in the closed-loop control of diabetes. The optimal insulin injection rate is obtained as a control signal for reference glucose tracking. Then, simulation results are observed to optimally track the patient's blood glucose in exchange for the optimal insulin injection rate. One of the advantages of the proposed method is high accuracy in the presence of physical disturbance such as food acceptance and uncertainty and its resistance.

Graphical Abstract

Control of delayed nonlinear model of type 1 diabetes using an improved sliding model strategy

Keywords


  1. Introduction

Various mathematical methods can be used in the wide variety of sciences like medical and engineering analysis [1, 2]. These studies embody, but are not limited to molecular studies, signal processing, different types of controllers, image processing and so on [3, 4, 5]. For example, researchers in [6], proposed a novel method for data compression and designed a compressor for wearable ECG recorder based on their CS algorithm. Diabetes is a metabolic disorder in which blood glucose levels rise above the normal range of 110-170 mg/dl in the fasting blood sugar test. In diabetics, beta cells, which are the source of insulin production in the body, are destroyed and the body cannot control blood glucose levels alone [7]. In this case, blood sugar should be regulated using insulin injection. The best way to regulate the states of a system is to design a proper closed controller for the system. Closed loop controllers get a feedback from output of the system and compare it with the required condition to reach to the desired input. There are wide variety of applications for these types of controllers such as intelligent traffic control [8, 9], biomedical applications [10] and etc. The closed loop control method works like an artificial pancreas. In this method, the patient's blood glucose concentration is measured immediately by a glucose sensor and the amount of insulin injected is determined by an control algorithm [11]. Finally, this amount of insulin is injected into the patient continuously using an insulin injection pump. Therefore, closed loop control has been considered as the most accurate control strategy in this paper. Blood glucose sensor, insulin injection pump and proper control algorithm are the main components of closed loop control strategy. Figure 1 shows an outline of a closed loop control system for a diabetic patient using an insulin injection pump. Numerous mathematical models for diabetes based on the interaction of insulin and glucose in the body have
been proposed.

                                   

Figure 1. Closed-loop control system of a diabetic patient using an insulin injection pump

 

Common and valid models for diabetes are DDE and ODE. Early models for diabetes could not model the time delay from when blood glucose levels rise to the release of insulin. In some models, in order to obtain fluctuations in insulin secretion, insulin had to be divided into two components, plasma insulin and cellular Ga, which is one of the disadvantages of the proposed model. In this paper, due to the nonlinear behavior of insulin-glucose interaction in the body of type 1 diabetic patient, the Palumbo delayed nonlinear model is used. One of the most important advantages of this model is compliance with the behavior of the patient population according to the IVGTT glucose intravenous injection test [12,13]. Numerous control methods have been used, including methods based on linear and nonlinear predictive  control, classical controllers with linearized equations around the equilibrium point, nonlinear slip-fuzzy model controller to regulate model glucose without delay in state variables for diet and some of which are mentioned in references [14,15]. Control methods based on nonlinear models are based on more knowledge of the physiological behavior of the patient community and provide the possibility of presenting different control theory methods in regulating blood glucose. Obviously, the closer the model's behavior is to the nature of the patient's body, the more accurate the resulting control law will be. Providing an appropriate solution for controlling delayed nonlinear models of diabetes is still an important issue. Lack of sensitivity to internal and external disturbances, final accuracy and pervasiveness, as well as limited time convergence, which are the main features of SMC with the human body, make it a good choice for control algorithms, where high accuracy is one of the special importance.

 

  1. Modeling of Type 1 Diabetes

Delayed nonlinear models have received more attention in recent years due to their consistency with the results of diagnostic tests on the patient population. The presence of delays in nonlinear relationships reduces the order of the model, but their analysis will be more complex. In 2007, Palumbo et al. developed a model that includes the τ_g delay and the K_xr index as the plasma insulin concentration index [18]. Plasma glucose concentration (in mM) Plasma insulin concentration (in pM) are considered as the intravenous insulin injection rate (in min / pM) and the controlled signal. Equation (1) represents this delayed nonlinear model.

 

           (1)

        

 

According to the presentation model, G_0 (τ), I_0 (τ) are the initial values of the patient’s glucose and insulin plasma, which are considered based on the values of  and . Introduction of palumbo DDE nonlinear model parameters are as follows:

  • is a constant and positive parameter in the pulverization of glucose in plasma
  • Insulin injection rate or control signal
  • Delay in pancreatic insulin secretion in return for increased blood glucose concentration [min]
  • Glucose Distribution Rate
  • Plasma Insulin Distribution Rate (L / KGBW)
  • The amount of glucose-dependent stored glucose produced in terms of [ ]
  • Determine the relationship between hepatic glucose and glucose found [ ]
  • Maximum rate of insulin secretion in the second phase [ ]

The nonlinear function , which represents the rate of insulin delivered, is expressed as Equation (2):

 

                           (2)

 

 Plasma glucose level is in the condition that insulin secretion has reached half of the maximum value [17,19].

 

 

  1. Improved Slip Mode Controller Design to Regulate Blood Glucose

The Bergman minimal model is used to design the object slip mode controller in various articles. Lack of sensitivity to internal and external disturbances, unlimited accuracy, resistance, as well as limited convergence time are the main features of the sliding mode that make SMC a suitable choice for algorithms related to the human body, in which accuracy is very important [16,20]. Traditional SMCs have several inherent problems such as rupture in the control signal, which is used to overcome the problems as well as higher accuracy of the higher order slip modeling method. HOSM puts the k-th slider mode variable at a suitable stable source for the controller design:

  • Design accuracy is improved using this method.
  • This method provides a continuous control.

 

3.1. HOSM controller design

The system introduced by Bergmann's minimal model is written in the form of state equations:

 

                                   (3)

 

Where ,  and  are blood plasma glucose concentration (mg / dl), the effect of insulin on blood glucose loss (1/min) and plasma insulin concentration (ml / Uµ), respectively. Stabilization of glucose concentration in the blood of a diabetic patient at the basic level of  is an issue of output tracking. Therefore, tracking error is defined as the difference between the concentration level and its essential amount in the blood of a diabetic patient as the following equation:

 

                                                          (4)

 

For the system introduced in Equation (2), the controller  must be designed in such a way that the error in the presence of uncertainties, parameter changes and perturbations, the absorption of food  is reduced to zero. At first, the relative degree of the system must be determined. Assuming , the relative degree is defined as the number of consecutive derivatives until the control appears in the equation. Thus, the relative degree r-th means that the controller  first appears in  of the final derivative . Using Equation (3), the control function appears in the equations after three derivations, i.e.:

 

                                                         (5)

 

Where:

                (6)

 

Since , , and , system (3) has a well-defined relative degree . This causes the controller designed for system (3) to satisfy the relation . For design, the slider variable controller is introduced as (7):

 

                                                 (7)

 

Which is obtained by deriving the following relation:

 

                                                     (8)

                         (9)

 

In the design process, it is assumed that Equation (9) is finite, i.e. . HOSM, which stabilizes  at the origin for a limited time, and is considered as the following equation:

 

                                                      (10)                                          

It is clear that with the introduction of virtual control, a derivative is added, which leads to an increase in the relative degree of the system from three to four. To calculate , the HOSM derivative is used. The general form for the order nth derivative of a uniform function  is as follows:

 

        (11)

 

Increasing the degree from three to four leads to a high frequency switch in the virtual control. While the local control u is continuous so that  also the fourth-order quasi-continuous control in Equation (11) can be used instead of the high-order control of Equation (10) [15,21]. The dynamic system is defined as follows:

 

 

 

                                                               (12)

 

R is a relative degree and the output of equation  is defined as follows:

 

                                                          (13)

 

If , then we have:

 

                                   (14)

 

In this case, the value of  from the input function  is changing in intervals as  and the changes  are checked for each moment, and according to the amount of insulin required, it changes the level of control input slip. In fact, in addition to commanding the slip surface, the proposed slip mode control also provides an optimal  based on insulin changes (ISMC) [22,23].

 

  1. Simulation using Improved Sliding Mode

To design a glucose detector system, the variable q(t) is added to the equations of state according to Equation (15).

 

 

           (15)

 

For quasi-linear representation of the above equations, the matrices  and  will be as follows:

 

 

 

                                                                  (16)

 

Assuming the initial conditions of the system are positive, the values  and  are positive for all times. Therefore, the quasi-linear matrix elements will be continuous and non-zero. On the other hand, the problem is justified by the degree of completeness of the controllability and visibility matrices depending on the situation. The parameters of the nonlinear model of diabetes for the studied patients are based on the fitting of the least squares generalized to the experimental data of the intravenous insulin injection test [24,28]. The values of  and  can be measured directly, some parameters such as  and  are fixed and known, and parameters such as ,  and  are estimated for each patient [29,30]. The parameters  and  are also parameters that are determined based on the conditions of physical stability of each patient based on algebraic relations. In the simulated example, the patient has a body mass index of 50 and  indicates a higher than normal level of the patient's blood glucose and the insulin resistance index is . These factors indicate an abnormal insulin secretion rate for a newly diagnosed diabetic patient whose factors such as obesity, inactivity, and genetics have led to a gradual decrease in the patient's insulin secretion rate. This patient also has symptoms of type 2 diabetes if left untreated. The parameters of the patient delayed nonlinear model of Equation (14) are as follows:

 

                                                                 (17)

 

The reference glucose signal as a reduction of the initial value of the patient's blood glucose  to the normal value of is considered as follows:

 

                                    (18)

 

Figures 2, 3, and 4 show the convergence of the diabetic patient's blood glucose to the glucose reference signal, the changes in insulin in the patient's blood plasma, and the rate of insulin injection as the control signal provided in the closed-loop improved control system.

 

Figure 2. Tracing of blood glucose reference with control actions for the patient

 

Figure 3. Patient insulin injection rate (control signal applied) (u) and patient blood insulin with normal control (I)

 

Figure 4. Reference blood glucose tracking by applying ISMC control to the patient

 

Below are the simulation results of the same example with the improved sliding mode control design of ISMC. The simulation results show the ability to properly track optimal glucose for patients based on the optimal insulin injection rate. The amplitude of the control signal or insulin injection rate for the patient is reduced compared to the feedback linearization method (the method used by palumbo), which indicates that the patient reaches the normal state of the body earlier in exchange for a lower injection rate and reduced costs. Also, by comparing Figures 2 and 4, tracking in the improved sliding mode has far better convergence results than in the traditional sliding mode

 

 

Figure 5. Patient insulin injection rate (control signal applied) (u) and Patient blood insulin with ISMC control

 

 

  1. Conclusion

In this paper, a sliding mode design for a delayed  nonlinear model of palumbo diabetes is considered. In this paper, the ISMC method was applied and a continuous and comprehensive controller was designed to maintain blood glucose at the basal level. The example provided by the ode23 function was simulated for delayed problems in MATLAB. In similar studies, only the sliding model controller has been used for the delay model. Simulation results also show that by changing the structure of the sliding mode and optimizing it, delayed biological systems have better controllability. Another advantage of this method is that the controller is resistant to changes in the slip of the object and the control is different for different patients, but eventually the glucose stabilizes to its basic level within a reasonable period of time. The ability to use the proposed method for other similar models and the use of the  control method for the studied model are among the suggestions for future projects.

 

 

Conflict of interest

The authors certify they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.

 

Acknowledgments

No applicable.

[1]    Taghibakhshi, Ali, Maryam Barisam, Mohammad Said Saidi, Navid Kashaninejad, and Nam-Trung Nguyen. "Three-dimensional modeling of avascular tumor growth in both static and dynamic culture platforms." Micromachines 10, no. 9 (2019): 580.  
[2]    Pishkenari, Hossein Nejat, Fartash Samie Yousefi, and Ali Taghibakhshi. "Determination of surface properties and elastic constants of FCC metals: a comparison among different EAM potentials in thin film and bulk scale." Materials Research Express 6, no. 1 (2018): 015020.  
[3]    R. Talebitooti, M. Morovati, “Study on TVD Parameters Sensitivity of a Crankshaft Using Multiple Scale and State Space Method Considering Quadratic and Cubic Non-Linearities”, Latin American Journal of Solids and Structures, Vol. 11, No. 14, pp. 2672-2695, 2014  
[4]    A. T. Adnani, A. Dokami, M. Morovati, “Fault detection in high speed helical gears considering signal processing method in real simulation”, Latin American Journal of Solids and Structures, Vol. 13, No. 11, pp. 2113-2140, 2016.
 [5]    Pishkenari, Hossein Nejat, Mohammad Reza Barzegar, and Ali Taghibakhshi. "Study and Simulation of Nanoparticle Translocation Through Cell Membrane." Iranian Journal of Science and Technology, Transactions of Mechanical Engineering (2019): 1-22.  
[6]    Izadi, Vahi, Pouria Karimi Shahri, and Hamed Ahani. "A compressed-sensing-based compressor for ECG." Biomedical engineering letters (2020): 1-9.  
[7]    Ahren, B., Taborsky Jr., G.J., 2002. " B-cell function and insulin secretion" Baron, A. (Eds.).In: Porte, D., Sherwin, R.S., Ellenberg and Rifkin’s Diabetes Mellitus, sixth ed. McGraw-Hill Professional, New York,pp. 43–65 (Chapter 4).
 [8]    Karimi Shahri, Pouria, Shubhankar Chintamani Shindgikar, Baisravan HomChaudhuri, and Amir H. Ghasemi. "Optimal Lane Management in Heterogeneous Traffic Network." In ASME 2019 Dynamic Systems and Control Conference. American Society of Mechanical Engineers Digital Collection, 2019.
 [9]    hahri, Pouria Karimi, Amir H. Ghasemi, and Vahid Izadi. Optimal Lane Management in Heterogeneous Traffic Network Using Extremum Seeking Approach. No. 2020-01-0086. SAE Technical Paper, 2020.
 [10]    Ahani, Hamed, Milad Familian, and Reza Ashtari. "Optimum Design of a Dynamic Positioning Controller for an Offshore Vessel." Journal of Soft Computing and Decision Support Systems 7, no. 1 (2020): 13-18.
 [11]    J. jaremco, O. Rorstad, "Advances toward the implantable artificial pancreas for treatment of diabetes", diabetes care, vol.21, pp.444-450, 1998.
 [12]    P. Palumbo, P. Pepe, S. Panunzi, A. De Gaetano, “Time-Delay Model Based Control of the Glucose–Insulin System, by Means of a State Observer,” European Journal of Control (2012)6:591–606.  
[13]    P. Palumbo, P. Pepe, S. Panunzi, A. De Gaetano, “Glucose control by subcutaneous insulin administration: a DDE modelling approach,” Preprints of the 18th IFAC World Congress, Milano (Italy) August 28 - September 2, (2011).
 [14]    L. Magni, D. M. Raimondo, C. Dalla Man, “Model predictive control of glucose concentration in type I diabetic patients: An in-silico trial,” Biomedical Signal Processing and Control 4 (2009) 338–346.  
[15]    Taghian Dinani, S.; Zekri, M.; Nazari, B., "Fuzzy high-order sliding mode control of blood glucose concentration," in Computer and Knowledge Engineering (ICCKE), 2013 3th International eConference on , vol., no., pp.40-45, Oct. 31 2013-Nov. 1 2013 doi: 10.1109/ICCKE.2013.6682818  
[16]    Edwards C, Spurgeon S. Sliding Mode Control: Theory and Applications. Taylor & Francis: London, 1998.  
[17]    Levant A. Quasi-continuous high order sliding mode controllers. IEEE Transactions on Automatic Control 2005; 50(11):1812–1816.  
[18]    Izadi V, Ghasemi AH. Modulation of control authority in adaptive haptic shared control paradigms. Mechatronics. 2021 Oct 1;78:102598.  
[19]    D.V. Giang, Y. Lenbury, A. De Gaetano, P. Palumbo, “Delay model of glucose-insulin systems: global stability and oscillated solutions conditional on delays,” J. Mathematical Analysis and Applications, 343 (2008), 996– 1006.  
[20]    S. Panunzi, P. Palumbo and A. De Gaetano, “A discrete single delay model for the intra venous glucose tolerance test, Theoretical Biology and Medical Modelling, (2007).  
[21]    Levant A. ''Higher-order sliding modes, differentiation and output feedback control'', International Journal of Control 2003; 76:924–941.
 [22]    Yousefiankalareh A, Najari A, Hosseynzadeh M. Tree-based Routing Protocol in Wireless Sensor Networks using Optimization Algorithm Batch Particles with a Mobile Sink. In2020 IEEE 17th International Conference on Smart Communities: Improving Quality of Life Using ICT, IoT and AI (HONET) 2020 Dec 14 (pp. 1-5). IEEE.
 [23]    Yousefiankalareh A, Manoochehri MA, Khoshnevis SA, Surakanti SR. Improve the Frequency Identification in SSVEP based BCI Systems with Moving Windows Algorithm. In2019 6th International Conference on Control, Instrumentation and Automation (ICCIA) 2019 Oct 30 (pp. 1-5). IEEE.
 [24]    Levant A.'' Universal SISO sliding mode controllers with finite-time convergence'', IEEE Transactions on Automatic Control 1998; 46(9):1447– 1451.  
[25]    P. Palumbo, P.Pepe, S. Panunzi, “Robust closed-loop loop control of plasma glycemia: a discrete-delay model approach,” Mathematcial biosciences and engineering, (2007).  
[26]    Izadi V, Abedi M, Bolandi H. Supervisory algorithm based on reaction wheel modelling and spectrum analysis for detection and classification of electromechanical faults. IET Science, Measurement & Technology. 2017 Nov 9;11(8):1085-93.
 [27]    Izadi V, Abedi M, Bolandi H. Verification of reaction wheel functional model in HIL test-bed. In2016 4th International Conference on Control, Instrumentation, and Automation (ICCIA) 2016 Jan 27 (pp. 155-160). IEEE.
 [28]    P. Palumbo, S. Panunzi, and A. De Gaetano, “Qualitative behavior of a family of delay differential models of the glucose-insulin system,” Discrete Contin. Dyn. Syst. Ser. B, 7 (2007), 399–424.
 [29]    Surakanti SR, Khoshnevis SA, Ahani H, Izadi V. Efficient Recovery of Structrual Health Monitoring Signal based on Kronecker Compressive Sensing. International Journal of Applied Engineering Research. 2019;14(23):4256-61.  
[30]    Tavakoli M, Ahani H. Removing EOG Artifacts from EEG Signals Using a Modified Wavelet-RLS Method. Journal of Bioengineering Research. 2020 Jun 1;2(2).
 [31]    Rafieipour H, Zadeh AA, Mirzaei M. Distributed Frequent Itemset Mining with Bitwise Method and Using the Gossip-Based Protocol. Journal of Soft Computing and Decision Support Systems. 2020 May 7;7(3):32-9.
 [32]    Rafieipour H, Abdollah Zadeh A, Moradan A, Salekshahrezaee Z. Study of Genes Associated With Parkinson Disease Using Feature Selection. Journal of Bioengineering Research. 2020 Dec 1;2(4):1-2.