We study equations of the form r(X) = A, where r is a rational function and A and X are square matrices of the same size. We develop two techniques for solving these equations by inverting (through a substitution strategy) two schemes for the evaluation of rational functions of matrices. For triangular matrices, the new methods yield the same computational cost as the evaluation schemes from which they are obtained. A general equation can be reduced to upper triangular form by exploiting the Schur decomposition of the given matrix. For real data, the algorithms rely on the real Schur decomposition in order to compute real solutions using only real arithmetic. Numerical experiments show that our implementations are faster than existing alternatives without sacrificing accuracy.

Substitution algorithms for rational matrix equations

Iannazzo, Bruno;Fasi, Massimiliano
2020

Abstract

We study equations of the form r(X) = A, where r is a rational function and A and X are square matrices of the same size. We develop two techniques for solving these equations by inverting (through a substitution strategy) two schemes for the evaluation of rational functions of matrices. For triangular matrices, the new methods yield the same computational cost as the evaluation schemes from which they are obtained. A general equation can be reduced to upper triangular form by exploiting the Schur decomposition of the given matrix. For real data, the algorithms rely on the real Schur decomposition in order to compute real solutions using only real arithmetic. Numerical experiments show that our implementations are faster than existing alternatives without sacrificing accuracy.
2020
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11391/1452384
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 2
  • ???jsp.display-item.citation.isi??? 2
social impact