Gradient-descent iterative algorithm for solving a class of linear matrix equations with applications to heat and Poisson equations

Advances in Difference Equations

article

Kittisopaporn, Adisorn¹  Chansangiam, Pattrawut¹ 

[1] Department of Mathematics, Faculty of Science, King Mongkut’s Institute of Technology Ladkrabang

关键词: Generalized Sylvester matrix equation;    Gradient descent;    Iterative method;    Matrix norms and conditioning;    Heat equation;    Poisson’s equation;   

DOI  :  10.1186/s13662-020-02785-9

学科分类：航空航天科学

来源: SpringerOpen

【 摘 要 】

In this paper, we introduce a new iterative algorithm for solving a generalized Sylvester matrix equation of the form$\sum_{t=1}^{p}A_{t}XB_{t}=C$ which includes a class of linear matrix equations. The objective of the algorithm is to minimize an error at each iteration by the idea of gradient-descent. We show that the proposed algorithm is widely applied to any problems with any initial matrices as long as such problem has a unique solution. The convergence rate and error estimates are given in terms of the condition number of the associated iteration matrix. Furthermore, we apply the proposed algorithm to sparse systems arising from discretizations of the one-dimensional heat equation and the two-dimensional Poisson’s equation. Numerical simulations illustrate the capability and effectiveness of the proposed algorithm comparing to well-known methods and recent methods.

【 授权许可】

CC BY