Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
metods / Теория принятия решений.pdf
Скачиваний:
109
Добавлен:
26.03.2015
Размер:
1.42 Mб
Скачать

3.2. Функционал, его свойства, необходимые и достаточные условия достижения экстремума

Функционал – это числовая функция, определенная на некотором множестве функций. Другими словами, функционал формализует закон, по которому каждой функции X из некоторого множества, или класса, функций ставится в соответствие число J:

t1

 

 

J (X)=

(3)

f0 (X, X,t)dt .

t0

 

 

Здесь X – возможно векторная функция, называемая также кривой или точкой функционального пространства, f0 – интегрант.

Решение задачи оптимизации функционала состоит в нахожде-

нии его экстремума, а также функции X из заданного множества (допустимой области C функционального пространства), достав-

ляющей этот экстремум: J (X)= min J (X) или

J (X)= max J (X).

x C

x C

Рассмотрим основные понятия и свойства функционала на примере одномерного случая X=(x1 )=x:

t1

f0

(x, x,t)dt .

(4)

J (x)=

 

 

 

 

t0

 

 

 

Расстояние между функциями определяется следующим образом:

ρ(x , x

)=

 

 

 

x

x

 

 

 

= max

 

 

x

(t)x (t)

 

.

 

 

 

 

 

 

 

1 2

 

 

 

 

1

2

 

 

 

t [t

0

; t

]

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Функционал называется линейным, если для любых функций x1 и x2 из области C выполняется условие J(x1 +x2 )=J(x1 )+J(x2 ) и

для любого числа α J(αx1 )=αJ(x1 ).

Функционал называется билинейным, если он является линейным относительно x1 при фиксированной x2 и линейным относительно x2 при фиксированной x1.

Билинейный функционал, в котором x1=x2, называется квадратичным.

Функционал называется дифференцируемым в точке x=x0 функционального пространства, если для любого допустимого приращения h функции x0 приращение функционала может быть представлено в виде

25

J (x0 ,h)= J (x0 ,h)J (x0 )= ϕ1(x0 ,h)(x0 ,h) h ,

где h – приращение x0 , ϕ1 (x0 ,h) – линейный относительно h

функционал, причем lim α(x0 ,h)= 0 .

h0

Функционал ϕ1 (x0 ,h) называется дифференциалом или первой вариацией функционала J(x) в точке x0: δJ=ϕ1 (x0 ,h).

Функционал называется дважды дифференцируемым в точке x=x0 функционального пространства, если для любого допустимого приращения h функции x0 приращение функционала может быть представлено в виде

J (x

,h)= ϕ

(x

,h)+

1

ϕ

 

(x

,h)(x

,h)

 

h

 

 

 

,

 

 

 

 

0

1

0

 

2

 

2

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функционал ϕ2 (x0 ,h) называется вторым дифференциалом

или

второй вариацией функционала

J(x) в то кеч x

0 :

δ2 J= ϕ2 (x0 ,h).

 

,h) –

где

