Algorithms for polynomial spectral factorization and bounded-real balanced state space representations

We illustrate an algorithm that starting from the image representation of a strictly bounded-real system computes a minimal balanced state variable, from which a minimal balanced state realization is readily obtained. The algorithm stems from an iterative procedure to compute a storage function, based on a technique to solve a generalization of the Nevanlinna interpolation problem