Abstract Numerical methods for developing port-Hamiltonian representations of general linear time-invariant systems are studied. The approach extends previous port-Hamiltonian characterizations to include the general non-minimal case and the case where the feedthrough term fails to have an invertible symmetric part. The resulting construction is able to identify infeasibility when the system fails to be port-Hamiltonian, and allows for the incorporation of perturbations in order to arrive at a nearby port-Hamiltonian system. Results are illustrated via numerical examples.
Beattie et al. (Wed,) studied this question.