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

Труды Карельского научного центра РАН №7. 2018. С. 62-67 DOI: 10.17076/mat768

УДК 519.179.4

О СТРУКТУРЕ КОНФИГУРАЦИОННОГО ГРАФА С НОРМАЛЬНО РАСПРЕДЕЛЕННЫМ ПАРАМЕТРОМ СТЕПЕННОГО РАСПРЕДЕЛЕНИЯ СТЕПЕНЕЙ ВЕРШИН

Ю. Л. Павлов

Институт прикладных математических исследований КарНЦ РАН, ФИЦ «Карельский научный центр РАН», Петрозаводск, Россия

Рассматриваются конфигурационные графы с N вершинами. Степени вершин являются независимыми одинаково распределенными случайными величинами, имеющими дискретное степенное распределение с положительным параметром т. Они равны числу занумерованных в произвольном порядке полуребер. Граф строится путем попарного равновероятного соединения полуребер для образования ребер. Изучается подмножество таких случайных графов при условии, что сумма степеней известна и равна п. Пусть т является случайной величиной, имеющей усеченное нормальное распределение на произвольном фиксированном конечном интервале. Для максимальной степени вершины и числа вершин заданной степени найдены предельные распределения в различных зонах стремления N и п к бесконечности.

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

Yu. L. Pavlov. THE STRUCTURE OF A CONFIGURATION GRAPH WITH A NORMALLY DISTRIBUTED PARAMETER OF THE POWER SERIES DISTRIBUTION OF VERTEX DEGREES

We consider configuration graphs with N vertices. The degrees of the vertices are independent random variables identically distributed according to the power law, with a positive parameter т. They are equal to the number of vertex semiedges that are numbered in an arbitrary order. The graph is constructed by joining all of the semiedges pairwise equiprobably to form edges. We study the subset of such random graphs under the condition that the sum of vertex degrees is known and it is equal to n. Let т be a random variable following a truncated normal distribution on an arbitrary fixed finite interval. We obtained the limit distributions of the maximum vertex degree and the number of vertices with a given degree for various zones of N and n tendency to infinity.

Keywords: random configuration graph; vertex degree; limit theorems.

Введение роко используются конфигурационные графы

со случайными степенями вершин (см., напри-Для моделирования современных сложных мер, [14]). В таких моделях предполагается, сетей коммуникаций в настоящее время ши- что степени вершин являются независимы-

ми одинаково распределенными случайными величинами. Многочисленные наблюдения за реальными сетями показали [14], что распределения этих случайных величин можно считать подчиняющимися дискретному степенному закону (аналог распределения Парето). Ребра таких графов образуются путем попарного равновероятного соединения друг с другом различимых полуребер, инцидентных вершинам. Существует несколько способов обеспечить, в случае необходимости, четность суммы степеней; эти способы, как замечено в [15], не влияют на асимптотические свойства графов. Обозначим £ случайную величину, равную степени любой вершины графа. В статье [15] предложена достаточно адекватная модель конфигурационного графа, в которой

Рк = Р{£ = к} = к-т - (к + 1)-

(1)

дели. Желательно, чтобы плотность распределения положительного параметра была симметричной. В настоящей статье впервые рассматривается усеченное нормальное распределение т на интервале [а,Ь], 0 < а < Ь < ж. Обозначим

М = Е£, а2 = Б£.

В качестве математического ожидания £ естественно выбрать середину интервала [а, Ь]:

М=

а + Ь 2 .

(2)

Среднеквадратическое отклонение а может принимать любое положительное значение. Далее нам будет удобно определить а на основе известного правила «трех сигм»:

Ьа

где к = 1, 2,..., т - положительный параметр распределения. В реально существующих сетях обычно т € (1, 2) [14]. Однако в наших последних исследованиях, посвященных моделированию лесных пожаров [12, 13], показано, что представляют интерес и модели, в которых т > 2.

В [9] впервые рассматривались условные конфигурационные графы при условии, что число ребер известно. Для таких графов в ряде работ исследовалось предельное поведение максимальной степени вершины и числа вершин заданной степени при различных соотношениях между стремящимися к бесконечности числе вершин и числе ребер [2, 3, 7]. Полученные результаты можно использовать как для описания сетей с ограничениями на число связей, так и для сетей общего вида. Более того, они могут быть полезны и для изучения поведения и других числовых характеристик, зависящих от степенной структуры графа.

