Научная статья на тему 'Об одной модификации условия регулярности постоянного ранга в задачах математического программирования'

Об одной модификации условия регулярности постоянного ранга в задачах математического программирования Текст научной статьи по специальности «Математика»

CC BY
47
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
нелинейное программирование / условие регулярности / оптимизация

Аннотация научной статьи по математике, автор научной работы — О Ф. Борисенко, Л И. Минченко, С М. Стаховский

Условия регулярности (constraint qualifications) играют важную роль в задачах математического программирования, поскольку позволяют гарантировать выполнение принципа Лагранжа в невырожденной форме. В то же время условия регулярности различаются между собой общностью, сравнительной простотой проверки и условиями применения. Ввиду этого значительный интерес вызывает поиск новых типов условий регулярности, применимых к более широким классам задач и более простых для проверки их выполнения. Целью данной заметки является обобщение широко известного в литературе условия регулярности постоянного ранга, а также сравнительный анализ некоторых типов условий регулярности в задачах математического программирования.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ON SOME MODIFICATION OF CONSTANT RANK CONSTRAINT QUALIFICATION IN MATHEMATICAL PROGRAMMING

Constraint qualification of constant rank is studied and generalized. It is proved that generalized constant rank constraint qualification implies metric regularity condition

Текст научной работы на тему «Об одной модификации условия регулярности постоянного ранга в задачах математического программирования»

2008

Доклады БГУИР

№ 6 (36)

УДК 517.977

ОБ ОДНОЙ МОДИФИКАЦИИ УСЛОВИЯ РЕГУЛЯРНОСТИ ПОСТОЯННОГО РАНГА В ЗАДАЧАХ МАТЕМАТИЧЕСКОГО ПРОГРАММИРОВАНИЯ

ОФ. БОРИСЕНКО, ЛИ. МИНЧЕНКО, С.М. СТАХОВСКИЙ

Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь

Поступила в редакцию 17 января 2008

Условия регулярности (constraint qualifications) играют важную роль в задачах математического программирования, поскольку позволяют гарантировать выполнение принципа Ла-гранжа в невырожденной форме. В то же время условия регулярности различаются между собой общностью, сравнительной простотой проверки и условиями применения. Ввиду этого значительный интерес вызывает поиск новых типов условий регулярности, применимых к более широким классам задач и более простых для проверки их выполнения. Целью данной заметки является обобщение широко известного в литературе условия регулярности постоянного ранга, а также сравнительный анализ некоторых типов условий регулярности в задачах математического программирования.

Ключевые слова: нелинейное программирование, условие регулярности, оптимизация.

Введение

Пусть hi (y), i=1,...,p — непрерывно дифференцируемые функции из Rm в R. Введем непустое множество допустимых точек

С= y<=Rm\ к,(у)<0, /е/, ht(y) - 0, /е/0 ,

где yeRm, I={l,...,s}, I0={s+\,...,p} или /о=0, и рассмотрим задачу (/') математического программирования f(y) —> min, у е С с непрерывно дифференцируемой целевой функцией/ Для задачи (Р) введем функцию Лагранжа

L{y,X} = f{y) + {\h{y)),mQX = (\,...,Xp),h = (hl,...,hp)

и множество множителей Лагранжа в точке y

А(у) = /I е | VxL(y, X) - О, Хг>0 и ХД(у) = 0, ге/ .

Условиями, существенно влияющими на возможность эффективного решения задачи (Р), являются условия регулярности. Одним из наиболее известных условий регулярности является условие (RLT) линейной независимости градиентов V//( (у0) / е /(}'„) UА, всех активных в точке j0 ограничений. Более общий характер носит обобщающее его условие (RMF) регу-

лярности Мангасаряна-Фромовица [1-8], требующее, чтобы система векторов У//( (у) / е/0 была линейно независимой и существовал вектор у0 такой, что

<УШ, Ро> = 0, /е/0, (Щ(У\ у0)< О, /е/(>').

Здесь /(^) = {/е/=

Пусть р(х,С) = ¡пГ |л" — , где — евклидова норма вектора, В — открытый единичный шар с центром в 0 в пространстве К".

Будем говорить, что множество С Л-регулярно в точке у0 [9-11], если найдутся числа а>0, Ь>0 такие, что

р(у,С)<атах 0, /г(у) / е /, \кг{у)\ /' е /0 , для всех у е у0+8В.

Известно [9], что выполнение условия регулярности Мангасаряна-Фромовица в точке у0 гарантирует Я-регулярность множества С в этой точке.

Непосредственным обобщением условия (ЯП) является также условие постоянного ранга [4-6]. Напомним [4] , что множество С удовлетворяет в точке у0 условию постоянного

ранга или СЯ-регулярно в этой точке, если для любого подмножества индексов 3 с: 1(у0)^10 система векторов \7/гг (_)'), / е./ имеет постоянный ранг в некоторой окрестности точки у0. Известно, что условия регулярности (ЯМЕ) и (СЯ) независимы друг от друга [4].

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

Очевидным недостатком условия постоянного ранга является трудность его проверки. Ниже предлагается более простая в данном отношении модификация этого условия.

Определение 1. Будем говорить, что множество С удовлетворяет в точке у0 е С модифицированному условию постоянного ранга, или МСЯ-регулярно в этой точке, если для любого подмножества индексов 3 = где К а I(у,, ). система векторов V// (_)'), / е./ имеет

постоянный ранг в некоторой окрестности точки у .

Покажем, что модифицированное условие постоянного ранга (МСК) обеспечивает выполнение принципа Лагранжа в невырожденной форме, т.е. является условием регулярности.

Пусть у0 е С . Рассмотрим касательный конус (нижний касательный конус) к множеству С в точке у0 е С :

Тс(Уо)= У е К™ I 3 функция о(/) такая, что _у + (у + о(/) е С \/t> О и введем множество

Гс(у0)= усКт\ (Щ(у0),у)<0, /е/ОО; (Щ^0\у) = О, ;е/0,

которое будем называть линеаризованным касательным конусом к С в точке у0 е С .

Лемма 1. Пусть множество С М( '/¿-регулярно в точке у0 еС . Тогда Г,. ( у) = Тс (_у0) . Доказательство. Пусть у е Гс (у0) и пусть ./ = -КУ ) -12 (у0, у) и , где 12(у0,у)= / е /(уи) (УИ^у^у) = О . Тогда для любой га-векторной функции г(1) такой, что г(7)// —> 0при / X 0. найдется число > 0 такое, что 1г(у0 +1у + !"(()) < 0 для всех

i el\l2(y0,y)и всех te(o,t0). Действительно, если iel\l(y0), то hj(_>',,) < 0 и, следовательно, /г Оо + (F + r(t)) = /г (>'0) + /(V/г (_у0 + 0(ty + /-(/ )),у) < 0 (где 0 < 0<1 ) при достаточно малых t >0.

Если 7 е 1(у0 ), но г «ё /2 О0, JQ, то

^ 0„ + ^+КО) - ^ Оо)+КЩ(у0\ у)+у(0 - КЩ (у0), у)+у(0,

где у(0 - (Щ(у0), г(0> + (Уh, (у0 + 9(tj + r(t)) - Щ(у0), ty + r(t)), 0 < 0<1.

В этом случае, поскольку (V^.(_y0), >>) < 0 и j(t)/t—>0,то кг.(_у0 + fy + r{t)) < 0 для всех достаточно малых положительных t.

Пусть ранг системы функций h (уп + fy + r) i е J относительно г в точке (/, г) (0. 0)

равен l. В силу условия (МСЕ) он сохранится и в некоторой окрестности (0, 0). Тогда (теорема 2 [12], с. 203) l функций системы (для определенности перенумеруем их так, чтобы это были h,...,h ) независимы, а остальные (если они есть) от них зависят, т.е. hl+l =<pl(hl,...hl),...,hl+q =q>q(hl,...hl), где ф1з...,ф — непрерывные функции с непрерывными

частными производными.

Рассмотрим в окрестности точки (0,0) систему уравнений

My0+fy + r) = 0

h,(y0 +fy + r) = 0 (1)

k,+q(yo+fy + r) = °

относительно переменных t, r.

Очевидно, она равносильна системе

А, С0+$7+ /■) = ()

(2)

с дополнительным условием

ViOo +^7+r) = ^(/2iOo +fy+rl~A(y0 +fy+r)) = о, Kq (Уо+$ + г) = Фч (К (y0+fy + r), . . .h, (y0+fy + r)) = 0.

Приэтом фД/гДУоХ-ЛОо)) = i = \~;<l и, следовательно, фг(0,...0) = 0, i = l,...,q. Если 1=т, то по теореме о неявной функции ([12], с. 188) система (2) определяет в окрестности (0,0) неявную непрерывно дифференцируемую функцию г = r(t) такую, что

г(0) = 0 И — (0) = lim^ = 0. dt ^ t

Если /< m, то, не ограничивая общности, можно предположить, что ранг системы (2) равен / относительно первых / координат вектора г. Положим в этом случае г — (г,г), где г ={г[,...,г1), ? = (rM,...,rm).

Тогда в силу теоремы о неявной функции система (2) определяет в окрестности точки (О, 0, 0) неявную непрерывно дифференцируемую функцию г — r(t, г ) , удовлетворяющую ВТ

условиям г (0,0) = 0, — (0,0) - 0 .

dt

Пусть г -0, положим г = r(t) = г(7,0) . Тогда функция г = r(t) = (г(t\0) удовлетворяет системе (2), а значит и (1). Кроме того, r(t)/t —> 0 при t X 0.

Таким образом, для у е Г,, (у0 ) существует функция г(1) такая, что

Ъ(у0+*Р + г(!)) = 0 i <е J , {y0+fy + r{t)) < 0 ieI\J,

для всех t g (о, t0 ), где t0 достаточно малое положительное число, и г(1)1] —» 0 при t X 0. Последнее означает, что y0+fy + ¡"(t ) е С при t е [0,t0] и, следовательно, у е Т(,(у(] ).

Теорема 1. Пусть множество С MCR-регулярно в точке у0 е С . Тогда для любой непрерывно дифференцируемой функции f, имеющей в точке у0 е С локальный минимум на множестве С, справедливо условие A( j'0 ) Ф 0 .

Доказательство. В силу леммы 1 для любого у е Г,. ( у) найдется векторная функция o(t) такая, что у0 + fy + o(t) еС, t > 0, и, следовательно,

f(y0+fy+o(t))-f(y0)>0,

откуда (У/(уа\у)> 0 для всех <Е Г,. (j'0 ). В таком случае по теореме двойственности линейного программирования [13] множество допустимых точек двойственной задачи, совпадающее с Л(_у0), не пусто.

Связь MCR-регулярности с ^-регулярностью

Установим связь MCR-регулярности с Я-регулярностью.

Пусть V е R'\ viC. Обозначим Пс(у) множество точек из С, ближайших к точке v . Очевидно, эти точки являются решениями задачи нелинейного программирования

/vW^min, уеС,

где fv(y) = \y~v\-

Обозначим

Lv{y,X) = fv{y) + {\h{y)\ hv{y)= VLv(y,X) = 0, >0 и Wj) = 0,/g/ .

Лемма 2. Пусть условие MCR-регулярности выполнено в точке у0 е frC . Тогда существуют числа M >0 и 5 > 0 такие, что при всех V £ С таких, что v е уп + дБ условие

р

V'/ (у) = {Я е Av (у) : < М} Ф 0 выполнено во всех точках y(v) е Пс (v) .

i=\

Доказательство. Допустим противное. В таком случае найдется последовательность ^ -> >'„ такая, что vkiC и р(0, AVt (yk)) оо при yk= y(vk) . При этом ук = y(vk) уи,

поскольку \vk - ук | <\vk - у01 .

Из условия МСЯ-регулярности в точке у0 е /гС следует, что это условие будет выполнено и в некоторой окрестности данной точки. В таком случае , не ограничивая общности, можно считать, что Л (ук ) Ф 0 . При этом ук удовлетворяет ограничениям

кг{ук) = 0 /е/0иад Ь(Ук)<0 ^(/0и/(Л))

при достаточно больших к . Поскольку I(ук) может принимать только конечное число значений, то из последовательности {ук} можно извлечь подпоследовательность (для простоты будем считать, что это сама {ук} ), для которой множество индексов /(у к) = I не зависит от к . Пусть J = I0 и/*. Тогда Иг(ук) = О, /' е 10 и/*; Иг(ук) <0, £ (70 и/*).

В силу условия МСЯ-регулярности система {V// (у) / е./{ не меняет ранг в окрестности точки у0. Пусть этот ранг равен /. Тогда система векторов {V// (у0), / е./{ имеет ранг I. Обозначим через В(у0) базисный минор этой системы. В силу непрерывности данного минора по у можно, не ограничивая общности, считать, что В(ук )#0 и, следовательно, в силу условия МСЯ-регулярности В( ук) остается базисным минором для системы {У^(ук\ / е./1 . Обозначим через У0 множество индексов из ./ . для которых соответствующие векторы являются базисными для базисного минора В (у0). Поскольку

р

Ач(ук)= ХеКр\^\УИ1(ук) + У/ч(ук) = 0, \>0 /е/Д=0 /£./},

р

г=\

то для всех к=1, 2, ..., система

/-1

имеет единственное решение Хк . При этом базисный минор не обращается в 0 при ук —> у0 и, следовательно, последовательность {Xк } ограничена, что противоречит предположению о том, что р(0, Л (ук)) —> оо . Полученное противоречие доказывает справедливость утверждения леммы.

Лемма 3 ([8]). Пусть у0 е /гС и существуют числа М > 0 и 80 > 0 такие, что при всех

V е у0 + 2-1805, V <£ С в точках у - у(у) е Пс (у) выполнено условие

р

г= 1

Тогда множество С Л-регулярно в точке у0.

Теорема 2. Пусть условие МСЯ-регулярности выполнено в точке у0 е /гС . Тогда

множество С Л-регулярно в данной точке.

Доказательство. Справедливость утверждения теоремы непосредственно получается последовательным применением лемм 2 и 3.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Рассмотрим ослабленный вариант условия из Определения 1. Будем говорить, что множество С удовлетворяет в точке у0 е С слабому условию (МСК), если для любого подмножества индексов 3 = К ^ /0, где К с 1(у0), система векторов V// (у), / е./ имеет постоянный ранг для всех у е С из некоторой окрестности точки у0.

Легко видеть, что данное условие тривиально выполняется в случае одноточечного множества С — {у0} и, как нетрудно убедиться, не является условием регулярности, так как не

гарантирует непустоту множества множителей Лагранжа в точках экстремума на множестве С.

Поскольку все условия регулярности обязаны гарантировать не пустоту множества множителей Лагранжа, то самый широкий характер, очевидно, носят так называемые условия регулярности по Лагранжу [3].

Определение 2. Множество С регулярно в точке J0 по Лагранжу, если при любой гладкой функции f, имеющей в точке локальный минимум на множестве С, множество А(_у0) не пусто.

Лемма 4. Пусть множество С регулярно по Лагранжу в точке у0 е frC и в ее окрестности на С, ив этой точке выполнено слабое условие (MCR). Тогда существуют числа М > О

и 8 > 0 такие, что при всех v £ С таких, что v е уп + 5В, условие

р

A^(y) = {XeAv(y): ^\\\<М}ф0

г=\

выполнено во всех точках у(у) е Пс (v) .

Доказательство данной леммы получается незначительной модификацией доказательства леммы 2.

Последовательное применение лемм 3 и 4 приводит к к следующему обобщению теоремы 2.

Теорема 3. Пусть множество С регулярно по Лагранжу в точке у0 е frC и в ее окрестности на C, и в этой точке выполнено слабое условие (MCR). Тогда множество С R-регулярно в данной точке.

ON SOME MODIFICATION OF CONSTANT RANK CONSTRAINT QUALIFICATION IN MATHEMATICAL PROGRAMMING

OF. BORISENKO, L.I. MINCHENKO, S.M. STAKHOVSKI

Abstract

Constraint qualification of constant rank is studied and generalized. It is proved that generalized constant rank constraint qualification implies metric regularity condition.

Литература

1. Bonnans J.F., Shapiro A. Perturbations analysis of optimization problems. N.-Y., 2000.

2. Luderer B., Minchenko L., Satsura T. Multivalued analysis and nonlinear programming problems with perturbations. Dordreht, 2002.

3. Березнев В.А., Завриев С.Н., Третьяков А.А. // Докл. АН СССР. 1988. Т. 300, № 6. С. 1289-1291.

4. Janin R. // Mathematical Programming Study. 1984. Vol. 21. P. 110-126.

5. Pang J.-S., Ralph D. //Mathematics of Operations Research. 1993. P. 102-154.

6. Ralph D., Dempe S. // Mathematical Programming. 1995. Vol. 70. P. 159-172.

7. Минченко Л.И., Тараканов А.Н. // Докл. НАН Беларуси. 2004. Т. 48, № 4. С. 24-28.

8. Минченко Л.И., Гвоздь Е.И. // Докл. НАН Беларуси. 2008. Т. 52, № 1. С. 7-12.

9. Robinson S.M. // SIAM J. Numer. Analysis. 1976. Vol. 13. P. 497-513.

10. Ioffe A.D. // Trans. Amer. Math. Soc. 1979. Vol. 251. P. 61-69.

11. Федоров В.В. Численные методы максимина. М., 1979.

12. Зорич В.А. Математический анализ. Ч. 1. М., 1981.

13. КармановВ.Г. Математическое программирование. М., 1978.

i Надоели баннеры? Вы всегда можете отключить рекламу.