Approximation of low rank solutions for linear quadratic control of partial differential equations

Algebraic Riccati equations (ARE) of large dimension arise when using approximations to design controllers for systems modeled by partial differential equations. We use a modified Newton method to solve the ARE that takes advantage of several special features of these problems. The modified Newton method leads to a right-hand side of rank equal to the number of inputs regardless of the weights. Thus, the resulting Lyapunov equation can be more efficiently solved. The Cholesky-ADI algorithm is used to solve the Lyapunov equation resulting at each step. The algorithm is straightforward to code. Performance is illustrated with a number of standard examples. An example on controlling the deflection of the Euler-Bernoulli beam indicates that for weakly damped problems a low rank solution to the ARE may not exist. Further analysis supports this point.

[1]  Alfredo Germani,et al.  Approximation of the algebraic Riccati equation in the Hilbert space of Hilbert-Schmidt operators , 1993 .

[2]  Kirsten Morris,et al.  SOLUTION OF ALGEBRAIC RICCATI EQUATIONS ARISING IN CONTROL OF PARTIAL DIFFERENTIAL EQUATIONS , 2007 .

[3]  K. Morris,et al.  Solving the linear quadratic optimal control problem for infinite-dimensional systems , 1996 .

[4]  Mathukumalli Vidyasagar,et al.  A Comparison of Different Models for Beam Vibrations From the Standpoint of Control Design , 1990 .

[5]  Eugene L. Wachspress,et al.  Alternating direction implicit iteration for systems with complex spectra , 1991 .

[6]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[7]  Thilo Penzl,et al.  A Cyclic Low-Rank Smith Method for Large Sparse Lyapunov Equations , 1998, SIAM J. Sci. Comput..

[8]  Kazufumi Ito,et al.  Approximation in LQR problems for infinite dimensional systems with unbounded input operators , 1994 .

[9]  Ruth F. Curtain,et al.  Model reduction for control design for distributed parameter systems , 2003 .

[10]  H. Banks,et al.  A numerical algorithm for optimal feedback gains in high dimensional linear quadratic regulator problems , 1991 .

[11]  John A. Burns,et al.  Numerical methods for approximating functional gains in LQR boundary control problems , 2001 .

[12]  J. Gibson.  Linear-Quadratic Optimal Control of Hereditary Differential Systems: Infinite Dimensional Riccati Equations and Numerical Approximations , 1983 .

[13]  Kirsten Morris,et al.  Introduction to Feedback Control , 2001 .

[14]  J. S. Gibson,et al.  The Riccati Integral Equations for Optimal Control Problems on Hilbert Spaces , 1979 .

[15]  Gene H. Golub,et al.  Matrix computations , 1983 .

[16]  Thilo Penzl.  Eigenvalue decay bounds for solutions of Lyapunov equations: the symmetric case , 2000 .

[17]  Paul Van Dooren,et al.  A collection of benchmark examples for model reduction of linear time invariant dynamical systems. , 2002 .

[18]  K. A. Morrisy.  Design of Finite-dimensional Controllers for Innnite-dimensional Systems by Approximation , 1994 .

[19]  Eugene L. Wachspress,et al.  Solution of lyapunov equations by alternating direction implicit iteration , 1991 .

[20]  Mansour Eslami,et al.  Mathematics of finite-dimensional control systems, theory and design , 1979, IEEE Transactions on Systems, Man, and Cybernetics.

[21]  Hans Zwart,et al.  An Introduction to Infinite-Dimensional Linear Systems Theory , 1995, Texts in Applied Mathematics.

[22]  Kazufumi Ito,et al.  An Approximation Theory of Solutions to Operator Riccati Equations for $H^\infty$ Control , 1998 .

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

[24]  Jacob K. White,et al.  Low-Rank Solution of Lyapunov Equations , 2004, SIAM Rev..

[25]  Karl Kunisch,et al.  The linear regulator problem for parabolic systems , 1984 .

[26]  Peter Benner,et al.  Efficient Algorithms for Large‐Scale Quadratic Matrix Equations , 2002 .

[27]  Kazufumi Ito,et al.  Strong convergence and convergence rates of approximating solutions for algebraic riccati equations in Hilbert spaces , 1987 .

[28]  K. A. Morris,et al.  A Variational Calculus Approach to the Modelling of Flexible Manipulators , 1996, SIAM Rev..

[29]  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.

[30]  J. Potter.  Matrix Quadratic Solutions , 1966 .

[31]  E. Wachspress.  Iterative solution of the Lyapunov matrix equation , 1988 .

[32]  R. A. Smith.  Matrix Equation $XA + BX = C$ , 1968 .

[33]  P. Benner,et al.  Solving large-scale control problems , 2004, IEEE Control Systems.

[34]  Y. Zhou,et al.  On the decay rate of Hankel singular values and related issues , 2002, Syst. Control. Lett..

[35]  Daniel J. Inman,et al.  On damping mechanisms in beams , 1991 .

[36]  K. Ito,et al.  AN APPROXIMATION THEORY OF SOLUTIONS TO OPERATOR RICCATI EQUATIONS FOR H∞ CONTROL∗ , 1994 .

[37]  Tongxing Lu,et al.  Solution of the matrix equation AX−XB=C , 2005, Computing.

[38]  I. Rosen,et al.  A multilevel technique for the approximate solution of operator Lyapunov and algebraic Riccati equations , 1995 .

[39]  R. Triggiani,et al.  Control Theory for Partial Differential Equations: Continuous and Approximation Theories , 2000 .