In this paper we consider Rellich's diagonalization theorem for analytic self-adjoint operator functions and investigate variational principles for their eigenfunctions and interlacing statements. As an application, we present a characterization for the eigenvalues of hyperbolic operator polynomials.