В теория управления, Формула Аккермана это система контроля расчетный метод решения распределение полюсов проблема для систем с инвариантным временем Юрген Акерманн.[1] Одной из основных проблем при проектировании систем управления является создание контроллеров, которые изменят динамику системы путем изменения собственных значений матрицы, представляющей динамику замкнутой системы.[2] Это эквивалентно изменению полюсов связанных функция передачи в случае, если нет отмены полюсов и нулей.
Контроль состояния обратной связи
Рассмотрим линейную инвариантную систему с непрерывным временем с представление в пространстве состояний


куда Икс - вектор состояния, ты - входной вектор, а А, B и C матрицы совместимых размеров, которые представляют динамику системы. Описание ввода-вывода этой системы дается функция передачи

Поскольку знаменатель правого уравнения задается характеристический многочлен из А, полюса грамм находятся собственные значения из А (обратите внимание, что обратное не обязательно верно, так как между членами числителя и знаменателя могут быть сокращения). Если система неустойчивый, или имеет медленный отклик или любую другую характеристику, которая не определяет критерии проектирования, может быть полезно внести в нее изменения. Матрицы А, B и Cоднако могут представлять физические параметры системы, которые нельзя изменить. Таким образом, одним из подходов к этой проблеме может быть создание петли обратной связи с коэффициентом усиления K который будет кормить переменную состояния Икс во вход ты.
Если система управляемый, всегда есть вход
такое, что любое государство
может быть переведен в любое другое состояние
. Имея это в виду, в систему можно добавить контур обратной связи с управляющим входом.
, так что новая динамика системы будет
![{ displaystyle { dot {x}} (t) = Ax (t) + B [r (t) -Kx (t)] = [A-BK] x (t) + Br (t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46b2f272a6d9f6d9961e14fb365a04cedef8bffb)

В этой новой реализации полюса будут зависеть от характеристического полинома
из
, то есть

Формула Аккермана
Вычисление характеристического полинома и выбор подходящей матрицы обратной связи может оказаться сложной задачей, особенно в больших системах. Один из способов упростить вычисления - использовать формулу Аккермана. Для простоты рассмотрим один входной вектор без ссылки параметра
, Такие как


куда
- вектор обратной связи совместимых размеров. Формула Аккермана утверждает, что процесс проектирования можно упростить, вычислив только следующее уравнение:
![{ Displaystyle к ^ {T} = left [0 0 cdots 0 1 right] { mathcal {C}} ^ {- 1} Delta _ { text {new}} (A) ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b8203c830fa477535079e2847a234c424dc3934)
в котором
- искомый характеристический многочлен, вычисленный на матрице
, и
это матрица управляемости системы.
Доказательство
Это доказательство основано на Энциклопедия систем жизнеобеспечения запись о контроле за размещением полюсов.[3] Предположим, что система управляемый. Характеристический многочлен
дан кем-то

Расчет степеней
приводит к
![{ displaystyle { begin {align} (A_ {CL}) ^ {0} & = (A-Bk ^ {T}) ^ {0} = I (A_ {CL}) ^ {1} & = (A-Bk ^ {T}) ^ {1} = A-Bk ^ {T} (A_ {CL}) ^ {2} & = (A-Bk ^ {T}) ^ {2} = A ^ {2} -ABk ^ {T} -Bk ^ {T} A + (Bk ^ {T}) ^ {2} = A ^ {2} -ABk ^ {T} - (Bk ^ {T}) [A -Bk ^ {T}] = A ^ {2} -ABk ^ {T} -Bk ^ {T} A_ {CL} vdots (A_ {CL}) ^ {n} & = (A- Bk ^ {T}) ^ {n} = A ^ {n} -A ^ {n-1} Bk ^ {T} -A ^ {n-2} Bk ^ {T} A_ {CL} - cdots - Bk ^ {T} A_ {CL} ^ {n-1} end {выровнено}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b7321dd8c36adc83c038a386ab6b99e36c8131e)
Замена предыдущих уравнений на
дает

Переписывая приведенное выше уравнение в виде матричного произведения и опуская члены, которые

не появляется единичных урожаев
![{ Displaystyle Delta (A_ {CL}) = Delta (A) - left [B AB cdots A ^ {n-1} B right] left [{ begin {array } {c} star vdots k ^ {T} end {array}} right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69da073f81563572410f590f3701a57968a459e5)
От Теорема Кэли – Гамильтона,
, таким образом
![{ displaystyle left [B AB cdots A ^ {n-1} B right] left [{ begin {array} {c} star vdots k ^ { T} end {array}} right] = Delta (A)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8df1b69bbb265e1bfe60d154e7906be65b21782)
Обратите внимание, что
это матрица управляемости системы. Поскольку система управляема,
обратимо. Таким образом,
![{ displaystyle left [{ begin {array} {c} star vdots k ^ {T} end {array}} right] = { mathcal {C}} ^ {- 1} Delta (A)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e15b78f4811def2fd01771b2e52f8a16357562b5)
Найти
, обе части можно умножить на вектор
давая
![{ displaystyle left [{ begin {array} {ccccc} 0 & 0 & 0 & cdots & 1 end {array}} right] left [{ begin {array} {c} star vdots k ^ {T} end {array}} right] = left [{ begin {array} {ccccc} 0 & 0 & 0 & cdots & 1 end {array}} right] { mathcal {C}} ^ {- 1} Delta (A)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0348e89a0003bdbbf755b54e7dc936fc375f854c)
Таким образом,
![{ displaystyle k ^ {T} = left [{ begin {array} {ccccc} 0 & 0 & 0 & cdots & 1 end {array}} right] { mathcal {C}} ^ {- 1} Delta (A )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a69c0cb248d6e52117000700a83285063dba4f27)
Пример
Учитывать[4]
![{ displaystyle { dot {x}} = left [{ begin {array} {cc} 1 & 1 1 & 2 end {array}} right] x + left [{ begin {array} {c} 1 0 конец {массив}} right] u}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f85b344485a422e2a500e66a550160ced77622d)
Мы знаем из характеристического полинома
что система нестабильна, так как
, матрица
будут иметь только положительные собственные значения. Таким образом, для стабилизации системы положим коэффициент обратной связи ![{ displaystyle K = left [{ begin {array} {cc} k_ {1} & k_ {2} end {array}} right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/729831fc20767ef795e56e0133dbb3d79923207f)
Из формулы Аккермана можно найти матрицу
это изменит систему так, что ее характеристическое уравнение будет равно желаемому многочлену. Предположим, мы хотим
.
Таким образом,
и вычисление матрицы управляемости дает
и ![{ displaystyle { mathcal {C}} ^ {- 1} = left [{ begin {array} {cc} 1 & -1 0 & 1 end {array}} right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c39f68769097cc9a6feecc819a09ec59dc22e063)
Также у нас есть это ![{ displaystyle A ^ {2} = left [{ begin {array} {cc} 2 & 3 3 & 5 end {array}} right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80eff75ad33eb017b8a69b51689ad22d3944cb76)
Наконец, из формулы Аккермана
![{ Displaystyle к ^ {T} = left [{ begin {array} {cc} 0 & 1 end {array}} right] left [{ begin {array} {cc} 1 & -1 0 & 1 end {array}} right] left [ left [{ begin {array} {cc} 2 & 3 3 & 5 end {array}} right] +11 left [{ begin {array} {cc} 1 и 1 1 и 2 end {array}} right] + 30I right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f5f9030a690179cdb94b8231ce357b2fd0bd6dcc)
![{ Displaystyle к ^ {T} = left [{ begin {array} {cc} 0 & 1 end {array}} right] left [{ begin {array} {cc} 1 & -1 0 & 1 end {array}} right] left [{ begin {array} {cc} 43 & 14 14 & 57 end {array}} right] = left [{ begin {array} {cc} 0 & 1 end { массив}} right] left [{ begin {array} {cc} 29 & -43 14 & 57 end {array}} right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a61e0d18412ce3873d86daeecfec5a73e0170537)
![{ displaystyle k ^ {T} = left [{ begin {array} {cc} 14 и 57 end {array}} right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26f77f5d5a2f0d870b5fd8f422744daeeea665d1)
Рекомендации
Смотрите также
внешняя ссылка