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

книги из ГПНТБ / Цирлин, А. М. Основы оптимального управления конспект лекций

.pdf
Скачиваний:
22
Добавлен:
19.10.2023
Размер:
6.17 Mб
Скачать

so

6.3. В невыпуклой

задаче;

 

 

 

1. Множество

сравнения

L выбирается

таким, чтобы на нем

максимизируемую

функцию и функции, определяющие

множество Л,

можно

было линеаризовать.

 

 

 

2.

Множество

вариаций,

приводящих к

росту J-0

, и множе­

ство допустимых вариаций для линеаризованной задачи представ-

яяют собой выпуклые конусы, имеющие обцую зернину

У

и не

имеющие общих внутренних точек. Значит должна существовать

разделяющая

их

гиперплоскость (см. 2.2). Условие

ее существо­

вания и формируется как условие существования i -

множителей

со свойствами

(6 . 7) .

 

 

Обычно

Ь

есть линь подмножество-JD. Это означает,

что

условия оптимальности выделяют много ''претендентов0 на реше­

ние, а

набор инозителей

Л

не

единственен.

 

 

Исключение составляет

случай,

когда множество

Л и макси­

мизируемая функция J-o

строго

выпуклы. Л ля этого нужно,

чтобв

функции

j / t '

были

линейны, а

ограничения (6.2)

выделяли

выпую

жое множество

( F j

- выпуклы).

 

 

 

6.*. Особый случай.

Впредыдущем параграфе мы останавливались на особых случа­ ях в задаче о максимуме функции при наличия связей. Эта особые сдучаз касались точек в пространстве X, в которых градиент функции ' j / j или максимизируемой функции обращался в нуль. Аналогичная ситуация возникает а в задаче с ограничениями.

Однако здесь возможны изолированные точки несколько иного рода,

Аименно, точки, лежащие на границе Л, и такие, что движение ив нзх внутрь допустимой области возможно по изолированному

направлению. Таков случай изображен на рис.6.3, где линии Р | я О и Р = О касаются з точке У * . Градиенты этих

SI

функций линейно

зависимы

 

, и уаловие

(6.8а) не выполнено ни

при

каких

конечных

значениях j\ .

 

 

 

 

 

 

 

 

Случай,

когда

в

У * градиент

равен нулю, совершенно

 

аналогичен соответствующему случаю в 5.5. Всюду

ниже мы будем,

часто не оговаривая этого, предполагать,

что условия

регуляр­

ности

(общности

положения)

выполнены.

 

 

 

 

 

 

 

 

6.5. Седловая точка функции Лагранжа

 

 

 

 

 

 

 

Будем рассматривать малую выпуклую окростносвь Е

точки X

 

такую,

для которой

функцию

/ „

и пересечение

L

множеств £ и

Л

можно было бы считать выпуклыми. Таким

образом,

задача об у с лоз­

ном максимуме

j.0

на Е имеет

единственное

решение. Покажем,что

функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R =

) -+ J. /

 

ft

 

(6.9)

 

имеет

в

X * седловую

точку.

Обозначим

черев

6** функции

 

 

Й » ( Л ) = , т у

