isStable

Класс: LagOp

Определите устойчивость полинома оператора задержки

Синтаксис

[indicator,eigenvalues] = isStable(A)

Описание

[indicator,eigenvalues] = isStable(A) берет объект A полинома оператора задержки и проверяет, стабильно ли это. Условие устойчивости требует, чтобы значения всех корней характеристического полинома были меньше чем 1 к в маленьком числовом допуске.

Входные параметры

`A`	Изолируйте объект полинома оператора, как произведено `LagOp`.

Выходные аргументы

indicator

Булево значение для теста устойчивости. true указывает, что A(L) стабилен и что значение всех собственных значений его характеристического полинома - меньше чем один; false указывает, что A(L) нестабилен и что значение по крайней мере одного из собственных значений его характеристического полинома больше, чем или равно одному.

eigenvalues

Собственные значения характеристического полинома сопоставлены с A(L). Длина eigenvalues является продуктом степени и размерностью A(L).

Примеры

развернуть все

Проверяйте полином оператора задержки на устойчивость

Скрипт Open Live Script

Разделите два объекта полинома Оператора Задержки и проверку, если получившийся полином стабилен:

A = LagOp({1 -0.6 0.08});
B = LagOp({1 -0.5});
[indicator,eigenvalues]=isStable(A\B)

indicator = logical
   1

eigenvalues = 4×1 complex

   0.3531 + 0.0000i
  -0.0723 + 0.3003i
  -0.0723 - 0.3003i
  -0.3086 + 0.0000i

Советы

Полиномы нулевой степени всегда стабильны.
Для полиномов степени, больше, чем нуль, присутствие NaN-ценных коэффициентов возвращает индикатор устойчивости false и вектор NaN s в eigenvalues.
При тестировании на устойчивость сравнение включает маленький числовой допуск. Индикатором является true, когда значения всех собственных значений являются меньше, чем 1-10*eps, где eps является точностью машины. Пользователи, которые хотят включить их собственный допуск (включая 0) могут просто проигнорировать indicator и определить устойчивость можно следующим образом:
```
[~,eigenvalues] = isStable(A);
indicator = all(abs(eigenvalues) < (1-tol));
```
для некоторого маленького, неотрицательного допуска tol.

Ссылки

[1] Гамильтон, J. D. Анализ timeseries. Принстон, NJ: Издательство Принстонского университета, 1994.

Смотрите также

mldivide | mrdivide

Темы

Документация Econometrics Toolbox

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.