Моделирование MPC

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

Модель объекта управления

Можно задать модель объекта управления в одном из следующих форматов линейного независимого от времени (LTI):

Числовые модели LTI — Передаточная функция (tf), пространство состояний (ss), нули, полюса и усиление (zpk)
Идентифицированные модели (требует System Identification Toolbox™) — idss (System Identification Toolbox), idtf (System Identification Toolbox), idproc (System Identification Toolbox), и idpoly (System Identification Toolbox)

Контроллер MPC выполняет все вычисления оценки и оптимизации с помощью дискретного времени, без задержек, система в пространстве состояний с безразмерными переменными ввода и вывода. Поэтому, когда вы задаете модель объекта управления в контроллере MPC, программное обеспечение выполняет следующее в случае необходимости:

Преобразование в пространство состояний — ss команда преобразует предоставленную модель в модель в пространстве состояний LTI.
Дискретизация или передискретизирующий — Если шаг расчета модели отличается от диспетчера MPC шаг расчета (заданный в Ts свойство), одно из следующего происходит:
- Если модель является непрерывным временем, c2d команда преобразует его в объект LTI дискретного времени с помощью шага расчета контроллера.
- Если модель является дискретным временем, d2d команда передискретизирует его, чтобы сгенерировать объект LTI дискретного времени с помощью шага расчета контроллера.
Задержите удаление — Если модель дискретного времени включает какой-либо вход, выход или внутренние задержки, absorbDelay команда заменяет их на соответствующее количество полюсов в z = 0, увеличивая общее число дискретных состояний. InputDelay, OutputDelay, и InternalDelay свойства получившейся модели в пространстве состояний являются всем нулем.
Преобразование в безразмерные переменные ввода и вывода — контроллер MPC позволяет вам задать масштабный коэффициент для каждой переменной ввода и вывода объекта. Если вы не задаете масштабные коэффициенты, они принимают значение по умолчанию к 1. Программное обеспечение преобразует переменные ввода и вывода объекта в безразмерную форму можно следующим образом:
$\begin{matrix} x_{p} (k + 1) = A_{p} x_{p} (k) + B S_{i} u_{p} (k) \\ y_{p} (k) = S_{o}^{- 1} C x_{p} (k) + S_{o}^{- 1} D S_{i} u_{p} (k) . \end{matrix}$
где _Ap, B, C и D являются постоянными матрицами пространства состояний нулевой задержки от шага 3, и:
- _Si является диагональной матрицей входных масштабных коэффициентов в технических модулях.
- _So является диагональной матрицей выходных масштабных коэффициентов в технических модулях.
- _xp является вектором состояния от шага 3 в технических модулях (включая любые поглощенные состояния задержки). Никакое масштабирование не выполняется на переменных состояния.
- _up является вектором из безразмерных входных переменных объекта, включая переменные, которыми управляют, измеренные воздействия и неизмеренные входные воздействия.
- _yp является вектором из безразмерных выходных переменных объекта.
Получившаяся модель объекта управления имеет следующую эквивалентную форму:
$\begin{matrix} x_{p} (k + 1) = A_{p} x_{p} (k) + B_{p u} u (k) + B_{p v} v (k) + B_{p d} d (k) \\ y_{p} (k) = C_{p} x_{p} (k) + D_{p u} u (k) + D_{p v} v (k) + D_{p d} d (k) . \end{matrix}$
Здесь, $C_{p} = S_{o}^{- 1} C$ , _Bpu, _Bpv и _Bpd являются соответствующими столбцами _BSi. Кроме того, _Dpu, _Dpv и _Dpd являются соответствующими столбцами $S_{o}^{- 1} D S_{i}$ . Наконец, u (k), v (k) и d (k) является безразмерными переменными, которыми управляют, измеренными воздействиями и неизмеренными входными воздействиями, соответственно.
Контроллер MPC осуществляет ограничение _Dpu = 0, что означает, что контроллер не позволяет прямое сквозное соединение ни от какой переменной, которой управляют, ни до какого объекта выход.

Введите возмущение

Если ваша модель объекта управления включает неизмеренные входные воздействия, d (k), входное возмущение задает тип сигнала и характеристики d (k). Смотрите Оценку состояния Контроллера для получения дополнительной информации о модели.

getindist команда обеспечивает доступ к модели в использовании.

Входное возмущение является ключевым фактором, который влияет на следующие атрибуты эффективности контроллера:

Динамический ответ на очевидные воздействия — символ ответа контроллера, когда измеренный объект вывел, отклоняется от своей предсказанной траектории, из-за неизвестного воздействия или ошибки моделирования.

