Optimal identification of discrete multivariable systems from noisy impulse response data