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

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

.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
533 Кб
Скачать

4. ДОСТАТОЧНЫЕ УСЛОВИЯ ОПТИМАЛЬНОСТИ

4.1.Постановка задачи. Основная теорема

Рассмотрим задачу оптимального управления

T

 

 

I (x,u) = f (t,x(t),u(t))dt min ,

(1)

0

 

 

x(t) = ϕ (t,x(t),u(t)),

t [0,T ] ,

(2)

x(0) = x0 ,

(3)

(x(t),u(t)) V t ,

t [0,T ].

(4)

В этой задаче условие (4) представляет собой обобщение условия (4) из предыдущего параграфа, множество V t R × R задано и при каждом t задает ограничение на положение точки (x(t),u(t)). Таким образом, ограничение

(4) – смешанное ограничение, связывающее фазовую переменную x и

управляющую переменную u . Ограничение u(t) U, t [0,T ] , очевидно,

можно записать в виде (4), если положить

V t = V = R ×U .

Для формулировки достаточных условий оптимальности введем вспомогательные функции W (t, x) и R(t, x,u) . Функция R строится по функции

W , при этом роль функции R близка к той роли, которую играет в принци-

пе максимума функция Гамильтона – Понтрягина H :

 

R(t,x,u) =

W +

W ϕ (t,x,u) f (t, x,u).

(5)

 

t

x

 

Заметим, что сумма первых двух слагаемых в правой части этого выражения представляет собой производную от функции W в силу системы (2):

d

(W (t, x(t))) =

W

(t, x)

 

x=x(t) +

W

(t, x)

dx

 

x=x(t) =

 

 

 

 

dt

t

 

 

x

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

W (t, x)

 

x=x(t) +

W (t, x)

 

x=x(t) ϕ

(t, x(t),u(t)),

 

 

 

 

 

 

 

 

t

 

 

 

 

x

 

 

 

 

 

 

 

 

где x(t) – траектория системы (2).

41

Допустимым процессом управления в задаче (1)–(4) называется пара

(x,u) KC1 [0,T ] × KC[0,T ] , удовлетворяющая условиям (2)–(4).

Основная теорема о достаточных условиях оптимальности. Пусть

функция W (t, x) и допустимый процесс управления (x,u ) таковы, что

 

(

t,

x

(

t

)

 

 

u

(

t

))

=

(x,u) Vt

 

(

t,x,u

)

 

t

[

0,T

]

 

а)R

 

 

 

 

,

 

 

 

max

R

 

 

,

 

 

,

б) W

 

(

T,

x

(T )

)

 

= minW (T, x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x VxT

 

 

 

 

 

 

 

 

 

 

 

где VxT

– проекция множества V T

на ось x.

 

 

 

 

 

Тогда (

x

,

u

) – решение задачи (1)–(4).

 

 

 

 

 

Доказательство. Введем вспомогательный функционал Λ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Λ(x,u)

= −R(t, x(t),u(t))dt +W (T,x(T ))W (0,x0 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

Этот функционал обладает следующими свойствами:

1. Пара (x,u ), очевидно в силу условий а) и б), доставляет минимум этому функционалу на множестве всех пар (x,u), удовлетворяющих условию (4). 2. На множестве допустимых процессов управления исходный функционал I

и функционал Λ совпадают.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Λ

(

x,u

)

 

= −T

d

W

(

t, x

(

t

dt + T

f

(

t, x

(

t

)

,u

(

t

))

dt +

 

 

 

 

dt

 

 

 

))

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+W (T, x(T ))W (0, x0 ) =

 

 

 

 

 

 

 

= −W (t, x(t))

 

T0 +W (T, x(T ))W (0, x0 )

T

 

 

(t, x

(t),u(t))dt =

 

 

 

 

+ f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

(t, x(t),u(t))dt = I (x,u).

 

 

 

 

 

 

 

 

 

= f

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

42

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

Заметим, что имеющуюся свободу в выборе функции W (t, x) можно

использовать, чтобы соответствующая функция R(t, x,u) допускала эффективное нахождение точек экстремума.

4.2. Примеры

Пример 1. Рассмотрим задачу.

I (x,u) = 5 (2x2u x)dt min ,

(6)

 

 

0

 

 

 

 

 

 

x = xu, t [0,5],

(7)

 

 

x(0) = 2 ,

 

(8)

x(5) = 4, u(t) U = [0,1], t [0,5].

(9)

Запишем условия (9) в виде (4). Для этого определим множество V t

следующим образом:

 

 

 

 

 

 

 

R×[0,1],

если t [0,5),

 

V t = {4}×[0,1],

если t = 5.

 

1. Запишем общий вид функции R:

 

 

 

 

R(t,x,u)

= W +

W xu (2x2u x) =

 

 

 

t

x

 

 

 

=

W

W

 

+ x .

(10)

t

+ xu

x

2x

 

 

 

 

 

2. Выберем W (t, x) так, чтобы R(t, x,u)

не зависела от u

(обратилось

в нуль выражение в круглых скобках):

W = −2x. x

43

Такой W (t, x) является

W (t, x) = −x2 .

При этом

 

R(t,x,u) = x

(11)

и для выполнения условия а) теоремы из всех возможных значений x следует выбрать наибольшее.

3. Дадим

описание множества всех допустимых значений. При каждом

t [0,5]

точка (t, x(t))должна лежать на допустимой траектории (на гра-

фике решения уравнения (7)). С учетом условий (8) и (9) нас интересуют

только траектории, выходящие из точки (0, 2) и заканчивающиеся в точке (5, 4) при всех допустимых управлениях u(t).

Все траектории, выходящие из точки (0,2) при управлениях u(t) [0,1] , образуют так называемую интегральную воронку, т.е. располагаются между «верхней» траекторией, соответствующей управлению u(t) = 1, и «нижней»

– соответствующей управлению u(t) = 0.

Для нахождения границ этой интегральной воронки следует решить две

задачи:

 

 

 

x(t) = x(t),

t [0,5],

x(0) = 2

(12)

и

 

 

 

x(t) = 0,

t [0,5],

x(0) = 2 .

(13)

Их решения: 2et и 2, соответственно.

 

 

Таким образом,

 

 

 

2 x(t) 2et , t [0,5].

(14)

Все траектории, заканчивающиеся в точке (5, 4), также образуют инте-

гральную воронку. Ее границы является решениями задач

 

x(t) = x(t),

t [0,5],

x(5) = 4

(15)

44

и

x(t) = 0, t [0,5], x(5) = 4.

(16)

Решая эти задачи (см. приложение 2), получаем

4et5 x(t) 4, t [0,5].

Так как допустимые траектории удовлетворяют условиям (8) и (9) одновременно, они заполняют множество, представляющее собой пересечение двух построенных интегральных воронок. Это множество Vxt показано на рис. 8.

Рис. 8

4. Учитывая, что функционал I рассматривается только на допустимых про-

цессах управления, множество V t

можно заменить более узким множеством

 

 

t =Vxt ×[0,1] = {(x,u): x Vxt , u [0,1]}, t [0,5].

V

5. Так как функция R(t,x,u) = x

(см. (11)) не зависит от u , то

(

max R(t,x,u) = maxR(t, x,u).

)

x

 

 

 

x,u V t

x V t

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

45

точку множества Vxt с максимально большой координатой x. Таким образом,

 

 

t

 

t [0, ln 2),

x

(t) =

2e

,

 

 

4,

 

t [ln 2, 5].

 

 

 

 

 

6. На рис. 8 эта траектория отмечена знаком «0».

Условие б) теоремы выполняется автоматически, так как множество VxT со-

7. По оптимальной траектории

 

 

x

(t) оптимальное управление

u

(t) можно

построить, используя уравнение (7), связывающее

x и u . Из него следует,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

(t) = 1,

t [0,ln 2),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

t [ln 2, 5].

 

 

 

 

 

 

 

 

 

8. Найдем Imin :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln2

 

 

 

 

2t

t

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Imin = I (x,u ) =

 

 

(4)dt

=

 

 

 

2(4e

 

)2e dt +

 

 

 

 

 

ln 2

=

{

0

 

 

 

 

ln 2

 

ln2

 

 

 

 

 

 

стоит из одной точки: VT

4 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 4 e2td (2t) 2 etdt 4(5 ln2) =

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 4(5 ln 2) = 4ln2 10.

 

 

 

 

= 4e2t

ln0

2 2et

 

ln0

 

 

 

 

 

 

 

Для сравнения приведем значение I1 , соответствующее допустимой траекто-

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

рии x1 (t) = 2 +

t и допустимому управлению u1 (t) =

5

 

 

 

:

 

 

 

 

2

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 +

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

I1 = I (x1,u1 ) = −4.

Пример 2. Рассмотрим задачу

46

 

 

 

 

4

 

 

(x2 8x) dt min,

 

 

I (x,u) = x2u + 2

(17)

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

x = 2x(1u),

 

 

 

 

(18)

 

 

 

 

 

 

x(0) = 1,

 

 

 

 

(19)

 

 

 

 

x(4) = 5,

u [0,1] .

 

 

(20)

1. Запишем общий вид функции R:

 

 

 

 

 

 

R(t,x,u) =

W +

W 2x(1u) x2u 2(x2 8x) =

 

 

 

 

t

 

x

 

 

 

 

 

 

=

W +

W 2x

W 2xu x2u 2(x2 8x) =

 

 

t

 

x

 

 

x

 

 

 

 

 

 

=

W

 

 

W

+ x

 

2

8x)+

W

2x .

 

t

xu

2

x

2(x

 

x

 

 

 

 

 

 

 

 

 

 

 

2. Выберем W (t, x) так, чтобы R(t, x,u) не зависела от u . Для этого в ка-

честве W (t, x) возьмем решение уравнения 2 W + x = 0:

x

W (t, x) = − x2 .

4

При этом окончательный вид функции R определяется равенством

1 def

R(t,x,u) = −2x2 +16x x2x = −3x2 +16x = r(x).

2

Эта функция принимает максимальное значение, если x является решением

уравнения

d

r(x) = 0: 6x +16 = 0, т.е. при

x = 2

2

. Это значение дос-

 

 

 

dx

3

 

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

x = 2 2 . 3

3. Построим, как и в примере 1, множество всех допустимых траекторий – пересечение двух интегральных воронок. Их границы находятся как решения

47

Vxt

следующих четырех задач (последние возникают после подстановки в уравнение (18) «граничных» управлений u(t) 1 и u(t) 0 с использованием начального и конечного условий):

x(t) = 0, t [0,4]; x(0) =1,

(21)

x(t) = 2x(t),

t [0,4];

x(0) = 1,

(22)

x(t) = 0,

t [0,4];

x(4) = 5,

(23)

x(t) = 2x(t),

t [0,4];

x(4) = 5.

(24)

Эти задачи имеют решения x(t) 1, x(t) e2t , x(t) 5,

x(t) 5e2(t4)

соответственно.

 

 

 

Множество Vxt для данного случая показано на рис. 9. На этом же рисунке штрихами показана прямая значений x, доставляющих максимальное

 

 

 

2

 

значение функции R: Rmax

= rmax

= r 2

 

.

3

 

 

 

 

Рис. 9

4. V t =Vxt ×[0,1] = {(x,u): x Vxt , u [0,1]}, t [0,4].

5. Построение оптимальной траектории (траектории, удовлетворяющей усло-

вию а) теоремы): фиксируем t = t , находим сечение , из всех точек x

этого множества выбираем точку, ближайшую к x = 2 2 .

3

48

Реализуя этот алгоритм, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2t ,

 

 

 

 

 

t [ 0,

1

ln2

2

),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

2

 

 

 

1

 

2

 

1

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

(t) = 2

 

,

t [

 

ln2

 

 

 

,

 

 

[8

+ ln(

 

 

2

 

) ] ],

 

3

2

3

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

5e

2(t4)

,

t (

1

[8

+ ln(

1

 

2

2

)],

4 ].

 

 

 

 

 

 

 

 

 

 

 

 

2

5

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Впервой строке этой формулы на указанном промежутке точки кривой x = e2t ближе к точкам прямой x = 2 2 , чем остальные допустимые значения.

3

Во второй строке на указанном промежутке значение x = 2

2

принадле-

 

 

 

3

 

жит множеству допустимых значений x.

 

В третьей строке на указанном промежутке точки кривой

x = 5e2(t4)

ближе к точкам прямой x = 2

2

, чем остальные допустимые значения.

 

3

 

 

 

6. Условие б) теоремы выполнено, так как V 4 ={5}.

 

 

 

x

 

7. Оптимальное управление u(t) находится по x(t) с использованием уравнения (18):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

0,

 

 

 

 

 

 

 

 

t [0,

 

ln 2

 

 

],

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

 

1

 

 

2

 

1

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

(t)= 1,

t (

 

 

ln2

 

 

,

 

[8

+ ln(

 

2

 

) ] ],

2

3

2

5

3

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

0,

 

t (

 

 

[8 + ln(

 

2

 

 

)], 4 ].

 

 

 

 

 

 

3

 

 

 

 

2

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

49

8. Найдем Imin :

4

 

 

 

 

 

 

 

 

t1

Imin = I(

x

,

u

) = [

x

2 (t)

u

(t) 2(

x

2 (t) 8

x

(t))]dt = −2(e4t 8e2t )dt +

0

0

t2

4

 

 

 

 

+[e4t 2(e4t

8e2t )]dt 2(25e4(t4) 8 5e2(t4) )dt =- 619182.

t1

 

 

 

 

t2

Пример 3. Воспользуемся достаточными условиями оптимальности для доказательства того, что найденный с помощью принципа максимума для за-

дачи (1.13)–(1.16)

о максимизации интегрального потребления процесс

(

x

,u

) (см. рис. 6)

является решением этой задачи (при условии T > 1).

1. Зададим функцию W (t,x) равенством

W (t,x) = (T t) x.

Отметим сразу, что

W(T ,x) 0,

ипоэтому условие б) теоремы о достаточных условиях оптимальности выполняется автоматически. Остается проверить условие а).

2.По выбранной W (t,x) построим функцию R:

R(t,x,u) = (T t)xu xu = (T 1t)xu.

3. Для каждого t [0,T 1) функция (T 1t)xu монотонно возрас-

тает по x и по u , поэтому пара (x (t),u (t)) = (x0 et ,1), где каждая компонента принимает максимальное из всех допустимых значений, удовлетворяет условию а).

4. Для каждого t (T 1,T ] функция (T 1t)xu монотонно убывает по x и по u , поэтому пара (x (t),u (t)) (x0 eT 1 ,0), где каждая компонента принимает минимальное из всех допустимых значений, удовлетворяет условию а).

5. Из 3 и 4 следует, что процесс (x (t),u (t)), t [0,T ].

50

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