用armax()在hydromad包里面
armax.sim(U, a_1 = 0, a_2 = 0, a_3 = 0,
b_0 = 1, b_1 = 0, b_2 = 0, b_3 = 0,
pars = NULL,
delay = 0, init = 0, na.action = na.pass,
epsilon = hydromad.getOption("sim.epsilon"),
return_components = FALSE)
ssg.armax(theta)
normalise.armax(theta) 还有更详细的,可以看这篇文章:http://hydromad.catchment.org/man/armax.html