model<-auto.arima(AQI,trace=T)model
就是 season 季节的啊,你不都是Arima(1, 0, 0)(1, 0, 0)乘积季节模型了嘛,Arima(1, 0, 0)就只是 ar1;(1-sar1B^12)(1-ar1B)X_t = Z_t,均值你自己加吧,就是这样