Maximum likelihood gradient identification for multivariate equation‐error moving average systems using the multi‐innovation theory