Научная статья на тему 'Внешняя оценка компакта лебеговым множеством выпуклой функции'

Внешняя оценка компакта лебеговым множеством выпуклой функции Текст научной статьи по специальности «Математика»

CC BY
48
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КАЛИБРОВОЧНАЯ ФУНКЦИЯ / ВНЕШНЯЯ ОЦЕНКА / СУБДИФФЕРЕНЦИАЛ / КВАЗИВЫПУКЛАЯ ФУНКЦИЯ / СИЛЬНО ВЫПУКЛОЕ МНОЖЕСТВО / СИЛЬНО ВЫПУКЛАЯ ФУНКЦИЯ / GAUGE FUNCTION / EXTERNAL ESTIMATE / SUBDIFFERENTIAL / QUASICONVEX FUNCTION / STRONGLY CONVEX SET / STRONGLY CONVEX FUNCTION

Аннотация научной статьи по математике, автор научной работы — Абрамова Вероника Валерьевна, Дудов Сергей Иванович, Осипцев Михаил Анатольевич

Рассматривается конечномерная задача о вложении заданного компакта D ⊂ Rp в нижнее лебегово множество G(α) = {y ∈ Rp : f(y) 6 α} выпуклой функции f(·) с наименьшим значением α за счет смещения D. Ее математическая формализация приводит к задаче минимизации функции φ(x) = max y∈D f(y x) на Rp. Исследованы свойства функции φ(x), получены необходимые и достаточные условия и условия единственности решения задачи. В качестве базового для приложений выделен случай, когда f(·) калибровочная функция Минковского некоторого выпуклого тела M. Показано, что если M многогранник, то задача сводится к задаче линейного программирования. Предложен подход к получению приближенного решения, в котором при построении последовательности приближений {xi}i=0,1,..., зная приближение xi, для получения xi+1 требуется решить более простую задачу вложения компакта D в лебегово множество калибровочной функции множества Mi = G(αi), где αi = φ(xi). Дается обоснование сходимости последовательности приближений к решению задачи.The finite-dimensional problem of embedding a given compact D ⊂ Rp into the lower Lebesgue set G(α) = {y ∈ Rp : f(y) 6 α} of the convex function f(·) with the smallest value of α due to the offset of D is considered. Its mathematical formalization leads to the problem of minimizing the function φ(x) = max y∈D f(y x) on Rp. The properties of the function φ(x) are researched, necessary and sufficient conditions and conditions for the uniqueness of the problem solution are obtained. As an important case for applications, the case when f(·) is the Minkowski gauge function of some convex body M is singled out. It is shown that if M is a polyhedron, then the problem reduces to a linear programming problem. The approach to get an approximate solution is proposed in which, having known the approximation of xi to obtain xi+1 it is necessary to solve the simpler problem of embedding the compact set D into the Lebesgue set of the gauge function of the set Mi = G(ai), where ai = f(xi). The rationale for the convergence for a sequence of approximations to the problem solution is given.

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

Похожие темы научных работ по математике , автор научной работы — Абрамова Вероника Валерьевна, Дудов Сергей Иванович, Осипцев Михаил Анатольевич

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

Текст научной работы на тему «Внешняя оценка компакта лебеговым множеством выпуклой функции»

МАТЕМАТИКА

УДК 519.853

НАУЧНЫЙ ОТДЕЛ

Внешняя оценка компакта лебеговым множеством выпуклой функции

B. B. Абрамова, С. И. Дудов, М. А. Осипцев

Абрамова Вероника Валерьевна, аспирант кафедры математической экономики, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, г. Саратов, ул. Астраханская, д. 83, Veronika0322@rambler.ru

Дудов Сергей Иванович, доктор физико-математических наук, заведующий кафедрой математической экономики, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, г. Саратов, ул. Астраханская, д. 83, DudovSI@info.sgu.ru

Осипцев Михаил Анатольевич, кандидат физико-математических наук, доцент кафедры математического анализа, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, г. Саратов, ул. Астраханская, д. 83, Osipcevm@gmail.com

Рассматривается конечномерная задача о вложении заданного компакта В с в нижнее лебегово множество О(а) = {у е : /(у) < а} выпуклой функции /(•) с наименьшим значением а за счет смещения В. Ее математическая формализация приводит к задаче минимизации функции ф(х) = тах /(у - х) на . Исследованы свойства функ-

уеБ

ции ф(х), получены необходимые и достаточные условия и условия единственности решения задачи. В качестве базового для приложений выделен случай, когда /( ) — калибровочная функция Минковского некоторого выпуклого тела М. Показано, что если М — многогранник, то задача сводится к задаче линейного программирования. Предложен подход к получению приближенного решения, в котором при построении последовательности приближений {хг}г=0д,..., зная приближение хг, для получения хг+1 требуется решить более простую задачу вложения компакта В в лебегово множество калибровочной функции множества Мг = О(аг), где аг = ф(хг). Дается обоснование сходимости последовательности приближений к решению задачи.

Ключевые слова: калибровочная функция, внешняя оценка, субдифференциал, квазивыпуклая функция, сильно выпуклое множество, сильно выпуклая функция.

Поступила в редакцию: 12.03.2019 / Принята: 05.06.2019 / Опубликована: 01.06.2020

Статья опубликована на условиях лицензии Creative Commons Attribution License (CC-BY4.0)

DOI: https://doi.org/10.18500/1816-9791-2020-20-2-142-153

1. ПОСТАНОВКА ЗАДАЧИ

Оценка сложных множеств множествами простой структуры — одно из актуальных направлений негладкого анализа [1-4]. В рамках этого направления находится и рассматриваемая здесь задача.

Пусть D — ограниченное замкнутое множество из конечномерного пространства Rp, а f (•) — выпуклая конечная на Rp функция. Требуется вложить множество D в нижнее лебегово множество G(a) = {y е Rp : f (y) ^ a} этой функции с наименьшим значением a за счет смещения D, т. е. решить задачу

a ^ min ,

(«,X)gRp+1 (1)

D - х С G(a).

Задача такого вида возникает, например, при расчете параметров проектируемого технического устройства [5,6]. При этом компакт D может отражать точностные возможности оборудования, используемого при изготовлении устройства, а функция f (•) — качественные характеристики проектируемого устройства.

Для корректности и нетривиальности задачи будем предполагать, что функция f(•) ограничена снизу на Rp, ее минимальное значение достигается, множество G(a0), где a0 = min f (x), является ограниченным и при этом

xERp

*

G(ao) - D = {z е Rp : z + D С G(ao)} = 0. Нетрудно видеть, что задача (1) эквивалентна следующей минимаксной задаче

ф(х) = max f (y — x) ^ min . (2)

rv y yeDJKU J XGRP v 7

При этом, если ф(х*) = minф(х), то пара (ф*,х*), где ф* = ф(х*), является решением

XERp

задачи (1).

Нетрудно заключить, что при сделанных предположениях решение задачи существует. Отметим также, что на задачу (2) можно смотреть как на обобщение задачи о чебышевском центре множества D (случай, когда f (•) — некоторая норма на Rp).

Цель работы — указать на основные свойства целевой функции задачи (2), получить критерий решения и условия единственности решения, а также предложить подход к получению приближенного решения задачи.

Далее используются следующие обозначения: A, intA, coA — соответственно замыкание, внутренность, выпуклая оболочка множества A;

||х|| — евклидова норма элемента х е Rp;

df (x) — субдифференциал выпуклой функции f (x) в точке x;

B(x,r) = {y G Rp : ||x — y|| ^ r} — шар с центром в точке x и радиусом r;

(x, y) — скалярное произведение элементов x и y из Rp;

s(-, A) = sup(^, a) — опорная функция множества A, 0p = (0,0,..., 0) G Rp.

aeA

2. СВОЙСТВА ЦЕЛЕВОЙ ФУНКЦИИ. КРИТЕРИЙ РЕШЕНИЯ

Основные свойства функции 0(x) отражает следующая теорема.

Теорема 1. Пусть f (x) является выпуклой конечной на Rp функцией. Тогда функция 0(x) является:

1) выпуклой конечной на Rp функцией, а ее субдифференциал в любой точке x G Rp может быть выражен в виде

d0(x) = —co{df (z — x) : z G Q*(x)}, (3)

где (x) = {z G D : f (z — x) = 0(x)};

2) строго выпуклой, если функция f (•) строго выпукла на Rp;

3) сильно выпуклой на Rp с константой C > 0, т. е. для любых точек x0, xi из Rp и a G [0,1] выполняется

