On the numerical solution of matrix polynomial equations

From MaRDI portal





The authors present certain direct and iterative methods for solving matrix polynomial equations (appearing in the problem of parameter identification for a stationary linear stochastic system) of the form \[ AX + AX^2 + \dots + AX^n = C. \] Here \(X\in\mathbb R^{n\times n}\) is the matrix which defines the dynamics of a stochastic system, \(A\in\mathbb R^{m\times n}\) is the matrix of observations of this system, and \(C = \sum_{i=1}^n C_i\), where \(C_i\in\mathbb R^{m\times n}\). The authors prove that, for the problem under consideration, there exists an equivalent representation which admits an effective and numerically stable solution under different ways of describing the original problem. A local convergence theorem of iterative algorithms is given, and the restrictions involved by this theorem are discussed. In the case of the problem of the so-called total rank, the authors expose an algorithm for constructing the set of all solutions to the equation. Estimations of the convergence rate of the iterative methods are established and, moreover, certain recommendations for their numerical realizations are presented.











This page was built for publication: On the numerical solution of matrix polynomial equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2729357)