Abstract

It is known that the solution of the semilinear matrix equation $$X - A\overline X B = C$$ can be reduced to solving the classical Stein equation. The normal case means that the coefficients on the left-hand side of the resulting equation are normal matrices. We propose a method for solving the original semilinear equation in the normal case that permits to almost halve the execution time for equations of order n = 3000 compared to the library function dlyap, which solves Stein equations in Matlab.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call