Solution of the matrix eigenvalue problem VA+A*V=µV with applications to the study of free linear dynamical systems

  • Authors:
  • L. Kohaupt

  • Affiliations:
  • TFH Berlin, FB II, Luxemburger Str. 10, D-13353 Berlin, Germany

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2008

Quantified Score

Hi-index 7.29

Visualization

Abstract

The new idea is to study the stability behavior of the solution x=x(t) of the initial value problem x@?=Ax,t=t"0,x(t"0)=x"0, with A@?C^n^x^n, in a weighted (semi-) norm @?.@?"R where R is taken as an appropriate solution of the matrix eigenvalue problemRA+A^*R=@rR, rather than as the solution of the algebraic Lyapunov matrix equationRA+A^*R=-S with given positive (semi-) definite matrix S. Substantially better results are obtained by the new method. For example, if A is diagonalizable and all eigenvalues @l"i(A) have negative real parts, i.e., Re@l"i(A)0(t-~), which is much more than the old result, which only states that lim"t"-"~x(t)=0. Especially, the semi-norms @?.@?"R"""i have a decoupling and filter effect on x(t). Further, new two-sided bounds (depending on x"0) for the asymptotic behavior can be derived. The best constants in the bounds are obtained by the differential calculus of norms. Applications are made to free linear dynamical systems, and computations underpin the theoretical findings.