Numeričko rješavanje linearnih matričnih jednadžbi

Abstract

U ovom smo radu dali pregled postojećih rezultata u području numeričkog rješavanja linearnih matričnih jednadžbi na primjeru Sylvesterove i Lyapunovljeve jednadžbe. Naveli smo zašto su od interesa u primijenjenoj matematici, posebno u analizi stabilnosti dinamičkih sustava. Razmotrili smo cesto korišten Bartels–Stewartov algoritam koji je pogodan za rješavanje jednadžbi malih dimenzija. Za matrične jednadžbe velikih dimenzija razmotrili smo iterativne načine rješavanja: metodu ADI iteracija i projekcijsku Krylovljevu metodu. Pritom smo uzeli u obzir strukture matrica koje se često pojavljuju u primjenama. Implementaciju ovih algoritama smo napravili u Pythonu.In this work we reviewed main computational methods for solving linear matrix equations, particularly Sylvester and Lyapunov equations. We explained their role in applied mathematics, especially in stability analysis of linear dynamic systems. We described Bartels–Stewart algorithm which is frequently used for solving matrix equations of small dimensions. For large-scale matrix equations, we considered special structure of matrices which are often encountered in applications and we described iterative methods for finding solution: ADI method and projection (Krylov) method. We implemented listed algorithms using Python