Идентификация неустойчивых систем

В этом небольшом примере мы строим очень простую неустойчивую линейную систему:

где  — управляющий вход, а  — возмущение. Сначала мы моделируем систему со стабилизирующей обратной связью и сохраняем данные. Затем мы снова моделируем ту же систему с сохраненными данными в качестве входных, но с немного иным начальным условием. Интуитивно может казаться, что эти две модели должны быть одинаковыми: система одна и та же и входные данные одни и те же, лишь начальное условие немного поменялось.

Результат показывает, что сначала вторая модель выглядит идентично, но затем начинает отклоняться по экспоненте. Чтобы понять причину, рассмотрим разницу между двумя моделями подробнее. В первом случае моделировалась система с замкнутым контуром.

Эта система является экспоненциально устойчивой, пока больше 1. Если мы собираем экспериментальные данные для неустойчивой системы со стабилизирующей обратной связью, мы по сути собираем данные из устойчивой системы.

Во втором случае моделировалась исходная система с разомкнутым контуром:

где данные управляющего входа поступали от контроллера стабилизации. Так почему же эта вторая модель оказалась неустойчивой? Чтобы понять это, посмотрим, что происходит при внесении в начальное условие небольшого возмущения . Разница между возмущенной и невозмущенной моделями выглядит так:

Это экспоненциально неустойчивая система, и, если не зависит от , разница будет возрастать по экспоненте.

Теперь понятно, почему происходит отклонение второй модели. Входные данные были собраны с контроллера стабилизации, воздействующего на , но небольшое отклонение состояния после сбора данных не корректируется контроллером и модель быстро отклоняется.

Как можно исправить эту ситуацию к лучшему? При чистом моделировании особо ничего не поделаешь. Однако при моделировании системы в целях оценки ее параметров есть очень эффективный прием: Решением является введение в модель также и обратной связи по измерению. Иными словами, вместо моделирования неустойчивой системы мы производим оценку ее состояния с использованием доступных данных измерений. Средство оценки состояния, например фильтр Калмана, имеет внутреннюю обратную связь по измерению, которая корректирует отклонения между состоянием системы и данными измерений. Для установившегося фильтра Калмана в дискретном времени эта коррекция выглядит так:

То есть мы производим линейную коррекцию на основе разницы между прогнозируемым выходом и измеренным выходом . Подобная коррекция должна быть вам знакома — сравните ее с уравнением замкнутой системы управления с опорным узлом :

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

Чтобы оценить параметры динамической системы с помощью средства оценки состояния, можно запустить средство оценки состояния, например фильтр Калмана, вперед во времени, сравнить оцененные состояния с измеренными данными и вычислить градиент погрешности относительно параметров. Именно так работает метод на основе погрешности измерения (PEM), реализованный в функции newpem. У этой функции есть именованный аргумент focus со значением по умолчанию focus = :prediction, которое, однако, можно изменить на focus = :simulation, чтобы отключить обратную связь по измерению. Однако для успешной идентификации неустойчивых систем обратная связь по измерению имеет ключевое значение и любой алгоритм моделирования без такой связи имеет тенденцию к отклонению.

Алгоритм PEM применяется для идентификации модели неустойчивой системы шара и балки в руководстве Шар и балка. Ближе к концу этого руководства эффективность оцениваемых моделей сравнивается с использованием все более длительных горизонтов прогнозирования, но уже при горизонте в 20 шагов эффективность прогнозирования существенно снижается вследствие неустойчивости.

У метода на основе погрешности прогнозирования есть и еще ряд привлекательных особенностей, о которых будет рассказано ниже.

По сути, метод PEM преобразует задачу минимизации потерь в плане производительности моделирования в задачу минимизации потерь в плане длительности прогнозирования ^{[1] [2]}. Это дает ряд преимуществ, два из которых показаны в следующем примере:

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

