Open Access
Translator Disclaimer
2014 Riemannian Gradient Algorithm for the Numerical Solution of Linear Matrix Equations
Xiaomin Duan, Huafei Sun, Xinyu Zhao
J. Appl. Math. 2014: 1-7 (2014). DOI: 10.1155/2014/507175


A Riemannian gradient algorithm based on geometric structures of a manifold consisting of all positive definite matrices is proposed to calculate the numerical solution of the linear matrix equation Q=X+i=1mAiTXAi. In this algorithm, the geodesic distance on the curved Riemannian manifold is taken as an objective function and the geodesic curve is treated as the convergence path. Also the optimal variable step sizes corresponding to the minimum value of the objective function are provided in order to improve the convergence speed. Furthermore, the convergence speed of the Riemannian gradient algorithm is compared with that of the traditional conjugate gradient method in two simulation examples. It is found that the convergence speed of the provided algorithm is faster than that of the conjugate gradient method.


Download Citation

Xiaomin Duan. Huafei Sun. Xinyu Zhao. "Riemannian Gradient Algorithm for the Numerical Solution of Linear Matrix Equations." J. Appl. Math. 2014 1 - 7, 2014.


Published: 2014
First available in Project Euclid: 26 March 2014

zbMATH: 07010658
MathSciNet: MR3166768
Digital Object Identifier: 10.1155/2014/507175

Rights: Copyright © 2014 Hindawi


Vol.2014 • 2014
Back to Top