Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Аналитическое конструирование оптимальных регуляторов. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Литература

1.Певзнер Л. Д., Чураков Е. П. Математические основы теории систем: учебное пособие. М.: Высшая школа, 2009. 503 с. : ил.

2.Афанасьев В. Н., Колмановский В. Б., Носов В. Р. Математическая теория конструирования систем управления. М.: Высшая школа, 1998. 574 с.

3.Ванько В. И., Ермошина О. В., Кувыркин Г. Н. Вариационное исчисление и оптимальное управление: учебник для вузов / под ред. В. С. Зарубина, А. П. Крищенко. 2-е изд. М.: Изд-во МГТУ, 2001. 488 с.

4.Куропаткин П. В. Оптимальные и адаптивные системы: учеб. пособие для вузов. М.: Высш. школа, 1980. 287 с.

11

Раздел 2. ВАРИАЦИОННЫЕ МЕТОДЫ ТЕОРИИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ

План

2.1.Основные понятия вариационного исчисления.

2.2.Вариационные задачи с фиксированными концами.

2.3.Условие трансверсальности.

2.4.Вариационные задачи с подвижными концами.

2.5.Решение задач вариационного исчисления на условный экстремум (синтез оптимального управления).

2.1. Основные понятия вариационного исчисления

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

вариационного исчисления.

Функционал является естественным обобщением понятия функция, используемым в математическом анализе. Поэтому кратко остановимся на свойствах экстремальных функций и функционалов. Каждому значению аргумента t из некоторой области изменения t соответствует значение x функции x(t). Аналогично каждой функции x(t) из некоторого класса функций x(t) соответствует значение J функционала J=J[x(t)]. Например, для

функционала

J = πsinωtdt

каждому

набору

функций

 

0

 

 

 

sin ω1t, sin ω2t, , sin ωnt

соответствуют числа J1 , J2 , , Jn .

Приращением аргумента функции является разностьt = t2 t1 ,

приращением или вариацией δx аргумента x(t) функционала J[x(t)] – разность между двумя функциями δx = x2(t) – x1(t). В данном примере вариация

x = sin ω2t – sin ω1t.

Таким образом, функционал можно рассматривать как функцию особого рода, в которой роль независимой переменной играет другая функция.

12

Примером функционала общего вида может быть выражение:

J = b f (t, x(t), x(t))dt .

a

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

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

Функция x(t) непрерывна, если малому изменению t соответствует малое изменение x(t). Функционал непрерывен, если малому изменению x(t) соответствует малое изменение J. Для функционалов, содержащих производные от функций, понятие непрерывности усложняется. Считается, что функционал непрерывен, если малому изменению функции и ее производных соответствует малое приращение функционала. Здесь необходимо уточнить, какие изменения функции x(t), являющейся аргументом функционала, называется малыми или, что то же самое, какие кривые x = x2(t) и x = x1(t) считаются мало отличающимися или близкими.

Уточним понятие «близость».

Пусть на координатной плоскости заданы кривые x0 (t) и x(t) .

Считается, что они имеют близость нулевого порядка, если выполняется условие:

x(t) x0 (t) ε .

Если для всех производных по времени до k-го порядка включительно выполняется условие:

di x(t)

di x (t)

ε,

i =1, ,k ,

dti

0

 

dti

 

 

то считается, что кривые x0 (t) и x(t) имеют близость порядка k.

Разность между близкими кривыми называется вариацией переменных x:

δ x(t) = x(t) x0 (t) . 13

Если функционал J[x(t)] достигает на кривой x0 (t) экстрему-

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

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

Аналогично понятию линейной функции линейным функционалом называется функционал J[x(t)], удовлетворяющий следующим условиям:

J [cx(t)]= cJ [x(t)] и J [x1 (t)+ x2 (t)]= J [x1 (t)]+ J [x2 (t)]

где с – произвольная постоянная.

Если в приращении функции x = x(t + ∆t) x(t) можно вы-

делить линейную по отношению к t часть приращения, то она называется дифференциалом функции и обозначается dx. Анало-

гично

выделение

в

приращении

функционала

J = J [x(t) +δ x]J [x(t)]

линейной по отношению к δx части

приращения дает нам вариацию функционала, которая обозначается δJ. Известно, что если некоторая функция в точке t* имеет экстремум, то dx = 0. Соответственно для экстремума функционала необходимо, чтобы его вариация равнялась нулю, т. е. δJ = 0.

Таким образом, решением задачи отыскания экстремума функционала служит функция x = x*(t), при подстановке которой в выражение функционала его величина принимает экстремальное

значение.

Например,

если

из

всех

функций

sinω1t, sinω2t, , sinωnt функция

sinωk t

доставляет

экстремум

функционалу J = πsinωk tdt , то любые вариации δ J 0 и только

0

та вариация δ J = 0 , для которой δx = sin ωкt – sin ωкt = 0. Функция, на которой вариация функционала δ J = 0 , называется экс-

тремаль.

В зависимости от того, придает функция x*(t) функционалу J[x(t)] максимальное или минимальное значение, она называется

максималью или минималью.

14

2.2. Вариационные задачи с фиксированными концами

Решение так называемого уравнения Эйлера является экстремалью. Рассмотрим упрощенный вывод этого уравнения. Исследуем на экстремум функционал

tк

 

J [x(t)]= f (t, x(t), x(t))dt

(2.1)

tн

Предполагаем, что искомая функция удовлетворяет краевым условиям X(tн)=Xн; X(tк)=Xк. Иными словами, искомая кривая проходит через заданные точки А(t0,x0) и В(t1,x1) (рис. 2.1).

Рис. 2.1. К вариационной задаче с закрепленными концами

Такая задача называется вариационной задачей с закрепленными концами или с закрепленными граничными точками. Пред-

положим, что экстремум достигается на кривой x(t) (кривая 1). Возьмем близкую к ней кривую x(t) (кривая 2), также проходя-

щую через точки А и В. Тогда вариация функции δ x = x(t) x(t) , для граничных точек δ x(tн ) =δ x(tk ) = 0 .

Проварьируем функции x(t) и x(t) , т. е. образуем функции x +δ x и x +δ x . В результате функционал J получит приращение, которое представим в виде ряда Тейлора:

15

tк

J = f (t, x +d x, x

tн

tк f (t, x, x) d xdt +

tн x

 

t

к

 

 

 

 

 

 

+d x)dt f (t, x, x)dt

 

 

 

(2.2)

tк

tн

tк

 

 

 

 

f (t, x, x)

 

2 f (t, x, x)

 

(d x)2

 

x

d xdt

+

x

2

 

 

dt + .

2

tн

 

tн

 

 

 

 

 

 

 

 

 

 

Поскольку значения δ x и δ x малы, то определяющими в правой части (2.2) являются линейные относительноδ x и δ x чле-

ны. Линейную часть приращения функционала называют егопервой вариацией δ J . Для функционала типа (2.1) первая вариация будет

tк

f (t, x, x)

 

 

f (t, x, x)

 

 

 

dJ =

 

d

x +

 

d

x

dt = 0 . (2.3)

x

x

tн

 

 

 

 

 

Необходимым условием существования экстремального зна-

чения интеграла типа (2.3) является равенство нулю его первой вариации, т. е.

tк

f (t, x, x)

d x +

f (t, x, x)

d x

 

(2.4)

x

x

dt = 0 .

tн

 

 

 

 

Интегрируя по частям второе слагаемое в выражении (2.4) получим

tк

f (t, x, x)

 

 

 

f (t, x, x)

 

 

 

tк

 

tk

d

f (t, x, x)

 

 

 

 

 

 

 

 

d xdt

=

 

 

 

 

 

d x

 

tн

 

 

 

 

d x dt .

x

 

 

x

 

 

dt

x

tн

 

 

 

 

 

 

 

 

 

 

tн

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно краевым условиям

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

δ x(tн ) =δ x(tк ) = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tк

 

 

 

 

 

 

 

 

 

 

 

 

f (t, x, x)

δ x

 

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

tн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (2.4) примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

tк

f (t, x, x)

 

 

 

tк

d f (t, x, x)

 

 

 

x

 

 

 

d xdt

 

 

 

 

 

x

 

d x dt =

 

 

 

 

 

 

 

 

 

 

 

 

tн

 

 

 

 

 

 

tн

dt

 

 

 

 

 

 

 

tк f

(t, x, x)

 

d

 

f (t, x, x)

 

 

 

=

 

 

x

 

 

 

 

 

 

 

 

x

 

d xdt = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tн

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

(2.5)

(2.6)

Так как величина δ x имеет произвольный характер, то на основании леммы Лагранжа последнее равенство выполнимо, если

f (t, x, x)

d

f (t, x, x)

= 0 .

(2.7)

x

 

 

x

 

 

 

dt

 

 

 

Уравнение (2.7) называют дифференциальным уравнением Эйлера, которое является необходимым условием экстремума интеграла (2.1) при фиксированных граничных условиях и отсутствии ограничений на координаты. Если граничные условия нефиксированны, то при решении уравнений Эйлера требуется учитывать условия трансверсальности.

Рассмотрим методику применения уравнения Эйлера. Пусть требуется определить экстремаль для функционала

 

 

π / 2

 

 

 

 

 

 

 

J[x(t)] =

[x

2

2

]dt,

x(0) =1; x(π / 2) = 0.

 

 

(x )

 

 

 

 

0

 

 

 

 

 

 

В этом случае уравнение Эйлера имеет вид x′′+ x = 0 , так как

f

f

 

. Его общим решением является x = c1cos t +

x = 2x;

x= −2x

 

c2cos t.

Используя граничные условия, получим с1 = 1, с2 = 0. Следовательно, экстремалью является лишь кривая x = cos t.

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

tк

 

J [x(t)]= f (t, x, x, x, x(n) )dt .

(2.8)

tн

Проведя для этого функционала аналогичные преобразования, как и для функционала (2.1), получим уравнение

f

d

f

+

d 2 f

− + (1)n

d (n)

 

f

= 0 .

(2.9)

 

 

 

 

 

 

dt2 x

dt(n) x(n)

x dt x

 

 

 

 

Это дифференциальное уравнение носит название уравнения Эйлера-Пуассона, а его решение называется экстремалью рассматриваемой вариационной задачи.

17

Докажем формулу (2.9), ограничившись для простоты выкладок случаем второй производной, т. е.

t к

 

J [x(t)]= f (t, x, x, x)dt .

(2.10)

tн

Прибавим к функции x(t) вариацию δ x такую, что при t = tн и t = tк δ x и δ x обращаются в нуль.

Главная, линейная, часть приращения функционала (его первая вариация) будет равна

 

tк

f

 

f

 

 

f

 

 

d J =

d x +

 

d x

+

 

 

(2.11)

 

x

x

x

d x dt .

 

tн

 

 

 

 

 

Преобразуем выражение (2.11), пользуясь интегрированием по частям:

 

 

 

tк

f

 

 

 

 

 

f

 

 

 

 

tк

 

tк

 

 

d f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d xdt

=

 

 

 

d x

t

 

 

 

 

 

 

 

d xdt,

 

 

 

 

x

x

н

 

dt x

 

 

 

 

tн

 

 

 

 

 

 

 

 

tн

 

 

 

 

 

 

tк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tк

 

 

 

 

 

 

 

f

 

 

f

 

tк

 

 

d

f

 

 

 

tк

 

 

d 2 f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d xdt =

 

 

 

d x

t

 

 

 

 

 

 

 

d x

t

 

 

+

 

 

 

 

 

d

 

x

x

н

dt

x

н

dt

2

x

 

tн

 

 

 

 

 

 

 

tн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так

как

в

 

 

задаче

 

 

 

с

 

 

 

закрепленными

(2.12)

xdt. (2.13)

концами

δx(tн )= δx(tk )= δx(tн )= δx(tk )= 0 , то все неинтегральные члены (2.12) и (2.13) обращаются в нуль, а вариация функционала будет

d J = tк

f

d

f

+

d 2

f d xdt .

(2.14)

 

 

 

 

 

x

 

dt x

 

dt

2

 

 

tн

 

 

 

x

 

Так как необходимым условием экстремума является равенство первой вариации нулю, то из (2.14) следует уравнение Эйлера – Пуассона вида

f

d f

+

d 2

f

= 0

,

(2.15)

 

 

 

 

 

x

dt x

dt2

x

 

 

 

 

 

которое в данном случае является дифференциальным уравнением четвертого порядка.

Определим экстремаль функционала.

J[x(t)] = π/ 2[(x′′)2 x2 ]dt ,

0

18

удовлетворяющего граничным условиям

 

 

 

π

 

π

x(0)

=1; x(0)

= 0; x(

2 ) = 0;

x(

2 ) = −1.

В этом случае уравнение Эйлера-Пуассона имеетследующий вид

x(IV ) x = 0 .

Его общим решением является выражение

x = c1et + c2et + c3 cost + c4 sin t .

Используя граничные условия, получаем

c1 = c2 = c4 = 0; c3 =1.

Таким образом, экстремалью является кривая x = cost.

2.3. Условие трансверсальности

Общее решение уравнения Эйлера содержит две произвольные постоянные, для определения которых необходимо иметь два условия. В задаче с неподвижными граничными точками такими условиями являются X(tн) = Xн; X(tk) = Xk. В задаче с подвижными граничными точками последние могут находиться в любых точках некоторых поверхностей или линий. В этом случае для определения постоянных интегрирования нужно иметь некоторые условия, которые получаются из основного необходимого условия экстремума – равенства нулю вариации функционала δJ = 0. Требуется найти экстремаль, которая соединила бы две линии ψ(t) и φ(t) по кратчайшему пути. Граничные точки подвижны и экстремаль x(t) может начинаться в любой точке на линии ψ (t) и оканчиваться

в любой точке на ϕ(t) . Точки эти находятся из следующих соот-

ношений, которые называются условиями трансверсальности

f (t, x, x)

 

x

 

 

f (t, x, x)

 

 

x

 

+ (ϕ x)

f (t, x, x)

= 0;

x

 

 

t=t1

(2.16)

 

f (t, x, x)

+ (ψ x)

= 0.

x

 

 

t=t0

 

Условия трансверсальности дают возможность определить постоянные интегрирования в уравнении Эйлера и тем самым получить конкретное решение.

19

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

Геометрический смысл трансверсальности виден из рисунка 2.2.

Рис. 2.2. Условие трансверсальности

2.4. Вариационные задачи с подвижными концами

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

Приращение функционала при переходе от x до x +δ x можно записать в виде:

tк +d tк

 

tк

 

J = J (x +d x) J (x) = f (t, x +δ x, x +δ x)dt f (t, x, x)dt =

 

tн +d tн

 

tн

 

tк

 

 

tк

(2.17)

f (t, x d x, x

d x)dt

f (t, x, x)dt +

tн

 

 

tн

 

tк +d tк

 

tн +d tн

 

+ f (t, x +d x, x +d x)dt

 

f (t, x +δ x, x +δ x)dt.

 

tк

 

 

tн

 

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]