Начнем с распространенной проблемы, возникающей при минимизации погрешности моделирования. Допустим, нужно оценить неизвестную длину маятника. Из-за небольшой погрешности в длине маятника изменяется частота колебаний. При достаточно большом горизонте два синусоидальных сигнала разной частоты становятся почти ортогональными. При использовании потерь квадратической ошибки в той или иной форме характер потерь будет в данном случае крайне невыпуклым, что и иллюстрируется ниже.

Другой случай, представляющий сложность для оценки погрешности моделирования, — неустойчивая или хаотичная система. Небольшая погрешность в начальном условии или параметрах может привести к расхождению погрешности моделирования и бессмысленности его градиента.

В обоих этих примерах можно использовать имеющиеся результаты измерения эволюции системы, чтобы предотвратить расхождение погрешности моделирования. Например, если измерить угол отклонения маятника, то с помощью результата этого измерения можно скорректировать угол во время моделирования так, чтобы он оставался близким к измеренному углу. Вместо чистого моделирования мы, так сказать, прогнозируем состояние на некоторое время вперед с учетом всех результатов измерений до текущего момента времени. Минимизируя погрешность прогнозирования вместо чистого моделирования, часто можно предотвратить расхождение погрешности модели даже при плохом начальном предположении.

Начнем с определения модели маятника. Модель принимает параметр , соответствующий длине маятника.

Мы предполагаем, что истинная длина маятника равна , и генерируем данные для этой системы.

Мы также определяем функции, которые моделируют систему и вычисляют потери для заданного параметра p, соответствующего длине.

Далее характер потерь рассматривается как функция от длины маятника:

Этот рисунок вызывает интерес: потери при истинном значении , конечно, равны 0, но для значений наклон в целом неверный! Более того, потери имеют колебательный характер, а это значит, что оптимизировать данную функцию крайне непросто и для схождения локального поиска к истинному значению потребуется очень хорошее начальное предположение. Обратите внимание: мы намеренно выбрали одномерный пример для наглядности, и, хотя одномерные задачи обычно решаются легко, тот же ход рассуждений применим и к более сложным задачам высокой размерности.

Теперь перейдем к определению модели предиктора. Наш предиктор будет очень простым: на каждом временном шаге будет вычисляться погрешность между смоделированным углом и измеренным углом . Доля этой погрешности будет использоваться для коррекции состояния маятника. Используемая коррекция является линейной и выглядит так: . Мы построили то, что обычно называется (линейным) наблюдателем. Фильтр Калмана — это линейный наблюдатель особого рода, где рассчитывается на основе статистической модели возмущений, действующих на систему. Чтобы не усложнять задачу, мы ограничимся простым наблюдателем с постоянным усилением.

Для подачи выборки данных в непрерывную модель мы используем интерполятор. Мы также определим новые функции: функцию predictor, которая содержит динамику маятника с коррекцией наблюдателя, функцию prediction, которая выполняет развертывание (во избежание путаницы мы не используем слово «моделирование»), и функцию потерь.

Потери снова рассматриваются как функция от параметра, но на этот раз ситуация выглядит гораздо лучше. Потери не являются выпуклыми, но градиент имеет правильное направление на гораздо большем интервале. Здесь мы произвольно задали коэффициент усиления наблюдателя . Однако при оценке методом PEM оптимизатор обычно определяет этот параметр самостоятельно, что и происходит внутри newpem.

Может возникнуть вопрос, почему мы использовали коррекцию в форме , а не задали угол в модели равным результату измерения. Причин две.

В этом примере был продемонстрирован базовый вариант использования метода на основе погрешности прогнозирования для оценки параметров. В нем имеющийся результат измерения напрямую соответствовал одному из состояний и построить подходящий наблюдатель или предиктор было нетрудно. В более сложных случаях можно прибегнуть к нелинейному наблюдателю, например расширенному фильтру Калмана (EKF), или спроектировать фильтр Калмана на основе линеаризации системы вокруг некоторой рабочей точки.

Ссылки: