УДК 519.6
Вестник СПбГУ. Прикладная математика... 2017. Т. 13. Вып. 1
В. В. Карелин , В. М. Буре, М. В. Свиркин
ОБОБЩЕННАЯ МОДЕЛЬ РАСПРОСТРАНЕНИЯ ИНФОРМАЦИИ В НЕПРЕРЫВНОМ ВРЕМЕНИ*
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
В данной статье предлагается модель, которая является обобщением классической модели Маки—Томпсона и стохастической модели Дэйли—Кендалла. Вводится новый параметр для оценки вероятности распространения слухов. Обобщенная модель формулируется и рассматривается в непрерывном времени. Процесс распространения слухов описывается системой линейных дифференциальных уравнений. Строится общее решение для динамики распространения слухов. Приводятся численные примеры. Библиогр. 13 назв. Ил. 2. Табл. 2.
Ключевые слова: информационные сети, дифференциальные уравнения, вероятность, распространение слухов.
V. V. Karelin, V. M. Bure, M. V. Svirkin
GENERALIZED MODEL OF INFORMATION SPREADING IN CONTINUOUS TIME
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
The diffusion of information or spreading of rumours is a social phenomenon and plays a significant role in a daily life. Rumours play a very important role in social life and have existed as a social fact since ancient times. The rumour model explains the spread of rumours and serves as a tool for understanding this social phenomenon. Rumours in economics have become more intensively discussed and investigated in recent decades. There are examples of rumour dynamics based on communication and exchange at auctions, in the stock markets and during trading. These backgrounds and motivations give the basis for a mathematical model of the diffusion of information or the spreading of rumours. We give a model that is a generalization of the classical Maki—Thompson model and the stochastic Daley—Kendall model of rumour spreading using the probability approach. A new parameter is suggested for the probability of rumour spreading. The generalized model is formulated and is considered in continuous time. The process of spreading rumours is described by a system of linear differential equations. The general solution for dynamics of spreading of rumours is constructed. Refs 13. Figs 2. Tables 2. Keywords: information networks, differential equations, probability, spreading rumors.
Карелин Владимир Витальевич — кандидат физико-математических наук, доцент; [email protected]
Буре Владимир Мансурович — доктор технических наук, профессор; [email protected] Свиркин Михаил Владимирович — кандидат физико-математических наук, доцент; smv01_01@ mail.ru
Karelin Vladimir Vitalievich — PhD of physical and mathematical sciences, associate professor; [email protected]
Bure Vladimir Mansurovich — doctor of technical sciences, professor; [email protected] Svirkin Mihail Vladimirovich — PhD of physical and mathematical sciences, associate professor; [email protected]
* Работа выполнена при финансовой поддержке Санкт-Петербургского государственного университета (грант № 9.38.205.2014).
© Санкт-Петербургский государственный университет, 2017
Введение. Интерес к проблематике, связанной с распространением информации в социальных сетях, с каждым годом возрастает, и связано это в первую очередь с тем, что они вовлекают в свою сферу огромное количество пользователей. В фундаментальной работе [1] исследуются математические модели социальных сетей и модели информационного влияния внутри этих сетей. Участники социальной сети имеют определенную репутацию, обмениваются информацией, объединяются в группы, решают задачи, связанные с достижением различных целей. В рамках указанной проблематики естественным образом возникают задачи исследования информационных конфликтов, в отдельных случаях речь может идти даже об "информационных войнах". По существу, близкий круг вопросов рассматривался ранее в проблематике, которая получила название "model of rumor spreading" — модели распространения слухов, например в [2-7]. Ученые из Гарвардского университета — Гордон У. Оллпорт и Лео Постман в 1947 г. в книге «Психология слухов» одними из первых разработали математические модели распространения слухов. По их определению, слух — это информационное сообщение, которое распространяется между людьми, как правило, в устной форме, без предоставления доказательств его достоверности. Главный признак слуха заключается в том, что он получает широкое распространение при отсутствии доказательств его достоверности. В этом состоит основное отличие слухов от новостей и фактов. Было отмечено, что люди порождают слухи в попытке осмыслить неоднозначные, неопределенные или запутанные ситуации [2-5]. Слухи играют очень важную роль в социальной жизни общества и существуют как социальный феномен с древних времен. Слухи в экономике стали более интенсивно обсуждаться и исследоваться в последние десятилетия [2]. В ряде экономических сфер слухи играют ключевую роль, например, динамика цен на аукционах, на фондовых рынках во многом определяется ими. Перечисленные причины дали основание для создания математических моделей распространения слухов [2, 6, 7]. Модель Маки—Томпсона [6] и стохастическая модель Дэйли—Кендалла [7] представляют собой классические модели распространения слухов. Модель Маки—Томпсона была рассмотрена в непрерывном времени и новое общее решение было получено для динамики распространения слухов [2]. В настоящей статье предложено обобщение модели Маки—Томпсона и стохастической модели Дэйли—Кендалла на основе использования вероятностного подхода. Она продолжает исследование нетрадиционных задач теории управления, опубликованных ранее в работах [8-11]. Проблематика исследования близка к проблемам, изучаемым в теории сетевых игр [12].
Основные результаты и обобщенная модель распространения слухов в непрерывном времени. Следуя подходу и обозначениям из [2], будем рассматривать три подгруппы населения: "ignorants" — лица, которые прежде не были проинформированы о слухе, "spreaders" — лица, активно распространяющие слух среди населения, "stiflers" — лица, которые были проинформированы о слухе, но в силу каких-то причин отказываются его распространять. В табл. 1 приведены схема распространения слуха и динамика изменения численности подгрупп населения при информационном общении их представителей: I — подгруппа "ignorants", S — подгруппа "spreaders", R — подгруппа "stiflers".
Обобщенная модель распространения слуха определяется по схеме вероятностей перехода. В табл. 2 присутствует ключевой для обобщенной схемы распространения слухов параметр: p — вероятность того, что в данном взаимодействии реализуется модель Маки—Томпсона, (1 — p) — вероятность того, что в конкретном взаимодействии
Таблица 1. Взаимодействие представителей подгрупп в процессе распространения слуха
S -> I Модель Дэйли—Кендалла Модель Маки—Томпсона S -> R
ik = ifc-1 - 1 ik = ik-1 ik = ik-1 ik = ik-1
Sk = Sfc-l + 1 Sk = Sk-1 - 2 Sk = Sk-1 - 1 Sk = Sk-1 ~ 1
Гк = Г к-1 rk = rk-1+2 rk = rk-i + 1 Гк = rk-l + 1
Aifc = ik ~ Ч-1 = -1 A ik = 0 A ik = 0 A ik = 0
Asfc = sk ~ Sfc-1 = 1 A sk = -2 A s к = -1 A sk = -1
Art = 0 A rk = 2 A rk = 1 Art = 1
реализуется схема Дэйли—Кендалла. Именно наличие введенного параметра p является главным отличием от работы [2], в которой был предложен подход к перенесению процесса распространения слухов на случай непрерывного времени, там же были получены дифференциальные уравнения описанных ниже типов, но рассмотрение ограничивалось только моделью Маки—Томпсона. Введение параметра p и усложнение схемы взаимодействия позволяют рассматривать более общую модель, включающую модели Маки—Томпсона и Дэйли—Кендалла как частные случаи: при p = 1 получаем модель Маки—Томпсона для непрерывного времени, а при p = 0 — соответственно модель Дэйли—Кендалла. Следует отметить, что при информационном общении двух представителей подгруппы "ignorants" изменение численности происходить не может. Схема информационного взаимодействия, содержащаяся в табл. 2, является непосредственным обобщением классических схем из моделей Маки—Томпсона и Дэйли—Кендалла для непрерывного времени.
Таблица 2. Схема вероятностей перехода слухов при информационном общении населения
Взаимодействия Переходы Переходы вероятности
S-> I S -> S S -> R (i - 1, s + l, r) {(¿, s _ 1, r + l),p } V {(¿. s - 2, r + 2), (1 - p)^} Модель Маки—Томпсона Модель Дэйли—Кендалла {г, s- 1,'Г+ 1) i/п (S - 1)/п г /п
Общая численность населения представляет собой сумму численностей подгрупп
i + s + r = n +1 = N. (1)
Введем начальные условия
i(0) = a • n, s(0) = в • n, r(0) = y • n, a, в, Y € (0,1), a + в + Y = 1, в • n > 1.
Условное ожидаемое изменение численности подгруппы "ignorants", в соответствии со схемой из табл. 1, оценивается как
л • / 14 ik-1 Ai^ N ifc-1
Ач = (-1)-Д t =>-— = (-1)-.
n At n
Переходя к пределам при At ^ 0, получим дифференциальное уравнение
^ , ,1 -Л1\ ^
— = (-1)— г * =Се
dt n
С = i(0) = а ■ n, i(t) = а ■ n ■ е~~.
Условное ожидаемое изменение численности подгруппы "spreaders" с учетом соотношения (1) для p € [0; 1] равно
Ask ik-i n - ik-i sk-i - 1 .
--(1 — p) --, At —> 0.
At
Тогда получим дифференциальное уравнение
ds i(t) 1 — p . . 1 — p — n -jT = 2—--- • s i +---,
dt n n n
ds 1 — p , ч t 1 — p — n
— +---s(t)= 2 -£-. (2)
dt n n
Условное ожидаемое изменение численности подгруппы "stiflers" с учетом соотношения (1) для p € [0; 1] оценивается как
Агк = (1 _p)sJ^±At + ~ д^о.
nn Следовательно, получим дифференциальное уравнение
dr J(t) , , s(t) - 1
"77= -1 — + 1-Р • —- +1.
dt n n
Таким образом, имеем систему уравнений
dj_ _ ) i(t)
dt ^ п '
& = 2 ^ •«(*) +(3) § = + +1-
Будем решать уравнение (2) методом Бернулли [13], сделаем замену переменной:
s = u • v, s = U • v + u • v.
У нас есть система уравнений
1p
v -\--- • v = 0,
n
о -i- . 1 -р-п
■и ■ v = 2 ■ а ■ е "И--.
n
Из первого уравнения находим функцию v
dV 1 - Р f dv (1 - P) f ,, , | | 1 - P , / (1 - P) • ^
v —>■ / — =--/ ai —In\v\ =---t^i/ = exp(--),
dt n v n n n
из второго получаем
du _p± 1 — p — n (i-p)-t
—— = 2 • a ■ e " H--- e
dt n
Отсюда
(1 — p — n) (i-p)-t 2 • a ■ n _pt_
и = —:--— . e "--- e " + C-
(1 - p ) p
Общее решение для s(t) имеет вид при p € (0,1)
^ (1 — р — п) 2 • а ■ п _± ^ iр (1 — р — п) ^ 2 • а ■ ?? Л (i-p)-t
(1-р) р V Р J
или
при р = 0 при р = 1
в(г) = (1 - п) - [2 • а • г - 1 - п(в - 1)] • е"
= ((3 + 2 • а) • п - 2 • се • п • е"- -
Общее решение для г (г) имеет вид при р Е (0,1)
г(г) = 7 • п + п •
(1 - р - п) 2 • а • (1 - р)
п
+ п •
0_(1-р-п) + 2_а (1 - р) • п р
п
(1 - р - п) 2 • а • (1 - р)
а
е п — п •
(1 —р - п) ^ 2 • а^ ^ (1 - р) • п р
(1 -У)-Ь 1 — р — п х е ----
при р = 0
г(г) = 2 • а • п • е
1 +
1 - п(в - 1)
а
п • е п —
1 - п(в - 1)
при р = 1
• п - 2 • а • п + 7 • п,
г(£) = £ + се • п • е п — се • п + 7 • п.
Рассмотрим несколько вариантов поведения 5 (г) в зависимости от р Е (0,1) (рис. 1). Точка Т, в которой в(Т) = 0, представляет собой момент времени, когда
а
р
п
р
г
п
п
а
п
заканчивается процесс распространения слуха. Этот процесс можно рассматривать на интервале времени (0,Т).
Динамика изменения численности подгруппы "ignorants" показана на рис. 2.
Заключение. В работе сформулирована и исследована модель распространения слухов в непрерывном времени, обобщающая классические модели распространения слухов Маки—Томпсона и Дэйли—Кендалла.
Литература
1. Губанов Д. А., Новиков Д. А., Чхартишвили А. Г. Социальные сети: модели информационного влияния, управления и противоборства. М.: Физматлит, 2010. 228 с.
2. Belen S., Kropat E., Weber G.-W. On the classical Maki—Thompson rumour model in continuous time // Central European Journal of Operations Researches. 2011. Vol. 19. P. 1-17.
3. Rosnow R. L. Psychology of rumor reconsidered // Psychol. Bull. 1980. Vol. 87, N 3. P. 578-591.
4. Rosnow R. L. Inside rumor: a personal journey // Amer. Psychol. 1991. Vol. 46, N 5. P. 484-496.
5. Difonzo N., Bordia P. Rumor and prediction: making sense (but losing dollars) in the stock market // Organ. Behav. Hum. Decis. Process. 1997. Vol. 71, N 3. P. 329-353.
6. Maki D. P., Thompson M. Mathematical models and applications. Prentice-Hall: Englewood Cliffs, 1973. 274 p.
7. Daley D. J., Kendall D. G. Stochastic rumours // J. Inst. Math. Appl. 1965. Vol. 1. P. 42-55.
8. Karelin V. V., Bure V. M., Polyakova L. N., Elfimov A. N. Control problem of a deterministic queuing system // Applied Mathematical Sciences. 2016. Vol. 10, N 21-24. P. 1023-1030.
9. Bure V. M., Karelin V. V., Polyakova L. N. Probabilistic model of terminal services // Applied Mathematical Sciences. 2016. Vol. 10, N 39. P. 1945-1952.
10. Bure V. M., Karelin V. V., Polyakova L. N. The problem of resource allocation between the protection system and constructing redundant components // Applied Mathematical Sciences. 2015. Vol. 9, N 93-96. P. 4771-4779.
11. Yakushev V. P., Karelin V. V., Bure V. M., Parilina E. M. Soil acidity adaptive control problem // Stochastic enviromental research and risk assessment. 2015. Vol. 29, N 6. P. 1671—1677 (see also URL: http://dx.doi.org/10.1007/s00477-014-0965-5, дата обращения: 15.10.2016).
12. Bure V., Parilina E., Sedakov A. Consensus in social networks with heterogeneous agents and two centers of influence // Intern. conference "Stability and control processes" in memory of V. I. Zubov (SCP) joined with 21st Intern. workshop on beam dynamics and optimization (BDO). Saint Petersburg, 05-09 Oct. 2015. Saint Petersburg, 2015. P. 233-235.
13. Kamke E. Differentialgleichungen: Losungsmethoden und Losungen, I, Gewohnliche Differentialgleichungen. Leipzig: Teubner, 1959. 668 p.
Для цитирования: Карелин В. В., Буре В. М., Свиркин М. В. Обобщенная модель распространения информации в непрерывном времени // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2017. Т. 13. Вып. 1. С. 74-80. DOI: 10.21638/11701/spbu10.2017.107
References
1. Gubanov D. A. , Novikov D. A., Chkhartishvili A. G. Socialnye seti: modeli informacionnogo vliyaniya upravleniya i protivoborstva [Social networks: models of informational influence, control and confrontation]. Moscow, Fizmatlit Publ., 2010, 228 p. (In Russian)
2. Belen S., Kropat E., Weber G.-W. On the classical Maki—Thompson rumour model in continuous time. Central European Journal of Operations Researches, 2011, vol. 19, pp. 1-17.
3. Rosnow R. L. Psychology of rumor reconsidered. Psychol. Bull., 1980, vol. 87, no. 3, pp. 578-591.
4. Rosnow R. L. Inside rumor: a personal journey. Amer. Psychol., 1991, vol. 46, no. 5, pp. 484-496.
5. Difonzo N., Bordia P. Rumor and prediction: making sense (but losing dollars) in the stock market. Organ. Behav. Hum. Decis. Process, 1997, vol. 71, no. 3, pp. 329-353.
6. Maki D. P., Thompson M. Mathematical models and applications. Prentice-Hall, Englewood Cliffs, 1973, 274 p.
7. Daley D. J., Kendall D. G. Stochastic rumours. J. Inst. Math. Appl., 1965, vol. 1, pp. 42-55.
8. Karelin V. V., Bure V. M., Polyakova L. N., Elfimov A. N. Control problem of a deterministic queuing system. Applied Mathematical Sciences, 2016, vol. 10, no. 21-24, pp. 1023-1030.
9. Bure V. M., Karelin V. V., Polyakova L. N. Probabilistic model of terminal services. Applied Mathematical Sciences, 2016, vol. 10, no. 39, pp. 1945-1952.
10. Bure V. M., Karelin V. V., Polyakova L. N. The problem of resource allocation between the protection system and constructing redundant components. Applied Mathematical Sciences, 2015, vol. 9, no. 93-96, pp. 4771-4779.
11. Yakushev V. P., Karelin V. V., Bure V. M., Parilina E. M. Soil acidity adaptive control problem. Stochastic enviromental research and risk assessment, 2015, vol. 29, no. 6, pp. 1671-1677. Available at: http://dx.doi.org/10.1007/s00477-014-0965-5 (accessed: 15.10.2016).
12. Bure V., Parilina E., Sedakov A. Consensus in social networks with heterogeneous agents and two centers of influence. Intern. conference "Stability and control processes" in memory of V. I. Zubov (SCP) joined with 21st intern. workshop on beam dynamics and optimization (BDO), Saint Petersburg (2015/05-09/10). Saint Petersburg, 2015, pp. 233-235. (In Germany)
13. Kamke E. Differentialgleichungen: Losungsmethoden und Losungen, I, Gewohnliche Differentialgleichungen. Leipzig, Teubner Publ., 1959, 668 p.
For citation: Karelin V. V., Bure V. M., Svirkin M. V. Generalized model of information spreading in continuous time. Vestnik of Saint Petersburg University. Applied mathematics. Computer science. Control processes, 2017, volume 13, issue 1, pp. 74-80. DOI: 10.21638/11701/spbu10.2017.107
Статья поступила в редакцию 1 ноября 2016 г. Статья принята к печати 19 января 2017 г.