The Matrix Square Root from a New Functional Perspective: Theoretical Results and Computational Issues

We give a new characterization of the matrix square root and a new algorithm for its computation. We show how the matrix square root is related to the constant block coefficient of the inverse of a suitable matrix Laurent polynomial. This fact, besides giving a new interpretation of the matrix square root, allows one to design an efficient algorithm for its computation. The algorithm, which is mathematically equivalent to Newton's method, is quadratically convergent and numerically insensitive to the ill-conditioning of the original matrix and works also in the special case where the original matrix is singular and has a square root.

[1]  Alan J. Laub,et al.  On the Iterative Solution of a Class of Nonsymmetric Algebraic Riccati Equations , 2000, SIAM J. Matrix Anal. Appl..

[2]  C. R. DePrima,et al.  The range of A−1A∗ in GL(n,C) , 1974 .

[3]  Nicholas J. Higham,et al.  Stable iterations for the matrix square root , 1997, Numerical Algorithms.

[4]  Dario Bini,et al.  Computations with infinite Toeplitz matrices and polynomials , 2002 .

[5]  Alan J. Laub,et al.  A hyperbolic tangent identity and the geometry of Padé sign function iterations , 2005, Numerical Algorithms.

[6]  C. Pearce,et al.  On the square-root method for continuous-time algebraic Riccati equations , 1999 .

[7]  P. Laasonen.  On the iterative solution of the matrix equation ²-=0 , 1958 .

[8]  E. Denman,et al.  The matrix sign function and computations in systems , 1976 .

[9]  N. Higham.  Computing the polar decomposition with applications , 1986 .

[10]  A. Laub,et al.  Rational iterative methods for the matrix sign function , 1991 .

[11]  Linzhang Lu,et al.  PROPERTIES OF A QUADRATIC MATRIX EQUATION AND THE SOLUTION OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION , 1995 .

[12]  B. Iannazzo,et al.  A note on computing the matrix square root , 2003 .

[13]  Å. Björck,et al.  A Schur method for the square root of a matrix , 1983 .

[14]  Charles R. Johnson,et al.  Topics in Matrix Analysis , 1991 .

[15]  N. Higham.  Computing real square roots of a real matrix , 1987 .

[16]  C. Loan.  A Symplectic Method for Approximating All the Eigenvalues of a Hamiltonian Matrix , 1982 .

[17]  Beatrice Meini,et al.  On the Solution of a Nonlinear Matrix Equation Arising in Queueing Problems , 1996, SIAM J. Matrix Anal. Appl..

[18]  N. Higham.  Newton's method for the matrix square root , 1986 .

[19]  N. Higham.  The matrix sign decomposition and its relation to the polar decomposition , 1994 .