Этот пример показывает, как оценить авторегрессивные интегрированные модели скользящего среднего значения (ARIMA).
Иногда требуются модели временных рядов, содержащие нестационарные тренды (сезонность). Одной из категорий таких моделей являются модели ARIMA. Эти модели содержат фиксированный интегратор в источнике шума. Таким образом, если управляющее уравнение модели ARMA выражено как A (q) y (t) = Ce (t), где A (q) представляет авторегрессивный член, а C (q) - скользящий средний член, соответствующая модель модели ARIMA выражается как
где термин представляет интегратор в дискретном времени. Точно так же можно сформулировать уравнения для моделей ARI и ARIX.
Использование команд оценки модели timeseries ar
, arx
и armax
можно ввести интеграторы в источник шума e(t)
. Вы делаете это при помощи IntegrateNoise
параметр в команде оценки.
Подход оценки не учитывает никаких постоянных смещений в данных timeseries. Способность вводить интегратор шума не ограничивается одними только данными timeseries. Вы можете сделать это также для моделей ввода-вывода, где нарушения порядка могут быть подвержены сезонности. Одним из примеров являются полиномиальные модели структуры ARIMAX:
Смотрите armax
страница с описанием для примеров.
Оцените модель ARI для скалярного timeseries с линейным трендом.
load iddata9 z9 Ts = z9.Ts; y = cumsum(z9.y); model = ar(y,4,'ls','Ts',Ts,'IntegrateNoise', true); % 5 step ahead prediction compare(y,model,5)
Оцените многомерную модель timeseries таким образом, чтобы шумовое интегрирование присутствовало только в одном из двух временных рядов.
load iddata9 z9 Ts = z9.Ts; y = z9.y; y2 = cumsum(y); % artificially construct a bivariate time series data = iddata([y, y2],[],Ts); na = [4 0; 0 4]; nc = [2;1]; model1 = armax(data, [na nc], 'IntegrateNoise',[false; true]); % Forecast the time series 100 steps into future yf = forecast(model1,data(1:100), 100); plot(data(1:100),yf)
Если выходы были связаны (na
не была диагональной матрицей), ситуация будет сложнее и просто добавление интегратора во второй канал шума не сработает.