Асимптотическое отклонение длительных воздействий — Если возмущение предсказывает длительное воздействие, корректировки контроллера, продолжается до объекта выход возвращается к своей желаемой траектории, эмулируя классический интегральный контроллер обратной связи.

Можно обеспечить входное возмущение как пространство состояний LTI (ssПередаточная функцияtf), или нули, полюса и усиление (zpk) объектное использование setindist. Контроллер MPC преобразует входное возмущение в дискретное время, без задержек, система в пространстве состояний LTI с помощью тех же шагов раньше преобразовывала модель объекта управления. Результат:

$\begin{matrix} x_{i d} (k + 1) = A_{i d} x_{i d} (k) + B_{i d} w_{i d} (k) \\ d (k) = C_{i d} x_{i d} (k) + D_{i d} w_{i d} (k) . \end{matrix}$

где _Aid, _Bid, _Cid и _Did являются постоянными матрицами пространства состояний, и:

_xid (k) является вектором из _nxid ≥ 0 входных состояний возмущения.
_dk (k) является вектором из _nd безразмерные неизмеренные входные воздействия.
_wid (k) является вектором из _nid ≥ 1 безразмерные белые шумовые входные параметры, принятые, чтобы иметь нулевое среднее значение и модульное отклонение.

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

Выведите возмущение

Выходное возмущение является особым случаем более общего входного возмущения. Его выход, _yod (k), непосредственно добавляется к объекту выход вместо того, чтобы влиять на состояния объекта. Выходное возмущение задает тип сигнала и характеристики _yod (k), и это часто используется на практике. Смотрите Оценку состояния Контроллера для получения дополнительной информации о модели.

getoutdist команда обеспечивает доступ к выходному возмущению в использовании.

Можно задать пользовательское выходное возмущение как пространство состояний LTI (ssПередаточная функцияtf), или нули, полюса и усиление (zpk) объектное использование setoutdist. Используя те же шаги что касается модели объекта управления, контроллер MPC преобразует заданное выходное возмущение в дискретное время, без задержек, система в пространстве состояний LTI. Результат:

$\begin{matrix} x_{o d} (k + 1) = A_{o d} x_{o d} (k) + B_{o d} w_{o d} (k) \\ y_{o d} (k) = C_{o d} x_{o d} (k) + D_{o d} w_{o d} (k) . \end{matrix}$

где _Aod, _Bod, _Cod и _Dod являются постоянными матрицами пространства состояний, и:

_xod (k) является вектором из _nxod ≥ 1 выходное состояние возмущения.
_yod (k) является вектором из _ny безразмерные выходные воздействия, которые будут добавлены к безразмерному объекту выходные параметры.
_wod (k) является вектором из _nod безразмерные белые шумовые входные параметры, принятые, чтобы иметь нулевое среднее значение и модульное отклонение.

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

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

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

Модель шума измерения

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

Используя те же шаги что касается модели объекта управления, контроллер MPC преобразует модель шума измерения в дискретное время, без задержек, система в пространстве состояний LTI. Результат:

$\begin{matrix} x_{n} (k + 1) = A_{n} x_{n} (k) + B_{n} w_{n} (k) \\ y_{n} (k) = C_{n} x_{n} (k) + D_{n} w_{n} (k) . \end{matrix}$

Здесь, _An, _Bn, _Cn и _Dn являются матрицами пробела постоянного состояния, и:

_xn (k) является вектором из _nxn ≥ 0 шумовых состояний модели.
_yn (k) является вектором из _nym безразмерные шумовые сигналы, которые будут добавлены к безразмерному измеренному объекту выходные параметры.
_wn (k) является вектором из _nn ≥ 1 безразмерные белые шумовые входные параметры, принятые, чтобы иметь нулевое среднее значение и модульное отклонение.

Если вы не предоставляете шумовую модель, значением по умолчанию является единица статическое усиление: _nxn = 0, _Dn является _nym-by-_nym единичная матрица, и _An, _Bn, и _Cn пуст.

Для mpc объект контроллера, MPCobj, свойство MPCobj.Model.Noise обеспечивает доступ к модели шума измерения.

Примечание

Если минимальное собственное значение $D_{n} D_{n}^{T}$ меньше ^1x10–8, контроллер MPC добавляет ^1x10–4 в каждый диагональный элемент _Dn. Эта корректировка делает успешное значение по умолчанию вычислением усиления Кальмана более вероятно.

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