exponenta event banner

Анализ моделей временных рядов

В этом примере показано, как анализировать модели временных рядов.

Модель временных рядов не имеет входных данных. Однако в таких моделях можно использовать множество команд вычисления отклика. Программное обеспечение обрабатывает (неявно) источник шума e(t) в качестве измеренного входного сигнала. Таким образом, step(sys) строит график ответа шага, предполагая, что входной сигнал шага был применен к каналу шума e(t).

Чтобы избежать неоднозначности в обращении программного обеспечения с моделью временных рядов, ее можно явно преобразовать в модель ввода-вывода с помощью noise2meas. Эта команда вызывает шум на входе e(t) обрабатывается как измеренный входной сигнал и преобразует линейную модель временных рядов в Ny выходные данные в модель «вход-выход» с Ny результаты и Ny входные данные. Результирующую модель можно использовать с такими командами, как, bode, nyquist, и iopzmap для изучения характеристик H передаточная функция.

Оценка модели временных рядов.

load iddata9
sys = ar(z9,4);

Преобразование модели временных рядов в модель ввода-вывода.

iosys = noise2meas(sys);

Постройте график ответа на шаг H.

step(iosys);

Figure contains an axes. The axes with title From: e@y1 To: y1 contains an object of type line. This object represents iosys.

Постройте график полюсов и нулей H.

iopzmap(iosys);

Figure contains an axes. The axes with title From: e@y1 To: y1 contains 2 objects of type line. This object represents iosys.

Вычислите и постройте график спектра временных рядов непосредственно без преобразования в модель «вход-выход».

spectrum(sys);

Figure contains an axes. The axes with title From: e@y1 To: y1 contains an object of type line. This object represents sys.

Команда готовит амплитуду спектра временного ряда Φ (ω) = ‖H (ω) ‖2.

См. также

|

Связанные темы