Научная статья на тему 'Аппроксимация многоугольников наилучшими наборами кругов'

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

CC BY
407
33
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЧЕБЫШЁВСКИЙ ЦЕНТР / НАИЛУЧШАЯ СЕТЬ / ХАУСДОРФОВО ОТКЛОНЕНИЕ / ВЫЧИСЛИТЕЛЬНАЯ ГЕОМЕТРИЯ / CHEBYSHEV CENTER / BEST NET / HAUSDORFF DISTANCE / COMPUTATIONAL GEOMETRY

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

Изучаются наилучшие аппроксимации многоугольников на плоскости наборами кругов. Основным компонентом их построения являются наилучшие сети, обобщение понятия чебышёвского центра и оптико-геометрический подход

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

Approximation of polygons with the best set of circles

The best approximations of flat polygons with circles are considered. The main component of their construction is the best net. It is the generalized case of the Chebyshev center. About the best segmentation based on the optics-geometrical approach.

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

Серия «Математика»

2013. Т. 6, № 3. С. 72-87

Онлайн-доступ к журналу: http://isu.ru/izvestia

УДК 514.174.2:3

Аппроксимация многоугольников наилучшими наборами кругов *

П. Д. Лебедев

Институт математики и механики УрО РАН

Д. С. Бухаров

Институт динамики систем и теории управления СО РАН

Аннотация. Изучаются наилучшие аппроксимации многоугольников на плоскости наборами кругов. Основным компонентом их построения являются наилучшие сети, обобщение понятия чебышёвского центра и оптико-геометрический подход.

Ключевые слова: чебышёвский центр; наилучшая сеть; хаусдорфово отклонение; вычислительная геометрия.

В теории управления [11] часто требуется провести аппроксимацию множеств со сложной геометрией более удобными в работе фигурами [19, 20]. Одним из самых простых и одновременно информативных способов подмена плоского множества набором из конечного числа кругов равного радиуса [14]. Схожие задачи об аппроксимации множеств эллипсами рассматривались в работах А.Б. Куржанского и его учеников [21, 9]. Вопросы существования и единственности оптимальных покрытий множеств шарами в различных пространствах, в том числе ив М2, были изучены А.Л. Гаркави [3, 4, 5, 7, 6] и Е.Н. Сосовым [16, 17, 18].

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

* Работа выполнена при частичной финансовой поддержке РФФИ: гранты 13-01-96055-р_урал_а, 12-07-13116_ОФИ_м_РЖД, 12-07-33045_мол_а_вед, 11-07-00245_а, 12-07-31080_мол_а, и проекта 12-С-1-1017/3 «Качественная теория и численные методы для задач динамики, управления и оптимизации», гранта президента РФ по поддержке ведущей научной школы №НШ-5927.2012.1, проекта молодых ученых и аспирантов УрО РАН «Оптимальные конструкции и аппроксимации в динамических игровых задачах»

Введение

проксимируют многоугольники. Важным практическим приложением данного исследования является решение задачи о наилучшем размещении логистических центров, и сегментации логистических зон обслуживания [10, 12]: размещение почтовых отделений на территории населенного пункта; размещение магазинов, предоставляющих товары первой необходимости; планирование мест размещения пунктов утилизации бытовых отходов; размещение центров захоронения радиоактивных (или токсичных) отходов; размещение спасательных пунктов (или пунктов скорой медицинской помощи) в местах, из которых обслуживание прилегающей территории производилось бы за время не превышающее минимально возможное.

С целью автоматизации решения подобных задач разработан модуль программного комплекса «ВИГОЛТ» [2], который применялся, например, при решении следующих задач: идентификация и сегментация зон утилизации старых автомобилей — определение границ зон утилизации, обслуживаемых определенным приемным пунктом; определение наилучшего места расположения социальных логистических объектов [12], к которым относятся аптечные пункты, школы, отделения банков, магазины и т.п.

1. Чебышёвский центр множества

Введем следующее определения.

