Octave TSA: как имитировать процесс AR(p) из параметров `rmle()`
Я использую Octave с пакетом TSA. У меня есть тест временных рядов (синусоида) и запускаю функцию rmle Я могу получить AR(4) (порядок 4 определяется из графика автокорреляции) параметры:
t = [0:0.1:2*pi];
x = sin(t);
pkg load tsa;
[a4, var4] = rmle(x', 4);
a4 =
1.00000 -3.88649 5.76185 -3.86093 0.98696
var4 = 1.2466e-07
Я хотел бы рассчитать смоделированные значения x, как в:
xar4 = f (x, a4, var4)
и затем нарисуйте результат. Что было бы лучшим способом достичь того, чего я хочу?