A Precise Integration Algorithm for Matrix Riccati Differential Equations

Wan-Xie Zhong,Jianping Zhu

DOI: https://doi.org/10.1007/3-540-45545-0_107

2001-01-01

Abstract:An efficient precise integration method for solving the matrix Riccati differential equation is described in this paper. The method is based on repeated combination of extremely small time intervals, which leads to solutions with an accuracy within the machine precision.

What problem does this paper attempt to address?

A Note on the Necessary Conditions for the Algebraic Riccati Equation.

Yongliang Zhu,Prabhakar R. Pagilla

DOI: https://doi.org/10.1093/imamci/dni017

2005-01-01

IMA Journal of Mathematical Control and Information

Abstract:We derive some useful and easily computable necessary conditions for the existence of a positive semi-definite solution to the algebraic Ricatti equation (ARE). A motivating example is given to highlight the usefulness of the conditions for controller and observer designs for nonlinear systems. Further, an upper bound on the trace of the solution to the ARE is also derived.
A Memory-Efficient Formulation of Precise-Integration Time-Domain Method With Riccati Matrix Differential Equations

Xiaojie Zhu,Xikui Ma,Jinghui Shao,Shuli Yin

DOI: https://doi.org/10.1109/TMAG.2019.2952152

IF: 1.848

2020-01-01

IEEE Transactions on Magnetics

Abstract:A novel numerical method, referred to as the Riccati precise-integration time-domain (Riccati-PITD) method, is proposed to reduce the memory requirement of the conventional PITD method. In the Riccati-PITD method, the electromagnetic field components of the transverse electric/transverse magnetic (TE/TM) wave are arranged in a matrix for constructing the Riccati matrix differential equations (RDEs) about Maxwell's curl equations. Then the precise integration technique is adopted for the numerical integration of the RDEs. Theoretical analyses about the memory requirements of the conventional and the Riccati-PITD methods are presented. It is found that the memory requirement can be more greatly reduced by solving the electromagnetic-field components in the form of the matrices in the Riccati-PITD method rather than in the form of the vectors in the conventional PITD method. Numerical examples are also given to confirm that the proposed method has significant advantages over the conventional PITD method with respect to the execution time and the memory requirement.
An Implicit Iterative Algorithm for Algebraic Riccati Matrix Equations

Ai-Guo Wu,Yu Wang,Ying Zhang

DOI: https://doi.org/10.23919/chicc.2017.8027333

2017-01-01

Abstract:In this paper, an implicit iterative algorithm is proposed to obtain the unique positive definite solution of the continuous algebraic Riccati matrix equation. In this proposed algorithm, there exists a tuning parameter which can be appropriately chosen such that the algorithm achieves better convergence performance. A sufficient condition is given for the convergence of the proposed algorithm. Moreover, an approach is also provided to choose an appropriate initial condition. Finally, a numerical example is employed to show the effectiveness of the proposed algorithm.
An Inversion-Free Iterative Algorithm for Riccati Matrix Equations in Discrete-Time Markov Jump Systems

Zhi Li,Ying Zhang,Ai-Guo Wu

DOI: https://doi.org/10.1109/tac.2022.3172266

2021-01-01

Abstract:In this article, a novel iterative algorithm is proposed to obtain the positive definite solutions of the discrete coupled Riccati matrix equations. In the presented algorithm, there are no operations of matrix inversion in all iterative steps. Another significant feature of this algorithm is that the latest updated information is used when the current estimation for each unknown matrix is conducted by taking the coupling structure of the discrete coupled Riccati matrix equations into consideration. The convergence performance of the algorithm is also analyzed. Finally, a numerical example is provided to demonstrate the effectiveness of the presented algorithm.
An implicit iterative algorithm for generalized algebraic Riccati matrix equations

Wu Aiguo,Wang Yu,Sun Huijie

DOI: https://doi.org/10.1109/ICICIP.2017.8113932

2017-01-01

