Рекурсивные алгоритмы оценки в System Identification Toolbox™ могут быть разделены на две категории:
Алгоритмы истории Бога — Эти алгоритмы стремятся минимизировать ошибку между наблюдаемыми и предсказанными выходными сигналами, навсегда продвигается с начала симуляции. System Identification Toolbox поддерживает оценку бесконечной истории:
Рекурсивные средства оценки командной строки для линейной регрессии наименьших квадратов, AR, ARX, ARMA, ARMAX, OE и структур модели BJ
Simulink® Recursive Least Squares Estimator и блоки Recursive Polynomial Model Estimator
Алгоритмы конечной истории — Эти алгоритмы стремятся минимизировать ошибку между наблюдаемыми и предсказанными выходными сигналами для конечного числа прошлых временных шагов. Тулбокс поддерживает оценку конечной истории для линейных в параметрах моделей:
Рекурсивные средства оценки командной строки для линейной регрессии наименьших квадратов, AR, ARX и структур модели OE
Блок Simulink Recursive Least Squares Estimator
Блок Simulink Recursive Polynomial Model Estimator, для AR, ARX и структур OE только
Алгоритмы конечной истории обычно легче настроить, чем алгоритмы бесконечной истории, когда параметры имеют быстрые и потенциально большие изменения в зависимости от времени.
Общая форма бесконечной истории рекурсивный алгоритм оценки следующие:
оценка параметра во время t. y (t), наблюдаемый выходной сигнал во время t, и предсказание y (t) на основе наблюдений до времени t-1. Усиление, K (t), определяет сколько текущей ошибки предсказания влияет на обновление оценки параметра. Алгоритмы оценки минимизируют остаточный член предсказания .
Усиление имеет следующую форму:
Рекурсивные алгоритмы, поддержанные продуктом System Identification Toolbox, отличаются на основе разных подходов для выбора формы Q (t) и вычисление . Здесь, представляет градиент предсказанного выхода модели относительно параметров .
Самый простой способ визуализировать роль градиента из параметров, должен рассмотреть модели с формой линейной регрессии:
В этом уравнении, вектор регрессии, который вычисляется на основе предыдущих значений измеренных вводов и выводов. представляет истинные параметры. e (t) является источником шума (инновации), который принят, чтобы быть белым шумом. Определенная форма зависит от структуры полиномиальной модели.
Для уравнений линейной регрессии предсказанный выход дан следующим уравнением:
Для моделей, которые не имеют формы линейной регрессии, не возможно вычислить точно предсказанный выход и градиент для текущей оценки параметра . Изучить, как можно вычислить приближение для и для общих структур модели смотрите раздел по рекурсивным методам ошибки предсказания в [1].
Программное обеспечение System Identification Toolbox предоставляет следующей бесконечной истории рекурсивные алгоритмы оценки для онлайновой оценки:
Фактор упущения и формулировки Фильтра Калмана более в вычислительном отношении интенсивны, чем градиент и ненормированные градиентные методы. Однако у них обычно есть лучшие свойства сходимости.
Упущение Фактора. Следующая система уравнений обобщает забывающий факторный алгоритм адаптации:
Программное обеспечение гарантирует, что P(t) является положительно-определенной матрицей при помощи алгоритма квадратного корня, чтобы обновить его [2]. Программное обеспечение вычисляет P
при предположении, что остаточные значения (различие между предполагаемыми и измеренными выходными параметрами) являются белым шумом, и отклонение этих остаточных значений равняется 1. R2 /2
P
приблизительно равно ковариационной матрице предполагаемых параметров, где R2 является истинным отклонением остаточных значений.
Q (t) получен путем минимизации следующей функции во время t:
Смотрите раздел 11.2 в [1] для деталей.
Этот подход обесценивает старые измерения, экспоненциально таким образом, что наблюдение, которое является старые переносы выборок вес, который равен времена вес нового наблюдения. представляет горизонт памяти этого алгоритма. Измерения, более старые, чем обычно несите вес, который меньше приблизительно 0,3.
называется фактором упущения и обычно имеет положительное значение между 0.98
и 0.995
Набор оценить независимые от времени (постоянные) параметры. Набор оценить изменяющиеся во времени параметры.
Забывающий факторный алгоритм для = 1 эквивалентно алгоритму Фильтра Калмана с R1=0 и R2=1. Для получения дополнительной информации об алгоритме Фильтра Калмана, смотрите Фильтр Калмана.
Фильтр Калмана. Следующая система уравнений обобщает алгоритм адаптации Фильтра Калмана:
Программное обеспечение гарантирует, что P(t) является положительно-определенной матрицей при помощи алгоритма квадратного корня, чтобы обновить его [2]. Программное обеспечение вычисляет P
при предположении, что остаточные значения (различие между предполагаемыми и измеренными выходными параметрами) являются белым шумом, и отклонение этих остаточных значений равняется 1. R2* P
ковариационная матрица предполагаемых параметров, и R1/R2 является ковариационной матрицей изменений параметра. Где, R1 является ковариационной матрицей изменений параметра, которые вы задаете.
Эта формулировка принимает форму линейной регрессии модели:
Q (t) вычисляется с помощью Фильтра Калмана.
Эта формулировка также принимает что истинные параметры описаны случайным обходом:
w (t) является Гауссовым белым шумом со следующей ковариационной матрицей или матрицей R1 дрейфа :
R2 является отклонением инноваций e (t) в следующем уравнении:
Алгоритм Фильтра Калмана полностью задан последовательностью данных y (t), градиент , R1, R2 и начальные условия (исходное предположение параметров) и (ковариационная матрица, которая указывает на ошибки параметров).
Это принято, что R1 и P (t = 0) матрицы масштабируются таким образом что R2 = 1. Это масштабирование не влияет на оценки параметра.
Нормированный и Ненормированный Градиент. В случае линейной регрессии градиентные методы также известны как методы наименьшее количество средних квадратичных (LMS).
Следующая система уравнений обобщает ненормированный градиент и нормированный алгоритм адаптации градиента:
В ненормированном подходе градиента Q (t) дают:
В нормированном подходе градиента Q (t) дают:
Нормированный алгоритм градиента масштабирует усиление адаптации, γ, на каждом шаге квадратом 2D нормы вектора градиента. Если градиент близко к нулю, это может вызвать скачки в предполагаемых параметрах. Чтобы предотвратить эти скачки, термин смещения вводится в масштабном коэффициенте.
Этот выбор Q (t) для алгоритмов градиента обновляет параметры в отрицательном направлении градиента, где градиент вычисляется относительно параметров. См. стр. 372 в [1] для деталей.
Методы оценки конечной истории находят, что параметр оценивает θ (t) путем минимизации
где y (k) является наблюдаемым выходным сигналом во время k, и предсказанный выход во время k. Этот подход также известен как оценку раздвижного окна. Подходы оценки конечной истории минимизируют ошибки предсказания для последних временных шагов N. В отличие от этого методы оценки бесконечной истории минимизируют ошибки предсказания, запускающиеся с начала симуляции.
System Identification Toolbox поддерживает оценку конечной истории для линейных в параметрах моделей (AR, и ARX), где предсказано выведенный имеет форму . Программное обеспечение создает и обеспечивает буфер регрессоров ψ (k) и наблюдало выходные параметры y (k) для k = t-N+1, t-N+2, …, t-2, t-1, t. Эти буферы содержат необходимые матрицы для базовой проблемы линейной регрессии минимизации по θ. Программное обеспечение решает эту задачу линейной регрессии с помощью QR, учитывающего с поворотом столбца.
[1] Ljung, L. System Identification: теория для пользователя. Верхний Сэддл-Ривер, NJ: PTR Prentice Hall, 1999.
[2] Карлсон, N.A. "Быстро треугольная формулировка фильтра квадратного корня". Журнал AIAA, Издание 11, Номер 9, 1973, стр 1259-1265.
[3] Чжан, Q. "Некоторые Аспекты Реализации Алгоритмов Наименьших квадратов Раздвижного окна". Продолжения IFAC. Издание 33, Выпуск 15, 2000, стр 763-768.
Recursive Least Squares Estimator | Recursive Polynomial Model Estimator | recursiveAR
| recursiveARMA
| recursiveARMAX
| recursiveARX
| recursiveBJ
| recursiveLS
| recursiveOE