• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Predicting the output error of the suboptimal state estimator to improve the performance of the MPC-based artificial pancreas

    2023-12-01 09:51:34MartinDodekEvaMiklovicov
    Control Theory and Technology 2023年4期

    Martin Dodek·Eva Mikloviˇcová

    Abstract The error of single step-ahead output prediction is the information traditionally used to correct the state estimate while exploiting the new measurement of the system output.However,its dynamics and statistical properties can be further studied and exploited in other ways.It is known that in the case of suboptimal state estimation,this output prediction error forms a correlated sequence,hence it can be effectively predicted in real time.Such a suboptimal scenario is typical in applications where the process noise model is not known or it is uncertain.Therefore, the paper deals with the problems of analytical and empirical modeling, identification, and prediction of the output error of the suboptimal state estimator for the sake of improving the output prediction accuracy and ultimately the performance of the model predictive control.The improvements are validated on an empirical model of type 1 diabetes within an in-silico experiment focused on glycemia prediction and implementation of the MPC-based artificial pancreas.

    Keywords Suboptimal state estimation · Prediction · Autoregressive model · Moving average model · Model predictive control·Diabetes mellitus·Artificial pancreas

    1 Introduction

    In many applications and problems, like the model predictive control of glycemia in diabetic subjects,there is a need to estimate the system state based solely on noisy measurements of the system output.Traditional recursive state estimators/observers,such as the Kalman filter[1],are algorithms based on the prediction and correction of the state estimate by the new output measurement.It means that the state estimate is corrected(innovated)according to the single step-ahead output prediction error produced by the system model.

    It is known from the theory of Kalman filtering,that for an optimal state estimator,the sequence of the output error of the state observer(OESO),which is also called the innovation sequence,has the properties of Gaussian white noise.However,inthecaseofsuboptimalstateestimation,thisinnovation sequence is correlated[2,3]and thus can be effectively predicted in real time.

    However,why should one even consider using the suboptimal state estimation instead of the optimal one?To answer this question,the suboptimal state estimation is basically not an option of choice,but rather an inevitable consequence that arises from the unknown or poorly estimated parameters of the process noise model.The assumption of suboptimality is because the design of Kalman filter and its optimality relies on the exact knowledge of the process noise model,parameters of which are highly uncertain in many applications,or even the entries of the process noise covariance matrix are often used simply as the tuning parameters.It can be claimed that if there is a mismatch between the noise model and the actual statistical properties of the system,the OESO forms a correlated sequence.

    The main practical motivation for studying the dynamics of the OESO is to try to predict it in real time and then correct the predictions of the output variable accordingly,yet the ultimate aim is to involve this prediction in the model predictive control.Therefore,the outlined strategy can be seen as a relatively feasible and cheap way to improve prediction and control performance by effectively compensating for the effect of suboptimal state estimation.

    Concerning the target application domain outlined in the title, performing highly accurate predictions yields better forecasting of severe hyper- and hypo-glycemia states, as these are the major risks linked to diabetes and its treatment [4, 5].Another application of the proposed strategy is possible within an implementation of the model predictive control-based artificial pancreas [6, 7] to control glycemia in subjects with type 1 diabetes by automating the insulin dosing.

    The rationale for choosing this application domain to demonstrate the effectiveness of the proposed strategy is supported by the fact that in most available studies,e.g.in[8–17],the crucial entries of the covariance matrix of the process noise are considered tuning parameters in the Kalman filter design.Therefore,it can be concluded that the state estimator works suboptimally in such scenario.

    It should be mentioned for completeness, that there also exist sophisticated and dedicated methods[18–22]for estimating the noise models, which can potentially eliminate the problem of suboptimality of the state estimation.Unfortunately, these covariance matrix estimation methods have limited practical applicability, as they often provide biased estimates and typically require very large datasets to be supplied,what can be infeasible in many applications.

    In this paper, we provide important insights into the dynamics of the OESO from the analytical point of view,but alsofromthepracticaldata-drivenperspectivewherereduced models are considered.The primary motivation for using the reduced structures,such as the autoregressive and the moving average model, rather than the full analytical model, is their feasible prediction and identifiability from the experimentally obtained OESO sequence.

    The paper has been divided into the following sections:Sect.2 introduces the basic preliminaries,formulation of the stochastic state-space model, and equations describing the conventional state observer.In Sect.3, the full analytical model of the OESO is derived to provide some theoretical background.Section4 considers the reduced autoregressive model and comprises the formulation of the corresponding identification problem and the predictive equation.In Sect.5, the moving average model is studied in a similar way.Section6 covers the idea of using the identified models to enhance the output prediction accuracy and their inclusion in the model predictive control.The setup of the experiment aimed at prediction and model predictive control of glycemia in subjects with type 1 diabetes is outlined in Sect.7,while its results are discussed in Sect.8.

    It should be marginally mentioned that the stochastic nature of glycemia dynamics does not necessarily have to be modeled in the state space by the process noise and the measurement noise, since basic input–output transfer function models like the autoregressive-exogenous, autoregressive moving average with exogenous inputs, and Box–Jenkins model can also be used as reported in[23–26].The stochastic part of these input–output models is usually estimated in one step together with the deterministic submodels as a result of the system identification procedure from the experimental data,so there is no need to worry about the optimality of the state estimation.

    An application of the unconstrained model predictive control with the state-space model along with the state estimation basedontheKalmanfilterwaspresentedinstudy[27].Asimilar control strategy was reported in[28,29],but due to the input–output problem formulation the state estimator was not required.As another example of similar artificial pancreas scheme, in [30] the Kalman filter was used to estimate the state for the linear model predictive control with disturbance rejection.

    In our recent work [31], a novel optimal state estimator was proposed as an alternative to the Kalman filter.However,this algorithm was not based on the traditional recursive correction of the state estimate by the OESO, but it used the generalized least squares formulation of the problem.Also in this case,the optimality of the state estimate depended on the exact knowledge of the process noise model.

    From the perspective of unique contributions of this paper,it is important to remark that in any of the above referenced studies or in the latest comprehensive survey papers[32–34],the strategy of predicting the OESO and compensating for the suboptimality of the state estimation was not proposed or discussed.It can also be concluded that most authors relied on suboptimal state estimation with ad hoc tuning of the process noise model,hence leaving a significant headroom for improving the control performance by targeting this problem.

    Unlike the conventional MPC-based artificial pancreas that typically uses the suboptimal state estimator,which normally results in a correlated OESO sequence,we propose a new strategy to effectively compensate for this effect.The proposed modification involves the prediction of the OESO to directly correct the prediction of the system free response.It is worth noting that this modification is easy to embed in already existing MPC schemes [32–34] while inducing little more computational cost yet requiring no additional hardware modifications.In other words,the proposed OESO models and the prediction/correction strategy can eventually be retrofitted to any advanced MPC-based artificial pancreas design that utilizes the state estimator at the expense of relatively straightforward structural modification.Note that the other features of the artificial pancreas, such as the safety algorithms and constraints,do not directly interact with the proposed modification.

    2 Model structure and preliminaries

    The general discrete-time stochastic state-space empirical model of glycemia dynamics in subjects with type 1 diabetes is postulated as[31]

    wherek∈N is the current sample, the outputy[mmol/l]stands for the deviation of glycemia from its steady-state value,the inputu[U/min]denotes the deviation of the insulin administration rate from the basal rate,andd[g/min]represents the carbohydrate intake rate input.The state vector of thisnth order model is denotedx[n×1],w[n×1] is the process noise vector and the zero-mean uncorrelated random processv~N(0,R)represents the measurement noise of the glucose monitoring device[35].

    The parameters of model (1) include the state-transition matrixA[n×n],the input matrixB[n×2],and the output vectorC[1×n].The state-transition matrixAconsists of the submatricesAu,Adandthezeromatrices0 oftheconforming dimensions as

    where matricesAu[nu×nu] andAd[nd×nd] are in the canonical form and comprise the model coefficientssuch that

    The input matrixBis simply

    where 0 are the zero vectors of conforming dimensions andBu[nu×1],Bd[nd×1]are equal to

    The output vectorCgets

    whereCu[1×nu]andCd[1×nd]will comprise the model coefficients as

    The state vectorxalso holds the canonical form

    wherenuandndare the orders of the corresponding submodels,so the overall model order isn=nu+nd.The state variablesxu,xdrepresent the partial effects of insulin administration and carbohydrate intake,respectively.

    Consider that the process noisewin model(1)represents the effect of input uncertainties,so one can write

    Since all stochastic terms in(1)were defined as zero-mean uncorrelated stationary random processes, the covariance matrixQ[n×n]of the process noise(9)and the varianceRof the measurement noise are equal to

    2.1 State observer

    The state vectorxof model(1)is usually estimated using the recursive state observer

    where ?x[n×1] is the estimated state andK[n×1] is the gain vector,which is the subject of the observer design.The design is usually based on the optimal Kalman approach[1,36]in the case of stochastic system assumption,or using the pole-placement method if the system is deterministic.

    The state estimate residuale[n×1]will be defined as

    Finally,the single step-ahead model output prediction error?,i.e.the OESO,gets

    3 Dynamics of the state observer output error

    The model of the state estimate residualecan be derived by substituting ?x(k+1)from (12) andx(k+1)according to (1a)into(13)while substituting the outputy(k+1)in the terms of 1b as

    Eq.(15) can be transformed by applying the forward timeshift operatorzase(k+1)=ze(k)obtaining

    The above equation can be reshaped to separate vectore(k)as

    whereIis the unit matrix of the conforming dimensions.

    According to(14)and(17),?(k)holds

    Eq.(18) implies that the dynamics of the OESO is represented by a stochastic system with multiple independent noise inputs.One may realize that termC(zI-A+KC)-1results in a row vector of rational functionswith the common denominators(z)as the characteristic polynomial of this system.This consideration yields the transfer function model

    where

    IfQis the covariance matrix of the process noise according to(10),thenw(k)can be written as

    where 0 [n×1]is the zero vector,η[n+1×1]is the vector of uncorrelated noise inputs with the unit variance, i.e.,cov(η,η)=I,

    whereRis the variance of the measurement noise according to(11).

    Finally,model(19)can be generalized as the sum ofn+1 autoregressive-moving-average (ARMA)models by substitutingw(k)in the terms of(21)andv(k)from(23)as

    Since the process noise(9)has only two nonzero components and diagonal covariance matrix(10),general model(24)can be reduced to

    However, model (24) or (25) cannot be used to predict the OESO,since the random input vectorηas well as the partial outputsare unmeasurable in practice.

    Another important paradox concerning this analytical model is that since the covariance matrix of the process noise is considered unknown in the case of suboptimal state estimation,analytical model(24)simply cannot be determined.On the contrary,if the covariance matrix of the process noise is known,what implies that the state estimator works optimally,then model (24) will be known, but it will be unnecessary since the OESO is uncorrelated and hence it cannot be predicted.

    Concerning the identification of full model (24) directly from experimental data, i.e.based on the OESO sequence,this would be hardly possible primarily due to its structure and the large number of parameters to be estimated.

    The aforementioned issues with predictability and identifiability are the main motivation for further considering two reduced model structures,particularly the autoregressive and the moving average model.

    3.1 Optimality test

    To test whether the state estimator works optimally or not,the sample autocorrelation function of the OESO sequence has to be analyzed.In the case of optimal state estimation,this autocorrelation function should show a character similar to that of Dirac delta function.

    Supposing a finite-length experiment withNsamples,the autocorrelation functionR??(nTs)is estimated as[38,39]

    wheren∈Z satisfiesn

    4 Autoregressive model

    In this section,the dynamics of the OESO will be approximated by the single-input single-output autoregressive model defined as

    The polynomialq(z)of thisnqth order model gets

    The equivalent difference equation of model(27)holds

    The parameter vectorqgets

    4.1 Identification strategy

    According to difference equation(29),the corresponding linear regression system consideringNavailable samples takes

    or using the shorthand notation

    whereqis the parameter vector (30),is the regression matrix and?[N×1],η[N×1]are vectors.

    The parameter vectorqcan be estimated as ?qin a straightforward way using the least squares method with the optimal parameter estimate determined analytically as[40]

    4.2 Predictive form

    For model(27),the explicit prediction formula can be derived based on difference equation (29).The future values of the whitenoiseinputareobviouslyunknown,soassumingthatits statistically unbiased prediction is zero,i.e.E■η(k+i)■=0,thepredictiveformconsideringthepredictionhorizonnpgets

    Sincethenoiseinput in(27)is unmeasurable,byreshaping equation(29),η(k)can be estimated as

    5 Moving average model

    In this section,the dynamics of the OESO will be approximated by the moving average model

    The polynomialg(z)in thength order model(40)gets

    The difference equation for model(40)can be written as

    5.1 Identification strategy

    It is well known that estimating the moving average processes is more difficult than estimating the autoregressive processes[41].Since the input noiseηis unmeasurable in practice,the straightforward approach based on the least squares minimization of the model single step-ahead prediction error cannot be directly applied in this case.

    Therefore, to estimate the coefficient vector (43) using the available OESO sequence,the two-step method of Durbin[41,42]will be adopted.The first step of this method consists of fitting an autoregressive model to the OESO sequence via the ordinary least squares method in the terms of Sect.4.In the second step,the identified autoregressive model is used to estimate the input noiseηby filtering the OESO sequence by the inverse of the estimated autoregressive model according to(39).

    The second step uses this estimated input noise sequence ?ηto create the regression system and to estimate the parameters of the moving average process in the least squares sense.

    The corresponding regression system then gets

    or using the shorthand notation,

    wheregis the parameter vector(43),is the regression matrix and?[N×1],η[N×1]are vectors.The optimal parameter vectorgcan be estimated as

    5.2 Predictive form

    Having the model parameters estimated, the OESO (36)can be predicted.Assuming that the statistically unbiased prediction of the input zero-mean white noise is zero, the predictive form of the moving average model (40) can be derived according to the difference equation(42)as

    In practice,the input noiseηcannot be measured,so it has to be estimated based on the inverse filtering of?according to the difference equation(42)as

    6 Prediction and model predictive control with the OESO compensation

    In this section,the algorithm of model predictive control will be adopted from[31],while an important modification that concerns the prediction of the OESO will be proposed.

    Prediction of the state vectorxand the outputycan be expressed considering(1)as

    wherek∈N is the current sample andi∈N getsi=1···np,while assumingnp∈N is the prediction horizon.

    Notice that in(52),the OESO ??,which was predicted by the identified autoregressive or the moving average model,was taken into account by correcting the output prediction ?y.This is the most important modification of the traditional prediction and predictive control algorithms as it allows us to effectively compensate for the suboptimality of the state estimation.

    The predictive control minimizes the quadratic cost function of the model-based predictions of chosen system variables over the prediction horizonnp.The corresponding quadratic form gets[43,44]

    whereΔu f[nc×1] is the vector of future changes of the manipulated variable,while assumingncis the control horizon.MatrixA[nc×nc], vectorb[nc×1] and scalarcare defined as

    The optimization problem(53)can be solved by quadratic programming if linear inequalities constraints are considered.For the sake of simplicity,the elements of the reference vectoryrare equal to an appropriately chosen constantGtrepresenting the target glycemia.Pursuing the receding horizon strategy, only the first element of the optimal solution

    Δu fis actually applied,so one can write

    Note that the manipulated variable has to be constrained,so the minimal insulin infusion rateumin=0 U/min, whileumaxwill be adopted from [27].The corresponding linear inequalities system with respect to the decision vectorΔu fcan be formed by involving matrixΨ(57)as in[45]

    Compared to[31],the crucial modification of the control algorithm is made here by adding the prediction of the OESO ??(k+i)to equation(52).

    Concerning the safety features of automated insulin therapy, various additional strategies could be considered to enhance the current configuration.To avoid the adverse and dangerous insulin stacking phenomenon,the insulin on board[46–48]representing the amount of insulin still active from the previously administered doses can be involved.The dynamics of insulin on board can be represented by simple linear models like in[49],[50]or[51],while this signal can be used to form a special quadratic penalty that is added to the cost function(53)of the MPC.Alternatively,hard constraints for the insulin on board signal can also be assumed by extending the linear inequalities system(59)of the quadratic program accordingly.Another type of safety feature is to hard constrain the controlled variable by modifying the inequalities system (59) to prevent the risk of hypoglycemia and hyperglycemia as proposed in [26].As this strategy can potentially induce control infeasibility, soft formulation of the controlled variable constraints should be preferred, as proposed in[52,53].However,further details are beyond the scope of this paper.For information on safety features,see,e.g.[46,48,54].Also note that these safety features do not directly interact with the proposed strategy of predicting the output error of the suboptimal state estimator, which is the main contribution of the paper.

    7 Experimental setup

    To validate the proposed strategy and assess its practical effectiveness in application to the problem of prediction and predictive control of glycemia in subjects with type 1 diabetes, a simulation-based experiment was designed and evaluated.

    The glycemia response for this experiment was obtained by in-silico approach, simulating the complex physiologybased nonlinear model that was described in[55,56]and the references therein.The basal state of this model was determined with respect to the basal glycemiaGb=6 mmol/l and the corresponding basal insulin administration ratevb=0.01 U/min.

    The orders of empirical model (1) were chosen asnu=nd=4,implying the overall ordern=8.Note that the theory related to estimation of the model parametersin(3)andin(7),is not within the scope of this paper,so we suppose that model(1)was identified with parameters(60).However,for more details on this topic, we refer an interested reader to our recent works[25,57].

    Thepredictionhorizonandthecontrolhorizonwereassumed asnp=15,nc=10, while the sample time was chosen asTs=10 min.

    The variance of the measurement noise 11 and the variances of the process noise 10 were empirically adjusted as

    Note that since the variances of the process noise were just empirically tuned to obtain acceptable performance of the Kalman filter while the actual values cannot be determined because they are not based on any particular physiological mechanisms or characteristics of a diabetic subject,the state estimator will perform only suboptimally in this case, and hence the OESO sequence will be correlated.

    The observer gain vectorKwas calculated according to the Kalman filter design[1,36]while considering the model parameters(60)and the noise model parameters(61)as

    Concerning the initial tuning of the proposed empirical models of the OESO,the order of autoregressive model(27)was set asnq= 4, whereas the order of moving average model(40)was chosen asng=12.

    Fig.1 Evolution of the OESO ?(k) acquired during the experiment

    The experiments were designed to mimic the insulin treatment of a subject with type 1 diabetes during the two-day period.The first investigated problem concerns the prediction of glycemia during standard insulin bolus therapy that wascarriedoutaccordingtotheboluscalculatorrule(see[58]and the references therein).The second deals with automated insulin dosing managed by the model predictive control algorithm of the artificial pancreas.Both algorithms were modified in the terms of Sect.6.

    8 Discussion

    In this section,the results of the outlined simulation experiment will be comprehensively analyzed and discussed.

    The sequence of the OESO obtained in the terms of equation(14)during regular insulin treatment while simultaneously performing the state estimation according to(12)by considering the observer gain (62) is plotted in Fig.1.This figure suggests that although the state estimate asymptotically converges to the actual state,the character of the OESO sequence is far from ideal uncorrelated noise.

    To prove this, the autocorrelation function ?R??of the OESO was estimated according to (26) by processing the sequence from Fig.1 and is plotted in Fig.2.Analyzing this autocorrelation function, one can conclude that the OESO sequence is correlated,what confirms that the state estimator works suboptimally due to the empirically adjusted variances of the process noise(61).

    Performing the estimation of both reduced models of the OESO by pursuing the strategies presented in Sections 4 and 5,the following coefficients were estimated.

    Parameter vectorqof the autoregressive model(27):

    Parameter vectorgof the moving average model(40):

    Fig.2 Estimated autocorrelation function ?R??(nTs) of the OESO sequence from Fig.1

    Fig.3 Estimate of the autocorrelation function ?R?η ?η(nTs) of the estimated noise input for the autoregressive model

    To validate the models by filtering the correlated OESO sequence?by the inverse of each of the identified empirical models, i.e.byq(z) and, respectively, the input noise sequence was estimated according to (39) for the autoregressive model and according to (50) for the moving average model.The autocorrelation functions ?R?η?η(nTs) of these sequences for both model structures are depicted in Figs.3 and 4, which show their Dirac delta function-like nature,proving the estimated models valid.

    Now follows the prediction of the OESO using the predictive form(34)of the autoregressive model and the predictive form (47) of the moving average model respectively.Relatively accurate predictions for randomly chosen starting points can be observed in Fig.5.Both model structures showed almost identical performances and could predict the future evolution of the correlated OESO with a satisfying accuracy considering the highly stochastic nature of this signal.

    Fig.4 Estimate of the autocorrelation function ?R?η ?η(nTs) of the estimated noise input for the moving average model

    Fig.5 Prediction of the OESO using the both autoregressive model)and the moving average model

    The next comparison concerns the practical impact of correcting the glycemia prediction by predicted OESO as the original contribution of the paper.In Fig.6,one can see the uncorrected prediction of glycemia(?G)as the conventional strategy,as well as the predictions that involved the corrections by OESO predicted using the autoregressive(?GAR)and the moving average(?GMA)model.By a basic visual assessment, one can observe an improvement of the prediction accuracy, while the improvements concerned primarily the peaks of the response.Keep in mind that such differences between the uncorrected and the corrected prediction can be critical in situations such as decision making with regard to the application of insulin therapy.

    In addition to the graphical assessment,the prediction performance will be quantified by the quadratic metric

    which will provide a better assessment of the prediction performance.

    The last part of the experiment is focused on the model predictive control of glycemia in the context of the artificial pancreas implementation,where a positive effect of the proposed predictors on control performance is anticipated.To demonstrate this,Fig.7 shows the closed loop glycemia response,where involving the predictions of the OESO visibly improved the control performance in terms of tighter control with respect to the reference value, and reduced maximal and minimal observed glycemia,what is especially significant to reduce the risk of hyperglycemia and hypoglycemia.It can be concluded that both predictors performed almost identically, but way better than in the original case without compensating for the OESO.

    Fig.6 Prediction of glycemia without the OESO compensation G(t)compared to using the autoregressive model GAR(t) and the moving average model GMA(t)

    Fig.7 Predictive control of glycemia without the OESO compensation G(t)compared to using the autoregressive model GAR(t)and the moving average model GMA(t)

    Keep in mind that since typical values of the OESO are relatively low compared to the magnitude of the controlled variable(see Fig.1),the proposed strategy can naturally yield a limited effect.It can also be claimed that the strength of the desired effect is directly related to the performance level of the state observer and thus to the degree of mismatch between the process noise model and the actual statistical properties of the system.Therefore,for systems with an empirically tuned covariance matrix of the process noise,the strategy proposed in this paper is highly recommended.

    The control performance will be quantified by the maximalGmaxand the minimalGminobserved glycemia,and by the quadratic metric

    Table 1 Comparison of the prediction and control performance metrics

    Table 2 Comparison of the prediction and control performance metrics for the autoregressive model

    Table 3 Comparison of the prediction and control performance metrics for the moving average model

    whereGtis the target glycemia.

    The summary of prediction and control performance metrics obtained during the experiment is documented in Table 1,which also confirm the observations from Figs.6 and 7.

    Toinvestigatetheeffectofthechoiceofthetunableparameters of the OESO models, particularly the order of the autoregressive model (27)nqand the order of the moving average model (40)ngon the resulting performance of the proposed strategy, the experimentation was repeated under various configurations, yielding the results summarized in Tables 2 and 3.

    It can be concluded that the moving average model of the OESO performs slightly better in both prediction and predictive control than the autoregressive model,while the choice of the corresponding model ordersnq,ngalso affected the overall performance,yet not consistently for all the metrics considered.However, all studied models performed better than the original uncompensated configuration.

    9 Conclusions

    This study stressed that the OESO sequence in the case of suboptimal state estimation is correlated, while it can be effectively predicted by the autoregressive or the moving average models.These two reduced models could provide good predictability and identifiability from the experimental data.The predicted OESO was then involved to correct the output variable prediction and hence ultimately improve the performance of the model predictive control.It can be concluded that the presented strategy allowed to effectively compensate for the suboptimality of the state estimation caused by the inaccuracy of the process noise model in a relatively inexpensive and feasible way.

    We also obtained theoretical results demonstrating that the actual dynamics of the OESO is analytically described as the sum of ARMA processes,while this full structure was approximated by the autoregressive and the moving average models in practice.

    A promising application of our results would be possible within an implementation of the artificial pancreas in subjects with type 1 diabetes,where maximizing the performance of the model predictive control is of highest priority.The presented simulation-based experiment demonstrated that the dynamics of the OESO can be effectively predicted by both proposed reduced models,while there were documented positive effects on the accuracy of the glycemia prediction and on the performance of the predictive control.Keep in mind that although qualitative improvements do not appear to be significant at first sight,any feasible improvement matters a lot when managing glycemia and can have significant longterm consequences for patient health.

    Compared with the conservative structure of the MPCbased artificial pancreas[6,32–34]that utilizes the Kalman filter state estimation with typically empirically tuned covariance matrix of the process noise,which normally results in its suboptimal performance,the strategy proposed in this paper additionally involves an easily identifiable prediction model of the OESO to effectively compensate for the adverse effect of suboptimality of the state estimation by correcting the system free response prediction.Moreover, this structural modification is not computationally demanding and it can be easily embedded to the existing modular MPC schemes[6, 32–34] without involving any hardware adjustments or directly interacting with other features of the artificial pancreas such as constraints.

    On the contrary, compared with the alternative solution based on methodological estimation of the covariance matrix of the process noise according to the methods presented in[18–22]to implicitly ensure the optimality of the state estimation,it can be claimed that these methods typically require very large experimental datasets(tens of thousands of samples)to provide reliable and unbiased estimates,whereas the statistical models of the OESO proposed in this paper can be identified from relatively limited experimental data(hundreds of samples).

    In a nutshell,the most significant contributions presented in this work include the derivation of the analytical stochastic ARMAmodeloftheOESOanditsapproximationswhichcan be used to improve the performance of the suboptimal state observer in applications when the process noise model is not exactly known or is uncertain.It can be concluded that the actual effect of the proposed strategy depends primarily on the degree of plant-model mismatch of the noise model.

    Funding Open access funding provided by The Ministry of Education,Science,Research and Sport of the Slovak Republic in cooperation with Centre for Scientific and Technical Information of the Slovak Republic Open access funding provided by The Ministry of Education,Science,Research and Sport of the Slovak Republic in cooperation with Centre for Scientific and Technical Information of the Slovak Republic.

    Declarations

    Conflict of interest The authors have no competing interests to declare that are relevant to the content of this article.

    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License,which permits use,sharing,adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitteduse,youwillneedtoobtainpermissiondirectlyfromthecopyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.

    国产精品久久久久久精品古装| 在线观看免费高清a一片| 国产精品一区二区在线不卡| 亚洲一区二区三区不卡视频| 少妇 在线观看| 免费黄频网站在线观看国产| 亚洲av电影在线进入| 母亲3免费完整高清在线观看| bbb黄色大片| 亚洲在线自拍视频| 亚洲熟女毛片儿| 亚洲片人在线观看| 高清黄色对白视频在线免费看| 99热网站在线观看| 啦啦啦在线免费观看视频4| 精品少妇一区二区三区视频日本电影| 搡老岳熟女国产| 国产欧美亚洲国产| 国产国语露脸激情在线看| 91老司机精品| 悠悠久久av| 国产一卡二卡三卡精品| 久久九九热精品免费| 免费av中文字幕在线| 9191精品国产免费久久| 天天躁狠狠躁夜夜躁狠狠躁| 满18在线观看网站| 久久草成人影院| 搡老熟女国产l中国老女人| 飞空精品影院首页| 成人亚洲精品一区在线观看| 亚洲免费av在线视频| 窝窝影院91人妻| 妹子高潮喷水视频| 黄色视频不卡| 中文欧美无线码| 国产人伦9x9x在线观看| 亚洲精品国产区一区二| 捣出白浆h1v1| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3 | 99久久人妻综合| 久久久水蜜桃国产精品网| 两性夫妻黄色片| 久久午夜综合久久蜜桃| 精品一区二区三区视频在线观看免费 | 免费一级毛片在线播放高清视频 | av免费在线观看网站| 国产真人三级小视频在线观看| 丝袜美腿诱惑在线| 18禁裸乳无遮挡动漫免费视频| 99久久99久久久精品蜜桃| 韩国精品一区二区三区| 99国产综合亚洲精品| 热99re8久久精品国产| 日日夜夜操网爽| 午夜精品国产一区二区电影| 免费高清在线观看日韩| 久久久久视频综合| www日本在线高清视频| 免费看a级黄色片| 老司机影院毛片| 国产成人一区二区三区免费视频网站| 亚洲va日本ⅴa欧美va伊人久久| 丁香欧美五月| 久久天躁狠狠躁夜夜2o2o| 国产精品欧美亚洲77777| 欧美日韩av久久| 在线观看免费视频网站a站| 99国产精品一区二区蜜桃av | 深夜精品福利| 操出白浆在线播放| 啦啦啦视频在线资源免费观看| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| av网站在线播放免费| 亚洲国产精品合色在线| 亚洲精品久久成人aⅴ小说| 国产欧美日韩一区二区精品| 精品少妇久久久久久888优播| 热re99久久精品国产66热6| 欧美激情 高清一区二区三区| 国产黄色免费在线视频| 国产精品欧美亚洲77777| 精品少妇一区二区三区视频日本电影| 高清视频免费观看一区二区| 欧美日韩一级在线毛片| 欧美人与性动交α欧美软件| av视频免费观看在线观看| 国产精品久久视频播放| 亚洲精品久久午夜乱码| 久久国产乱子伦精品免费另类| 国产精品二区激情视频| 视频在线观看一区二区三区| 久久国产精品人妻蜜桃| 欧美成人午夜精品| 日本黄色日本黄色录像| 无人区码免费观看不卡| netflix在线观看网站| 999精品在线视频| 精品一区二区三卡| a级片在线免费高清观看视频| 欧美午夜高清在线| 午夜福利一区二区在线看| 亚洲人成77777在线视频| 在线看a的网站| tocl精华| av网站在线播放免费| 亚洲 欧美一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩国产mv在线观看视频| 咕卡用的链子| 看黄色毛片网站| 欧美精品av麻豆av| 免费av中文字幕在线| 脱女人内裤的视频| 久久精品国产综合久久久| 欧美日韩精品网址| 一级a爱视频在线免费观看| 午夜福利视频在线观看免费| 一本一本久久a久久精品综合妖精| 午夜福利乱码中文字幕| 国产高清激情床上av| 人妻 亚洲 视频| 久热这里只有精品99| 国产主播在线观看一区二区| 国产三级黄色录像| 亚洲性夜色夜夜综合| 国产一区二区三区在线臀色熟女 | 国产精品影院久久| 99国产极品粉嫩在线观看| 9191精品国产免费久久| 国产亚洲精品一区二区www | 国产精品九九99| 黄片小视频在线播放| 国产免费av片在线观看野外av| av欧美777| 精品一区二区三卡| 国产av一区二区精品久久| 久久精品熟女亚洲av麻豆精品| 国产精品香港三级国产av潘金莲| 国产又色又爽无遮挡免费看| 国产精品av久久久久免费| 亚洲七黄色美女视频| 欧美成狂野欧美在线观看| 国产成人啪精品午夜网站| 老司机亚洲免费影院| 午夜免费成人在线视频| videosex国产| 久99久视频精品免费| 久久国产亚洲av麻豆专区| 18禁观看日本| 亚洲精品一二三| 欧美激情 高清一区二区三区| 高清视频免费观看一区二区| 免费看十八禁软件| 国产精品 欧美亚洲| 热99国产精品久久久久久7| 婷婷成人精品国产| 亚洲精品美女久久久久99蜜臀| 视频在线观看一区二区三区| 大码成人一级视频| 国产在视频线精品| 极品少妇高潮喷水抽搐| 久久精品国产a三级三级三级| 国产野战对白在线观看| 国产精品久久久av美女十八| 久久精品熟女亚洲av麻豆精品| 亚洲视频免费观看视频| 中文字幕人妻丝袜一区二区| 日本黄色日本黄色录像| 一个人免费在线观看的高清视频| 91av网站免费观看| 欧美精品人与动牲交sv欧美| 久久精品成人免费网站| 欧美人与性动交α欧美精品济南到| 三级毛片av免费| 亚洲一区高清亚洲精品| www.精华液| 黄片播放在线免费| 国产成人精品久久二区二区免费| 好看av亚洲va欧美ⅴa在| 免费日韩欧美在线观看| 97人妻天天添夜夜摸| 在线观看免费视频日本深夜| 亚洲欧美日韩另类电影网站| 每晚都被弄得嗷嗷叫到高潮| 免费黄频网站在线观看国产| 成在线人永久免费视频| 久久精品国产清高在天天线| 国产一区二区三区综合在线观看| 大香蕉久久网| 欧美激情高清一区二区三区| 国产成人精品久久二区二区91| 免费观看a级毛片全部| 中出人妻视频一区二区| 夜夜夜夜夜久久久久| 国产极品粉嫩免费观看在线| 两性夫妻黄色片| 亚洲国产欧美一区二区综合| 国产区一区二久久| www.精华液| 欧美日韩国产mv在线观看视频| 女同久久另类99精品国产91| 国产高清视频在线播放一区| 啦啦啦 在线观看视频| 久久久水蜜桃国产精品网| 亚洲av日韩精品久久久久久密| 成在线人永久免费视频| 男女之事视频高清在线观看| 母亲3免费完整高清在线观看| 亚洲精品自拍成人| 男人的好看免费观看在线视频 | 午夜91福利影院| 成人永久免费在线观看视频| 岛国在线观看网站| 国产精品一区二区精品视频观看| 一边摸一边做爽爽视频免费| 久久婷婷成人综合色麻豆| 午夜精品国产一区二区电影| 久久精品人人爽人人爽视色| 午夜亚洲福利在线播放| 中文字幕另类日韩欧美亚洲嫩草| 午夜亚洲福利在线播放| av超薄肉色丝袜交足视频| 久久精品亚洲av国产电影网| 国产熟女午夜一区二区三区| 亚洲人成伊人成综合网2020| 最新在线观看一区二区三区| 欧美成人免费av一区二区三区 | 成年人午夜在线观看视频| 国产精品免费视频内射| 国产亚洲一区二区精品| 国产深夜福利视频在线观看| 欧美在线黄色| 18禁国产床啪视频网站| 日本黄色日本黄色录像| 国产精品av久久久久免费| 啪啪无遮挡十八禁网站| 如日韩欧美国产精品一区二区三区| 久久久精品国产亚洲av高清涩受| 国产精品一区二区免费欧美| 国产精品一区二区在线不卡| 国产激情欧美一区二区| 少妇 在线观看| 国产91精品成人一区二区三区| 久久精品亚洲av国产电影网| 看免费av毛片| 日韩免费av在线播放| 宅男免费午夜| 免费看十八禁软件| 女人被狂操c到高潮| 国产精品免费大片| 久久热在线av| 咕卡用的链子| 超碰97精品在线观看| 精品一品国产午夜福利视频| 18禁国产床啪视频网站| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区91| a在线观看视频网站| 精品一区二区三区四区五区乱码| 欧美中文综合在线视频| 一级毛片精品| 欧美精品啪啪一区二区三区| 欧美日韩黄片免| 91成人精品电影| 国产精品国产av在线观看| videosex国产| 亚洲一区二区三区欧美精品| bbb黄色大片| av一本久久久久| 黑人操中国人逼视频| 国产精品亚洲一级av第二区| 日本五十路高清| 欧美激情 高清一区二区三区| 精品国产一区二区三区久久久樱花| 丝袜人妻中文字幕| 十分钟在线观看高清视频www| av一本久久久久| 婷婷成人精品国产| 日韩中文字幕欧美一区二区| 大型黄色视频在线免费观看| 国产免费av片在线观看野外av| 国产精品一区二区在线观看99| av网站免费在线观看视频| 三级毛片av免费| 搡老乐熟女国产| 好看av亚洲va欧美ⅴa在| 热re99久久精品国产66热6| 亚洲国产精品sss在线观看 | 一二三四在线观看免费中文在| 亚洲 国产 在线| 国产一卡二卡三卡精品| 欧美激情高清一区二区三区| 后天国语完整版免费观看| 国产精品久久电影中文字幕 | 飞空精品影院首页| 日韩精品免费视频一区二区三区| 久久香蕉国产精品| 丝袜人妻中文字幕| 成年女人毛片免费观看观看9 | 中文亚洲av片在线观看爽 | 在线观看免费午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 天天躁日日躁夜夜躁夜夜| 超碰97精品在线观看| 亚洲精品久久成人aⅴ小说| 老司机影院毛片| 后天国语完整版免费观看| 国产激情久久老熟女| 搡老熟女国产l中国老女人| 成人特级黄色片久久久久久久| 久久国产精品男人的天堂亚洲| 99在线人妻在线中文字幕 | 国产精品久久久久成人av| 亚洲国产欧美日韩在线播放| av在线播放免费不卡| 99精国产麻豆久久婷婷| 欧美人与性动交α欧美精品济南到| 午夜成年电影在线免费观看| 一本综合久久免费| 纯流量卡能插随身wifi吗| 91麻豆精品激情在线观看国产 | 亚洲九九香蕉| 一级作爱视频免费观看| 成人影院久久| 欧美精品啪啪一区二区三区| 国产精品影院久久| 成年人免费黄色播放视频| 黄片大片在线免费观看| 午夜久久久在线观看| 欧美人与性动交α欧美软件| 乱人伦中国视频| www.999成人在线观看| a级毛片在线看网站| 在线av久久热| 色综合婷婷激情| 老汉色∧v一级毛片| 黄色女人牲交| 人成视频在线观看免费观看| 十分钟在线观看高清视频www| 国产1区2区3区精品| 身体一侧抽搐| 91麻豆精品激情在线观看国产 | 亚洲精品中文字幕一二三四区| x7x7x7水蜜桃| 久久国产亚洲av麻豆专区| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 欧美性长视频在线观看| 男女午夜视频在线观看| 日韩熟女老妇一区二区性免费视频| 桃红色精品国产亚洲av| 国产激情久久老熟女| 岛国在线观看网站| 丁香欧美五月| 黄频高清免费视频| 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 999久久久国产精品视频| 啦啦啦 在线观看视频| 制服诱惑二区| 夜夜躁狠狠躁天天躁| 午夜两性在线视频| 国产一卡二卡三卡精品| 18在线观看网站| 欧美精品av麻豆av| 国产精品免费视频内射| 国产高清激情床上av| 亚洲国产欧美一区二区综合| 午夜久久久在线观看| 欧美日韩精品网址| 日韩欧美一区视频在线观看| 亚洲黑人精品在线| 欧美成人免费av一区二区三区 | 免费看a级黄色片| 黄色视频,在线免费观看| 国产国语露脸激情在线看| 成人精品一区二区免费| av线在线观看网站| 国产一区二区三区在线臀色熟女 | 在线观看免费午夜福利视频| 国精品久久久久久国模美| 亚洲专区中文字幕在线| 久久精品人人爽人人爽视色| 精品国产一区二区久久| 看片在线看免费视频| 成人国产一区最新在线观看| 久久久国产精品麻豆| 国产单亲对白刺激| 亚洲国产欧美网| 亚洲精品国产区一区二| 99热只有精品国产| 麻豆乱淫一区二区| 国产精品一区二区在线观看99| 黄色丝袜av网址大全| 男女午夜视频在线观看| 精品免费久久久久久久清纯 | 午夜视频精品福利| 19禁男女啪啪无遮挡网站| 99久久国产精品久久久| 9191精品国产免费久久| 亚洲美女黄片视频| 高清黄色对白视频在线免费看| 熟女少妇亚洲综合色aaa.| av免费在线观看网站| 精品国产美女av久久久久小说| 18在线观看网站| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 免费高清在线观看日韩| 久久久久视频综合| 9热在线视频观看99| 国产91精品成人一区二区三区| 欧美精品一区二区免费开放| 国产99久久九九免费精品| 嫩草影视91久久| 黄频高清免费视频| 久久ye,这里只有精品| 国产无遮挡羞羞视频在线观看| 亚洲精品国产精品久久久不卡| 91麻豆av在线| 97人妻天天添夜夜摸| 一本一本久久a久久精品综合妖精| 亚洲色图 男人天堂 中文字幕| 成人精品一区二区免费| 久久国产乱子伦精品免费另类| 交换朋友夫妻互换小说| 最近最新中文字幕大全电影3 | 中文字幕av电影在线播放| 欧美激情极品国产一区二区三区| x7x7x7水蜜桃| 久久久久国产精品人妻aⅴ院 | 无限看片的www在线观看| 丝瓜视频免费看黄片| 亚洲欧美一区二区三区黑人| 丰满人妻熟妇乱又伦精品不卡| 日本黄色视频三级网站网址 | 国产高清国产精品国产三级| 久久久精品区二区三区| 午夜福利免费观看在线| 99久久综合精品五月天人人| 国产亚洲av高清不卡| 国产精品国产av在线观看| 69精品国产乱码久久久| 久久国产精品影院| 精品福利永久在线观看| 高潮久久久久久久久久久不卡| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕大全免费视频| 中文字幕高清在线视频| 久久天堂一区二区三区四区| 热99re8久久精品国产| av免费在线观看网站| 午夜亚洲福利在线播放| 国产1区2区3区精品| 19禁男女啪啪无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 免费不卡黄色视频| 一进一出抽搐gif免费好疼 | 一级作爱视频免费观看| 国产精品国产高清国产av | 美女福利国产在线| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 国产人伦9x9x在线观看| 国产一卡二卡三卡精品| 正在播放国产对白刺激| 亚洲国产看品久久| 欧美黄色淫秽网站| 老熟女久久久| 日本wwww免费看| 日本黄色视频三级网站网址 | 一二三四在线观看免费中文在| 亚洲av熟女| 成年人午夜在线观看视频| 99精国产麻豆久久婷婷| 精品国产国语对白av| 午夜福利免费观看在线| 美女福利国产在线| 成人亚洲精品一区在线观看| 99国产精品99久久久久| 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 一本一本久久a久久精品综合妖精| 国产人伦9x9x在线观看| 久久天堂一区二区三区四区| √禁漫天堂资源中文www| 在线观看一区二区三区激情| 18禁裸乳无遮挡动漫免费视频| 每晚都被弄得嗷嗷叫到高潮| 十八禁人妻一区二区| 国产精品免费视频内射| 狠狠婷婷综合久久久久久88av| 国产在线精品亚洲第一网站| 亚洲精品一二三| 国产在线精品亚洲第一网站| 天堂动漫精品| av国产精品久久久久影院| 国产精品九九99| 国产精品免费大片| 一级黄色大片毛片| 国产精品九九99| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩另类电影网站| 大香蕉久久网| 亚洲欧洲精品一区二区精品久久久| 亚洲少妇的诱惑av| 国产色视频综合| 曰老女人黄片| 国产免费现黄频在线看| 国产aⅴ精品一区二区三区波| 在线视频色国产色| 操出白浆在线播放| 欧美成人午夜精品| av天堂久久9| 欧美精品高潮呻吟av久久| 波多野结衣一区麻豆| 9热在线视频观看99| 国产精品自产拍在线观看55亚洲 | 午夜福利欧美成人| 很黄的视频免费| 久久久久久免费高清国产稀缺| 精品福利永久在线观看| 午夜91福利影院| 国产欧美亚洲国产| 成人影院久久| 久久狼人影院| 日韩一卡2卡3卡4卡2021年| 黑人巨大精品欧美一区二区mp4| 天堂中文最新版在线下载| 国精品久久久久久国模美| 18禁国产床啪视频网站| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 久久九九热精品免费| 亚洲欧美精品综合一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦在线免费观看视频4| 欧美激情 高清一区二区三区| 亚洲国产中文字幕在线视频| 午夜影院日韩av| 视频区图区小说| 老司机靠b影院| 超色免费av| 国产精品免费视频内射| 男女免费视频国产| av电影中文网址| 90打野战视频偷拍视频| 亚洲精品久久午夜乱码| 久久精品国产清高在天天线| 国产精品久久电影中文字幕 | 一个人免费在线观看的高清视频| 午夜免费鲁丝| 黄片小视频在线播放| 色尼玛亚洲综合影院| 久久精品熟女亚洲av麻豆精品| 女人被狂操c到高潮| 十八禁高潮呻吟视频| 精品熟女少妇八av免费久了| 国产精品久久久久久精品古装| 高清毛片免费观看视频网站 | 熟女少妇亚洲综合色aaa.| 欧美精品高潮呻吟av久久| 色婷婷av一区二区三区视频| 好男人电影高清在线观看| 午夜精品在线福利| 天堂俺去俺来也www色官网| 1024视频免费在线观看| 欧美精品亚洲一区二区| 99国产精品一区二区蜜桃av | 久久香蕉国产精品| av国产精品久久久久影院| 高清视频免费观看一区二区| 国产日韩一区二区三区精品不卡| 国产主播在线观看一区二区| 校园春色视频在线观看| 搡老岳熟女国产| 老司机午夜福利在线观看视频| tube8黄色片| 亚洲一区高清亚洲精品| 亚洲成人免费av在线播放| 国产区一区二久久| 精品久久久久久,| 视频区图区小说| 日韩制服丝袜自拍偷拍| 国产成人一区二区三区免费视频网站| 三上悠亚av全集在线观看| av福利片在线| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人看| 99久久人妻综合| 免费不卡黄色视频| 国产精品影院久久| 两人在一起打扑克的视频| 午夜影院日韩av| 亚洲精品乱久久久久久| 亚洲欧美色中文字幕在线| 中出人妻视频一区二区| 美女午夜性视频免费| 操出白浆在线播放| 91字幕亚洲| 性少妇av在线| 天天躁狠狠躁夜夜躁狠狠躁| 99久久综合精品五月天人人| 一a级毛片在线观看| 国产精品永久免费网站| 1024香蕉在线观看| 91老司机精品| 久久香蕉国产精品| 91在线观看av|