1. Введение 4
2. Модель случайного графа с фиксированным распределением степеней вершин 7
2.1. Производящие функции 7
2.2. Величина компонент связности 11
2.3. Ориентированные графы 13
3. Модель кросс-холдинга Гая-Кападия 18
3.1. Общая структура 18
3.2. Производящие функции 21
3.3. Фазовый переход 23
4. Алгоритм моделирования случайного графа и запуск дефолтов 26
4.1. Общая схема 26
4.2. Моделирование дискретной случайной величины 27
4.3. Преобразование набора вершин в граф 29
4.4. Моделирование кросс-холдинга 30
4.5. Моделирование дефолтов на полученном графе 31
4.6. Программа 31
5. Результаты 33
Список литературы 41
Существуют системы, в которых элементы тесно и сложно связаны между собой. Примером такой системы может являться замкнутая группа людей, в которой один участник заразился некоторой инфекцией. Тогда распространение инфекции в группе зависит от частоты контактов заболевшего с остальными членами группы, силы иммунитета и прочих факторов. Аналогичным примером может служить сеть связанных компьютеров, в которой распространяется вирус, или социальные сети и распространение через них различного вирусного контента. Подобные эффекты можно увидеть и в финансовой сфере. Объектом исследования данной работы являются кросс-холдинги. Кросс-холдинг — это система, состоящая из различных организаций (например, банков), в которой участники обладают друг перед другом рядом обязательств (взаимные кредиты в случае банков). Аналогом заболевания в таком случае будет служить дефолт банка, т.е. состояние неплатежеспособности и невозможность отдавать долги. Таким образом, дефолты могут распространяться от банка к банку, а неплатежеспособность единственного банка может привести к краху всей системы. Подробно системы кросс-холдингов были описаны в [1].
В качестве инструмента, который используется для исследования кросс-холдингов, были выбраны случайные графы с фиксированным распределением степеней вершин. Дело в том, что графы сами по себе удобны для представления и изучения различных систем, представимых в виде сетей, а теория случайных графов позволяет изучать некоторые типичные свойства подобных систем. Конкретная модель случайного графа — с фиксированным распределением степеней вершин — также была выбрана неслучайно. Несмотря на то, что, пожалуй, самой известной моделью является модель случайного графа Эрдёш-Реньи (описана в [2]), она мало подходит для исследования кросс-холлдингов из-за Пуассоновского распределения степеней вершин, которое нехарактерно для изучаемых моделей. Поэтому, как утверждается в [4], для подобных систем следует использовать схему, которая позволяет строить граф, основываясь на желаемом распределении степеней его вершин. Для этих целей была выбрана модель, предложенная Ньюменом, Строгатцем и Ватсом (мы будем ссылаться на неё по имени первого автора).
Наложение финансовой составляющей на модель случайного графа Ньюмена приводит нас к модели кросс-холдингов, которая описана в статье Гая-Кападии [3]. В ней аппарат производящих функций объединяется с понятиям уязвимости банка, и проводятся исследования соответствующих свойств полученных систем.
С практической точки зрения интересен алгоритм программного моделирования случайного графа с заданным законом распределения. Для его реализации требуется решить несколько проблем: моделирование случайных величин, процесс соединения вершин графа, отсеивание плохих случаев. Общая схема была взята из статьи [4]. Задача моделирования дискретной случайной величины была решена классическим методом Уолкера (описан в [5]).
Работа состоит из 5 параграфов. Во втором описываются различные модели случайных графов и аппарат производящих функций, который позволяет их исследовать. В третьем параграфе говорится о модели кросс-холдинга и о применении в ней теории случайных графов. Некоторые теоретические результаты, полученные в этих двух параграфах, в дальнейшем сравниваются с данными, полученными в результате моделирования. В четвертом параграфе идёт речь о непосредственно моделировании случайных графов и экономических процессов на них. В заключительном, пятом, параграфе анализируются результаты проведенного моделирования.
Ниже представлены графики, иллюстрирующие некоторые результаты моделирования. Для значения каждого параметра была взята выборка из 50 графов. В начале рассмотрим влияние изменения параметров распределения непосредственно на структуру графа, а именно на число ребер:
График 1
График 1 иллюстрирует явную линейную зависимость между числом вершин и числом ребер в графе при заданном в (46) распределении.
Далее представлены графики зависимости числа ребер от смещения параметров Tin, Tout и максимальной полустепени в графе:
График 2
График 3
График 4
Далее рассмотрим влияние различных параметров распределения на количество банков, которые стали неплатежеспособными после дефолта случайно выбранного в графе банка. Это количество выражается в процентном соотношении относительно размера всей системы:
График 5
Как видно из графика 5, процент банков, объявивших дефолт при таких параметрах распределения, примерно одинаков для всех значений N и находится в районе 42 процентов. Причем, если увеличить минимальное значение полустепени с 4 до 5, то при тех же параметрах моделирования дефолт исходного банка никак не влияет на систему.
График 6
График 6 иллюстрирует тот факт, что с увеличением Tin количество дефолтов становится чрезвычайно большим, что позволяет сделать вывод о желательном понижении этого параметра.
График 7
Небольшие ограничения максимальной степени (подразумевается как полустепень захода, так и полустепень выхода) графа не приводят к существенным изменениям в проценте неплатежеспосбных банков:
График 8
График 9
График 9 показывает, что независимо от количества вершин, мат. ожидание и среднее значение полустепени захода вершины при данном распределении почти совпадают.
Однако значения мат. ожидания и среднего значения полустепени выхода, напротив, стабильно различаются на достаточно большую величину:
График 10
В следующих графиках размер выборки был увеличен до 100. Можно посчитать значение Vj, которое было выведено в (32), как отношение числа банков системы с полустепенью захода j, для которых выполнено (31), к общему числу банков с полустепенью захода j. В графиках представлены зависимости среднего числа различных полустепеней захода j, для которых Vj = 1. Причем при q =1 это число стабильно равно 2, а значения j при этом равны 4 и 5, т. е. это две минимально возможные полустени захода. Аналогичные эффекты наблюдаются и для остальных значений q. Эти данные лишний раз иллюстрируют тот факт, что надежность системы укрепляется вместе с ростом полустепеней банков.
q = 0.95:
График 11
q = 0.9:
График 12
q = 0.5:
График 13
Кроме количества банков, которые объявляют дефолт, интересно посмотреть и среднюю глубину волны дефолтов (под глубиной подразумевается наибольшее расстояние от первого банка, объявившего дефолт, до некоторого другого банка, до которого дошла волна дефолтов):
График 14
Максимальные по выборке значения глубины дефолта:
График 15.
[1] M. Elliott, B. Golub, M.O. Jackson. Financial Networks and Contagion. American economic review, vol. 104, no. 10, 2014, pp. 3115-53.
[2] Райгородский А. М. Модели случайных графов. М .: МЦНМО, 2011.
[3] S. Kapadia, P. Gai. Contagion in financial networks. Proceedings of Royal Society Series A, 466(2120),2010, pp. 2401-2423.
[4] M. E. J. Newman, S. H. Strogatz, and D. J. Watts. Random graphs with arbitrary degree distributions and their applications. Physical review E 64 (2), 026118, 2001, arXiv:cond-mat/0007235v2.
[5] И. В. Романовский. Моделирование случайных распределений. Конспект лекций по курсу "Дискретный анализ". http://www.math.spbu.ru/ru/Archive/Courses/jvr/DA_html/ _lec_1_07.html
[6] A.V. Leonidov and E.L. Rumyantsev. "Russian interbank networks: main characteristics and stability with respect to contagion", Proc. "Instabilities and Control of Excitable Networks: from macro- to nano- systems MIPT, 2012, arXiv: 1210.3814.