C

ф((1 — a)xo + axi) ^ (1 — a)0(xo) + a0(xi) — — a(1 — a)||xo — xi||2,

если функция f (•) сильно выпукла на Rp с константой C;

4) строго квазивыпуклой, т. е. для любых точек x0 = x1 из Rp и a G (0,1) выполняется

ф((1 — a)x0 + ax1 ) < max{0(xo), 0(x1 )}, если функция f (•) строго квазивыпукла на Rp.

Доказательство. 1) Функция 0(x), как функция максимума от семейства выпуклых конечных функций f (y — x) по параметру y на компакте D, также является выпуклой конечной на Rp. Кроме того, функция ^(x,y ) = f (y — x) является непрерывной всюду по x и по y. Тогда, по известной в выпуклом анализе теореме о субдифференциале функции супремума от семейства выпуклых функций по параметру (см. [1, гл. 2, § 1], [7, гл. 4, § 4.2] или [8, гл. 1, § 1.17]), получаем

d0(x) = —cô{df (z — x) : z G Q*(x)}.

Знак замыкания здесь можно снять. Действительно, непрерывность выпуклой конечной функции f(•) и замкнутость множества D дают замкнутость множества Q^(x). Кроме того, субдифференциал df(-), как многозначное отображение, обладает полунепрерывностью сверху. Отсюда следует замкнутость множества в фигурной скобке.

Свойства 2)-4) легко вытекают из дополнительных условий, накладываемых на функцию f(•). □

Теорема 2. 1. Для того чтобы точка x* была точкой минимума функции 0(x) на Rp, необходимо и достаточно, чтобы

0p G co{df (z — x*): z G Q^(x*)}. (4)

2. Если существует ö > 0 такое, что

B(0p, ö) С co{df (z - x*) : z G Яф(x*)}, (5)

то x* является единственной точкой минимума функции ф(х) на Rp, причем

ф(х) ^ ф(х*)+ ö||x - x*У, V x G Rp. (6)

Доказательство. 1. Первое утверждение теоремы вытекает из формулы субдифференциала (3) и известного в выпуклом анализе критерия достижения выпуклой функции своего минимального значения (например, [1, гл. 4, § 2])

0(x*) = min 0(x) ^ 0p G d^(x*).

2. Из (3)-(5), в соответствии с определением субдифференциала [1-3]

Ö0(x) = {v G Rp : ф(у) - 0(x) ^ (v, y - x>, V y G Rp},

следует

0(x) - 0(x*) ^ (v, x - x*>, V v G B(0p, ö), x G Rp.

Отсюда, положив v = -—ö—-(x - x*), получаем (6). □

||x - x* |

Замечание. Очевидно, единственность решения задачи (2), кроме соотношения (5), обеспечивает также строгая квазивыпуклость выпуклой функции f(•) и тем более ее строгая или сильная выпуклость. Это вытекает из теоремы 1.

3. БАЗОВЫЙ СЛУЧАЙ

Базовым для предлагаемого далее подхода к получению приближенного решения задачи (2) является случай, когда функция f (•) является калибром (калибровочной функцией Минковского) некоторого выпуклого телесного компакта M, причем 0p G int M, т. е.

f (x) = k(x, M) = inf{а ^ 0 : x G aM}. (7)

Напомним основные свойства калибра (например, [1,9]):

а) если Л > 0, то k(Ax, M) = Ak(x, M), V x G Rp;

б) если x G M, то k(x, M) ^ 1, а если x / M, то k(x, M) > 1;

в) k(x + y, M) ^ k(x, M) + k(y, M), V x, y G Rp.

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

3.1. Пусть выпуклое тело M является многогранником, заданным в виде

M = {y G Rp : (Л, y) ^ 1, i = 1, m},

где Ai G Rp, i = 1, m, 0p G intcojA, : i = 1,m}.

Из свойств калибра следует: если max (A« , x) = 1, то точка x является граничной

i=1,m

для M и k(x,M) = 1. А поскольку калибр является положительно однородной функцией, мы получаем для него формулу

k(x, M) = max(Ai,x), Vx G Rp.

i=1,m

Тогда задача (2) принимает вид

max max (A, y — x) ^ min

i=1,m

max{a^ — (A^, x)} ^ min, (8)

