Interpolation by polynomials with symmetries on the imaginary axis

We here specialize the standard matrix-valued polynomial interpolation to the case where on the imaginary axis the interpolating polynomials admit various symmetries: Positive semidefinite, Skew-Hermitian, $J$-Hermitian, Hamiltonian and others. The procedure is comprized of three stages, illustrated through the case where on $i\R$ the interpolating polynomials are to be positive semidefinite. We first, on the expense of doubling the degree, obtain a minimal degree interpolating polynomial $P(s)$ which on $i\R$ is Hermitian. Then we find all polynomials $\Psi(s)$, vanishing at the interpolation points which are positive semidefinite on $i\R$. Finally, using the fact that the set of positive semidefinite matrices is a convex subcone of Hermitian matrices, one can compute the minimal scalar $\hat{\beta}\geq 0$ so that $P(s)+\beta\Psi(s)$ satisfies all interpolation constraints for all $\beta\geq\hat{\beta}$. This approach is then adapted to cases when the family of interpolating polynomials is not convex. Whenever convex, we parameterize all minimal degree interpolating polynomials.