
книги из ГПНТБ / Цирлин, А. М. Основы оптимального управления конспект лекций
.pdfso
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.- ^ о).