й ( х , Л

 

 

 

 

 

 

( 6 . 1 0 )

'

Решение задачи (6.103 зависит от значений

Ji-множителей

и не

 

всегда

является

допустимым

по условиям

(6,1) н (6.2),

то еоть

 

X*

( j f )

не всегда

принадлежит

L ~ В Л13

 

 

 

 

 

 

_Соглаоно теореме Куна-Таккера

найдется

такое

значение

 

J

= j(

 

Для которого

X *

( Jf )

 

£

L

. При двбоы j /

 

 

 

 

 

 

 

 

хьу

 

 

 

 

 

 

 

 

(б . п )

 

так

как

L

£

Е

. При j l

=

ji

 

это неравенство превращается

в равенство,

причем

всюду

i

L ,

$

( Л ) не завиоит

от с//

 

и равен

^Д-Х ^„ (X j так как два последних

слагаемых в

(6.9)

 

обращаются в__нуль._Таким образом}

К

( «/»

) достигает

мини­

мума при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

$*(IV

 

=•

m

^

5 ? * *

Q№

 

7)

 

 

 

(6-I2)

 

SI

 

j

Здесь множество вначений

J , на котором имеется ыинямун,

определено условиями (6.7).

Говорят,

что функция (?

 

имеет седловую точку (рис.6.4), если

для любых Х€Е

и j/

справедливы неравенства

 

 

 

G(k}3f*)±QC**J*)*

R(**,J).

(б.п,

6.6. до оих пор речь кжа о седловой точки функции R в неко

торой, вообще говоря, малой окреотности Е предполагаемого реше­ ния. Однако мы можем раооирить ату окрестность до' всего множе­ ства X или любого его подмножества. Действительно, пусть можяо указать некоторое достаточно широкое мноаеотво V заведомо охватывающее Д. Тогда

л о в

/

 

Rt*J)^tp*xR(z,J)-™**Mv

X£V*

_

_

***

(6.14)

При

j l

<//

это неравенство обращается в равеяотво,

с точке макоиыума еоть

искомое решение. Причем решение, обес­

печивающее нелокальный, а абсолютный максимум

J-0 за J .

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

 

 

 

 

Иначе говоря, волн найдется такой вектор

 

jj' «*

,

450 каксимуи

функции &

п о

I

на V

достигается

в точ-

к®

€ 7} ) 1

0 8 - 2 8

*Ьчка

является

сэдловой,

и функцня

 

&*С?) ' тЯ*

Л) иишшальиа 3

ней по J .

 

 

Для малого

нноаества

сравнения

и

и

_малой окрестности

Е

Ередполагаеного

решения такое

значение

^

должно найтись

s регулярном

случае.

 

 

 

 

 

 

 

 

 

 

 

 

 

„ 3 3

 

 

 

 

 

 

Senepb жо, когда

речь

идах об абсолютном максимуме, ньт

квкаких

гарантий того, что

Jf

 

отыщется. Тон не менее, неравенство

(6.14) может

оказаться

полезным. Еолн найти достаточно

проотуп

з не олишком отличную от J) область V ,

в которой легко

опреде­

лить абсолютный иакоимум

Q

по <У ,

при некотором фиксирован­

ной

значении

<у/

,

то полученное

число даот верхнюю оценку

ножо­

вого

решения.

 

 

 

 

 

 

 

 

 

 

 

Теперь

с помощь» определенной

процедуры можно так назначать

J j

чтобы этот

накоимун уменьшался, о эь-ачит верхняя

оценка

становилась

вое ближе и олиже к решению.

 

 

 

В некоторых задачах

получение верхней оценки яэдязтон доста­

точным. Так, в задачах

выбора

оптимального режима действующих

аппаратов, где некоторый допустимый (н обычно не самый плохой)

режим уже найден, бливооть верхней оценки ревення к фактнчеом

достигнутому

вначешш

целевое функции говорят о нецедеоообраа-

нооти детального

решения

задачи.

 

 

 

 

 

Эанетим,

что функции

R

можно записать в неокодьао

бонее

общей фошо

 

 

 

 

 

 

 

 

 

 

 

 

где

jj- ( X ) в

(X ) - функция, завноящяе от одной али нескожь*,

КИЕ

составляющих

вектора X и

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

неравепотваи

 