Abstract:In this paper, an implicit iterative algorithm is developed to obtain the unique positive definite solution of the generalized algebraic Riccati matrix equation. For this proposed algorithm, there exisits a tuning parameter which can be chosen such that this algorithm achieves better convergence performance. Some convergence results are given for the proposed algorithm. Moreover, an approach is also provided to choose an appropriate initial condition for this algorithm. Finally, a numerical example is presented to verify the effectiveness of this algorithm.
A Novel Iterative Algorithm for Solving Coupled Riccati Equations

Ai-Guo Wu,Hui-Jie Sun,Ying Zhang

DOI: https://doi.org/10.1016/j.amc.2019.124645

IF: 4.397

2019-01-01

Applied Mathematics and Computation

Abstract:In this paper, a novel iterative algorithm is developed for solving the coupled algebraic Riccati equation arising from the quadratic optimal control problem for continuous-time Markovian jump linear systems. First, two existing iterative algorithms to solve the coupled Riccati matrix equation are reviewed. Next, based on analysis for these two algorithms, a new iterative algorithm that combines both the information in the current iterative step and the information in the last iterative step is proposed. It is shown that the proposed algorithm with proper initial conditions can monotonically converge to the unique positive definite solution of the coupled Riccati matrix equation if the associated Markovian jump system is stochastically stabilizable. Also, numerical examples show that the presented algorithm is faster than some previous algorithms when the weighted parameter is appropriately selected.
Accurate solutions of M -matrix algebraic Riccati equations

Jungong Xue,Shufang Xu,Ren-Cang Li

DOI: https://doi.org/10.1007/s00211-011-0421-0

IF: 2.5

2011-01-01

Numerische Mathematik

Abstract:This paper is concerned with the relative perturbation theory and its entrywise relatively accurate numerical solutions of an M -matrix Algebraic Riccati Equations (MARE) XDX - AX - XB + C = 0 by which we mean the following conformally partitioned matrix ( [ B -D; -C A ]) is a nonsingular or an irreducible singular M -matrix. It is known that such an MARE has a unique minimal nonnegative solution Φ . It is proved that small relative perturbations to the entries of A , B , C , and D introduce small relative changes to the entries of the nonnegative solution Φ . Thus the smaller entries Φ do not suffer bigger relative errors than its larger entries, unlike the existing perturbation theory for (general) Algebraic Riccati Equations. We then discuss some minor but crucial implementation changes to three existing numerical methods so that they can be used to compute Φ as accurately as the input data deserve. Current study is based on a previous paper of the authors’ on M -matrix Sylvester equation for which D = 0.
Iterative Algorithms with the Latest Update for Riccati Matrix Equations in Itô Markov Jump Systems

Ma Chuang,Li Zhi,Wu WanQi,Zhang Ying

DOI: https://doi.org/10.1007/s11431-020-1668-4

2020-01-01

Abstract:This study is concerned with the problem to solve the continuous coupled Riccati matrix equations in Itô Markov jump systems. A new iterative algorithm is developed by using the latest estimation information and introducing a tuning parameter. The iterative solution obtained by the proposed algorithm with zero initial conditions converges to the unique positive definite solution of the considered equations. The convergence rate of the algorithm is dependent on the adjustable parameter. Furthermore, a numerical example is provided to show the effectiveness of the presented algorithms.
Backward differentiation formula method and random forest method to solve continuous‐time differential Riccati equations

Juan Zhang,Wenwen Zou,Chenglin Sui

DOI: https://doi.org/10.1002/asjc.3494

IF: 2.4

2024-09-21

Asian Journal of Control

Abstract:In this paper, we explore the utilization of machine learning techniques for solving the numerical solutions of continuous‐time differential Riccati equations. Specifically, we focus on generating a reduction matrix capable of transforming a high‐order matrix into a low‐order matrix. Additionally, we address the issue of differential terms in the continuous‐time differential Riccati equation and incorporate the backward differentiation formula of the matrix to improve stability and accuracy. Finally, by training samples through neural networks and machine learning methods, we could predict the solutions for high‐order matrix equations.

automation & control systems
The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

Zhen-Chen Guo,Xin Liang

DOI: https://doi.org/10.1007/s11075-022-01413-9

IF: 2.37

2022-01-01

Numerical Algorithms

Abstract:In this paper, we derive a Toeplitz-structured closed form of the unique positive semi-definite stabilizing solution for the discrete-time algebraic Riccati equations, especially for the case that the state matrix is not stable. Based on the found form and fast Fourier transform, we propose a new algorithm for solving both discrete-time and continuous-time large-scale algebraic Riccati equations with low-rank structure. It works without unnecessary assumptions, complicated shift selection strategies, or matrix calculations of the cubic order with respect to the problem scale. Numerical examples are given to illustrate its features. Besides, we show that it is theoretically equivalent to several algorithms existing in the literature in the sense that they all produce the same sequence under the same parameter setting.
Low-rank plus diagonal approximations for Riccati-like matrix differential equations

Silvère Bonnabel,Marc Lambert,Francis Bach

2024-07-01

Abstract:We consider the problem of computing tractable approximations of time-dependent d x d large positive semi-definite (PSD) matrices defined as solutions of a matrix differential equation. We propose to use "low-rank plus diagonal" PSD matrices as approximations that can be stored with a memory cost being linear in the high dimension d. To constrain the solution of the differential equation to remain in that subset, we project the derivative at all times onto the tangent space to the subset, following the methodology of dynamical low-rank approximation. We derive a closed-form formula for the projection, and show that after some manipulations it can be computed with a numerical cost being linear in d, allowing for tractable implementation. Contrary to previous approaches based on pure low-rank approximations, the addition of the diagonal term allows for our approximations to be invertible matrices, that can moreover be inverted with linear cost in d. We apply the technique to Riccati-like equations, then to two particular problems. Firstly a low-rank approximation to our recent Wasserstein gradient flow for Gaussian approximation of posterior distributions in approximate Bayesian inference, and secondly a novel low-rank approximation of the Kalman filter for high-dimensional systems. Numerical simulations illustrate the results.

Numerical Analysis
An Iterative Algorithm for Coupled Riccati Equations in Continuous-Time Markovian Jump Linear Systems

Ai-Guo Wu,Hui-Jie Sun,Wanquan Liu

DOI: https://doi.org/10.1080/00207721.2020.1800130

IF: 2.648

2020-01-01

International Journal of Systems Science

Abstract:In this paper, a novel implicit iterative algorithm with some tuning parameters is developed to solve the coupled algebraic Riccati matrix equation arising in the continuous-time Markovian jump linear systems. By introducing some tuning parameters in the proposed iterative algorithm, the current estimation for unknown variables is updated by using the information not only in the last step but also in the current iterative step and previous iterative steps. These tuning parameters can be appropriately chosen such that the proposed algorithm has faster convergence performance than some previous algorithms. It is shown that the proposed algorithm with zero initial conditions can monotonically converge to the unique positive semidefinite solution of the coupled Riccati matrix equation if the corresponding Markovian jump system is stabilisable. Finally, an example is provided to show the effectiveness of the developed algorithm.
Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

Changli Liu,Wei-Guo Wang,Jungong Xue,Ren-Cang Li

DOI: https://doi.org/10.1016/j.cam.2021.113614

IF: 2.872

2021-01-01

Journal of Computational and Applied Mathematics

