СТОХАСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗАМКНУТЫХ КРИВЫХ НА ПЛОСКОСТИ Текст научной статьи по специальности «Математика»

УДК 519.245

БСТ: 10.14529/ mmp210103

СТОХАСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЗАМКНУТЫХ КРИВЫХ НА ПЛОСКОСТИ

М.В. Куркина1, В.В. Славский1

1 Югорский государственный университет, г. Ханты-Мансийск, Российская Федерация

Наиболее универсальный метод имитационного моделирования - стохастическое моделирование. Первоначально Энрико Ферми в 1930-х годах в Италии, а затем Джон фон Нейман и Станислав Улам в 1940-х в Лос-Аламосе предложили использовать стохастический подход для аппроксимации многомерных интегралов в уравнениях переноса, возникших в связи с задачей о движении нейтрона в изотропной среде. После начала использования компьютеров произошел большой прорыв, и этот метод стал применяться в самых разных задачах, для решения которых стохастический подход оказался более эффективным, чем другие математические методы. В данной работе изучается форма случайного выпуклого овала на плоскости и более общая задача форма случайной замкнутой кривой на плоскости, исследуется изопериметрическое отношение - отношение квадрата длины кривой к площади ограниченной кривой. Величина этого отношения в силу изопериметрического неравенства ограниченна и характеризует отклонение кривой от окружности. Определяется конечномерное многообразие замкнутых регулярных кривых на плоскости и его бесконечномерный аналог. Изучается вероятностные распределения изопериметрического отношения на них. Основной результат состоит в установлении аналитического закона вероятностного распределения отношения - как распределения Фреше являющиеся частным случаем обобщенного распределения экстремальных значений. Основным используемым методом является разложение Фурье опорной функции множества на плоскости и применение математических пакетов МаШетайса и МаЙаЪ при стохастическом моделировании.

Ключевые слова: изопериметрическое отношения; 'распределения экстремальных значений.

Введение

Геометрические вероятности - один из самых старых и востребованных объектов исследований в приложениях теории вероятностей, достаточно упомянуть монографию [1]. В последние десятилетия интерес к предмету значительно возрос под влиянием приложений в физике, астрономии, биологии, кристаллографии, петрографии и т.д.

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

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

Рис. 1. Кривая с опорной функцией Л,(ф) = 3/2 + 1/2сов(ф) + 1 сов(2ф)

соответствующую статью в 1927 году, дальнейшая работа была проделана Фишером и Типпеттом в 1928 году и Гамбелем в 1958 году.

Основным используемым методом является разложение Фурье опорной функции множества на плоскости и применение математических пакетов Ма^еша^са и Ма^аЬ при стохастическом моделировании. Отметим определенную связь между теорией распределения экстремальных значений в статистике и тропической математикой.

1. Представление опорной функции рядом Фурье

Обозначим вслед [2] за множество всех компактных выпуклых подмножеств евклидовой арифметической плоскости R2 через CSC. В отличии от работы [2] мы воспользуемся хорошо разработанным аппаратом разложения в ряд Фурье [3-5].

Определение 1. Пусть K £ CSC, п(ф) = [cos(0), sin(0)]; ф £ [0, 2п] - единичный вектор на плоскости R2. Опорной функцией hK(ф) множества K называется функция:

hK(ф) = тах(п(ф), r), ф £ [0, 2п],

r€K

где (п(ф),г) - скалярное произведение вектора п(ф) и радиус вектора точки r £ K. По функции h(^), £ [0, 2п] можно найти границу множества K, как огибающую семейства опорных прямых:

x cos(ф) + y sin(ф) = h(ф), ^ £ [0, 2п].

Уравнение огибающей находится из системы уравнений [6]:

x cos^) + y sin(ф) = Л,(ф), —x sin(ф) + y cos^) = h/(ф).

В параметрической форме оно имеет вид;

x = Л,(ф) cos(ф) — sin^)h'(ф),

г(ф) =

y = cos^)h'^) + h(ф)sin(ф). (1)

Уравнение (1) мы будем рассматривать также, независимо от К, просто как удобный способ задания замкнутой кривой с помощью периодической функции Л-(ф), которую будем называть опорной, как и в случае выпуклой кривой. (Пример такой кривой на рис. 1).

