This paper considers the identification problem for multi-input multi-output nonlinear systems. The difficulty of the parameter identification of such systems is that the information vector in the identification model contains unknown variables. The solution is using the auxiliary model identification idea to overcome the difficulty. An auxiliary model based multi-innovation extended stochastic gradient algorithm is presented by expanding the innovation vector to an innovation matrix. The proposed algorithm uses not only the current innovation but also the past innovations at each recursion and thus the parameter estimation accuracy can be improved. The numerical example shows that the proposed algorithm is effective.