A numerically stable algorithm is derived to compute orthonormal bases for any deflating subspace of a regular pencil hB-A. The method is based on an update of the OZ-algorithm, in order to obtain any desired ordering of eigenvalues in the quasitriangular forms constructed by this algorithm. As applications we discuss a new approach to solve Riccati equations arising in linear system theory. The computation of deflating subspaces with specified spectrum is shown to be of crucial importance here.
We construct optimally robust port-Hamiltonian realizations of a given rational transfer function that represents a passive system. We show that the realization with a maximal passivity radius is a normalized port-Hamiltonian one. Its computation is linked to a particular solution of a linear matrix inequality that defines passivity of the transfer function, and we provide an algorithm to construct this optimal solution. We also consider the problem of finding the nearest passive system to a given non-passive one and provide a simple but suboptimal solution.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.