Определение 1. Чебышёвским центром ограниченного множества M С R2 называется такая точка c(M), что

sup ||x — c(M)|| = sup inf ||x — y||. (1.1)

xeM xeM yeR2

Здесь и далее норма вектора x = (x, у) понимается в евклидовой метрике ||х|| = ||(ж, у)|| = л/ж2 + у2, расстояние между точками х и у считается равным ||x — у||. В [4, 5] показано, что при этом чебышёвский центр всегда существует и единствен.

Определение 2. Чебышёвским радиусом R(M) компактного множества M С R2 называется величина (1.1).

Определение 3. Хаусдорфовым отклонением множества A от множества B называется число

h(A, B) = sup inf ||x — y||. xeA yeB

Чебышёвский центр множества является той точкой на плоскости, хаусдорфово отклонение h(M, {x}) от которой множества M минимально. Для любого компактного множества он существует и является единственным [3].

Определение 4. В пространстве М2 п-сетью называется непустое множество, состоящее не более чем из п точек в М2 (см. [17, с. 14]).

Обозначим через Хп множество всех п-сетей пространства М2.

Определение 5. Наилучшей п-сетью [5, 16] компактного множества М С М2 называется такая сеть Б*, для которой выполняется

Н(М,Б*) = М Н(М,Б). (1.2)

Задача 1. Пусть задано ограниченное замкнутое множество М С М2 с кусочно-гладкой границей. Рассматриваются наборы кругов 0(х1,г), 0(х2, г),0(хп, г) одинакового радиуса г с центрами в точках XI, х2, ..., хп соответственно (будем обозначать далее такие наборы через X). Требуется найти такой набор кругов 2(Х, г) = У 0(х^,г), который

г=1 ,п

при наименьшем г удовлетворяет включению

М С Н(Х, г).

Набор кругов 2, отвечающий условиям задачи 1 будем называть наилучшим для множества М при заданном п. В дальнейшем полагаем, что в задаче 1 целевое множество есть многоугольник с вершинами а1, а2,..., ат. Вообще под многоугольником понимаем замкнутое односвязное множество на плоскости, граница которого состоит из конечного числа отрезков. В работах [3, 4] показано, что наилучшая сеть в евклидовом пространстве существует для любого компактного множества, но не всегда единственна.

Соответственно в общем случае решение задачи 1 сводится к построению наилучшей п-сети Б многоугольника М с границей дМ = [а1, а2] и [а2, а3] и ... и [ат_1, ат] и [ат, а1]. Точки х^ е Б совпадают с центрами кругов из набора 2, а их радиус г равен хаусдорфовому отклонению множества М от Б. В случае если наилучшая п-сеть содержит менее чем п* < п точек, то положение (п — п*) кругов может быть произвольным. Мы ограничимся рассмотрением задачи 1 при п равном двум и трём. Заметим, что её решение при п = 1 сводится к нахождению чебышёвского центра и чебышёвского радиуса множества.

Задача 1 имеет важное практическое применение. Пусть для заданной области М С М2 требуется найти оптимальное расположение п логистических центров. А именно расположение, при котором потребитель, находящийся в любой точке множества М, мог бы проехав минимально возможное расстояние, оказаться в ближайшем из них. Тогда точки х; е 1, г = 1,п являются местами оптимального расположения логистических центров. Зоной обслуживания г-го центра (т.е. той частью области М, из которой ближе всего добираться именно до него) является 0(х, г) П М.

2. Аналитические выражения для наилучших сетей

Рассмотрим сначала самые простые многоугольники, т.е. треугольники.

Если М — остроугольный треугольник, то с(М) есть центр описанной вокруг него окружности. Чебышёвский центр тупоугольного треугольника совпадает с серединой наибольшей из его сторон. Если треугольник прямоугольный, то центр описанной вокруг него окружности совпадает с серединой его гипотенузы. Подробнее об отыскании чебышёвского центра для многоугольников см. в [14].

