Abstract and Applied Analysis

A Numerical Method for Computing the Principal Square Root of a Matrix

Abstract

It is shown how the mid-point iterative method with cubical rate of convergence can be applied for finding the principal matrix square root. Using an identity between matrix sign function and matrix square root, we construct a variant of mid-point method which is asymptotically stable in the neighborhood of the solution. Finally, application of the presented approach is illustrated in solving a matrix differential equation.

Article information

Source
Abstr. Appl. Anal., Volume 2014, Special Issue (2014), Article ID 525087, 7 pages.

Dates
First available in Project Euclid: 6 October 2014

https://projecteuclid.org/euclid.aaa/1412606413

Digital Object Identifier
doi:10.1155/2014/525087

Mathematical Reviews number (MathSciNet)
MR3256251

Zentralblatt MATH identifier
07022554

Citation

Soleymani, F.; Shateyi, S.; Khaksar Haghani, F. A Numerical Method for Computing the Principal Square Root of a Matrix. Abstr. Appl. Anal. 2014, Special Issue (2014), Article ID 525087, 7 pages. doi:10.1155/2014/525087. https://projecteuclid.org/euclid.aaa/1412606413

References

• N. J. Higham, Functions of Matrices: Theory and Computation, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 2008.
• C. Kanzow and C. Nagel, “Semidefinite programs: new search directions, smoothing-type methods, and numerical results,” SIAM Journal on Optimization, vol. 13, no. 1, pp. 1–23, 2002.
• G. W. Cross and P. Lancaster, “Square roots of complex matrices,” Linear and Multilinear Algebra, vol. 1, pp. 289–293, 19774.
• N. J. Higham, “Newton's method for the matrix square root,” Mathematics of Computation, vol. 46, no. 174, pp. 537–549, 1986.
• Å. Björck and S. Hammarling, “A Schur method for the square root of a matrix,” Linear Algebra and its Applications, vol. 52/53, pp. 127–140, 1983.
• E. D. Denman and J. Beavers, “The matrix sign function and computations in systems,” Applied Mathematics and Computation, vol. 2, no. 1, pp. 63–94, 1976.
• J. M. Gutierrez, M. A. Hernandez, and M. A. Salanova, “Calculus of nth roots and third order iterative methods,” Nonlinear Analysis: Theory, Methods & Applications, vol. 47, no. 4, pp. 2875–2880, 2001.
• B. Iannazzo, “A note on computing the matrix square root,” Calcolo, vol. 40, no. 4, pp. 273–283, 2003.
• M. Frontini and E. Sormani, “Some variant of Newton's method with third-order convergence,” Applied Mathematics and Computation, vol. 140, no. 2-3, pp. 419–426, 2003.
• S. Lakić and M. S. Petković, “On the matrix square root,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 78, no. 3, pp. 173–182, 1998.
• F. A. Potra, “On $Q$-order and $R$-order of convergence,” Journal of Optimization Theory and Applications, vol. 63, no. 3, pp. 415–431, 1989.
• B. Iannazzo, “On the Newton method for the matrix $p$th root,” SIAM Journal on Matrix Analysis and Applications, vol. 28, no. 2, pp. 503–523, 2006.
• W. Gander, “On Halley's iteration method,” The American Mathematical Monthly, vol. 92, no. 2, pp. 131–134, 1985.
• C. H. Guo, “On Newton's method and Halley's method for the principal $p$th root of a matrix,” Linear Algebra and its Applications, vol. 432, no. 8, pp. 1905–1922, 2010.
• C. S. Kenney and A. J. Laub, “Rational iterative methods for the matrix sign function,” SIAM Journal on Matrix Analysis and Applications, vol. 12, no. 2, pp. 273–291, 1991.
• F. K. Haghani and F. Soleymani, “On a fourth-order matrix method for computing polar decomposition,” Computational and Applied Mathematics, 2014.
• F. Soleymani, E. Tohidi, S. Shateyi, and F. K. Haghani, “Some matrix iterations for computin g matrix sign function,” Journal of Applied Mathematics, vol. 2014, Article ID 425654, 9 pages, 2014.
• C. H. Guo and N. J. Higham, “A Schur-Newton method for the matrix $p$th root and its inverse,” SIAM Journal on Matrix Analysis and Applications, vol. 28, no. 3, pp. 788–804, 2006.
• J. Hoste, Mathematica Demystified, McGraw-Hill, New York, NY, USA, 2009.
• B. Meini, “The matrix square root from a new functional perspective: theoretical results and computational issues,” Tech. Rep. 1455, Dipartimento di Matematica, Universit di Pisa, Pisa, Italy, 2003. \endinput