(6.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

на Э, как в прежде,

расишренная функция R

равна

J.e (X), следовательно, неравенство (6.14) справедливо

н идя этой

функции

 

 

_

 

 

 

 

 

 

 

 

 

 

rnay

£

 

 

 

 

=? прах

X? ^rnojc

J

.

,

5Напомним, что мы не оговаривали здесь никаких условий нек- •дшишпшн i чау та jnt» ч/i

оерывности или гладкости 1Г н f7'-</4 .

Si

Вовсе не осязательно учитывать нов условия задачи, вводя соответствующие множители в функцию Лаграняа.

Множеотво Д представляв? сооой пересеченно множеств V- ,

выделяемых каждый из условий. Цусть одна группа условие выделяет HHOseosBoVg. Тогда ножво оорвэовать функцию Леграну

В

• J. + Js

?r(*J

и искать

максимум этой

функцяа ва 7 j . Dps фиксированной Jr

пожучим вврхвош оценку

решошя.

0"с. 6./

Рис. 6.2.

Ь*0

Рос. 6. 3

ss

§ 7. Алгоритмы численного решения

задачи об условном

маноим.7И8

УНКЦИИ

 

 

Чи о ленные алгоритмы

нахождения условного максимума можно рвя-

 

Ф

 

 

бить на две группы: поисковые методы,

использующие для

поотрое-

ляя улучшающей последовательности информация о поведении

функции

э локальной окрестности точки поиска, и методы глобальные, осно­ ванные на получении оценок реиенил в последовательном их ужучм-

нян.

 

Методы

поиска условного

максимума

 

 

 

Рассмотрим

первоначально

задачу

поиска

максимуыа

функции

J-0QL) при

наличии

условие

типа

равенотв

 

 

 

 

 

 

 

 

 

J-i(X)*0

 

 

 

 

 

 

(7.1)

 

 

 

 

 

 

 

 

 

 

1=

i , 2 . . .

m.

я ограничениях,

наложенных

линь

на

отдельные составляющие J( .

Все методы решения этой задачи сводятся к движении э направ­

ленна

условного

градиента

{

1

).

 

 

 

 

7 . 1 . Метод проектирования градиента

{

49 3

 

 

Пусть

начальная

точка поиска

J(0

удовлетворяет

уоювияа

(7 . 1) . Б общем случае, чтобы ее вайтн, нужио ревить вопомота-

тедьную

задачу

минимизации

функция

 

 

 

 

 

ят прн

фиксированных значениях

(Г\ - Тг\ ) составляющих X оп­

ределить

остальные

ТА составляющих

иа

системы

(7 . 1) .

5

точке

XQ

 

находится направление

условного

градиента

j ^ (

X ) ,

т . е .

направление

наиоолеэ быстрого роста

целевой

функции в линейном подпространстве, касательном к поверхноств,

определяемой связями (7.1).

При достаточно малом жаге я найденном1

направлении

в точке

не слипком сильно нарушены уравнения

свявей. Поэтому

и з н е ё о помощью,

например,

процедуры минимизаций

56

(7.2) можно вернуться > допустимую область, найдя точку I j , удовлетворяю^» связям, ш сделать следующий шаг поиска. Поо. ледовательность поиска приведет к локальному условному мавд|

«уму

 

 

/ о ( Х ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Остановимся неоколько подробнее на определении направле­

ния уоловного градиента. Линейное подпространство,

касатель­

ное к поверхности (7.1), обладает тем овойотзом, что любой

вектор

 

£

,

лежащий

в нем,

нормален

к любому из

градиентов

функций

J.^

в

точке

 

Х0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

-

i,

2, . . .

v n

 

 

 

Будем

считать градиенты

 

^ .

линейно

независимыми, так

что

число уравнений

( 7 . 3 )

 

равно

Тп .

 

 

 

 

 

 

 

Так как нао интересует только направление наискорейшего

роста

j.0

,

 

нормируем вектор

£

потребовав,

чтобы

 

 

 

 

 

 

 

 

ft

"

'

1

 

 

 

 

 

 

 

 

 

 

Л )

Величина

 

пронвводной функции

Jc

 

по

направлению

t?

оп­

ределяется как скалярное произведение градиента

 

/ 0

 

на

век­

тор

i

 

 

 

 

 

 

 

 

 

 

 

 

_

 

 

 

 

 

 

 

Таким

образом,

необходимо

найти

такой

У\ -мерный вектор &%

который

доставляет

максимум

 

( 7 . 5 ) при

(W- )-OM условиях

( 7 . 3 )

и

( 7 . 4 ) .

Зто

задача

условного

экстремума,

облегченная

 

 

 

 

 

 

 

Z

 

 

 

 

 

тем

обстоятельством,

что выражения

( 7 . 5 ) и

1 7 . 3 )

линейны, а

( 7 . 4 )

квадратично

зависят

от

 

 

.

Таким

образом,

решение

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

s?

7.2. Локальный метод исключения зависимых составляющих

Разобьем

П. составляющих вектора

У

на

гл.

)

свободных

У с

и

УЛ - зависимых

У 9 .

'Гак что

при

задан­

ных переменных первой группы,вторые могут быть определены

ив уравнений

(7 . 1) .

 

 

 

 

 

Запишем

функцию

Лагранжа

 

 

 

 

R - 1Л*) +24-A-W,

 

 

'

(7.6)

 

 

 

 

 

I -

'

 

 

 

 

 

 

 

 

 

содержащую

ТУ]