В некоторых работах [10, 11] отмечается, что в процессе роста сетей распределение степеней вершин может меняться вместе с числом вершин и даже быть случайным. В [8] рассматривались графы, в которых степени вершин имеют распределение (1), с изменяющимся параметром т. В статьях [4, 5] этот параметр является случайной величиной, равномерно распределенной на интервале [а, Ь], 0 < а < Ь < ж. Однако более естественным является предположение о том, что случайная величина т должна иметь унимодальное распределение. Поэтому в [1] в качестве закона распределения т было предложено использовать гамма-распределение. Известно, что коэффициент асимметрии этого распределения не равен нулю, что представляется недостатком мо-

а=

6

(3)

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

Таким образом, в работе рассматриваются условные конфигурационные графы с N вершинами, при условии, что сумма степеней вершин равна п, сами степени являются независимыми случайными величинами, имеющими общее распределение (1), где случайный параметр т подчиняется усеченному нормальному закону, заданному на интервале [а, Ь], 0 < а < Ь < ж и имеющему параметры (2) и (3). Для таких графов при N, п ^ ж найдены предельные распределения максимальной степени вершины ) и числа вершин ^г, имеющих степень г. В следующем разделе формулируются полученные результаты в виде теорем 15, а в последнем разделе изложено доказательство этих результатов.

РЕЗУЛЬТАТЫ

Введем вспомогательную случайную величину п такую, что

Рк (А) = Р{п = к} = А крк/В (А), (4) где к = 1,2,..., 0 < А < 1 и

В(А) = £ Ак рк.

к=1

Из (4) и (5) следует, что

те

т(А) = Еп = В-1(А) ^ кАк

Рк,

к=1

а2(Л) = Dn = В-1(Л) J] k2Afcpk - т2(Л). fc=i

Пусть Л является единственным решением уравнения т(Л) = n/N. Обозначим F(ж) функцию нормального распределения с параметрами (M, а), а символами C, Ci, C2 - некоторые положительные постоянные.

Справедливы следующие результаты.

Теорема 1. Пусть N, n ^ ж так, что n/N ^ 1, (n - N)3/N2 ^ C > 0, а r = r(N, n) - наименьшие натуральные числа такие, что NЛrpr+1/p1 ^ Y, где y - некоторая положительная постоянная. Тогда

P{£(N) = r} ^ e-Y, P{) = r+1} ^ 1-e-Y.

Теорема 2. Пусть N, n ^ ж так, что 1 < C1 ^ n/N ^ C2 < M, а r = r(N, n) выбраны так, что

6a^r+1

(b - a)e9(F(b) - F(a))B(Л)га+1 ln r

где 7 - некоторая положительная постоянная. Тогда для любого фиксированного к = 0, ±1, ±2,...

) < г+к} = ехр{-7Лк(1-Л)-1}(1+о(1)).

Теорема 3. Пусть Ж, п ^ те так, что п/Ж ^ М и выполнено одно из следующих условий:

1. а > 2;

2. 0 < а ^ 2, и существует такое 5 > 0,

что N(1 - Л)2+г ^ те.

Тогда

P{| ln Л|£(^) - u < z} ^ e-

где -ж < z < ж, а u = u(N, n) выбраны так,

что

6a| ln л|ae-(u+9)

(b - a)(F(b) - F(а))В(Л)иа+1 ln(u/| ln Л|)

1.

Теорема 4. Пусть N, n ^ ж так, что n/N ^ 1,n - N ^ ж. Тогда для любого фиксированного натурального r ^ 3 и для целых неотрицательных k равномерно относительно (k - Npr(Л))/^Npr(Л) в любом фиксированном конечном интервале

PK = k} = (Npkpk e-NPr (л)(1 + 0(1)).

Теорема 5. Пусть N, n ^ ж так, что 1 < C1 ^ n/N ^ C2 < M, или выполнены условия теоремы 3. Тогда для любого натурального r равномерно относительно ur = (k -Npr (Л))/(а rr VN) в любом фиксированном конечном интервале

P{ßr = k}

1+ o(1) e-u?/2

arr

где

а^г(Л) = Pr(ЛИ 1 - pr(Л) -

(т(Л) - r)2 а2(Л)

Pr (Л)

Легко видеть, что условия теорем 1-5 не исчерпывают все возможные случаи поведения N и п, поэтому получение аналогичных результатов при других условиях предполагается осуществить в следующих работах.

Доказательство теорем

Обозначим / (х) плотность распределения параметра т. Тогда

( (х-М)2

е 2ст2 , х е [а, о];

f (ж) = ^ (b)-F(а))

0,

Отсюда и из (1) получаем, что

P{C = k} = (b)-

ж [a, b].

F (а))

(6)

(x-M)2

X Ja l kX - (kW 1 e "

fb 1

dx.

Полагая х = М + сту, легко получить соотношение:

Ь ь-м

Л к- к-.]у

Проведем в этом равенстве еще одну замену переменной интегрирования: £ = (у + ст 1п к)/л/2, тогда

e 2ст2 dx = , v ^ e 2 x

Ja kx km 2

(7)

x (erf (b-M+^) - erf

е

2 fz

erf (z) = e-t2/2dt

V 7 vWo

a-M+ff2 In fc

(8

- интеграл вероятности. Известно, что при больших z > 0

erf (z) ~ 1 -

(9)

1 + (-1)k

nz

k=1

1 ■ 3 ... (2k - 1)

2k z2k

64

x

2

z

Пусть к ^ ж. Из (7)-(9) находим, что

„ ( а - М + а21п к \ ег/ I -^- I = 1-

л/2<

а

- ехр< -<<" - М)/а + а 1П к)2 > х

2

/2

а

л/П(а - М + а21п к)

1

аМ

2

+ а 1п к + О, 4 а ) \ 1п4 к,

1

Отсюда следует, что

„ .'а - М + а21пк , ег/ I -^- 1=1 -

х ехр < --

1

лДа

а2 1п2 к

/2

X

2

а - М

а2 1п к

1 +

+ О

а - М

а2 1п к 1

л/Па 1п к 2

' х (10)

1п2 к

/2

ехр < --

а2 1п2 к Ь М

л/Па 1п к ~ г 1 2 Ь-М

1+

а2 1п к

1

+ О

1

Ь^к

а2 1п к

Из (2) находим, что (а - М)2 = (Ь - М)2, поэтому

а21п2 к 2

1+

- 1 +

Ь - М\' а2 1п к/

= (Ь - а) 1п к

а - М%

а2 1п к

и из (11) следует, что при к ^ ж ( Ь - М + а2 1п к \

егЧ у-

а - М + а21п к

- ег/

/2.

а

/2

л/Па 1п к

ехр -

а2 1п2 к

2

1+

Ь - М^ а2 1п к

Аналогично

х ехр

уДа

а2 1п2 к

/2

л/Па 1п к

1

2

ЬМ

а2 1п к

1+ + О

Ь - М"

а2 1п к 1

1п2 к

Отсюда и из (10) видим, что

( Ь - М + а2 1п к \

егЧ ;

а - М + а2 1п к

- ег/

/2.

а

/2

а2 1п2 к

л/Па 1п к

ехр

2

1+

Ь - М" а2 1п к

х ехр

а2 1п2 к

2

Ь -

а2 1п к

(11)

х ехр -

а2 1п2 к а М

2

1+

а2 1п к

а - М 1

х( 1 - а-Мк+О( л

к

Отсюда

Ь—а

1

а - М

а2 1п к

+О

1п2 к

'Ь - М + а21пк\

ег/1 у

а - М + а21п к

- ег/ /2

ехр

л/Па 1п к 2

а — М

л/2а

а2 1п2 к (Ь - а)2 4а2

к

1

а2 1п к

+О

1

1п2 к

Тогда из (7) вытекает, что

1 (х —М)2

ка 1п к Аналогично

—е 2ст2 ¿ж = кх

1 (Ь-д)2 / а - М / 1 е —^ ( 1--^^ + О

а2 1п к 1п2 к

1 (х-М)2

е 2^2 =

(к + 1) 1

(к + 1)а 1п(к + 1)

(Ь — а)2

е 4ст2 х

а - М 1

х (1 - „ . „, + О

а2 1п(к + 1)

1п2 к

х

х

X

X

1

X

X

X

Ь

X

X

Ь

X

Проводя элементарные преобразования, из последних двух равенств выводим, что

1

kx

1

(k + 1)x

(x-M )2

dx =

0(1 + 0(1)) (b a)2

= —-—— e

ka+1 ln k

Учитывая (3), отсюда и из (6) получаем, что при k ^ ж

Pie = k}~

6а

(12)

(b - a)e9\/2n(F(b) - F(a))ka+1 ln k'

В статье [6] рассматривалось предельное поведение ) и в условных конфигурационных графах при условии, что сумма степеней вершин равна п, а случайная величина £, равная степени любой вершины графа, при к ^ ж обладает свойством:

pie = k}

d(1 + o(1)) kg(ln k)h '

(13)

где 0,д ^ 1,# + Л, > 1. В этой же статье показано, что примеры таких графов возникают в случаях, когда параметр т распределения (1) равномерно распределен на любом конечном фиксированном интервале или подчиняется гамма-распределению. Из соотношения (12) следует, что в случае усеченного нормального распределения т условия теорем работы [6] выполнены, если в (13) положить у = а + 1, Л, = 1 и

d =

6а

(b - a)e^V^(F(b) - F(а))'

Таким образом, рассмотренный нами случай можно считать еще одним примером случайного конфигурационного графа, удовлетворяющего условиям статьи [6]. Это значит, что утверждения теорем 1-5 непосредственно следуют из соответствующих теорем этой статьи.

Финансовое обеспечение исследований осуществлялось из средств федерального бюджета на выполнение государственного задания КарНЦ РАН (Институт прикладных математических исследований КарНЦ РАН) и при финансовой поддержке РФФИ (проект № 16-01-0005а).

ЛИТЕРАТУРА

1. Лери М. М, Павлов Ю . Л. Об устойчивости конфигурационных графов в случайной среде // Информатика и ее применения (в печати).

2. Павлов Ю. Л. О предельных распределениях степеней вершин в условных Интернет-графах // Дискретная математика. 2009. Т. 21, вып. 3. С. 14-23. doi: 10.4213/dm1057

3. Павлов Ю. Л. Об условных Интернет-графах, степени вершин которых не имеют математического ожидания // Дискретная математика. 2010. Т. 22, вып. 3. С. 20-33. doi: 10.4213/dm1104

4. Павлов Ю. Л. Об условных конфигурационных графах со случайным распределением степеней вершин // Труды КарНЦ РАН. 2016. № 8. С. 62-72. doi: 10.17076/mat313

5. Павлов Ю. Л. Один случай предельного поведения степеней вершин в условных конфигурационных графах// Труды КарНЦ РАН. 2017. № 8. С. 66-75. doi: 10.17076/mat613

6. Павлов Ю. Л. Условные конфигурационные графы со случайным параметром степенного распределения степеней // Математический сборник. 2018. Т. 209, № 2. С. 120-137. doi: 10.4213/sm8832

7. Павлов Ю. Л., Дертишникова Е. Н. О предельном распределении максимальной степени вершины в случайном графе Интернет-типа // Труды КарНЦ РАН. 2010. № 3. С. 59-65.

8. Павлов Ю. Л., Феклистова Е. В. О предельном поведении максимальной степени вершины условного конфигурационного графа вблизи критических точек // Дискретная математика. 2016. Т. 28, вып. 2. С. 58-70. doi: 10.4213/dm1369

9. Павлов Ю. Л., Чеплюкова И. А. Случайные графы Интернет-типа и обобщенная схема размещения // Дискретная математика. 2008. Т. 20, вып. 3. С. 3-18. doi: 10.4213/dm1008

10. Самосват Е. А. Моделирование Интернета с помощью случайных графов: дис. ... канд. физ.-мат. наук. М., 2014. 98 с.

11. Bianconi G, Barabasi A.-L. Bose-Einstein condensation in complex networks // Physical Review Letters. 2001. Vol. 86. P. 5632-5635. doi: 10.1103/PhysRevLett.86.5632

12. Leri M., Pavlov Yu. Power-law random graph's robustness: link saving and forest fire model // Austrian Journal of Statistics. 2014. Vol. 43, no. 4. P. 229-236.

13. Leri M., Pavlov Yu. Forest Fire Models on Configuration Random Graphs // Fundamenta Informaticae. 2016. Vol. 145, no. 3. P. 313-322.

14. Hofstad R. Random Graphs and Complex Networks. Volume One. Cambridge University Press, 2017. 337 p.

15. Reittu H., Norros I. On the power-law random graph model of massive data networks // Performance Evaluation. 2004. Vol. 55, no. 4. P. 3-23. doi: 10.1016/S0166-53/6(3)00097-x

Поступила в редакцию 27.01.2018

66

b

2

e

r^j

References

1. Leri M. M, Pavlov Yu. L. Ob ustoichivosti konfiguratsionnykh grafov v sluchainoi srede [On robustness of configuration graphs in random environment]. Informatika i ee primenenija [Informatics and its applications] (in print).

2. Pavlov Yu. L. On the limit distributions of the vertex degrees of conditional Internet graphs. Discrete Mathematics and Applications. 2009. Vol. 19, iss. 4. P. 349-360. doi: 10.1515/DMA.2009.023

3. Pavlov Yu. L. On conditional Internet graphs whose vertex degrees have no mathematical expectation. Discrete Mathematics and Applications. 2010. Vol. 20, iss. 5-6. P. 509524. doi: 10.1515/dma.2010.031

4. Pavlov Yu. L. Ob uslovnykh konfiguratsionnykh grafakh so sluchainym raspredeleniem stepenei vershin [On conditional configuration graphs with random distribution of vertex degrees]. Trudy KarNTs RAN [Trans. KarRC RAS]. 2016. No. 8. P. 62-72. doi: 10.17076/mat313

5. Pavlov Yu. L. Odin sluchai predel'nogo povedenija stepenei vershin v uslovnykh konfiguratsionnykh grafakh [A case of limit behaviour of vertex degrees in conditional configuration graphs]. Trudy KarNTs RAN. [Trans. KarRC RAS]. 2017. No. 8. P. 66-75. doi: 10.17076/mat613

6. Pavlov Yu. L. Conditional configuration graphs with random parameter of the power-law degree distribution. SB MATH. 2018 (in print). doi: 10.1070/SM8832

7. Pavlov Yu. L., Dertishnikova E. N. O predel'nom raspredelenii maksimal'noj stepeni vershiny v sluchajnom grafe Internet-tipa [On limit distribution of maximum vertex degree in a random graph of Internet type]. Trudy KarNTs

СВЕДЕНИЯ ОБ АВТОРЕ:

Павлов Юрий Леонидович

главный научный сотрудник, д. ф.-м. н., проф. Институт прикладных математических исследований КарНЦ РАН, Федеральный исследовательский центр «Карельский научный центр РАН» ул. Пушкинская, 11, Петрозаводск, Республика Карелия, Россия, 185910 эл. почта: [email protected] тел.: (8142) 781218

RAN [Trans. KarRC RAS]. 2010. No. 3, iss. 1. P. 59-65. doi: 10.17076/mat313

8. Pavlov Yu. L., Feklistova E. V. On limit behaviour of maximum vertex degree in a conditional configuration graph near critical points. Discrete Mathematics and Applications. 2017. Vol. 27, iss. 4. P. 213-222. doi: 10.1515/dma-2017-0023

9. Pavlov Yu. L., Cheplyukova I. A. Random graphs of Internet type and the generalised allocation scheme. Discrete Mathematics and Applications. 2008. Vol. 18, iss. 5. P. 447-464. doi: 10.1515/DMA.2008.033

10. Samosvat E. A. Modelirovanie Interneta s pomosh'ju sluchainykh grafov [Internet modeling with help of random graphs]: DSe (Cand. of Phys.-Math) thesis. Moscow, 2014. 98 p.

11. Bianconi G, Barabasi A.-L. Bose-Einstein condensation in complex networks. Physical Review Letters. 2001. Vol. 86. P. 5632-5635. doi: 10.1103/PhysRevLett.86.5632

12. Leri M., Pavlov Yu. Power-law random graph's robustness: link saving and forest fire model. Austrian Journal of Statistics. 2014. Vol. 43, no. 4. P. 229-236.

13. Leri M., Pavlov Yu. Forest Fire Models on Configuration Random Graphs. Fundamenta Informaticae. 2016. Vol. 145, no. 3. P. 313-322.

14. Hofstad R. Random Graphs and Complex Networks. Volume One. Cambridge University Press, 2017. 337 p.

15. Reittu H., Norros I. On the power-law random graph model of massive data networks. Performance Evaluation. 2004. Vol. 55, no. 4. P. 3-23. doi: 10.1016/S0166-53/6(3)00097-x

Received January 27, 2018

CONTRIBUTOR:

Pavlov, Yury

Institute of Applied Mathematical Research,

Karelian Research Centre,

Russian Academy of Sciences

11 Pushkinskaya St., 185910 Petrozavodsk,

Karelia, Russia

e-mail: [email protected]

tel.: (8142) 781218