i=1,m

или

где a = max(Ai5 y).

yeD

Задача (8) известным приемом [10] сводится к задаче линейного программирования. А именно справедлива

Теорема 3. Задача (8) эквивалентна следующей задаче линейного программирования:

x(p+1) ^ min ,

(ж,ж(Р+1) )eRP+1 (9)

(A, x) + x(p+1) - a ^ о, г =

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

Причем, если X0 = (x0 , x0P+1)) — решение задачи (9), то x0 — решение задачи (8). А если x0 — решение задачи (8), то x0 = (x0,x0p+1)), где x0p+1) = max{a — (Ai5x0)}, решение задачи (9).

i=1,m

3.2. Теперь предположим, что M является сильно выпуклым множеством. Посмотрим, как это отразится на свойствах функции f (•) вида (7) и, соответственно, на функции ф(-).

Сначала напомним некоторые понятия и факты из параметрически выпуклого анализа [8,11].

Определение 1. Множество A с Rp называется r-сильно выпуклым, если оно представимо в виде пересечения евклидовых шаров радиуса r.

Наряду с давно введенными понятиями сильно и слабо выпуклых (вогнутых) функций (см. [8,11-13]) имеют место эквивалентные им альтернативные определения, которыми мы и будем пользоваться.

Определение 2. Пусть задано выпуклое множество X с Rp и функция f : X ^ R. Функция f (•) называется сильно (слабо) выпуклой с константой C > 0 на множестве X, если функция f (x) — -||x||2 (f (x) + —||x||2) выпукла на X.

Функция f(•) называется сильно (слабо) вогнутой с константой C > 0 на множестве X, если функция —f (•) сильно (слабо) выпукла на X с константой C > 0.

Лемма 1. ([11, гл. 2, § 2.8]). Если функция f : Rp ^ R слабо вогнута с константой C > 0, то ее сопряженная функция

f *(y ) = sup {(x,y) — f (x)}

xeRp

сильно выпукла с константой -.

Лемма 2. ([11, гл. 3, § 3.7]). Пусть M с Rp является r-сильно выпуклым множеством и задана дифференцируемая, слабо вогнутая с константой C > 0 функция F : R ^ R такая, что |F'(t)| ^ L|t| для любого t g R. Тогда функция f(x) = F(s(x, M)) слабо вогнута на Rp с константой CMC + LCMr, где CM = sup ||x||.

xeM

Теперь докажем важное свойство калибра.

Лемма 3. Пусть M с Rp является r-сильно выпуклым множеством и См = sup ||x||. Тогда функция Ф(-) = k2(^,M) является сильно выпуклой на Rp

xeM

2

с константой, равной

См (См + г)-

Доказательство. Для сопряженной к Ф(-) функции получаем

Ф*(y) = sup {(y,x) — k2(x, M)} = sup sup {(y, Ax) — A2k2(x,M)} =

xeRp A^0 x:k(x,M)=1