Для дальнейших рассуждений приведем одну оценку для наилучших п-сетей множества.

Лемма 1. Пусть множество М содержит п+1 точек = 1,п + 1,

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

V* = 1,п + 1, V? = 1,п + 1 {гф ^ => (||^ - ^-|| ^ К). (2.1)

Тогда для хаусдорфова отклонения М от любой п-сети Б выполняется оценка

Н(М,Б) ^ К/2. (2.2)

Доказательство. Допустим (2.2) не выполняется. найдётся п-сеть 5*, для которой Ь(М,Б*) = г < К/2. Поскольку точки = 1,п + 1

лежат в М, то для каждой из них найдётся точка х* е Б, такая что ||^ — х* || ^ г. Из (2.1) следует, что одна точка не может находиться от двух различных элементов из набора vn+l1 на расстоянии, не превышающем г < К/2. Следовательно сеть Б * должна содержать минимум п+1 точек (каждая из которых лежит в круге (Э(^,г), г = 1,п + 1). Получилось противоречие. □

Из леммы 1 следует простой критерий определения наилучшей сети. Пусть для множества М найдены п+1 точек, для которых выполняется оценка (2.2). Если для некоторой п-сети Б имеет место

Н(М,Б) < К/2, (2.3)

то эта сеть является наилучшей для М.

В общем случае задача о построении наилучшей п-сети даже при п = 2 решается только численно. Однако для одного класса треугольников можно указать аналитическое решение.

Теорема 1. Пусть М — треугольник с вершинами а1, а2, а3 и длинами сторон А1 = ||а2 — а3||,А2 = |а1 — а3||,А3 = |а1 — а2||. Если выполняется оценки

А2 ^ А2 + А2 (2.4)

и

А2 ^ Аз, (2.5)

то

Б = {хг }2г=1 = {Ы + аг) / 2} 2=1 (2.6)

является наилучшей 2-сетью для М. Здесь § — пересечение прямой, содержащей отрезок [а1, а3], со срединным перпендикуляром 2 к отрезку [а1, а2].

Доказательство. Заметим, что условие (2.4) означает, что треугольник тупоугольный либо прямоугольный, и угол при вершине а1 не меньше п/2. Условие (2.5) означает, что сторона [а1, аз] не меньше, чем [а1, а2], а следовательно срединный перпендикуляр Л к отрезку [а1, а2] пересекает сторону [а1, а3] в некоторой точке % (см. рис. 1). В случае равенства сторон А1 и А2 точка % совпадает с А3.

Рассмотрим треугольник Да^а2. Поскольку точка § лежит на срединном перпендикуляре к отрезку [а1, а2], то две его стороны, прилегающие к ней, равны. Обозначим К = ||а1 — §|| = ||а2 — §||.

Покажем, что отклонение множества М от любой 2-сети Б* не может быть меньше К/2. Для этого воспользуемся леммой 1. Расстояние между вершинами аг ,г = 1, 2 треугольника М и точкой § не меньше К. Действительно, по построение в треугольнике Да1§а2 две стороны равны К, а третья, [а1, а2], лежит против тупого (или прямого) угла, а значит не меньше их. Значит их можно взять в качестве точек } в условии (2.1). Поэтому для любой 2-сети Б * выполняется неравенство (2.2).

Следовательно для того чтобы доказать, что 2-сеть Б, заданная выражением (2.6) является наилучшей, достаточно показать, что выполняется оценка (2.3).

Рис. 1

Действительно, можно представить М в виде объединения трёх треугольников

М = Да^шз и Да2§ш3 и Да2§а3, (2.7)

где т3 — середина стороны [а1, а2]. Треугольники Да1 §т3 и Да2§т3 по построению является прямоугольными с длиной гипотенузы равной К. Точки х1 и х2 лежат на серединах гипотенуз, а значит совпадают с чебышёвскими центрам треугольников. Поэтому выполняются вложения

Да^тз С 0(хьК/2), Да2§тз С 0(х2,К/2). (2.8)

Если А2 > Аз, то треугольник Да2§аз является либо тупоугольным, либо прямоугольным, при это его наибольшая сторона совпадает с отрезком [§, а2], равным по длине К. Если А2 = Аз, то точки § и а2 совпадают, и треугольник вырождается в отрезок [§, а2]. В обоих случаях точка х2 является чебышёвским центром и для Да2§аз, а значит

Н(Да2§аз, {х2}) = К/2. (2.9)

Из представления (2.7) и вложений (2.8) и равенства (2.9) следует, что для расположения треугольника М и сети Б имеет место соотношение

М С 0(х1, К/2) и 0(х2, К/2) = Б + 0(0, К/2).

Из него следует оценка (2.3). Значит 2-сеть, заданная в (2.6), — наилучшая для М (при этом Н(М, Б) = ||а — %\\/2, г = 1, 2). □

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

Задачу о построении наилучшей 3-сети треугольника естественно попытаться свести к более простой задаче о построении его чебышёвского центра. В некоторых случаях это возможно

Теорема 2. Пусть М — треугольник с вершинами а1, а2, а3 и углами при них а1,а2,а3. Если выполняется оценки

шт{а1 ,а2,а3} ^ п/6 (2.10)

и

тах{а1,а2,а3} ^ п/2, (2.11)

то

Б = {хг }3=1 = {(аг + с)/2}3=1, (2.12)

где с — центр описанной вокруг треугольника окружности, является наилучшей 3-сетью для множества М.

Доказательство. Покажем, что отклонение множества М от любой 3-сети Б * не может быть меньше К/2, где К — радиус описанной вокруг М окружности. Для этого воспользуемся, как и при доказательстве

предыдущей теоремы, леммой 1. Из условия (2.11) следует что треугольник M остроугольный (либо прямоугольный), а значит точка с принадлежит самому треугольнику (см. рис. 2). При этом

||а* — с|| = R, г = 1Д (2.13)

Из неравенств (2.10), (2.11) и того факта, что функция sin а монотонно возрастает на отрезке [0,п/2] следует, что для синусов углов выполняется оценка

sinoii £ [sin(7r/6), sin(7r/2)] = [0.5,1], г = 1,3. (2-14)

Исходя из теоремы синусов для длин сторон треугольника выполняются равенства

Ai = 2R$mcti, г = 1,3. (2.15)

Из (2.14) и (2.15) вытекают оценки

Ai^R, г =173. (2.16)

Равенства (2.13) и неравенства (2.15) означают, что набор {ai, a2, a3, с} можно взять в качестве точек {vi} в условии (2.1). Поэтому для любой 3-сети S* выполняется неравенство (2.2).

Рис. 2

Покажем, что выполняется оценка (2.3) для 3-сети (2.12).

Поскольку центр описанной окружности принадлежит треугольнику М, то его можно представить в виде объединения шести прямоугольных треугольников

М = Да1СШ2иДа1СШзиДа2ст1 иДа2ст3иДа3ст1 иДа3ст2. (2.17)

Точка Х1 является серединой гипотенузы Х1 с треугольников Да1ст2 и Да1стз. Поэтому она является центром окружности, описанной вокруг этих треугольников. Поскольку длина гипотенузы [х1, с] равна

|х1 с| = Я, радиус описанной окружности равен Я/2, а значит выполняются вложения

Да1ст2 С 0(х1,Я/2), Да1ст3 С 0(х1,Я/2). (2.18)

Аналогичные рассуждения для точки Х2 и Хз показывают, что имеет место

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

Да2стз С 0(х2,Я/2), Да2стз С 0(х2,Я/2), (2.19)

Да3ст1 С 0(хз,Я/2), Дазст2 С 0(хз,Я/2). (2.20)

соответственно.

Из (2.17)—(2.20) следует, что для расположения треугольника М и сети 5 имеет место соотношение

М С 0(х1, Я/2) и 0(х2, Я/2) и 0(хз, Я/2) = 5 + 0(0, Я/2).

Из него в свою очередь следует оценка (2.3). Следовательно сеть , заданная в (2.12) является наилучшей. При этом Н(М,5) = ||с — аг||/2, г = 1Д □

С помощью теоремы 2 можно строить наилучшие 3-сети для ряда широко встречающихся в различных геометрических задачах треугольников. Под её условия подходит правильный треугольник и египетский треугольник (со сторонами равными 5г, 4г, 3г, г € [0, +то)). Углы при вершинах последнего равны а1 = п/2, а2 ^ 0.927, аз ^ 0.644, и все лежат на отрезке [п/6,п/2].

3. Численные алгоритмы построения наилучших сетей

В общем случае задачу о построении наилучшей п-сети нельзя решить аналитически. Её можно свести к более простой задаче о разбиении множества М на подмножества [8, 1] и отыскании их чебышёвских центров. В случае небольших п данный метод оптимален с точки зрения вычислительных операций.

Лемма 2. Пусть хаусдорфово отклонение множества М от его наилучшей п-сети равно г. Тогда найдутся такие т множеств М\, г =

1,т, т ^ п, что

М = У Мг, (3.1)

г=1,т

и сеть

5 = {с(Мг )}™1, (3.2)

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

Доказательство. Рассмотрим некоторую наилучшую п-сеть 5* = {х*}^=1 (1 ^ к ^ п) множества М. Обозначим

Рі = {х є М2: У.] = 1 ,к (||х-х*|| ^ ||х — х*||)}

множество точек плоскости, которые лежат от точки х* не дальше, чем от других точек сети 5*. Далее обозначим через I множество тех номеров і, для которых

Рі П М = 0.

Покажем, что имеет место вложение

м с у Мі. (3.3)

іЄІ

Действительно, рассмотрим произвольную точку т є М. К ней среди точек сети 5* найдётся минимум одна ближайшая Х[. Следовательно т є Рі, а значит Рі П М = 0 и І є I. Поскольку любая точка т є М лежит в У Рі, то выполняется вложение (3.3). іЄІ

Пусть через т число номеров в I, через I(і). Обозначим

Mi = Pj П М, 1 = /(г), г = 1, т.

Набор {Мг}т=1 в силу (3.3) удовлетворяет условию (3.1). Из построения множества Мг и наилучшей сети 5* видно, что найдётся точка х* € 5*, для которой

Н(Мг, {х*}) < г.

Из определения чебышёвского центра с(Мг) следует, что хаусдорфово отклонение Мг от него не больше, чем от любой точки на плоскости, в том числе и х* :

Н(Мг, {с(Мг)}) ^ г. (3.4)

Получаем оценку из вложения (3.3) и неравенства (3.4).

h(Mi,{S}) ^ maк_h(Mi,{c(Mi)}) ^ г. (3.5)

г=1,т

Выражение (3.3) означает, что для набора множеств {Мг}г=1 выполняется условие (3.1). А оценка хаусдорфова отклонения (3.5) доказывает, что сеть 5, заданная выражением (3.2), является наилучшей п-сетью для множества М. □

Рис. 3

4. Примеры решения задачи 1

На базе методов, описанных в пункте §3, разработаны численные методы построения наилучших п-сетей 5 при п = 2 и п = 3. Первоначальные оценки сетей получается при переборе составляющих многоугольники треугольников с использованием теорем 1 и 2. Строится решение задачи 1 в набора кругов 2, центры которых с точками сети 5.

Пример 1. Задан выпуклый шестиугольник М с вершинами:

{аг}6=1 = -1)> (0, -1-5), (1, -1), (4.5,1), (-2,1.5), (-4,1)}. (4Л)

Требуется построить для М решение задачи 1 при п = 2 .

Решение задачи сводится к нахождению наилучшей 2-сети 5 = {х1, Х2} множества М и его хаусдорфова отклонения от 5. Построенный численными методами результат: Х1 ~ (-1.857, -0.178), Х2 ~ (2.202,

0.163), г = Н(М, 5) & 2.45. Множество М, система кругов 2 = 0(х1, г)и 0(х2,г) и наилучшая 2-сеть 5 = {х1, Х2} показаны на рис. 3.

Пример 2. Пусть в задаче 1 множество М — выпуклый шестиугольник с набором вершин (4.1), количество кругов задано п = 3.

х

Рис. 4

Наилучшая 3-сети М, построенная численными методами, состоит из точек: х1 & (-2.192, 0.854), х2 & (0.084, -0.351), х3 & (2.691,1.139), г = Н(М,5) & 1.81. Множество М, система кругов 2 = 0(х1,г) и 0(х2, г) и 0(х3,г) и наилучшая 3-сеть 5 = {х1, х2, х3} показаны на рис. 4.

Рис. 5

Пример 3. Задан невыпуклый семиугольник М с вершинами:

{а}=1 = {(-1.5, -2.5), (1, -2), (3, 2), (0,1.5), (-2, 2), (-3,1), (-2, 0.5)}.

(4.2)

Требуется для М найти решение задачи 1 при п = 2 .

Решение задачи сводится к нахождению наилучшей 2-сети 5 = {х1, х2} множества М и хаусдорфова отклонения М от 5. Построенный чис-

ленными методами результат: хі & (—0.937, —0.16), х2 & (1.741, —0.052), г = Н(М,Б) & 2.41. Множество М, система кругов 2 = 0(х1 ,г) и 0(х2,г) и наилучшая 2-сеть 5 = {хі, х2} показаны на рис. 5.

Пример 4. Пусть в задаче 1 множество М — выпуклый шестиугольник с набором вершин (4.1), количество кругов задано п = 3.

Наилучшая 3-сети М, построенная численными методами, состоит из точек: х1 & (—1.38, 0.356), х2 & (1.6, 0.937), х3 & (0.045, —1.66), г = Н(М, Б) & 1.76. Множество М, система кругов 2 = 0(х1,г) и 0(х2, г) и 0(х3,г) и наилучшая 3-сеть Б = {х1,х2,х3} показаны на рис 6.

Как видно из примеров расположение кругов из набора 2(Х, г), содержащего три элемента может быть существенно различным. Для шестиугольника с набором вершин (4.1) двое кругов не имеет общих точек между собой (но пересекаются с третьим). Для шестиугольника с набором вершин (4.2) все три круга имеют общее множество.

X

Рис. 6

5. О наилучшей сегментации на основе оптико-геометрического подхода

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

Задача 2. Пусть задана ограниченная область М С М2 с кусочногладкой границей, заданы функция среды f (ж, у), определяющая мгновенную скорость «света» с(х,у) = 1 //(ж, у) в точке (ж,у) € М, и т объектов Ак{хк,Ук)-, к = 1,т, расположение которых заранее неизвестно. Требуется найти такое расположение объектов Ак, сегментировав область М на т зон, чтобы суммарное время достижения каждой точки области М «световыми» волнами было минимально возможным.

Известным способом определяется минимальное время перемещения из точки К (ж, у) € М в точку Ак

Требуется найти оптимальные расположения объектов (х*к, у£), к =

1 ,т, и такое разбиение области М на т сегментов, чтобы достигался минимум целевой функции

В основу алгоритма решения вышеописанной задачи (5.1) - (5.2) положена оптико-геометрическая аналогия, основанная на фундаментальных вариационных принципах механики, идея которой заключается в следующем: на исследуемой ограниченной области М С К2 задаются места расположения двух точек.

Одна из них определяется как источник «света» и из нее производится распространение световой волны, т.е. строятся фронты волны через некоторые промежутки времени Ді. Каждая точка (х, у) Є М, достигнутая фронтом волны, принимается за вторичный источник света, что соответствует принципу Гюйгенса. Построение фронтов продолжается до тех пор, пока не будет достигнута вторая заданная точка.

Так как для каждой точки известно время ее достижения из источника света (исходного или вторичного), то можно, двигаясь в обратном направлении по времени, восстановить кривую, которая согласно принципу Ферма является оптимальным по времени передвижения маршрутом (что на плоскости соответствует прямой линии). Поскольку область М ограниченная, то за конечное число итераций будет либо найдено решение, либо показано, что решения нет (область М «заполнена» фронтами волны, но заданная конечная точка не достигнута).

Для решения задачи о наилучшем размещении логистических центров и сегментации логистических зон обслуживания используется алгоритм, идея которого заключается в следующем:

(5.1)

т

(5.2)

1. Генерация начального расположения. Определяются начальные координаты объектов Ад. (к = 1 ,т) методом случайной генерации положений, так чтобы Ak Е M.

2. Сегментация области. Производится разбиение области M на m сегментов, например методом, представленным в работе [10], или строим диаграмму Вороного [15]. Для каждого сегмента А /е М (к = 1, т) определяется граница.

3. Переопределение координат центров. Определяется новое расположение объектов Ak относительно выделенных сегментов Mk (к = 1,m). Производится конструирование фронтов волны с границы сегментов Mk внутрь области и определяется точка, в которой суммарное время ее достижения минимально. Данная точка и определяет новое расположение центра Ak. Если координаты Ak изменяются, то переходим к сегментации. В противном случае определяется процесс сегментирования завершен.

Если значение функции среды f (x,y) = const, то вышеописанный алгоритм позволяет построить решение задачи о наилучшем размещении заданного количества непересекающихся равных кругов с максимально возможным радиусом (задача об упаковке, рис. 7).

Рис. 7

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

Использование многократной генерации необходимо для выявления глобального минимума в задаче (5.1) - (5.2), без применения которой, как правило, возможно нахождение только локального минимума, при этом вписанные круги будут иметь меньший радиус, чем при оптимальной конфигурации.

Список литературы

1. Болтянский В. Г. Разбиение фигур на меньшие части / В. Г. Болтянский, И. Ц. Гохберг. - М. : Наука, 1971. - 88 с. - (Популярные лекции по математике ; вып. 50).

2. Бухаров Д. С. Программная система «ВИГОЛТ» для решения задач оптимизации, возникающих в транспортной логистике / Д. С. Бухаров, А. Л. Казаков // Вычисл. методы и программирование. - 2012. - Разд. 2. - С. 65—74.

3. Гаркави А. Л. О существовании наилучшей сети и наилучшего поперечника множества в банаховом пространстве / А. Л. Гаркави // Успехи мат. наук. -1960. - Т. 15, вып. 2. - С. 210-211.

4. Гаркави А. Л. О наилучшей сети и наилучшем сечении множеств в нормированном пространстве / А. Л. Гаркави // Изв. АН СССР. Сер. мат. - 1962. -Т. 26, №1. - С. 87-106.

5. Гаркави А. Л. О чебышёвском центре и выпуклой оболочке множества / А. Л. Гаркави // Успехи мат. наук. - 1964. - Т. 19, вып. 6. - С. 139-145.

6. Гаркави А. Л. О методе циклического спуска в задаче наилучшего приближения / А. Л. Гаркави // Мат. заметки. - 1980. - Т. 27, №4. - С. 549—558.

7. Гаркави А. Л. Об условном чебышёвском центре компактного множества непрерывных функций / А. Л. Гаркави // Мат. заметки. - 1973. - Т. 14, № 4. - С. 469-478.

8. Гервер М. Л. О разбиении множеств на части меньшего диаметра: теоремы и контрпримеры // М. Л. Гервер // Математическое просвещение. Сер. 3. -1999. - Вып. 3. - С. 168-183.

9. Гусев М. И. Оценки множеств достижимости многомерных управляемых систем с нелинейными перекрёстными связями / М. И.Гусев // Тр. ИММ УрО РАН. - 2009. - Т. 15, № 4. - С. 82-94.

10. Казаков А. Л. Об одном подходе к решению задач оптимизации, возникающих в транспортной логистике / А. Л. Казаков, А. А. Лемперт // Автоматика и телемеханика. - 2011. - № 7. - С. 50-57.

11. Красовский Н. Н. Позиционные дифференциальные игры / Н. Н. Красовский ,

А. И. Субботин. - М. : Наука, 1974. - 456 с.

12. Лемперт А. А. Математическая модель и программная система для решения задачи размещения логистических объектов / А. А. Лемперт, А. Л. Казаков, Д. С. Бухаров // Управление большими системами. - 2013. - Вып. 41. - С. 270284.

13. Лейхтвейс К. Выпуклые множества / К. Лейхтвейс. - М. : Наука, 1985. - 335 с.

14. Лебедев П. Д. Аппроксимация множеств на плоскости оптимальными наборами кругов / П. Д. Лебедев В. Н. Ушаков // Автоматика и телемеханика. - 2012. -№ 3. - С. 79-90.

15. Препарата Ф. Вычислительная геометрия / Ф. Препарата, М. Шеймос. - М. : Мир, 1989. - 478 с.

16. Сосов Е. Н. Об аппроксимативных свойствах множеств в специальном метрическом пространстве / Е. Н. Сосов // Изв. вузов. Математика. - 1999. - №6. -С. 81-84.

17. Сосов Е. Н. Введение в метрическую геометрию. Ч. 2 : учеб. пособие / Е. Н. Сосов. - Казань : Казан. гос. ун-т, 2008. - 29 с.

18. Сосов Е. Н. Метрическое пространство всех Ж-сетей геодезического пространства / Е. Н. Сосов // Учён. зап. Казан. гос. ун-та. Сер. Физ.-мат. науки. -2009. - Т. 15, вып. 4. - С. 136-149.

19. Ушаков В. Н. Об одном дополнении к свойству стабильности в дифференциальных играх / В. Н. Ушаков, А. А. Успенский // Тр. Математического ин-та им. В. А. Стеклова. - 2010. - Т. 271. - С. 299-318.

20. Ушаков В. Н. Дефект стабильности в игровой задаче о сближении в момент /

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

В. Н. Ушаков, А. Р. Матвийчук, П. Д. Лебедев // Вестн. Удмурт. ун-та. Математика, механика, компьютер. науки. - 2010. - Вып. 3. - С. 87-103.

21. Kurzhanski A. B. Ellipsoidal calculus for estimation and control / A. B. Kurzhanski, I. Valyi. - Boston : Birkhauser, 1997. - 220 p.

P. D. Lebedev, D. S. Bukharov

Approximation of polygons with the best set of circles

Abstract. The best approximations of flat polygons with circles are considered. The main component of their construction is the best net. It is the generalized case of the Chebyshev center. About the best segmentation based on the optics-geometrical approach.

Keywords: Chebyshev center; best net; Hausdorff distance; computational geometry.

Лебедев Павел Дмитриевич, кандидат физико-математических наук, Институт математики и механики УрО РАН, 620990, Екатеринбург, ул. Софьи Ковалевской, д. 16, тел.: (343)3753438 (pleb@yandex.ru).

Бухаров Дмитрий Сергеевич, программист, Институт динамики систем и теории управления СО РАН, 664033, Иркутск, ул. Лермонтова, д. 134, тел.: (3952)453030, (bukharovds@gmail.com).

Lebedev Pavel, Institute of Mathematics and Mechanics of Ural Branch of the Russian Academy of Sciences, 16, S.Kovalevskaja st., Ekaterinburg, 620219, researcher, Phone: (343)3753438 (pleb@yandex.ru).

Bukharov Dmitry, Institute for System Dynamics and Control Theory SB RAS, 134, Lermontov st., Irkutsk, Phone: (3952)453030, programmer, (bukharovds@gmail.com).

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