добавочных переменных

^//- .

 

 

 

Пусть

в

начальной

точке

У0

 

, как

и раньие,

удовлет­

ворялись

условия (7 . 1), а

значит

Q

к

)

=

/ 0

0).

Вдоль

поверхности

(7.1)

второе

слагаемое

в (7.6)

равно

кулю, а значит

градиент

V

R

вдоль

этой

поверхности равен

уоловноыу

градиенту

 

. Выберем

теперь

 

jf«

так,

чтобы

функция

R.

в

точке

^

не

зависела

 

от составляющих J<^

о Изо

 

 

 

-7,

 

Ъ**о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ = 1. 2 , . . . l r w

вычисляем

составляющие

градиента

 

 

по

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

^

5 5

 

I ,

2 . . . .

( Г \ . - Щ ) ) иQ

изменяем

каждую И8

них.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* г ^ У « / Х . »

 

 

 

( 7 ' 8 )

где £ - величина шага. Найдя новые значения свободных

составляющих,. определим зависимые составляющие из уравнений

зованы в окрестности д(„

, и переменные d(s найдены из сис­

темы линейных уравнений

(19).

Отметим, что при поиске условного максиму!-», основанном на использовании функции Лагранжа, все переменные входят в R сим­ метрично, поэтому им удобнее пользоваться, когда некоторые из переменных содержатся только в уравнениях связей и отсутствуют

вцелевой функции.

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

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

зависимые

составляющие

затрудняет решение.

 

В этом случае удобен алгоритм, изложенный ниже, где все

составляющие вектора -X

равноправны.

 

7;3.

Поиск седаовой

точки функции Лагранжа

[ 22 ]

Представление о решении задачи условного максимума как о

седаовой

точке функции

Лагранжа полезно не только потому, что

приводна

к возможности

вычислять оценку искомого

решения.

Седдовуи точку можно находить поисковыми методами. Процедура поиска сходится к абсолютному максимуму J-0 на и лишь при вы­ полнении условий выпуклости.В остальных случаях, как и во всех поисковых нетодах, шкно лишь гарантировать, что полученное ре ­ шение не хуже других, допустимых и лежащих пососедству.

 

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

функции Лагранжа

I?

по

JC

черев

,

а

составляющую

градиента этой

функции

no

J

через

£ j

.

При поиске

седловой точки мы должны

двигаться

в направлении роста

/2

 

при изменении

У

и

уменьшении

Q

при изменении

J

,

что приводит к

оледующе -

цу

алгоритму

 

 

 

 

 

 

 

 

где величина шага £ положительна.

Возможны и различные модификации этого алгоритма, обеспе­ чивающие не окольно более быструю сходимость. Так, еоли в ок­ рестности максимума зависимость функции Лагранжа ст.Х близ­

ка к квадратичной,

то может быть использован алгоритм поиска

максимума по JC ,

приводящий в случае квадратичной зависимо­

сти в точку максимума за один шаг

~

? f e x (te)]4QM>X)

C7.II)

Здесь Q.xx -

квадратичная форма (аналог второй

производ­

ной для функции одной переменной). Множители Лагранжа меня7

ЕТСЯ

согласно зависимости

(7.10).

 

В задачах, где нахождение безусловного максимума

1^ по

У

на множестве V ,

определяемом ограничениями,

реали­

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

рому

на каждом

шаге

Х^ (

находится из условия абсолют­

ного

максимума

Q

, а

по

jj

проводятся итерации согласно

(7.10) с ограничением

на

 

 

J; ( jf.- ^ о).