= sup{A sup (y, x) — A2} = sup{As(y, M) — A2} = ^(y' M)

A^0 x:k (x,M)=1 A^0

Таким образом, сопряженная функция имеет вид

4

t2

Ф* (y) = F (s(y, M)), F (t) = -.

Функция F(•) является дифференцируемой, |F'(t)| ^ L|t| для L = 1/2 и, в соответствии с определением 2, слабо вогнута с константой 1/2. Поэтому по лемме 2 функция Ф*(y) слабо вогнута с константой 2См(См + r). Тогда по лемме 1 функция Ф**(•) = (k2)**(^,M) будет сильно выпуклой с константой 2(См(См + r))-1. Остается заметить, что функция k2(-,M) является выпуклой, как квадрат выпуклой функции с неотрицательными значениями, и непрерывной на Rp. Поэтому по теореме Фенхеля-Моро [7-9,11] выполняется (k2 )** (-,M ) = k2 (-,M). □

Очевидным следствием леммы 3 является

Теорема 4. Если множество M является r-сильно выпуклым и См = sup ||x||,

xeM

то функция ф2 (x) = max k2 (y — x,M) является сильно выпуклой на Rp с

yeD

2

константой п , ^.

CM (CM +r)

Замечание. Сама функция ф(^), как и f (•) вида (7), не является сильно выпуклой независимо от свойств множества M. Минимизация функции ф(-) при f (•), имеющей вид (7), эквивалентна минимизации функции 02(x). А применение к минимизации сильно выпуклой функции численных методов (например, субградиентного типа, см. [2,14]) позволяет рассчитывать на сходимость последовательности приближений к решению со скоростью геометрической прогрессии.

4. ПОДХОД К ПРИБЛИЖЕННОМУ РЕШЕНИЮ

Поскольку ф(-) является выпуклой функцией, известна формула ее субдифференциала, то для получения приближенного решения задачи минимизации этой функции можно применять хорошо известные численные методы выпуклого программирования (см., например, [2, 14]). Мы обозначим здесь еще один возможный подход к получению приближенного решения задачи (2) с произвольной выпуклой функцией f (x) и компактом D, удовлетворяющих заявленным в п. 1 условиям.

Принципиальная схема метода

*

Пусть f(y0) = minf(y). По исходному предположению G(a0) — D = 0, где

yeRp

a0 = f(y0). Поэтому в силу непрерывности выпуклой конечной функции f(•) выполняется y0 G intG(a) для любого a > a0 (или 0p G int(G(a) — y0)). Последовательность приближений {x^}i=0)1)... строится следующим образом. Точка начального приближения x0 выбирается произвольно. Предположим, уже получена точка x«. Возьмем в качестве M« = G(a^) — y0, где a« = ). После этого решаем задачу

Ф(x) = max k(y — x, Mj) ^ min (10)

yeD-yo ^ ^ xeRp V 7

и в качестве xi+1 берем одно из ее решений, т.е. ф(xi+1) = minф^(x). Если при

xeRp

этом окажется, что ф^) = 0(xi+1), то, как будет доказано, это будет означать,

что ) = min 0(x) и итерационный процесс на этом заканчивается. В противном

xeRp

случае, т.е. при ф^) = 0(xi+1), процесс построения последовательности продолжается.

Обоснование сходимости

Без потери общности будем далее считать y0 = 0p. Из предполагаемой ограниченности множества G(a0) и выпуклости функции f(•) следует ограниченность множества G(a) при любом a > a0 [2, гл. 1, § 9]. Напомним также, что выпуклая конечная функция удовлетворяет условию Липшица на любом ограниченном множестве [2,8,9].

Докажем следующий вспомогательный факт.

Лемма 4. Пусть a > ф* = min 0(x), L — константа Липшица для функции f (•) на множестве G(a) и C ^ sup ||x||. Если точка x*(a) такова, что

xeG(a)

max k(y — x*(a), G(a)) = min max k(y — x, G(a)), (11)

yeD xeRp yeD

то

ф(**(«)) < a — ^O—C+S?• d2)

Доказательство. Из липшицевости функции f (•) следует

£(ф*)+ B(0Р, А) с G(a), (13)

a — ф* I

где А = —-—. С другой стороны, поскольку G(a*) с G(a) и — G(a) с B(0p, I), то LC

) с С(ф*) + B(0p, Д). (14) Если x* G Arg min ф(х), то D — x* с G(0*). Поэтому, учитывая (13)—( 14), имеем

C

D — x* С G(a). (15)

148 Научный отдел

Из (11), (15) и определения калибра вытекает

C

max k(y — x*(a), G(a)) ^ max k(y — x*, G(a)) ^ .

yGD yGD C + A

Рассмотрим функцию

Ф(х, a) = a0 + (a — a0)k(x, G(a)). Нетрудно видеть, что в силу выпуклости функции f (•) и свойств калибра

f (x) < Ф(х,а), Vx G G(a). Очевидно D — x*(a) С G(a). Тогда, используя (16)—(18), получаем

0(x*(a)) = maxf (y — x*(a)) ^ max^(y — x*(a), a) =

yED yeD

C

= a0 + (a — a0) max k(y — x*(a), G(a)) ^ a0 + (a — a0) —-—.

yeD C + A

После очевидного преобразования этого неравенства и подстановки A получаем (12).

(16)

(17)

(18)

a — ф:

L

Обозначим через p(x, S*) = inf ||x — y || — расстояние от точки x до S* = {y G : ф(у) = min ф(х) — множества решений задачи (2).

Теорема 5. Если последовательность {х^}^=0,1,... строится в соответствии с указанной схемой, то возможны два варианта:

1) на некотором конечном шаге i0 + 1 будет выполняться равенство 0(xio) = 0(xio+1), что будет означать 0(xio) = min ф(х);

2) последовательность {ф^)}i=0)1... является строго убывающей и бесконечной, причем

ф(хт) ^ ф(хЛ —

(ф(х) — ao )(ф(х^) — ф*)

LC + ф(х*) — ф*

(19)

где L — константа Липшица для функции f(•) на множестве G(a0), а C ^ sup ||x||. Кроме того, имеет место

ireG(ao)

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

lim p(xi*) = 0.

(20)

Доказательство. В соответствии со схемой метода точка x^+i является одним из решений задачи (10). Если a = ф^) > ф*, то по лемме 4, при а = а и x*(а^) = xi+1, получаем 0(xi+1) < ). Следовательно, если на конечном шаге i0 + 1 выполняется равенство 0(xio) = 0(xio+1), то это будет означать 0(xio) = min0(x) = ф*.

В противном случае, т. е. если такое равенство не достигается на конечном шаге, последовательность {ф^)}i=0)1)..., подчиняясь, по лемме 4, неравенству (19), будет строго убывающей и бесконечной.

Предположим, (20) неверно. Тогда существует подпоследовательность {x^}=0д... такая, что для некоторого а > 0 выполняется

р(х?1 *) ^ a, j = 0,1,...

(21)

Последовательность {хг}г=0)1)... является ограниченной, поскольку

D — X С С(ф(х,)) с С(ФЫ),

а любое множество G(a), как уже отмечалось, в силу исходных требований функции к f (•) и множеству D является ограниченным. Поэтому без потери общности можно считать подпоследовательность {х^.}j=0)1... сходящейся: х^. ^ х*, j ^ то. Известно, что функция расстояния является непрерывной. Поэтому, переходя к пределу по по j ^ то в неравенстве (21), получаем р(х*, S * ) ^ a и, следовательно, х* / S * и ф(х*) > ф*.

С другой стороны, и функция ф(х), как выпуклая и конечная на Rp, является непрерывной. Поэтому из монотонности последовательности {ф(хг)}г=0)1)... следует lim ф(хг) = lim ф(хг.) = ф(х*). Тогда, переходя в (19) к пределу по i ^ то, получаем

г^ж j ^ж j

неравенство

«х- ) < «х- ) - '«%-лm -' ■

которое является противоречивым ввиду полученного выше неравенства ф(х *) > ф *. Теорема доказана полностью. □

Авторы еще раз подчеркивают, что лишь обозначили альтернативный применению известных методов выпуклого программирования подход к получению приближенного решения рассматриваемой задачи. Практическая реализация предполагает умение решать на каждом шаге некоторые нетривиальные (но более простые по отношению к исходной) вспомогательные задачи. В связи с этим приведем ряд комментариев, которые будут полезны при реализации данного подхода. При этом сразу отметим, что компакт D будем считать выпуклым, поскольку в случае невыпуклости его замена в задаче (2) на coD дает новую задачу, эквивалентную прежней.

1. Изначально можно аппроксимировать выпуклый компакт D выпуклым многогранником De с некоторой погрешностью 6 > 0. Тогда подсчет приближенного к ф(х) значения фе(х) = max f (y — х) будет сведен к вычислению значений f (y — х)

V^Ds

для y из конечного множества вершин многогранника De и выбору максимального из них.

2. Множество Мг = ^(аг) — y0 на каждом шаге также можно аппроксимировать многогранником. Отметим наличие обширного спектра методов полиэдральной аппроксимации (см., например, обзор [15]). В частности, в нашем случае из точки y0, которая является внутренней для G(a) при a > a0, можно провести лучи до пересечения с границей Мг ив этих граничных точках построить опорные к Мг гиперплоскости, которые и образуют многогранник М/, аппроксимирующий многогранник Мг. Направление этих лучей можно задать сеткой точек «мелкости» ö > 0 на единичной сфере (см. [8, § 4.9]).

3. В качестве приближенного значения k(х,Мг) можно брать k(х, М/). Если многогранник М/ выражен в форме М/ = {х G Rp : (Aj, х) ^ 1,j = 1,m}, то, как отмечалось в п. 3.1, справедлива формула

