The aim of this paper is to design a band-limited optimal input with power constraints for identifying a linear multi-input multi-output system, with nominal parameter values specified. Using spectral decomposition theorem, the power spectrum is written as φu(jω) = 1/2H(jω)H*(jω). The matrix H(jω) is expressed in terms of a truncated basis for L2([-ωc, ωc]), where ωc is the cut-off frequency. The elements of the Fisher Information Matrix and the power constraints become homogeneous quadratics in basis coefficients. The optimality criterion used are D-optimality, A-optimality, T-optimality and ε-optimality. This optimization problem is not known to be convex. A bi-linear formulation gives a lower bound on the optimum, while an upper bound is obtained through a convex relaxation. These bounds can be computed efficiently. The lower bound is used as a suboptimal solution, its sub-optimality determined by the difference between the bounds. Simulations reveal that the bounds match in many instances, implying global optimality. © 2018, © 2018 Informa UK Limited, trading as Taylor & Francis Group.