|
)+ [cɺ(tf |
|
)] |
∂F |
|
|
|
|
|
∂F |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||
F(tf |
)− xɺ(tf |
|
|
|
|
δtf |
− η(t0 ) |
∂ xɺ |
|
|
= 0 . |
(4.2.23) |
||
∂ xɺ |
|
|
|
|
||||||||||
|
|
|
|
|
t f |
|
|
|
|
t0 |
|
|
||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из первого условия, как и ранее, следует уравнение Эйлера-Лагранжа (4.2.16).
Второе слагаемое в условии (4.2.23) равно нулю, если величина ∂F
∂ xɺ t0
ограничена, поэтому равно нулю должно быть выражение в фигурных скобках соотношения (4.2.23). Отсюда следует уравнение
∂F |
|
|
F(tf )= [xɺ(tf )− cɺ(tf )]∂ xɺ |
. |
(4.2.24) |
|
t f |
|
Соотношение (4.2.24) называется условием трансверсальности. Нелишне отметить, что уравнение (4.2.16) в начальный момент времени (то есть при t = t0 ) может и не выполняться, поскольку η(t0 )= 0 и тогда условие (4.2.22) справедливо. Важно лишь, чтобы в начальный мо-
мент времени x(t) была бы непрерывна, а |
∂F – конечна. |
|
∂ xɺ |
Если речь идёт о системе регулирования, когда осуществляется переход от одного существенного режима работы к другому, то есть от одного значения x к другому значению, конечная кривая x = c(t) является горизонтальной прямой и cɺ(tf )= 0 . В этом случае условие трансверсальности (4.2.24) упрощается и принимает вид
221
∂F |
|
|
F(tf )= xɺ(tf )∂ xɺ |
. |
(4.2.25) |
|
t f |
|
Можно показать, что в случае, когда функция F явно не зависит от времени, соотношения (4.2.16) и (4.2.25) сводятся к одному, а именно:
F(x, xɺ)= xɺ |
∂F |
|
∂ xɺ . |
(4.2.26) |
Действительно, дифференцируя обе части последнего выражения по t как сложную функцию, после упрощения получаем
|
∂F |
|
d ∂F |
|
||||
xɺ |
|
|
− |
|
|
|
|
= 0 . |
|
|
|
||||||
|
∂ x |
|
dt |
|
ɺ |
|
|
|
|
|
|
|
∂ x |
|
|
||
|
|
|
|
|
|
|
|
|
Обычно в процессе управления производная xɺ не обращается в нуль, поэтому нулю должно равняться выражение в квадратных скобках, а это левая часть уравнения Эйлера-Лагранжа (4.2.16).
Таким образом, оптимальную систему регулирования, осуществляющую переход от одного рабочего режима к другому за время (t0 ,tf ),
определяет единственное условие (4.2.26).
4 . 2 . 3 . Оптим ально е упр авление по м иним ум у интегр ально го кр итер ия
Теперь рассмотрим применение вариационного исчисления к задаче оптимального управления с интегральным критерием.
Пусть динамика объекта задана уравнением 222
xɺ(t)= f (x,u), |
(4.2.27) |
где x – выход объекта, u – управляющее воздействие при заданном начальном значении
x(t0 )= x0 . |
(4.2.28) |
Требуется сформировать управление u(t), минимизирующее интегральный критерий
t f |
|
I = ∫ F(x,u)dt . |
(4.2.29) |
t0 |
|
Функция F в критерии (4.2.29) может описывать отклонение от желаемого состояния, потребление энергии, стоимость САУ и т.п.
Обозначая функции, доставляющие минимум критерия (4.2.29) через x(t) и u(t), а некоторые близкие к ним – через x(t)и u(t) соответственно, можно записать
x = x + εη, |
(4.2.30) |
u = u + εξ, |
где η и ξ – в достаточной мере произвольные функции времени, определённые на интервале (t0 ,tf ), а ε – малый параметр.
Подставляя выражения (4.2.30) в критерий (4.2.29) вместо x(t) и u(t), получим
223
I(ε)= ∫f F(x + εη,u + εξ)dt .
t0
Разложив в последней формуле подынтегральную функцию в ряд Тейлора вблизи оптимальной траектории, получим
|
f |
|
|
η |
∂F |
+ ξ |
∂F |
|
+ членыс ε2 |
, ε3 |
|
I(ε)= ∫ |
F(x,u)+ ε |
|
|
|
,... dt . (4.2.31) |
||||||
t |
0 |
|
|
|
∂ x |
|
∂u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Необходимое условие экстремума функции (4.2.31) – равенство
∂I(ε) |
|
|
|
= 0 . |
|
∂ε |
|
ε=0 |
|
||
Взяв производную по ε от правой части выражения (4.2.31) и поменяв местами порядок дифференцирования по ε и интегрирования по t, приравняем полученное выражение нулю при ε=0
t f |
|
∂F |
|
∂F |
|
∫ |
|
|
|
|
|
|
η |
|
+ ξ |
dt . |
(4.2.32) |
t0 |
|
∂ x |
|
∂ u |
|
Вернёмся к уравнению динамики объекта (4.2.27) и заменим x(t) и u(t) на x(t)и u(t) соответственно. Получим
xɺ(t)+ εηɺ = f (x + εηɺ,u + εξ).
224
Разложив правую часть последнего выражения в ряд Тейлора вблизи оптимальной траектории по малым значениям ε, и ограничиваясь линейными членами, получим
x(t)+ εη = |
|
η |
∂f |
+ ξ |
∂f |
|
|
f (x,u)+ ε |
|
|
. |
||||
ɺ |
ɺ |
|
|
∂x |
|
∂u |
|
|
|
|
|
||||
С учётом уравнения (4.2.27) последнее соотношение превратится в уравнение
ɺ |
∂f |
+ ξ |
∂f |
(4.2.33) |
η = η |
∂x |
∂u |
с начальными условиями η(t0 )= 0 .
Решая уравнение (4.2.33) относительно ξ, подставим это значение в уравнение (4.2.32)
|
|
|
|
|
t f |
|
∂F |
|
η − η |
∫ |
η |
|
+ |
|
|
|
∂x |
|
ɺ |
t0 |
|
|
|
|
|
|
|
|
|
∂F
∂f ∂u dt = 0 . (4.2.34) ∂x ∂f
∂u
Слагаемое, связанное с производной ηɺ в выражении (4.2.34), найдем, интегрируя по частям
t f |
∂F |
|
|
∂F |
|
t f |
|||
|
|
|
|||||||
|
|
|
|
|
|
|
|
||
∂u |
|
|
∂u |
|
|
|
|||
ɺ |
|
|
|
|
|
|
|||
∫η |
∂f |
dt = η |
|
∂f |
|
|
|
||
t0 |
|
|
|
|
|||||
∂u |
|
|
∂u |
|
t0 |
||||
|
|
|
|
|
|||||
|
|
|
|
|
|||||
|
|
∂F |
|
||||
t f |
d |
|
∂u |
|
|
|
|
− ∫η |
|
|
|
|
|
dt . |
(4.2.35) |
|
|
∂f |
|
||||
t0 |
dt |
|
|
|
|||
∂u |
|
|
|||||
|
|
|
|
|
|
||
225