Stability optimization of uncertain timedelay systems
MATLAB package
 Author
 Corresponding articles
 Download
 All scientific publications, for which the Matlab package has been used, must refer to the above publications.
 Description of the method
 Acknowledgments
The software was written by Fenzi Luca under supervision of Wim Michiels.
L. Fenzi and W. Michiels, Robust stability optimization for linear delay systems in a probabilistic framework, Linear Algebra and its Application, 526:126 (2017).
(Preprint.)L. Fenzi and W. Michiels, Robust stability optimization for linear delay systems in a probabilistic framework,
Internal report TW 671, Department of Computer Science, KU Leuven, August 2016.
Download the software.
All the codes are commented and there are several examples, including a model of an experimental heatexchanger. The software is released under the GNU GPL v3.0 license.
This eigenvaluebased stabilization method tunes the controller parameters (static or dynamic feedbacks) in order to improve the stability properties of a timedelay system, which can be affected by uncertainties, modeled by the realizations of a random vector. The closedloop system is described by a delay differential algebraic equations (DDAE) of retarded type, in this way we can take into account integral controller and distributed terms. The timedelay system may nonlinearly depend on the controller and uncertain parameters.
In order to take into account the uncertainty, the stabilization minimizes an objective function consisting of the mean of the spectral abscissa with a variance penalty. For every realization of the uncertainties, the spectral abscissa (real part of the rightmost eigenvalue) is computed with the Infinitesimal Generator Approach, and then corrected with Newton's method. The objective function and its gradient are numerically evaluated by computing integrals using quasiMonte Carlo methods. The minimization of the objective function relies on the software HANSO (Hybrid Algorithm for Non Smooth Optimization).
Special thanks goes to Dan Pilbauer for his careful tests and constructive comments.
Funding: This work was supported by the project C14/17/072 of the KU Leuven Research Council, by the project G0A5317N of the Research FoundationFlanders (FWO  Vlaanderen), and by the project UCoCoS, funded by the European Unions Horizon 2020 research and innovation program under the Marie SklodowskaCurie Grant Agreement No 675080.
Polynomial (chaos) approximation of maximum eigenvalue functions
MATLAB tutorial
 Authors
 Corresponding articles
 Download
 All scientific publications, for which the tutorial has been used, must refer to the tutorial, and the above publications.
 Abstract of the Tutorial

Funding: This work was supported by the project C14/17/072 of the KU Leuven Research Council, by the project G0A5317N of the Research FoundationFlanders (FWO  Vlaanderen), and by the project UCoCoS, funded by the European Unions Horizon 2020 research and innovation program under the Marie SklodowskaCurie Grant Agreement No 675080.
The software was written by Fenzi Luca under supervision of Wim Michiels.
L. Fenzi and W. Michiels.
Polynomial (chaos) approximation of maximum eigenvalue functions: efficiency and limitations, Numerical Algorithms, 82:1143–1169 (2019).
(Preprint) L. Fenzi and W. Michiels.
Polynomial (chaos) approximation of maximum eigenvalue functions: efficiency and limitations, arXiv: 1804.03881 (2018).
Download the tutorial. This tutorial is also available at TW internal reports KU Leuven.
Download the Matlab script. This script generates the latex tutorial by the Matlab code publish.
This tutorial reviews the numerical experiments contained in the article, Fenzi & Michiels (2019) "Polynomial (chaos) approximation of maximum eigenvalue functions: efficiency and limitations", providing a template that can be modified for explorations of your own.
The tutorial explores the polynomial approximation of smooth, nondifferentiable and nonLipschitz continuous functions in the univariate and bivariate cases. The analyzed functions arise from parameter eigenvalue problems; in particular, they are the real part of the rightmost eigenvalue (the socalled spectral abscissa).
The polynomial approximations are obtained by Galerkin and collocation approaches. In the Galerkin approach, the numerical approximation of the coefficients in the univariate case is achieved by extended (or composite) Trapezoidal and Simpson's rules or by Gauss and ClenshawCurtis quadrature rules. For the bivariate case, the coefficients are approximated by tensorial and nontensorial ClenshawCurtis cubature rules, based on tensorproduct Chebyshev grid and Padua points, respectively. The collocation approach interpolates the function on Chebyshev points in the univariate case, while for the bivariate case the interpolant nodes are given by tensorproduct Chebyshev grid and Padua points.
Stability and Stabilization of Periodic TimeDelay Systems
MATLAB package
 Authors
 Corresponding articles
 Download
 All scientific publications, for which the Matlab package has been used, must refer to the above publications.
 Description of the software

Funding: This work was supported by the project C14/17/072 of the KU Leuven Research Council, by the project G0A5317N of the Research FoundationFlanders (FWO  Vlaanderen), and by the project UCoCoS, funded by the European Unions Horizon 2020 research and innovation program under the Marie SklodowskaCurie Grant Agreement No 675080.
The software was written by Fenzi Luca under supervision of Wim Michiels.
W. Michiels and L. Fenzi.
Spectrumbased stability analysis and stabilization of a class of timeperiodic time delay systems,
arXiv: 1804.03881 (2019).
Download the software. All the codes are commented and there are several examples.
The software is released under the GNU GPL v3.0 license.
The software deals with spectrumbased stability assessment and stabilization methods for periodic linear timedelay systems, where the delays and the period of the system matrices are commensurate. The software relies on the dual interpretation of the Floquet multipliers, which can be characterized as nonzero eigenvalues either of the monodromy operor or of a finitedimensional nonlinear eigenvalue problem, whose solution is obtained by solving an initial value problem.
The stability assessment of a periodic timedelay system is determined by the largest in modulus Floquet multiplier. This Floquet multiplier is first approximated by a spectral collocation of the monodromy operator eigenproblem. Hence, its precision is refined up to machine precision by Broyden's method, considering the finitedimensional nonlinear eigenvalue problem.
The derivate of the Floquet multiplier with respect to system parameter can also be computed. The derivative computation is used in the stabilization method, which designs targeted parameters of the timedelay systems in order to minimize the modulus of the largest Floquet multipliers. The stabilization method requires the usage of the software HANSO (Hybrid Algorithm for Non Smooth Optimization).
The package presents additional routines for the stability and stabilization of largescale periodic timedelay systems. These largescale routines solve with Arnoldi's method the discretized aigenvalue problem arising from the monodromy operator. Moreover, they consider a dual timesystems for the Floquet multiplier derivative computation.