Касательный вектор к кривой при данной параметризации равен

% = т =

— sin^) (Л,"(ф) + h(ф)) cos(ф) (h//(ф) + h(ф))

Элемент длины кривой (1) принимает вид ¿8 = |^г| = |к''(ф) + к(ф)|^ф, коэффициент к"(ф) + к(ф) = щ равен величине обратной к кривизне кривой (1). Знак к(ф) + к"(ф) дает условие локальной выпуклости или вогнутости кривой, относительно направления обхода кривой и естественной ориентации плоскости. Формула для длины кривой имеет вид:

2п

Ь = I |к(ф) + Н"(ф)1йф. 0

В случае если к(ф) + к''(ф) > 0, формула длины упростится:

2п 2п

ь = I |к(ф) + к"(ф)^ф = ! к(ф)сф, 00

за счет периодичности функции к(ф). Площадь охваченная кривой (1) на плоскости (с учетом знака и кратности) равна:

2п 2п

3 = Ц Цф) (Цф) + к"(ф)) ¿Ф = Ц (Цф)2 - к'(ф)2) 00

Рассмотрим разложение функции к(ф) в ряд Фурье:

Нф) = ^ + ^,Ак<хж(кф + ек), (2)

к=1

где а0 > 0, Ак > 0 - амплитуды, 9к Е [0, 2п] - фазы гармоник разложения. Первое слагаемое - А1 сов(ф + в1) в силу уравнения (1) отвечает параллельному сдвигу кривой на вектор [А1 сов(01),А1 вт^)], поэтому при исследовании свойств кривой инвариантных относительно относительно сдвига его можно отбросить.

Замечание 1. Условие регулярности кривой (1) имеет вид:

2 к=2

условие (3) будет автоматически выполнено, если потребовать более сильное условие:

| - 2 - 1 )Ак > 0, (4)

при этом к(ф) + к''(ф) > 0.

Определение 2. Обозначим через С0 бесконечномерное многообразие всех регулярных замкнутых выпуклых кривых на плоскости коэффициенты Фурье опорных функций к(ф) которых удовлетворяют условию (4), аналогично через СП С С0 обозначим многообразие кривых для которых дополнительно Ак = 0 при к > п.

Замечание 2. Из условия (4) следует ограниченность (А2 — 1)Ак < С отсюда Ак = о (^з), поэтому функция 1г(ф) и кривая К Е Со будут кусочно непрерывно дифференцируемые [5], а если К £ СП, то опорная функция есть тригонометрический полином и кривая (1) аналитична.

Теорема 1. Площадь охваченная кривой К € С0 с учетом кратности обхода, а также длина кривой выражаются через коэффициенты Фурье функции Л,(ф) и равны:

S

п 4

ап -

2^>2 - , L

па0.

(5)

k=2

Доказательство. Проверяется непосредственно.

□

Определение 3. Изопериметрическим отношением замкнутой кривой (1) называется число

Z

L2 4тг

(6)

Известно [7], что Z > 1. Число Z =1, если и только если кривая есть окружность, величину Z часто используют как меру отклонения кривой от окружности.

Теорема 2. Для регулярных кривых К € С0 справедливо неравенство:

z<®.

~ 5

(7)

Доказательство. Для кривой K G C0 справедлива формула:

Z =

ао

- 2E(fc2 - 1)Ak

k=2

введем для кривых K G C0 обозначение bk =

Ak(k2 - 1)

ао

k > 1 - приведенные коэф-

фициенты Фурье функции Л,(ф). В этих обозначениях:

< И(ф) = а0

k=2

1 ь

2 k=2 k - 1

й(ф) + й"(ф) = ао

Очевидно, что:

о ~ cos (кф + Ок)

k=2

> 0, Z

1

,2

1 Z fc2_l k=2

bk

k=2

k2 1

k=2

k2 1

2

1

г2

22 1

1 б'

отсюда: Z < —.

5

□

2

а

о

ж

Рис. 2. График функции к(ф) к примеру 1: Я = 1, Ь = 2, 5, Z = 1, 24435

Пример 1. Пусть опорная функция кривой (рис. 3) равна:

Г Я + Ь сов(ф), если ф Е [-п/2,п/2], = ( Я, если фЕ [-п/2,п/2],

тогда Z 1 + >/6

(2пЯ + 2Ь)2 Ь

———-^ггт- гешая неравенство А > — получим, что при — >

4(пЯ2 + 2ЯЬ) 5 Я

5

п ~ 2,16738 кривая К Е С0. Опорная функция к(ф) примера 1 непрерывная и кусочно дифференцируемая, разложение ее в ряд Фурье имеет вид:

и АЛ 2(Ь + 7ГЯ) 1 (лЛ/2Ь (| СО8(20) - ^ соз(40) + ^ соз(60) - ^ со8(80) + ...

п\ф) =--\--Ьс08{ф)-\---—-—--

п 2 п

Рис. 3. Кривая к примеру 1: Я =1, Ь = 2, 5, Z = 1, 24435

6

Замечание 3. Как следует из доказательства теоремы 2 равенство Z = — для

_ 5 _

"1 1

кривой К Е Со означает, что && = 0 при к > 3 и Ь2 = = 1

Неравенство (7) достигается, примером служит двугольник:

х = — (9сов(ф) — сов(Зф)), 12

у = —г-

2 + g cos (20)

Рис. 4. Рисунок к замечанию 4: Ь2 = М0) = 1 • + | соз(2ф)] Определение 4. Пусть ш = [а0, {Ьк}П=2,№}П=2] система чисел таких, что:

k=2

=2'

1

2;

а0 >0, Ък > 0, bk < -, 0 < 9k < 2тг, А: = 2,... , /г,

и пусть Пп(1) совокупность всех таких ш при а0 = 1. Сопоставим набору чисел ш € Пга(1) опорную функцию Нш (ф) - тригонометрический полином:

К (ф)

1 n b

k=2

и соответствующую кривую Кш € СП определенную с точностью до движения.

Следствие 1. Имеется взаимнооднозначное соответствие между кривыми К € СП определенными с точностью до гомотетии и точками подмножества Пга (1) С

R+ 1 х [0, 2п]п 1 евклидова пространства.

Определение 5. Множество Пга(1) вместе с метрикой из евклидова пространства дает естественную параметризацию многообразия замкнутых регулярных кривых СП. Данное определение мотивировано тем, что приведенные коэффициенты Фурье {Ьк}П=2 опорной функции Л,(ф) инвариантны относительно гомотетий замкнутой выпуклой кривой на плоскости, а фазы {0к}П=2 определены с точностью до слагаемого. Аналогичным образом определим (1) - параметризацию бесконечномерного многообразие С0 - замкнутых регулярных выпуклых кривых.

n

2. Распределение значений изопериметрического отношения Z

на пространстве C(

Определение 6. Обозначим через ПП(1) проекцию (или фактор пространство) Пп (1):

' 1

{bk}eRn+~l: гдеЪк> 0 < -

k=2

П'п-

Пространство Пп представляет собой симплекс

П :

п

Функция Z задана на ПП формулой:

' 1

{bk} G : гдеЪк> 0 ]>>fc < -

k=2

Z(bk) = — 1

п 1,2 ' 1 _ 2 V к к=2

6

и принимает значения 1 < А < -. Пусть с € [1,6/5] рассмотрим подмножество

5

ПП(с) = [{Ьк}П=2 : Z(Ьк) < 1 + с] . Функция распределения значений Z равна отношению объемов:

pm та '

— PRK"TTXTTTYYR ппт^РА/Г РХТ

2n-1n!

где = п_1—г - евклидов объем симплекса На-

пример 2. Рассмотрим подробнее случай п = 2. Множество П'2(c) С Я2 на плоскости задано в виде:

О < Ь2, Ъ3, Ъ2 + Ъ3 < --72^<1 + С

1 - 2 ( +

х 3 8

или

1 62, b2 . c

0<ь2, ъ3-ъ2 + ъ3<~- +

2 3 8 2(1 + с) На рис. 5 изображена область П'2(с) при различных значениях с. Площадь (-¡^

2V / 1 1 | ' ''1 | '< 1 J.l 1-11 I 1 I }1 Л J I 1 V I 1 1-1 /1 л о. 1 1 .1 W I I I./I „I, I } I 2 у | -Т

равна:

/ 1 \ 192V3 - 18\/39 — 12V3 - 24\/13 - - 24\/l3 + 4^ + 18\/39 + 12^3 2 V17У 4396 +

242Убтг + 121V6sin"1 (Зл/З - 4)) - 121^6sin"1 (4 + Зу/Щ) + 4356 '

Вестник ^ЭУрГУ. Серия «Математическое моделирование 45

и программирование» (Вестник ЮУрГУ ММП). 2021. Т. 14, № 1. С. 39-49 5

Рис. 5. Рисунок к примеру 2: с = с = рр с = |

Замечание 4. Видно, что уже в случае n = 2 формула для функции распределения F2 (c) имеет составной, сложный вид. В принципе еще остается возможность того, что формула упростится, если перейти к пределу:

F» (c) = lim Fn (c). (8)

В следующем пункте будет дано представление эмпирической функции F» (с) распределения значений Z путем стохастического моделирования.

3. Стохастическое моделирование кривой на плоскости

Используя математические пакеты Matlab и Mathemaica были проведены стохастические моделирование замкнутых регулярных СП кривых на плоскости с целью получения эмпирической функции распределения значения изопериметрического отношения Z.

При помощи генератора случайных чисел в Matlab строились коэффициенты разложения Фурье функции h(0). Затем вычислялись уравнение кривых (1). Отсеивались нерегулярные кривые и вычислялись L, S, Z с помощью формул (4) - (7). Результатом этих вычислений на Matlab были серии длины 600 - 800 случайных тригонометрических полиномов h(0) степени n < 18, по которым вычислялись распределения функции Z. Далее, эти серии исследовались с помощью удобной процедуры из Mathematica:

FindDisiribuiion[dai].

Которая подбирает наиболее близкое вероятностное распределений из 23 стандартных, имеющихся в пакете Mathematica и смеси их (MixtureDistribution). На первое место чаще всего выходило трех параметрическое распределение Freche:

(x-ß\~a / \ -а-1

ае'У—) ' ^

FrechetDistribution{a, ß, ¡jl,x) = ^ ----? х > /х;

0, True.

В таблице приведены результаты для значений параметров {а, в, Y} распределения Freche в ходе 4 испытаний. На рис. 6 и 7 приведены графики соответствующих эмпирических распределений величины Z.

Таблица

Результаты моделирования

N а /3 ¡ь P-Value

1 4,01350 0,011038 0,994077 0,996642

2 6,60121 0,019339 0,986328 0,982874

3 4,46088 0,012403 0,992845 0,971323

4 5,29214 0,014452 0,990937 0,998749

150т

1 1.005 1.01

Рис. 6. Графики распределений Freche из таблицы

Рис. 7. Распределение значений Z по результатам стохастического моделирования Заключение

Основываясь на экспериментальных результатах стохастического моделирования и следуя философским идеям В.И. Арнольда [8] об <экспериментальной математике» можно выдвинуть гипотезу, что существуют такие значения параметров {а, в, Y}, что распределение значений (8) функции Z = на бесконечномерном многообразии CSC - всех плоских замкнутых выпуклых кривых (не обязательно регулярных, не обязательно из класса Co) подчиняется закону распределения Фреше с этими параметрами.

Статья посвящается 65-летнему юбилею профессора С.Г. Пяткова, руководителя научной школы <Математическое моделирование в решении задач естество-

знания и социально-экономической сферы» Югорского государственного университета.

Работа проводилась при финансовой поддержке Российского фонда фундаментальных исследований (код проектов 18-47-860016, 18-01-00620), при поддержке Научного Фонда ЮГУ № 05.5/19-ЮГУ-105.

Литература

1. Сантало, Л. Интегральная геометрия и геометрические вероятности / Л. Сантало. - М.: Наука, 1983.

2. Marckert, J.F. Compact Convex Sets of the Plane and Probability Theory / J.F. Marckert, D. Renault // ESAIM: Probability and Statistics. - 2014. - V. 18. - P. 854-880.

3. Зигмунд, А. Тригонометрические ряды. Т. 1,2 / А. Зигмунд. - М.: Мир, 1965.

4. Бари, Н.К. Тригонометрические ряды / Н.К. Бари. - М.: Физматлит, 1961.

5. Бесов, О.В. Тригономерические ряды Фурье / О.В. Бесов. - М.: МФТИ, 2004.

6. Залгаллер, В.А. Теория огибающих / В.А. Залгаллер. - М.: Наука, 1975.

7. Полиа, Г. Изопериметрические неравенства в математической физике / Г. Полиа, Г. Се-ге. - М.: Физматлит, 1962.

8. Арнольд, В.И. Математика: границы и перспективы. - М.: ФАЗИС, 2005.

Мария Викторовна Куркина, кандидат физико-математических наук, доцент, институт «Цифровая экономика», Югорский государственный университет (г. Ханты-Мансийск, Российская Федерация), mavi@inbox.ru.

Виктор Владимирович Славский, доктор физико-математических наук, профессор, институт «Цифровая экономика», Югорский государственный университет (г. Ханты-Мансийск, Российская Федерация), slavsky2004@mail.ru.

Поступила в редакцию 10 октября 2020 г.

MSC 00А71 DOI: 10.14529/mmp210103

STOCHASTIC MODELLING OF CLOSED CURVES IN THE PLANE M.V. Kurkina1, V.V. Slavsky1

:Ugra State University, Khanty-Mansiysk, Russian Federation E-mails: mavi@inbox.ru, slavsky2004@mail.ru

The most versatile simulation method is stochastic simulation. Initially, Enrico Fermi in the 1930s in Italy, and then John von Neumann and Stanislav Ulam in the 1940s at Los Alamos, proposed using a stochastic approach to approximate multidimensional integrals

in transport equations arising in connection with the problem of the motion of a neutron in an isotropic medium. After the start of the use of computers, there was a big breakthrough, and this method began to be applied in a wide variety of problems, for which the stochastic approach proved to be more effective than other mathematical methods. In this paper, we study the shape of a random convex oval in a plane and a more general problem, the shape of a random closed curve in a plane, investigate the isoperimetric ratio - the ratio of the squared length of a curve to the area of bounded curve. The value of this ratio, due to the isoperimetric inequality, is limited and characterizes the deviation of the curve from the circle. A finite-dimensional manifold of closed regular curves in the plane and its infinite- dimensional analog are defined. The probability distributions of the isoperimetric ratio on them are studied. The main result is to establish an analytical law for the probability distribution of the ratio - as Frechet distributions, which are a particular case of the generalized distribution of extreme values. The main method used is the Fourier expansion of the support set function on the plane and the use of mathematical packages Mathematica and Matlab for stochastic modeling.

Keywords: isoperimetric relationship; distribution of extreme values.

References

1. Santalo, L. Integralnaya geometria i geometricheskie verojatnosti [Integral Geometry and Geometric Probabilities]. Moscow, Nauka, 1983. (in Russian)

2. Marckert, J.F. Compact Convex Sets of the Plane and Probability Theory. ESAIM: Probability and Statistics, 2014, vol. 18, pp. 854-880. DOI: 10.1051/ps/2014008

3. Sigmund, A. Trigonometricheskie ryady [Trigonometric Series]. Moscow, Mir, 1965. (in Russian)

4. Bari, N.K. Trigonometricheskie ryady [Trigonometric Series]. Moscow, Fizmatlit, 1961. (in Russian)

5. Besov, O.V. Trigonomericheskie ryady Furie [Trigonomeric Fourier Series]. Moscow, MIPT, 2004. (in Russian)

6. Zalgaller, V.A. Teorija ogibajushhih [The Theory of Envelopes]. Moscow, Nauka, 1975. (in Russian)

7. Polya, G., Szege, G. Izoperimetricheskie neravenstva v matematicheskoj fizike [Isoperimetric Inequalities in Mathematical Physics]. Moscow, Fizmatlit, 1962. (in Russian)

8. Arnold, V.I. Matematika: granicy i perspektivy [Mathematics: Boundaries and Perspectives]. Moscow, FAZIS, 2005. (in Russian)

Received October 10, 2020