3.6.2. Модель объекта управления и решение задачи оптимизации
Рассмотрим сначала алгоритм ОПУ для систем с одним входом и одним выходом. При выводе алгоритма ОПУ Кларк и др. ввели модель УАРПСС (управляемая авторегрессия — проинтегрированное скользящее среднее); эта модель является более подходящей для использования в управлении процессами с целью устранения смещения, вызываемого возмущениями. Эти возмущения можно смоделировать как дрейф (броуновское
движение) или как ступенчатую случайную величину, изменяющуюся в случайные моменты времени. Модель УАРПСС с временной задержкой
равной 1, может быть записана в следующей форме:
Возмущающее воздействие может быть как детерминированным, так и стохастическим, однако за счет оператора А его математическое ожидание принимается равным нулю. В подходе на основе ОПУ процесс
предполагается стохастическим, а полином
— устойчивым и равным 1.
Для решения задачи оптимизации (3.6.2) требуется вычислить прогноз на
шагов вперед для выходного процесса
при
Вычисления выполняются на основе информации, известной в момент времени
а также информации о будущих значениях инкремента управления, которые требуется выбирать таким образом, чтобы оптимизировать критерий (3.6.2). Прогноз выходного процесса на
шагов можно описать в следующей форме:
где
и
получены путем рекурсивного решения диофантова уравнения:
где
Вычитая уравнение (3.6.5) из уравнения (3.6.6), получим
Таким образом,
Полагая
получим
Из уравнения (3.6.7) получим:
Так как слагаемое
содержит члены
значит, первый
член выражения
нулевой. Отсюда следует, что
Таким образом, имеем
и
Приравнивая члены при одинаковых степенях
получим:
где индекс! изменяется от 0 до степени полинома.
Прогноз выходного значения можно разбить на две компоненты:
которая известна в момент
и
— компонента, соответствующая управляющему воздействию. Таким образом,
где
Здесь Г — нижняя треугольная
-матрица. Ее элементы представляют собой величины реакции объекта управления на соответствующих шагах.
Отличительная особенность алгоритма ОПУ состоит в том, что выбор параметров
, является гибким и не влияет на устойчивость управляемой системы. Однако правильный выбор величин
обеспечивает сокращение времени вычислений.
Квадратичной показатель качества (3.6.2) теперь можно представить в виде
где
— вектор заранее заданных уставок. Получение инкрементного вектора управления
для будущих моментов времени достигается посредством дифференцирования показателя качества по
Таким образом, получим:
Так как требуется только первое управляющее воздействие, имеем
где
— первая строка матрицы
Чтобы сделать алгоритм самонастраивающимся,
выполняется оценка параметров объекта управления на основе рекурсивного метода наименьших квадратов, рассмотренного в разделе 3.2. Векторы данных и параметров для модели объекта управления записываются в следующем виде: