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

Дегтярёв Оптимальное управление

.pdf
Скачиваний:
63
Добавлен:
22.03.2016
Размер:
1.68 Mб
Скачать

 

11

 

 

J 0

(1.11)

Для доказательства предположим, что кривая

y(t) доставляет минимум

(максимум) функционалу J [ y(t)]. Это означает,

что при любых y(t)

J 0

( J 0) . Если теперь предположить, что условие (1.11) неверно, то в

силу линейности

J от

y(t) первая вариация J , а значит и J будут менять

знак при смене

знака

y , что противоречит условию J 0

( J 0) и

доказывает справедливость условия (1.11).

Условие (1.11) еще называют условием стационарности. Поскольку оно является необходимым условием экстремума, ему должны удовлетворять все кривые, доставляющие и минимум, и максимум функционалу J [ y(t)], т.е. это условие позволяет выделить все множество функций, среди которых находится и искомое, если оно существует. Используя условие (1.11), можно получить уравнения, которым с необходимостью должны удовлетворять искомые экстремали.

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

Лемма. Пусть (t) - заданная непрерывная функция, удовлетворяющая уравнению

 

t1

 

 

 

 

(t)h(t)dt 0 ,

 

(1.12)

 

t0

 

 

 

где h(t) - любая непрерывная функция, то (t) 0 .

 

Доказательство. Предположим, что условие (t) 0 не выполняется и

существует некоторая точка t* , в

 

 

которой (t* ) 0 . Для определенности

Ф h

 

 

 

 

 

Ф(t)

положим, что (t* ) 0 . Тогда в силу

 

h(t)

 

 

 

 

непрерывности

функции

(t)

 

 

существует в

окрестности

точки t*

t0 t

t t t1

 

 

 

Рис. 5

12

некоторый интервал (t* , t* ) , где (t) 0 .

Поскольку h(t) - любая непрерывная функция, выберем ее следующим

образом

 

при

0

h(t)

 

 

при

0

В этом случае

t (t* , t* ), t (t* , t* ).

t

t *

(t)h(t)dt

(t)h(t)dt 0 ,

t0

t *

что противоречит условию (1.12) и доказывает лемму.

 

 

 

 

2.4. Простейшая задача вариационного исчисления, уравнение Эйлера.

 

Рассмотрим простейшую задачу вариационного исчисления. Среди

непрерывно-дифференцируемых

скалярных

функций

сравнения

y(t) ,

удовлетворяющих

 

 

 

 

условиям

 

 

 

 

y(t0 ) y0 , y(t1) y1

(проходящих

через

y

 

 

 

 

B(y1,t1)

 

 

 

 

 

 

 

 

 

 

 

точки

A(t0 , y0 и

B(t1, y1)), найти

такую,

 

 

 

 

которая доставляет минимум функционалу

A(y0,t0)

 

 

 

y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.13)

 

 

 

 

 

J F (t, y(t), y (t))dt ,

t0

t1

t

 

 

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(t, y(t), y (t)) - достаточно гладкая

Рис. 6

 

 

 

где

 

 

 

 

функция своих аргументов.

 

 

 

 

 

 

 

Для решения задачи используем условие (1.11). Для этого приращение

функционала (1.13)

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J F (t, y(t) y(t), y (t) y (t)) F (t, y(t), y (t)) dt ,

 

 

 

t0

 

 

 

 

 

 

 

 

 

 

воспользовавшись формулой Тейлора, представим в следующим виде

 

 

 

 

t1

 

F y

F

 

 

 

(1.14)

 

J

 

y dt 0( y) .

 

 

 

t

0

 

y

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

где 0( y) - содержит величины более высокого порядка малости, чем

.

Из выражения (1.14) следует, что

t1

 

F

y

F

 

J

 

y dt.

t

0

 

y

 

y

 

 

 

 

 

 

 

y и y'

(1.15)

Поскольку решение мы ищем в классе непрерывно-дифференцируемых функций, то и вариации y(t) тоже будет непрерывно-дифференцируемыми.

Нетрудно показать, что справедливо равенство

d y

dy

 

 

 

.

 

 

dt

dt

Учитывая (1.16), представим

F

y

d

 

F

 

 

d

 

F

 

 

y

 

 

y .

 

 

 

 

 

y

 

 

 

y

 

 

 

 

 

 

 

dt

 

 

dt

y

Тогда, подставляя (1.17) в (1.15) и интегрируя, получим

 

t1

F

 

d

F

 

F

 

t1

 

 

 

J

 

 

 

 

 

 

y dt

 

 

y

.

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

t0 y

 

dt y

 

 

 

t0

 

 

 

 

 

(1.16)

(1.17)

Так как все функции

сравнения

проходят через заданные точки

A(t0 ,y0 ) и B(t1, y1) , то y(t0 ) y(t1) 0 и

 

t1

 

F

 

d

 

F

 

J

 

 

y

 

y d t .

 

 

 

 

 

 

 

 

 

 

 

t0

 

y

 

dt

 

y

Теперь в силу произвольности y(t) из условия J 0, применяя

основную лемму вариационного исчисления, получим уравнение экстремали

F

 

d

 

F

 

 

 

 

 

 

 

0

,

(1.18)

 

 

 

y

 

 

 

 

 

 

 

 

 

dt

y

 

 

 

которое называется уравнением Эйлера.

Уравнению Эйлера удовлетворяют кривые, доставляющие как относительный минимум, так и относительный максимум функционалу (1.13),

14

поэтому для идентификации экстремума необходимы дополнительные

исследования.

Учитывая, что F F(t, y, y'), уравнение Эйлера можно представить в виде

F

2 F

 

2 F

 

2 F

 

 

y t y

 

y 2 y

0 .

y y y

 

То есть уравнение Эйлера – это дифференциальное уравнение второго порядка относительно y(t) , его общее решение содержит две произвольные постоянные

c1

и c2 ( y(t) y(t,c1,c2 )) ,

для нахождения которых используются

условия

прохождения искомой

экстремали

через

заданные

точки

A и

B,

( y(t0 ,c1,c2 ) y0 , y(t1,c1,c2 ) y1) .

Таким

образом,

решение

задачи

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

Рассмотрим некоторые частные случаи уравнения Эйлера.

 

 

1. Функция F

не содержит

y , т.е. F F (t, y') . При этом

F

0 , и

 

 

 

 

 

y

 

уравнение Эйлера запишется в виде

 

 

 

 

 

d

F(t,y') 0 .

 

 

 

 

 

 

 

 

 

dt

y'

 

 

Откуда следует

F (t,y') c , где c - произвольная постоянная. Разрешая

 

y'

1

 

1

 

 

 

 

 

 

 

 

это уравнение относительно y' , получим:

y' f (x,c1) .

t

Интегрируя его, найдем y f (t ,c1) dt c2 ) , c2 - произвольная постоянная.

0

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

15

2. В функцию F не входят x и

y .

Тогда

F

 

F

0 , и уравнение

 

 

 

 

 

t

 

y

 

Эйлера преобразуется к виду

F( y') c ,

или

y' f (c ) const . В этом случае

 

y'

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

y f (c1)t c2 . Экстремали представляют прямые линии.

Пример 1. Среди всех гладких кривых, соединяющих две заданные точки

A(t0 , y0 ) и B(t1, y1) найти такую, которая имеет минимальную длину. В этом случае

t1

dy

2

 

 

 

 

J L

 

 

 

 

 

1

 

 

dt .

 

 

 

 

 

 

 

 

t0

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подынтегральная функция

F

1 y'2 не содержит

t

и y ,

а искомая

экстремаль по только что

доказанному представляет

прямую

y c1t c2 ,

проходящую через заданные точки. Постоянные c1 и c2

находятся из системы

уравнений

 

 

 

 

 

 

 

 

 

y0 c1t0 c2 ,

y1 c1t1 c2 ...

Пример 2. Рассмотрим обтекание тонкого профиля (рис. 7) в

линеаризованном сверхзвуковом потоке [2]. Пусть передняя кромка находится в начале координат (t0 0, y0 0) , а задняя – в точке (t1 a, y1 h), y(t) - ордината точки профиля. Тогда волновое сопротивление профиля можно представить в виде [2]

a

J K y'2 d t , где K - коэффициент пропорциональности, зависящий от

0

плотности и скорости набегающего потока.

Так как F (t, y, y') не зависит от t и y , уравнение экстремали имеет вид y(t) c1t c2 . Удовлетворяя граничным условиям, найдем y h a t .

16

y

Это отрезок прямой линии. Но какой

 

экстремум реализует эта прямая остается

h

неизвестно. Необходимы дополнительные

v

исследования.

a

t

Замечание. При

выводе

уравнения

 

 

Рис. 7

 

Эйлера предполагали,

что y(t)

- скалярная

 

 

 

 

функция. В том случае, когда

y(t) - есть

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

функция

имеет вид

F F (t, y1,... yn , y1',... yn ') ,

а граничные точки A и B задаются в (n 1) - мерном

пространстве yi (t0 ) yi0 , yi (t1) yi1 . В этом случае

 

 

 

t1

n

 

F

 

 

n

F

 

 

J

 

yi

yi d t .

 

 

yi

 

 

 

t0

i 1

yi

i 1

 

 

Используя необходимое условие экстремума J 0,

и основную лемму

вариационного исчисления, получим систему n уравнений второго порядка

 

F

 

d

 

 

 

 

 

 

 

 

 

 

 

F

0

 

i 1,....... n.

(1.19)

 

 

 

 

 

 

 

 

 

yi

 

d t

 

yi

 

 

 

 

Общее решение этой системы содержит 2n констант, для определения которых имеет 2n условий.

2.5.Уравнение Эйлера-Пуассона.

Рассмотрим задачу поиска экстремума функционала

 

t1

 

 

 

 

I F (t, y, y', y'') dt ,

 

 

t0

 

 

 

y(t0 ) y0 ;

 

;

y(t1) y1;

 

y'(t0 ) y0

y'(t1) y1 .

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

(граничным) условиям.

17

Первую вариацию

t1 F

I

y

t0

учитывая формулы

F y'y'

y

d

dt

F

y'

Fy' y

 

F

 

y'

 

y''

 

 

y''

 

 

 

 

 

d

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

y'

 

d t ,

y ,

 

F

 

 

d

 

F

 

 

 

 

 

 

d

 

 

d

 

 

F

 

 

d

2

 

F

 

 

y''

 

y'

 

 

y

 

 

 

 

 

 

 

y ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y''

 

 

 

 

y''

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

dt

2

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

dt

y''

 

 

 

y''

и y y' 0 при

t t0

и t t1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

преобразуем к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

F

 

d F

 

d 2

F

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y dt 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt2

 

 

 

 

 

 

 

 

 

 

 

a

y

 

dt y'

 

y''

 

 

 

 

 

 

 

 

Применяя основную лемму вариационного исчисления, получим Уравнение Эйлера – Пуассона:

F

 

d F

 

d 2

 

F

 

 

 

 

 

 

y

dt y'

dt2

 

y''

 

 

 

Это уравнение четвертого порядка.

констант, которые находятся из условия заданные граничные точки.

0 .

Решения его зависит уже от 4-х

прохождения экстремали через

2.6.Вариационные задачи с нефиксированными границами. Условие трансверсальности.

При рассмотрении простейшей задачи вариационного исчисления предполагалось, что граничные точки A и B являются заданными. В

прикладных задачах эти точки могут быть не фиксированы, т.е. t0 , y0 ,t1, y1

могут быть свободны все или частично. Можно выделить, например,

следующие случаи:

1. Координаты точек A и B свободны (t0 , y0 , t1, y1-не заданы).

18

2.Точка A – задана, точка B свободна ( t0 , y0 - заданы, t1, y1 - не

заданы).

 

 

 

3.

Точка A задана, координаты т. B

связаны зависимостью f (t1, y1) 0.

4.

Координаты точек

A и B

связаны совокупностью условий

 

f j (t0 , y0 ,t1, y1) 0

j 3

 

Могут рассматриваться и другие варианты.

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

Как и ранее, найдем приращение функционала (1.13). При нефиксированных граничных условиях J

 

t1

t1

t1

 

J

F t, y y, y y d t F t, y, y d t

 

t0 t0

t0

 

представим в следующем виде

 

 

 

t1

F t, y y, y y F t, y, y d t

J

 

t0

 

 

 

t1

t1

t0

t0

 

 

F

t, y y, y y d t

F

t, y y, y y d t ,

 

t1

 

t0

 

где t0 , t1 - независимые вариации параметров t0 и t1 .

Воспользовавшись формулой Тейлора и теоремой о среднем, получим

t1

 

F

y F

 

 

 

 

 

 

 

 

J

 

y d t F (t, y y, y y

 

t

 

 

t

 

t1

 

 

y

y

 

 

1

 

1

 

1

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

 

d F

 

 

F (t, y y, y y

 

t0 1

 

F

 

y d t Fy

 

 

 

 

 

 

 

 

t0 2 t0

 

t

0

y

dt y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

F (t, y, y )

 

t

t1

2

F (t, y, y )

 

t

 

t0 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

0

 

 

 

 

 

 

 

y

 

t1

 

(1.20)

 

 

 

t0

 

 

 

 

 

 

где 0 i 1, 1, 2 , 3 - величины более высокого порядка малость чем

t0 , t1, y(t), y'(t) .

 

 

 

 

 

 

 

 

 

 

 

19

 

 

 

 

 

 

 

 

 

 

Поскольку искомая экстремаль при нефиксированных граничных точках

должна быть экстремалью и среди кривых сравнения, имеющих с нею общие

граничные точки, она должна удовлетворять уравнению Эйлера и должно

выполняться уравнение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

d F 0

 

 

 

 

 

 

 

 

 

 

(1.21)

 

 

 

 

 

 

y

 

dt y'

 

 

 

 

 

 

 

 

 

 

 

Тогда из (1.20) следует, что

 

 

 

 

 

 

 

 

 

 

 

J

F

y(t1)

F

 

y(t0 )

 

 

 

 

 

 

 

 

t0

(1.22)

 

 

 

 

F (t, y, y )

 

t1 F (t, y, y )

t0

 

 

 

y t

 

 

 

y t

 

 

 

 

t1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В выражение (1.22),

кроме независимых вариаций t0 , t1

входят, так

называемые

изохронные

вариации

y(t1), y(t0 ) ,

которые

 

не

являются

независимыми. Выразим изохронные вариации через приращения конечных

значений y0

 

и y1 . Поскольку эта зависимость нелинейная, то линеаризируя

ее, получим (рис. 8).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

y(t1) y

t1

t1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0

y(t0 ) y

t

t0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

 

 

 

 

y(t1) y1 y t1

t1,

 

 

 

 

y(t1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(t0 ) y0 y

t

t0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

t1

 

 

 

t1 t1

t

 

 

 

Используя

 

 

 

полученные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 8

 

 

 

 

 

соотношения, найдем

из

(1.22) и

(1.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

F

y1

 

F

 

 

 

 

 

 

F

 

 

 

 

 

 

 

y t

y t

 

y0 F (t, y, y ) y

y

t1

 

 

 

 

 

 

 

 

 

0

 

 

 

 

t

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

(1.23)

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t0 0.

 

 

 

 

 

 

 

 

 

 

 

F (t, y, y ) y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y t

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

Полученное условие может быть использовано для определения оптимальных значений координат граничных точек. Оно называется условием

трансверсальности и используется совместно с уравнением Эйлера.

Рассмотрим некоторые частные условия трансверсальности.

1. Точки A и B заданы. В этом случае

t0 t1 y0 y1 0 , условие

(1.23) выполняется тождественно и дополнительной информации не дает.

2. Координаты t0 и t1 , заданы, а

y0

и

y1 произвольны. В этом случае

t0 t1 0 и из

условия

трансверсальности получаем

дополнительно два

условия для определения y0

и y1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

0

F

0

 

 

 

 

 

 

 

 

 

y

 

 

y

 

 

 

 

 

 

 

 

 

t

0

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3. Точка A

задана,

а координаты

точки

B связаны

зависимостью

f (t1, y1) 0 . В этом случае

t0 , y0 0 ,

а

вариации

t1 и

y1 связаны

зависимостью

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

f

 

 

 

 

 

 

 

 

f

t

 

 

 

t

 

y 0

 

 

y

 

 

1

t

,

 

 

 

 

 

 

 

 

 

t1

1

 

y1

1

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

f

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

f '

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (t,y,y') y' F

 

 

F

 

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'

t

y'

 

f ' y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В том случае, когда y(t) – n – мерный вектор достаточно гладких функций, искомая экстремаль удовлетворяет системе уравнений Эйлера

F

 

d

 

F

0,

i 1, 2,..., n

(1.24)

yi

dt yi

 

 

 

 

а условие трансверсальности имеет вид: