×

A fast Newton’s method for a nonsymmetric algebraic Riccati equation. (English) Zbl 1166.15300

Summary: A special instance of the algebraic Riccati equation \(XCX-XE-AX+B=0\) where the \(n\times n\) matrix coefficients \(A,B,C,E\) are rank structured matrices is considered. Relying on the structural properties of Cauchy-like matrices, an algorithm is designed for performing the customary Newton iteration in \(O(n^2)\) arithmetic operations (ops). The same technique is used to reduce the cost of the algorithm proposed by L.-Z. Lu in [Numer. Linear Algebra Appl. 12, No. 2–3, 191–200 (2005; Zbl 1164.65386)] from \(O(n^3)\) to \(O(n^2)\) ops while still preserving quadratic convergence in the generic case. As a byproduct we show that the latter algorithm is closely related to the customary Newton method by simple formal relations. In critical cases where the Jacobian at the required solution is singular and quadratic convergence turns to linear, we provide an adaptation of the shift technique in order to get rid of the singularity. The original equation is transformed into an equivalent Riccati equation where the singularity is removed while the matrix coefficients maintain the same structure as in the original equation. This leads to a quadratically convergent algorithm with complexity \(O(n^2)\) which provides approximations with full precision. Numerical experiments and comparisons which confirm the effectiveness of the new approach are reported.

MSC:

15A24 Matrix equations and identities
65F05 Direct numerical methods for linear systems and matrix inversion
65H10 Numerical computation of solutions to systems of equations

Citations:

Zbl 1164.65386
Full Text: DOI