ϕ1 (x0 ,h) – первая вариация функционала J(x), ϕ2 (x0

квадратичный относительно h функционал,

lim α(x0 ,h)= 0 .

 

 

 

h 0

 

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

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

Функционал достигает локального экстремума на кривой x0 (t), если разность J(x0 )–J(x) сохраняет знак в некоторой окрестности

x0 , т.е. при x x0 < ε , где x C – области определения функ-

ционала.

Если J (x0 )J (x) для всех x C , то x0 доставляет функционалу абсолютный максимум, если J (x0 )J (x) – абсолютный мини-

мум.

Первое необходимое условие достижения локального экстре-

мума: δJ=0.

условие: δ2 J 0 для минимума или

Второе необходимое

δ2 J 0 для максимума.

Выполнение этого условия проверяется

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

временном выполнении первого и второго необходимых условий при строгом неравенстве во втором.

26

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

3.3. Вариационные задачи на безусловный экстремум

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

Простейшая задача вариационного исчисления имеет следующую постановку: требуется

найти функцию x(t), определенную на интервале

t [t0 ; t1]

при заданных

 

 

 

 

 

значениях

x(t0 )=x0

,

 

 

 

 

 

x(t1 )=x1 , доставляющую

 

 

 

 

 

экстремум

функционалу

 

 

 

 

 

вида (4). Иногда говорят,

 

 

 

 

 

что требуется найти кри-

 

 

 

 

 

вую x(t), концы которой за-

 

 

 

 

 

креплены в точках x(t0 )=x0

 

 

 

 

 

и x(t1 )=x1 (рис. 8).

 

 

 

 

 

Рис. 8

 

Первое необходимое условие достижения экстремума в про-

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

 

 

f0

 

d

 

f0

 

 

 

 

 

 

 

= 0 .

(5)

 

 

 

x

 

 

 

 

 

 

dt

x

 

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера сводится к дифференциальному уравнению второго порядка относительно x(t). Его решение называют экстремалью. Далее с учетом граничных условий x(t0 )=x0 , x(t1 )=x1 находят допустимую экстремаль (одну или несколько).

Второе необходимое условие достижения экстремума прини-

мает вид условия Лежандра:

2 f0

0

для минимума или

 

 

xx

 

 

27

2 f0

0

для максимума вдоль всей допустимой экстремали. Ус-

 

xx

 

 

ловия Лежандра в форме строгих неравенств называются усиленными. Усиленное условие Лежандра вместе с (5) образует достаточное условие минимума или максимума функционала (4).

Пример 14. Обеспечить экстремум функционала J (x)= 1x2dt

0

при граничных условиях x(0)=0, x(1)=1 (рис. 9).

Найдем частные производные и составим уравнение Эйлера:

fx0 = 0 , fx0 = 2x ; dtd (2x)= 0 , 2x = 0 , x = 0 .

Уравнение экстремали получаем двойным интегрированием:

Рис. 9

Пример

J (x)= 1t2 x2dt

0

x = c1 , x=c1 t+c2 .

С учетом граничных условий находим допустимую экстремаль:

x(0)=c2 =0, x(1)=c1 =1, x = t .

Проверим выполнение условия

Лежандра: 2 f0 = 2 > 0 .

xx

Таким образом, на кривой x = t (кривая b на рис. 9) обеспечивается

1

минимум величинойSmin = 12 dt =1.

0

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

при граничных условиях x(0)=0, x(1)=1.

Найдем частные производные и составим уравнение Эйлера:

fx0 = 0 , fx0 = 2t2 x; dtd (2t2 x)= 0 .

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

тегрированию: 2t2 x = c1 , x = 2ct12 , x = − c21t +c2 .

28

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

Пример 16. Определить характеристический полином системы автоматического управления (САУ), для которой функционал

 

 

 

 

 

 

J (x)= (a0 x

2

 

2

)dt

минимален при

a0 , a1 >0 и граничных ус-

 

+ a1x

 

0

 

 

 

 

 

 

ловиях x(0)=x0 , x()= 0 .

Найдем частные производные и составим уравнение Эйлера:

fx0 = 2a0 x , fx0 = 2a1x ; 2a0 x 2a1x = 0 .

Решение полученного уравнения имеет вид: x=c1e-αt+c2eαt, где

α= a0 / a1 > 0 .

Сучетом граничных условий находим c1 =x0 , c2 =0 и допустимую экстремаль в виде x = x0e−αt .

Проверим выполнение условия Лежандра: 2 f0 = 2a1 > 0 .

xx

Найденная допустимая экстремаль обеспечивает минимум рассматриваемого функционала и является решением однородного дифференциального уравнения x + αx = 0 , левая часть которого определяет характеристический полином САУ [2]. Здесь он примет вид D(s)=s+α.

В векторном варианте простейшей задачи при X=(x1 ,x2 ,…,xn ) и функционале (3) задаются граничные условия xi (t0 )=xi 0 ,

xi (t1 )=xi 1 , i=1,2,...,n.

Первое необходимое условие достижения экстремума здесь принимает вид системы уравнений Эйлера-Лагранжа:

f0

 

d

 

f0

 

 

 

 

 

 

= 0 ,

i=1,2,....,n.

(6)

 

 

 

xi

 

 

 

 

 

dt

xi

 

 

 

 

 

 

 

 

 

 

 

 

Второе условие (условие Лежандра) предусматривает анализ матрицы вторых частных производных

29

 

 

 

2 f0

 

 

 

2 f0

...

 

2 f0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1x1

 

 

 

x1x2

 

 

 

 

x1xn

 

 

F..

=

 

2 f0

 

 

 

2 f0

 

...

 

2 f0

 

(7)

 

x

x

 

 

 

x

x

 

 

x

x

 

xx

 

2

1

2

2

 

 

 

 

2

n

 

 

 

...

...

...

 

...

 

 

 

 

 

 

 

 

 

 

2 f0

 

 

 

2 f0

 

 

...

 

2 f0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

x1

 

 

 

xn

x2

 

 

 

 

xn

xn

 

 

и по критерию Сильвестра полностью аналогично анализу матрицы Гессе в задаче обеспечения экстремума функции.

Условие Лежандра (необходимое условие): вдоль всей допустимой экстремали матрица F.x.x должна быть положительно полу-

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

Усиленное условие Лежандра (достаточное условие): вдоль всей допустимой экстремали матрица F.x.x должна быть положи-

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

Пример 17. Обеспечить экстремум функционала

1

2

+ 2x1 )dt

 

2

при граничных условиях x1 (0)=1,

J (x1, x2 )= (x1

+ x2

0

 

 

 

x2 (0)=0, x1 (1)=3/2, x2 (1)=1.

Найдем частные производные и составим систему уравнений Эйлера-Лагранжа:

 

f0

= 2

,

 

f0

 

= 0 ,

 

f0

= 2x1 ,

f0

= 2x2 ;

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

x2

 

x1

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

2

2x1 =

0 , 2x2 = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения экстремали:

x = t

2 / 2 + c t + c

2

, x2 =c3 t+c4 .

 

 

 

 

 

 

1

 

 

1

 

 

С учетом граничных условий находим допустимую экстремаль:

x1 (0)=c2 =1, x1 (1)=1/2+c1 +c2 =3/2, c1 =0,

 

/ 2 +1;

x = t2

x2 (0)=c4 =0, x2 (1)=c3 +c4 =1, c3 =1,

1

 

x2 = t .

 

 

 

 

Составим матрицу (7) и проверим выполнение условия Лежандра:

30

F.. =

 

2

0

 

 

 

, 1 =2>0, 2 =4>0.

 

 

 

 

 

 

 

 

 

xx

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, на кривой

 

 

обеспечивается ми-

 

x = (t2 / 2 +1; t)

 

1

 

 

 

 

 

 

2

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нимум величиной Smin =

(t2 +12

+t2 + 2)dt =

 

t3

+3t

 

0

= 3

 

.

3

3

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В общем случае граничные условия могут быть заданы неполностью. В этом случае первое необходимое условие достижения экстремума (5) или (6) дополняется условиями трансверсальности. Рассмотрим возможные варианты.

1. Задача с подвижными концами: значения всех или некоторых составляющих векторной функции X при t=t0 или t=t1 не заданы (концы частично не закреплены). Пример для одномерного случая, когда левый конец закреплен, правый подвижен, представлен на рис. 10. В таких задачах первое необходимое условие достижения экстремума вместо каждого недостающего граничного условия дополняется условием трансверсальности вида

f0

t = t0

= 0 или

f0

= 0 .

(8)

 

 

xi

 

xi

t = t1

 

2. Задача со свободными концами: не задана также граница временного интервала t0 или t1 (рис. 11 для одномерного случая, правый конец свободен). Вместо неопределенной границы первое необходимое условие достижения экстремума дополняется условием трансверсальности:

 

n

f0

 

 

 

 

 

n

f0

 

 

 

 

 

f0

 

xi

 

= 0 или

f0

 

xi

 

= 0 .

(9)

 

 

 

 

i=1 xi

 

 

t = t0

 

 

i=1 xi

 

 

t = t1

 

 

 

 

 

 

 

3. Задача со скользящими концами: вместо значений t0 и x0i за-

даны уравнения xi (t0 )=ϕi (t0 ), i=1,2,…,n, и аналогично для правого конца (пример для одномерного случая на рис. 12). Вместо условий (8), (9) вводится условие трансверсальности:

 

 

 

n f0

 

 

 

 

 

 

 

 

 

 

f0

 

 

 

(xi − ϕi )

 

 

= 0

 

 

 

 

 

 

 

 

i=1 xi

 

 

 

 

 

 

t = t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

f0

 

 

 

 

 

 

 

 

или

f0

 

 

(xi − ϕi )

 

= 0 .

(10)

 

 

 

 

 

i=1 xi

 

 

 

 

 

t = t1

 

 

 

 

 

 

 

 

31

Рис. 10

Рис. 11

Рис. 12

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

уравнения xi (t0 )=ϕi (t0 ) также рассматриваются как граничные условия и используются для нахождения коэффициентов экстремали вместе с прочими граничными условиями и условиями (10).

Пример 18. Обеспечить экстремум функционала

J (x)= 1(x2 + x)dt при граничном условии x(1)=0 (рис. 13).

0

В данной задаче левый конец подвижен. Поэтому потребуется составить и учесть одно условие трансверсальности (8).

Найдем частные производные и составим уравнение Эйлера и условие трансверсальности:

32

 

 

f0

=1 ,

 

f0

= 2x ;

 

 

 

 

 

 

 

 

 

x

 

 

 

 

x

 

 

 

 

 

 

 

 

1dt

(2x)= 0 , 1 2x = 0 ;

 

d

 

 

 

 

 

 

 

f0

 

 

 

 

 

 

 

 

 

= 2x(0)= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

x

t = 0

 

 

 

 

 

 

 

 

Получим

 

 

уравнение

экстремали:

x = t / 2 + c1 ,

 

x

= t

2

/ 4 + c1t + c2 .

 

 

 

 

 

 

 

 

 

 

 

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

играничного условия находим допус-

тимую

экстремаль:

x(0)= с1 = 0 ,

x(1)=1/ 4 + с2 = 0 ,

 

с2 = −1/ 4 ,

 

 

Рис. 13

x = (t2 1) / 4 .

Проверим выполнение условия Лежандра: 2 f0 = 2 > 0 .

xx

Таким образом, на кривой x = (t2 1) / 4 (кривая b на рис. 13) обеспечивается минимум величиной

 

 

1

t2

 

t2 1

 

1

 

 

 

 

Smin =

 

 

+

 

 

 

.

 

 

 

 

4

4

dt = −

12

 

 

 

0

 

 

 

 

 

Пример

19.

Обеспечить

экстремум

функционала

J (x)=T(x2 x +1)dt при T>0 и граничном условии x(0)=0.

0

В данной задаче правый конец свободен. Поэтому потребуется составить и учесть для правого конца два условия трансверсальности: одно (8) и одно (9).

Найдем частные производные и составим уравнение Эйлера и условия трансверсальности:

fx0 = −1, fx0 = 2x ; 1dtd (2x)= 0 , 12x = 0 ;

f0

 

= 2x(T)= 0 ;

 

 

x

t =T

 

 

33

 

f0

 

 

 

 

 

 

2

x(T)+1

 

2

 

f0

 

x

 

= x(T)

2x(T) =1 x(T)= 0 (с учетом

 

x

 

 

t =T

 

 

 

 

 

 

 

 

x(T)= 0 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнение

экстремали:

x = −t / 2 + c1 ,

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = −t2 / 4

+ c t + c

2

. Далее кроме неизвестных констант c1 и c2 на

 

 

 

1

 

 

 

 

 

 

 

основе граничного условия и условий трансверсальности находим

также значение T: x(0)=c2 =1,

 

x(T)= −T / 2 + с1 = 0 , с1 =T / 2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x(T)=1 +T 2 / 4 T / 2 T = 0 , T=2, c1 =1,

 

 

 

x = t t2 / 4 .

Проверим выполнение условия Лежандра:

 

2 f0

= 2 > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xx

 

Таким образом, на кривой

 

 

 

 

 

 

на интервале [0; T]

 

x = t t2 / 4

обеспечивается минимум величиной

 

 

 

 

 

 

 

 

 

2

 

t

 

2

 

t

2

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

.

 

Smin = 1

 

 

 

t +

 

 

+1 dt =

 

 

2

 

4

3

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Пример

20.

Обеспечить

 

экстремум

 

функционала

J (x)=T(x2 + x)dt при граничном условии x(T)=T, T>0 (рис. 14).

0

В данной задаче левый конец подвижен, правый конец – скользящий по кривой ϕ(T)=T. Поэтому потребуется составить и учесть для левого конца условие трансверсальности (8), для пра-

вого (10).

Найдем частные производные и составим уравнение Эйлера и

условия

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

f0

=1

,

f0

= 2x ;

1

d

(2x)= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

=

1

;

f0

 

 

= 2x(0)= 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

x

t = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x − ϕ)

t =T

= x(T)

 

+ x(T)+ 2x(T)(1 x(T))=

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −x(T)

 

+ x(T)+ 2x(T)= 0.

 

 

 

 

 

 

 

 

 

 

 

34