Numerical integration of the differential matrix Riccati equation

Two new Bernoulli substitution methods for solving the Riccati differential equation are tested numerically against direct integration of the Riccati equation, the Chandrasekhar algorithm, and the Davison-Maki method on a large set of problems taken from the literature. The first of these new methods was developed for the time-invariant case and uses the matrix analog of completing the square to transform the problem to a bisymmetric Riecati equation whose solution can be given explicitly in terms of a matrix exponential of order n . This method is fast and accurate when the extremal solutions of the associated algebraic Riccati equation are well separated. The second new method was developed as a means of eliminating the instabilities associated with the Davison-Maki algorithm. By using reinitialization at each time step the Davison-Maki algorithm can be recast as a recursion which is over three times faster than the original method and is easily shown to be stable for both time-invariant and time-dependent problems. From the results of our study we conclude that the modified Davison-Maki method gives superior performance except for those problems where the number of observers and controllers is small relative to the number of states in which ease the Chandrasekhar algorithm is better.

[1]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[2]  I. Bar-Itzhack,et al.  Eigenfactor solution of the matrix Riccati equation - A continuous square root algorithm , 1984, The 23rd IEEE Conference on Decision and Control.

[3]  Some New Analytic and Computational Results for Operator Riccati Equations , 1975 .

[4]  Alan J. Laub,et al.  Schur techniques for riccati differential equations , 1982 .

[5]  Comments on "The numerical solution of the matrix Ricatti differential equation" , 1974 .

[6]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[7]  R. Ward.  Numerical Computation of the Matrix Exponential with Accuracy Estimate , 1977 .

[8]  Roy B. Leipnik.  A canonical form and solution for the matrix Riccati differential equation , 1985, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

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

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

[11]  Edward J. Davison,et al.  The numerical solution of the matrix Riccati differential equation , 1973 .

[12]  W. A. Coppel.  Matrix quadratic equations , 1974, Bulletin of the Australian Mathematical Society.

[13]  Demetrios Lainiotis.  Generalized Chandrasekhar algorithms: Time-varying models , 1975, 1975 IEEE Conference on Decision and Control including the 14th Symposium on Adaptive Processes.

[14]  Demetrios G. Lainiotis,et al.  Partitioned Ricatti solutions and integration-free doubling algorithms , 1976 .

[15]  Thomas Kailath,et al.  Some new algorithms for recursive estimation in constant linear systems , 1973, IEEE Trans. Inf. Theory.

[16]  Numerical experiments in linear control theory using generalized X - Y equations , 1976 .