StateSpaceTransform

StateSpaceTransform[sys, {p, q}]

StateSpaceTransform[sys, {{x1 -> p1[z], …}, {z1 -> q1[x], …}}]