k (х,М/ ) = max (Aj ,х).

j=1,m

Следовательно, в качестве приближенного значения для фг(х) можно брать

ф7е' (х) = max k (y — х, М/ ) = max {a. — (A., х)},

^ ta? — (A? V^Ds j=1,m

где aj = max(Aj, y).

vGDs

4. Решение задачи (10) можно заменить решением «приближенной» задачи

(x) ^ min,

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

5. Конечно, такой подход потребует исследования устойчивости и чувствительности решения задачи к погрешности полиэдрального приближения выпуклого компакта D и множеств M¿.

Благодарности. Авторы выражают признательность рецензенту за указанные погрешности и конструктивную критику.

Библиографический список

1. Пшеничный Б. Н. Выпуклый анализ и экстремальные задачи. М. : Наука, 1980. 320 с.

2. Демьянов В. Ф., Васильев Л. В. Недифференцируемая оптимизация. М. : Наука, 1981. 384 с.

3. Демьянов В. Ф., Рубинов А. М. Основы негладкого анализа и квазидифференциальное исчисление. М. : Наука, 1990. 431 с.

4. Кларк Ф. Оптимизация и негладкий анализ. М. : Наука, 1988. 280 с.

5. Брейтон Р. К., Хэчтел Г. Д., Санджованни-Винчентелли А. Л. Обзор методов оптимального проектирования интегральных схем // ТИИЭР. 1981. Т. 69, № 10. С. 180-215.

6. Дудов С. И., Мещанов В. П. Параметрическая оптимизация проектируемых устройств по критериям стоимости и качества // Обзоры по электронике. Сер. 1. Электроника СВЧ. М. : ЦНИИ «Электроника», 1990. Вып. 1 (1512). 40 с.

7. Иоффе А. Д., Тихомиров В. М. Теория экстремальных задач. М. : Наука, 1974. 481 с.

8. Половинкин Е. С., Балашов М. В. Элементы выпуклого и сильно выпуклого анализа. М. : Физматлит, 2004. 240 с.

9. Рокафеллар Р. Т. Выпуклый анализ. М. : Мир, 1973. 340 с.

10. Зуховицкий С. И., Авдеева Л. И. Линейное и выпуклое программирование. М. : Наука, 1964. 460 с.

11. Иванов Г. Е. Слабо выпуклые множества и функции. М. : Физматлит, 2006. 352 с.

12. Поляк Б. Т. Введение в оптимизацию. М. : Наука, 1983. 384 с.

13. Vial J.-P. Strong and weak convexity of set and funtions // Math. Ops. Res. 1983. Vol. 8, № 2. P. 231-259.

14. Васильев Ф. П. Методы оптимизации. М. : МЦНМО, 2011. 433 с.

15. Бронштейн E. M. Аппроксимация выпуклых множеств многогранниками // СМФН. 2007. Т. 22. С. 5-37.

Образец для цитирования:

Абрамова В. В., Дудов С. И., Осипцев М. А. Внешняя оценка компакта лебеговым множеством выпуклой функции // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2020. Т. 20, вып. 2. С. 142-153. 001: https://doi.org/10.18500/1816-9791-2020-20-2-142-153

The External Estimate of the Compact Set by Lebesgue Set

of the Convex Function

V. V. Abramova, S. I. Dudov, M. A. Osiptsev

Veronika V. Abramova, https://orcid.org/0000-0002-7926-1347, Saratov State University, 83 As-trakhanskaya St., Saratov 410012, Russia, Veronika0322@rambler.ru

Sergei I. Dudov, https://orcid.org/0000-0002-7926-1347, Saratov State University, 83 Astrakhan-skaya St., Saratov 410012, Russia, DudovSI@info.sgu.ru

Mikhail A. Osiptsev, https://orcid.org/0000-0002-7926-1347, Saratov State University, 83 As-trakhanskaya St., Saratov 410012, Russia, Osipcevm@gmail.com

The finite-dimensional problem of embedding a given compact D c Rp into the lower Lebesgue set G(a) = {y 6 Rp : f (y) < a} of the convex function f (•) with the smallest value of a due to the offset of D is considered. Its mathematical formalization leads to the problem of minimizing the function 0(x) = max f (y - x) on Rp. The properties of the function 0(x) are researched, necessary and sufficient conditions and conditions for the uniqueness of the problem solution are obtained. As an important case for applications, the case when f (•) is the Minkowski gauge function of some convex body M is singled out. It is shown that if M is a polyhedron, then the problem reduces to a linear programming problem. The approach to get an approximate solution is proposed in which, having known the approximation of x* to obtain xi+1 it is necessary to solve the simpler problem of embedding the compact set D into the Lebesgue set of the gauge function of the set Mi = G(a*), where a* = f (xi). The rationale for the convergence for a sequence of approximations to the problem solution is given.

Keywords: gauge function, external estimate, subdifferential, quasiconvex function, strongly convex set, strongly convex function.

Received: 12.03.2019 / Accepted: 05.06.2019 / Published: 01.06.2020

This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0)

