SDP-based approximation of stabilising solutions for periodic matrix Riccati differential equations

ABSTRACT Numerically finding stabilising feedback control laws for linear systems of periodic differential equations is a nontrivial task with no known reliable solutions. The most successful method requires solving matrix differential Riccati equations with periodic coefficients. All previously proposed techniques for solving such equations involve numerical integration of unstable differential equations and consequently fail whenever the period is too large or the coefficients vary too much. Here, a new method for numerical computation of stabilising solutions for matrix differential Riccati equations with periodic coefficients is proposed. Our approach does not involve numerical solution of any differential equations. The approximation for a stabilising solution is found in the form of a trigonometric polynomial, matrix coefficients of which are found solving a specially constructed finite-dimensional semidefinite programming (SDP) problem. This problem is obtained using maximality property of the stabilising solution of the Riccati equation for the associated Riccati inequality and sampling technique. Our previously published numerical comparisons with other methods shows that for a class of problems only this technique provides a working solution. Asymptotic convergence of the computed approximations to the stabilising solution is proved below under the assumption that certain combinations of the key parameters are sufficiently large. Although the rate of convergence is not analysed, it appeared to be exponential in our numerical studies.

[1]  J. Willems Least squares stationary optimal control and the algebraic Riccati equation , 1971 .

[2]  Vladimir B. Larin,et al.  Algorithm for Solving the Discrete Periodic Riccati Equation , 2002 .

[3]  Patrizio Colaneri,et al.  A note on the maximal solution of the periodic Riccati equation , 1989 .

[4]  Leonid B. Freidovich,et al.  Virtual-Holonomic-Constraints-Based Design of Stable Oscillations of Furuta Pendulum: Theory and Experiments , 2007, IEEE Transactions on Robotics.

[5]  Leonid B. Freidovich,et al.  Transverse Linearization for Controlled Mechanical Systems With Several Passive Degrees of Freedom , 2010, IEEE Transactions on Automatic Control.

[6]  H. Walker,et al.  Numerical Linear Algebra with Applications , 1994 .

[7]  Zhi-Gang Wu,et al.  Fourier expansion based recursive algorithms for periodic Riccati and Lyapunov matrix differential equations , 2011, J. Comput. Appl. Math..

[8]  E. Yaz Linear Matrix Inequalities In System And Control Theory , 1998, Proceedings of the IEEE.

[9]  Vlad Ionescu,et al.  The periodic Riccati equation , 2003 .

[10]  Sergei V. Gusev Kalman-Yakubovich-Popov Lemma for matrix frequency domain inequality , 2009, Syst. Control. Lett..

[11]  Carlos Canudas-de-Wit,et al.  Constructive tool for orbital stabilization of underactuated nonlinear systems: virtual constraints approach , 2005, IEEE Transactions on Automatic Control.

[12]  Mark A. Shayman On the Phase Portrait of the Matrix Riccati Equation Arising from the Periodic Control Problem , 1985 .

[13]  András Varga,et al.  On solving periodic Riccati equations , 2008, Numer. Linear Algebra Appl..

[14]  Theodore Kilgore An elementary simultaneous approximation theorem , 1993 .

[15]  Andreas Varga,et al.  Computational issues for linear periodic systems: paradigms, algorithms, open problems , 2013, Int. J. Control.

[16]  D. Kleinman On an iterative technique for Riccati equation computations , 1968 .

[17]  Kim-Chuan Toh,et al.  Solving semidefinite-quadratic-linear programs using SDPT3 , 2003, Math. Program..

[18]  R. Johansson,et al.  Periodic motions of the Pendubot via virtual holonomic constraints: Theory and experiments , 2008, Autom..

[19]  Leonid B. Freidovich,et al.  New approach for swinging up the Furuta pendulum : theory and experiments , 2009 .

[20]  Brian D. O. Anderson,et al.  A New Iterative Algorithm to Solve Periodic Riccati Differential Equations With Sign Indefinite Quadratic Terms , 2008, IEEE Transactions on Automatic Control.

[21]  Leonid B. Freidovich,et al.  LMI approach for solving periodic matrix Riccati equation , 2007, PSYCO.

[22]  A. Laub A schur method for solving algebraic Riccati equations , 1978, 1978 IEEE Conference on Decision and Control including the 17th Symposium on Adaptive Processes.

[23]  V. Yakubovich A linear-quadratic optimization problem and the frequency theorem for nonperiodic systems. I , 1986 .

[24]  R. S. Bucy,et al.  Periodic equilibria for matrix riccati equations , 1979 .

[25]  J. Lofberg,et al.  YALMIP : a toolbox for modeling and optimization in MATLAB , 2004, 2004 IEEE International Conference on Robotics and Automation (IEEE Cat. No.04CH37508).

[26]  Anders Lindquist,et al.  From Finite Covariance Windows to Modeling Filters: A Convex Optimization Approach , 2001, SIAM Rev..

[27]  Bo Kågström,et al.  A numerical evaluation of solvers for the periodic Riccati differential equation , 2010 .

[28]  Y. Ershov Mapping a restriction of spaces of orderings of fields , 1986 .

[29]  Yang Zhou Chen,et al.  Comparison and uniqueness theorems for periodic Riccati differential equations , 1998 .

[30]  P. Parrilo On the numerical solution of LMIs derived from the KYP lemma , 1999, Proceedings of the 38th IEEE Conference on Decision and Control (Cat. No.99CH36304).

[31]  Arkadi Nemirovski,et al.  Lmi Control Toolbox For Use With Matlab , 2014 .

[32]  R. E. Kalman,et al.  Contributions to the Theory of Optimal Control , 1960 .