Abstract:This paper is concerned with a M-matrix algebraic Riccati equation (MARE) XDX - AX - XB + C = O for which A is block-diagonal and its defining matrix W = [GRAPHICS] is a nonsingular or irreducible singular M-matrix. Such an MARE can be decomposec into many coupled algebraic Riccati equations (AREs) that can be solved by the Jacobior Gauss-Seidel-like iteration updating scheme at the outer-loop while by a doubling algorithm in the inner loop for each coupled ARE, as first proposed by Meini (2013). The goals of this paper are two-fold. One is to resolve a critical technical detail in Meini's algorithm that was not addressed. It is about whether each ARE in the inner loop has a minimal nonnegative solution. It is proved that the defining matrix of each coupled ARE during a doubling iteration is indeed a nonsingular or irreducible singular M-matrix and, as a result, they do have minimal nonnegative solutions and a doubling algorithm is an efficient way to compute them. The other goal is to design a highly accurate implementation of the doubling algorithm for the inner loop so that all entries of the minimal nonnegative solution to the original MARE are calculated with high entrywise relative accuracies, regardless of their magnitudes. This is made possible by a novel way of constructing triplet representations for the coupled ARE during doubling iterations. Numerical examples are presented to demonstrate that the resulting algorithm can indeed deliver an entrywise relatively accurate solution. (C) 2021 Elsevier B.V. All rights reserved.
A Hamiltonian approach to the factorization of the matrix Riccati equation

Faris A. Badawi,Anders Lindquist

DOI: https://doi.org/10.1007/BFb0120970

1982-01-01

Abstract:In this note the theory of Hamiltonian systems and an idea due to L. E. Zachrisson is used to obtain the factorizations of the matrix Riccati difference equation on which the derivations [8, 9, 11, 12, 13] of fast (non-Riccati) algorithms are based. Although we are primarily interested in discrete-time Riccati equations, for comparison, the corresponding, continuous-time result is briefly discussed.
Highly accurate decoupled doubling algorithm for large-scale M-matrix algebraic Riccati equations

Zhen-Chen Guo,Eric King-wah Chu,Xin Liang

2020-01-01

Abstract: We consider the numerical solution of large-scale M-matrix algebraic Riccati equations with low-rank structures. We derive a new doubling iteration, decoupling the four original iteration formulae in the alternating-directional doubling algorithm. We prove that the kernels in the decoupled algorithm are small M-matrices. Illumined by the highly accurate algorithm proposed by Xue and Li in 2017, we construct the triplet representations for the small M-matrix kernels in a highly accurate doubling algorithm. Illustrative numerical examples will be presented on the efficiency of our algorithm.
An Inversion-Free Iterative Algorithm with a Scalar Tuning Parameter for Coupled Riccati Matrix Equations Arising in LQ Optimal Control of Markov Jump Systems

Kaiwen Jiang,Zhi Li,Ying Zhang

DOI: https://doi.org/10.1109/tac.2024.3476993

IF: 6.549

2024-01-01

IEEE Transactions on Automatic Control

Abstract:In this paper, iterative algorithms are investigated to solve the Riccati algebraic matrix equations arising in the context of linear quadratic optimal control of discrete-time Markov jump systems. By using the matrix inversion lemma, an equivalent form for the considered coupled Riccati matrix equation is constructed, and then two intermediate variables are introduced to further construct a new equivalent form in order to avoid the operation of matrix inversion. With the aid of the obtained equivalent form, an inversion-free iterative algorithm is developed to obtain the positive definite solutions of the discrete coupled Riccati matrix equation via the principle of fixed-points. To improve the convergence performance, the latest updated information is exerted in the presented algorithm. Furthermore, the convergence analysis is conducted by mathematical induction for the obtained inversion-free iterative algorithm. At the end, numerical examples are given to demonstrate the effectiveness of the proposed algorithm by making comparisons with existing results.
Precise Transfer Matrix Method for Solving Differential Equation Systems

X. G. Zhang,Q. Zhang,J. P. Sun,T. Wang,Z. P. Song,J. J. Wang

DOI: https://doi.org/10.1063/1.5075644

2018-01-01

Abstract:The precise transfer matrix method is going to be used to solve a p system of linear ordinary differential equations. The complete problem is extended by adding index calculation involved. It is chosen an index calculation scheme to obtain transfer expressions. A recurrence strategy is used in these equations and permits one to relate different points in the domain where boundary equations are defined. Finally a p algebraic system of transfer matrix equations is noted and solved. To show the efficiency and accuracy, the method is applied to determine the structural behavior of a bending beam with different supports and to solve a differential equation of second degree with different boundary conditions.
On the explicit Hermitian solutions of the continuous‐time algebraic Riccati matrix equation for controllable systems