Acknowledgements: The authors are grateful to the referee for the indicated errors and constructive criticism.

References

1. Pschenichnyi B. N. Vypuklyi analiz i ekstremal'nye zadachi [Convex Analysis and Extremal Problems]. Moscow, Nauka, 1980. 320 p. (in Russian).

2. Demyanov V. F., Vasiliev L. V. Nondifferentiable Optimization. New York, Springer-Optimization Software, 1986. 452 p. (Russ. ed.: Moscow, Nauka, 1981. 384 p.).

3. Demyanov V. F., Rubinov A. V. Osnovy negladkogo analiza i kvazidifferentsial'noe is-chislenie [Elements of nonsmooth analysis and quasidifferential calculus]. Moscow, Nauka, 1990. 431 p. (in Russian).

4. Clarke F. H. Optimization and Nonsmooth Analysis. SIAM, 1990. 308 p. (Russ. ed.: Moscow, Nauka, 1988. 280 p.).

5. Brayton R. K., Hachtel G. D., Sangiovanni-Vincentelli A. L. A survey of optimization techniques for integrated-circuit design. Proceedings of the IEEE, 1981, vol. 69, iss. 10, pp. 1334-1362.

6. Dudov S. I., Meshchanov V. P. Parametric optimization of the designed devices according to the criteria of cost and quality. Obzory po elektronike. Ser. 1. Elektronika SVCh [Reviews on Electronics. Ser. 1. Electronics Microwave]. Moscow, TsNII "Electronika", 1990, iss. 1 (1512). 40 p. (in Russian).

7. Ioffe A. D., Tikhomirov V. M. Theory of extremal problems. Amsterdam, Nord Holland PC, 1979. 460 p. (Russ. ed.: Moscow, Nauka, 1974. 481 p.).

8. Polovinkin E. S., Balashov M. V. Elementy vypuklogo i sil'no vypuklogo analiza [Elements of convex and strongly convex analysis]. Moscow, Fizmatlit, 2004. 240 p. (in Russian).

9. Rockafellar R. T. Convex Analysis. Princeton, Princeton Univ. Press, 1970. 472 p. (Russ. ed.: Moscow, Mir, 1973. 340 p.).

10. Zukhovitski S. I., Avdeeva L. I. Linejnoe i vypukloe programmirovanie [Linear and Convex Programming]. Moscow, Fizmatlit, 1967. 460 p. (in Russian).

11. Ivanov G. E. Slabo vypuklye funktsii i mnozhestva [Weakly convex sets and functions]. Moscow, Fizmatlit, 2004. 352 p. (in Russian).

12. Poliak B. T. Introduction to optimization. New York, Optimization Software, 1987. 460 p. (Russ. ed.: Moscow, Nauka, 1983. 384 p.).

13. Vial J.-P. Strong and weak convexity of set and funtions. Math. Ops. Res., 1983, vol. 8, no. 2, pp. 231-259.

14. Vasiliev F. P. Metody optimizatsii [Optimization methods]. Moscow, MTsNMO, 2011. 433 p. (in Russian).

15. Bronstein E. M. Approximation of convex sets by polytopes. Journal of Mathematical Sciences, 2008, vol. 153, iss. 6, pp. 727-762.

Cite this article as:

Abramova V. V., Dudov S. I., Osiptsev M. A. The External Estimate of the Compact Set by Lebesgue Set of the Convex Function. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2020, vol. 20, iss. 2, pp. 142-153 (in Russian). DOI: https://doi.org/10.18500/1816-9791-2020-20-2-142-153

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