Least squares and stochastic gradient parameter estimation for multivariable nonlinear Box‐Jenkins models based on the auxiliary model and the multi‐innovation identification theory