Liangyin Zhang,Michael Z. Q. Chen,Zhiwei Gao,Lifeng Ma

DOI: https://doi.org/10.1049/cth2.12618

IF: 2.67

2024-01-11

IET Control Theory and Applications

Abstract:The explicit solutions of the algebraic Riccati equation for systems in controllable canonical form are obtained using the entries of the system matrix, the closed‐loop system matrix, and the weighting matrix. The unknown entries of the closed‐loop system matrix are solved by scalar quadratic equations. If the system is near singularity, the explicit solution has higher numerical precision compared with the solution computed by numerical algorithms. This paper proposes explicit solutions for the algebraic Riccati matrix equation. For single‐input systems in controllable canonical form, the explicit Hermitian solutions of the non‐homogeneous Riccati equation are obtained using the entries of the system matrix, the closed‐loop system matrix, and the weighting matrix. The unknown entries of the closed‐loop system matrix are solved by scalar quadratic equations. For a homogeneous Riccati equation with a zero weighting matrix, the explicit solutions are proposed analytically in terms of the system eigenvalues. The advantages of the explicit solutions are threefold: first, if the system is controllable, the solution is directly given and the invariant subspaces of the Hamiltonian matrix are not required; second, if the system is near singularity, the explicit solution has higher numerical precision compared with the solution computed by numerical algorithms; third, for a real system in the controllable canonical form, the non‐negativity can be analysed for the explicit almost stabilizing solution.

automation & control systems,engineering, electrical & electronic,instruments & instrumentation
Efficient computation of the sinc matrix function for the integration of second-order differential equations

Lidia Aceto,Fabio Durastante

2024-06-12

Abstract:This work deals with the numerical solution of systems of oscillatory second-order differential equations which often arise from the semi-discretization in space of partial differential equations. Since these differential equations exhibit (pronounced or highly) oscillatory behavior, standard numerical methods are known to perform poorly. Our approach consists in directly discretizing the problem by means of Gautschi-type integrators based on $\operatorname{sinc}$ matrix functions. The novelty contained here is that of using a suitable rational approximation formula for the $\operatorname{sinc}$ matrix function to apply a rational Krylov-like approximation method with suitable choices of poles. In particular, we discuss the application of the whole strategy to a finite element discretization of the wave equation.

Numerical Analysis

A Precise Integration Algorithm for Matrix Riccati Differential Equations

A Note on the Necessary Conditions for the Algebraic Riccati Equation.

A Memory-Efficient Formulation of Precise-Integration Time-Domain Method With Riccati Matrix Differential Equations

An Implicit Iterative Algorithm for Algebraic Riccati Matrix Equations

An Inversion-Free Iterative Algorithm for Riccati Matrix Equations in Discrete-Time Markov Jump Systems

An implicit iterative algorithm for generalized algebraic Riccati matrix equations

A Novel Iterative Algorithm for Solving Coupled Riccati Equations

Accurate solutions of M -matrix algebraic Riccati equations

Iterative Algorithms with the Latest Update for Riccati Matrix Equations in Itô Markov Jump Systems

Backward differentiation formula method and random forest method to solve continuous‐time differential Riccati equations

The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

Low-rank plus diagonal approximations for Riccati-like matrix differential equations

An Iterative Algorithm for Coupled Riccati Equations in Continuous-Time Markovian Jump Linear Systems

Accurate Numerical Solution for Structured M-matrix Algebraic Riccati Equations

A Hamiltonian approach to the factorization of the matrix Riccati equation

Highly accurate decoupled doubling algorithm for large-scale M-matrix algebraic Riccati equations

An Inversion-Free Iterative Algorithm with a Scalar Tuning Parameter for Coupled Riccati Matrix Equations Arising in LQ Optimal Control of Markov Jump Systems

Precise Transfer Matrix Method for Solving Differential Equation Systems

On the explicit Hermitian solutions of the continuous‐time algebraic Riccati matrix equation for controllable systems

Efficient computation of the sinc matrix function for the integration of second-order differential equations