Sample Complexity of the Robust LQG Regulator with Coprime Factors Uncertainty

This paper addresses the end-to-end sample complexity bound for learning the H2 optimal controller (the Linear Quadratic Gaussian (LQG) problem) with unknown dynamics, for potentially unstable Linear Time Invariant (LTI) systems. The robust LQG synthesis procedure is performed by considering bounded additive model uncertainty on the coprime factors of the plant. The closed-loop identification of the nominal model of the true plant is performed by constructing a Hankel-like matrix from a single time-series of noisy finite length input-output data, using the ordinary least squares algorithm from Sarkar et al. (2020). Next, an H-infinity bound on the estimated model error is provided and the robust controller is designed via convex optimization, much in the spirit of Boczar et al. (2018) and Zheng et al. (2020a), while allowing for bounded additive uncertainty on the coprime factors of the model. Our conclusions are consistent with previous results on learning the LQG and LQR controllers.Comment: Minor Edits on closed loop identification, 30 pages, 2 figures, 3 algorithm

Approximation methodologies for explicit model predictive control of complex systems

This thesis concerns the development of complexity reduction methodologies for the application of multi-parametric/explicit model predictive (mp-MPC) control to complex high fidelity models. The main advantage of mp-MPC is the offline relocation of the optimization task and the associated computational expense through the use of multi-parametric programming. This allows for the application of MPC to fast sampling systems or systems for which it is not possible to perform online optimization due to cycle time requirements. The application of mp-MPC to complex nonlinear systems is of critical importance and is the subject of the thesis. The first part is concerned with the adaptation and development of model order reduction (MOR) techniques for application in combination to mp-MPC algorithms. This first part includes the mp-MPC oriented use of existing MOR techniques as well as the development of new ones. The use of MOR for multi-parametric moving horizon estimation is also investigated. The second part of the thesis introduces a framework for the ‘equation free’ surrogate-model based design of explicit controllers as a possible alternative to multi-parametric based methods. The methodology relies upon the use of advanced data-classification approaches and surrogate modelling techniques, and is illustrated with different numerical examples.Open Acces

Advanced multiparametric optimization and control studies for anaesthesia

Anaesthesia is a reversible pharmacological state of the patient where hypnosis, analgesia and muscle relaxation are guaranteed and maintained throughout the surgery. Analgesics block the sensation of pain; hypnotics produce unconsciousness, while muscle relaxants prevent unwanted movement of muscle tone. Controlling the depth of anaesthesia is a very challenging task, as one has to deal with nonlinearity, inter- and intra-patient variability, multivariable characteristics, variable time delays, dynamics dependent on the hypnotic agent, model analysis variability, agent and stability issues. The modelling and automatic control of anaesthesia is believed to (i) benefit the safety of the patient undergoing surgery as side-effects may be reduced by optimizing the drug infusion rates, and (ii) support anaesthetists during critical situations by automating the drug delivery systems. In this work we have developed several advanced explicit/multi-parametric model predictive (mp-MPC) control strategies for the control of depth of anaesthesia. State estimation techniques are developed and used simultaneously with mp-MPC strategies to estimate the state of each individual patient, in an attempt to overcome the challenges of inter- and intra- patient variability, and deal with possible unmeasurable noisy outputs. Strategies to deal with the nonlinearity have been also developed including local linearization, exact linearization as well as a piece-wise linearization of the Hill curve leading to a hybrid formulation of the patient model and thereby the development of multiparametric hybrid model predictive control methodology. To deal with the inter- and intra- patient variability, as well as the noise on the process output, several robust techniques and a multiparametric moving horizon estimation technique have been design and implemented. All the studies described in the thesis are performed on clinical data for a set of 12 patients who underwent general anaesthesia.Open Acces

Numerical Techniques for Optimization Problems with PDE Constraints

The development, analysis and implementation of efficient and robust numerical techniques for optimization problems associated with partial differential equations (PDEs) is of utmost importance for the optimal control of processes and the optimal design of structures and systems in modern technology. The successful realization of such techniques invokes a wide variety of challenging mathematical tasks and thus requires the application of adequate methodologies from various mathematical disciplines. During recent years, significant progress has been made in PDE constrained optimization both concerning optimization in function space according to the paradigm ’Optimize first, then discretize’ and with regard to the fast and reliable solution of the large-scale problems that typically arise from discretizations of the optimality conditions. The contributions at this Oberwolfach workshop impressively reflected the progress made in the field. In particular, new insights have been gained in the analysis of optimal control problems for PDEs that have led to vastly improved numerical solution methods. Likewise, breakthroughs have been made in the optimal design of structures and systems, for instance, by the socalled ’all-at-once’ approach featuring simultaneous optimization and solution of the underlying PDEs. Finally, new methodologies have been developed for the design of innovative materials and the identification of parameters in multi-scale physical and physiological processes

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties

Structure-Preserving Model Reduction of Physical Network Systems

This paper considers physical network systems where the energy storage is naturally associated to the nodes of the graph, while the edges of the graph correspond to static couplings. The first sections deal with the linear case, covering examples such as mass-damper and hydraulic systems, which have a structure that is similar to symmetric consensus dynamics. The last section is concerned with a specific class of nonlinear physical network systems; namely detailed-balanced chemical reaction networks governed by mass action kinetics. In both cases, linear and nonlinear, the structure of the dynamics is similar, and is based on a weighted Laplacian matrix, together with an energy function capturing the energy storage at the nodes. We discuss two methods for structure-preserving model reduction. The first one is clustering; aggregating the nodes of the underlying graph to obtain a reduced graph. The second approach is based on neglecting the energy storage at some of the nodes, and subsequently eliminating those nodes (called Kron reduction).</p

ON THE STATE-OPTIMIZATION APPROACH TO SYSTEM PROBLEMS: CLOSED-LOOP THINKING SOLUTIONS

State-optimization approach has been proposed to treating various different system problems in optimal projection equations (OPEQ). While the OPEQ for problems of open-loop thinking is found consisting of two modified Lyapunov equations, excepting the rank conditions whereas required in system identification and its related robust problems, the one for closed-loop thinking consists of two modified either Reccatti or Lyapunov equations, excepting conditions for compensating system happened to be in a problem like that of order reduction for controller. Apart from addditonally constrained-conditions and simplicity in the solution form have been obtainable for each problem, it has been found the system identification problem switching over to computing the solution of OPEQ and the physical nature of medeled states possibly retaining in optimal order reduction problem