Изучение нейронных сетей является активно развивающейся в настоящее время областью прикладной математики. Искусственная нейронная сеть - математическая модель, построенная по принципу организации и функционирования сетей нервных клеток живых организмов. В качестве нейронной сети можно рассматривать любую модель, содержащую узлы и связи между ними. Нейронные сети находят применение во многих областях: социологии, экономике, биохимии, технических науках.
Топология связей в сетях обычно бывает либо полностью регулярной, либо полностью произвольной. Но, поскольку многие реальные сети находятся между двумя этими крайностями, была выделена группа сетей, обладающая высокой степенью кластеризации, как сети с решетчатой топологией, и имеющая малую среднюю длину пути, как сети со случайной (произвольной) топологией. Такие сети носят название «small-world» и были описаны Дунканом Ваттсом и Стивеном Строгацем в 1998 году [15]. «Small-world» определяется как сеть, в которой расстояние Lмежду двумя произвольно выбранными вершинами растет пропорционально логарифму от числа вершин Nв сети, то есть L х log2N. За этой работой по следовало множество других исследований проявлений свойств «small world» [1, 5, 9, 10, 13]. Сети «small world» описывают различные объекты, в том числе, различные отделы головного мозга [7, 11], а в сфере информационных технологий они проявляют себя в файлообменных сервисах [8], пиринговых сетях.
Такие сети имеют много вершин с редкими связями, но не настолько редкими, чтобы граф мог оказаться несвязным, поэтому принимаем N >>к >>lnN >> 1, где к >>ln N гарантирует связность графа
[3] . Также «small-world» сети имеют тенденцию содержать в себе клики (термин пришел из работы [6]), это следует из определения такого свойства сети, как коэффициент кластеризации - степень, с которой вершины склонны группироваться вместе.
Алгоритм построения сетей «small-world» (описанный далее в разделе 2.1) позволяет «настроить» граф на промежуток между регулярностью и произвольностью и исследовать область 0 <р < 1. В работе [15] показано, что есть большой интервал р, на котором длина пути L(p)почти также мала как в произвольном графе Lrandomпри коэффициенте кластеризации много большем, чем в произвольном графе, С (р) >>Crandom (см. рис. 1). Данные собирались для графов с числом вершин N = 1000 и соседей к = 10 и были нормированы с помощью значений L(0) и С(0) для регулярного графа. Логарифмическая горизонтальная шкала использована, чтобы показать быстрое падение длины пути L(p),что соответствует «small-world». При этом коэффициент кластеризации С(р) почти не изменяет своего значения, соответствующего регулярному графу. Это значит, что переход сети к виду «small-world» почти невозможно заметить на локальном уровне.
Согласно графикам, сеть типа «small-world» получается при значении вероятности от 0.01 до 0.1. То есть, чтобы отнести моделируемую сеть к типу «small-world», нужно проанализировать полученные значения длины кратчайшего пути в среднем и коэффициента кластеризации. Но для исследования областей устойчивости нейронных сетей будет интересна любая сеть, моделируемая со значением вероятности из области 0 <р < 1, упомянутой в оригинальной работе [15].
Одним из предметов исследования в работе [17] является важное свойство - устойчивость нейронной сети. Устойчивость - способность продолжать стабильно функционировать при возникновении колебаний при передаче сигналов между нейронами, например, вследствие запаздывания. Достаточно разработана тема устойчивости непрерывных моделей нейронных сетей стандартных конфигураций (кольцевой и линейной) [21]. Особенностью данной работы является то, что модель нейронной сети дискретна, то есть для изучения проблемы ее устойчивости необходимо применение теории разностных матричных уравнений с запаздываниями [16].
Эта проблема разработана гораздо меньше, чем аналогичная для непрерывных моделей. В работе [17] был разработан метод анализа устойчивости дискретных моделей нейронных сетей - метод конусов устойчивости, для которого был реализован программный продукт [18] в программной среде MATLAB.
Цель и задачи исследования
Целью данной работы является разработка и реализация параллельных алгоритмов для построения моделей нейронных сетей Ваттса-Строгаца, расчета их параметров, генерации входных данных для программного продукта [18], осуществляющего анализ устойчивости нейронных сетей.
Для достижения поставленной цели нужно выполнить следующие задачи:
1) рассмотреть существующие варианты последовательного алгоритма по строения моделей нейронных сетей Ваттса-Строгаца и его реализации, а также алгоритмов для расчета его свойств;
2) изучить методы разработки параллельных алгоритмов для сопроцессора Intel Xeon Phi;
3) на основе существующих последовательных алгоритмов спроектировать и реализовать параллельные алгоритмы;
4) провести эксперименты по оценке эффективности разработанных алгоритмов.
Структура и объем работы
Работа состоит из введения, 4 основных разделов, заключения и
библиографии. Объем работы составляет 34 страницы, объем библиографии - 21 наименование.
Содержание работы
В первом разделе проведен обзор существующих алгоритмов и их реализаций для построения модели нейронной сети Ваттса-Строгаца и расчета ее коэффициентов.
Во втором разделе приведены описания последовательных алгоритмов.
Третий раздел содержит описание параллельных алгоритмов на основе технологии OpenMP, а также приведено краткое описание архитектуры сопроцессора Intel Xeon Phi и методов разработки для этого сопроцесора.
В четвертом разделе представлены результаты экспериментов и анализ эффективности разработанных алгоритмов.
В заключении подведены итоги проделанной работы.
Данная работа посвящена разработке и реализации параллельного алгоритма для построения моделей нейронных сетей Ваттса-Строгаца, а также алгоритмов для расчета их параметров и генерации матриц, которые являются входными данными для продукта [18], производящего расчет областей устойчивости нейронных сетей. Были получены следующие основные результаты:
1) рассмотрены существующие варианты последовательного алгоритма по строения моделей нейронных сетей Ваттса-Строгаца и его реализации, а также алгоритмов для расчета его свойств;
2) изучены методы разработки параллельных алгоритмов для сопроцессора Intel Xeon Phi;
3) на основе существующих последовательных алгоритмов спроектированы и реализованы параллельные алгоритмы с помощью технологии OpenMP;
4) проведены эксперименты по оценке эффективности разработанных алгоритмов на процессоре и сопроцессоре.
Генерируемые программой данные пригодны не только для исследования дискретных моделей нейронных сетей, но и для аналогичных непрерывных моделей. Результаты экспериментов показывают, что реализованные алгоритмы пригодны для достаточно быстрой обработки нейронных сетей больших размерностей.
1. Barrat A., Weight M. On the properties of small-world network models. // The European Physical Journal, 2000. - No. 13. - P 547-560.
2. Boccaletti S. Complex Networks: Structure and Dynamics /
S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, D.U. Hwang // Physics Reports, 2006. - Vol. 424. - No. 4. - P 175-308.
3. Bollobas B. Random graphs // London Mathematical Society Monographs, Academic Press, London, 1985. - P 447.
4. Building a Native Application for Intel®Xeon PhiTMCoprocessors | Intel® Software. [Электронный ресурс] URL: https://software.intel.com/en- us/articles/building-a-native-application-for-intel-xeon-phi-coprocessors/(дата обращения: 25.03.2016).
5. Dorogovtsev S.N., Mendes J.F.F. Evolution of networks. // Advances in Physics, 2002. - No. 13. - P 1079-1187.
6. Duncan L.R., Perry A.D. A method of matrix analysis of group structure. // Psychometrika, 1949. - No. 14. - P 95-116.
7. Gray R.T., Fung C.K.C., Robinson P.A. Stability of small-world networks of neural populations. // Neurocomputing, 2009. - Vol. 72. - No. 7-9. - P 1565-1574.
8. lamnitchi A., Ripeanu M., Foster I. Small-world file-sharing communities. // INFOCOM 2004. Twenty-third AnnualJoint Conference of the IEEE Computer and Communications Societies, 7-11 March, 2004. - Vol. 2. - P. 952-963.
9. Ivanov S., Kipnis M. On the stability of a neural network with links based on the Watts-Strogatz model. // International Journal of Pure and Applied Mathematics, 2015. - Vol. 105. - No. 3. - P 431-438.
10. Kleinberg J. The small-world phenomenon: an algorithmic perspective // STOC ’00 Proceedings of the thirty-second annual ACM symposium on Theory of computing, New York, NY, USA, 2000. - P. 163-170.
11. Netoff T.I. Epilepsy in Small-World Networks / T.I. Netoff,
R. Clewley, S. Arno, T. Keck, J.A. White // The Journal of Neuroscience.
2004. - Vol. 24. - No. 37. - P. 8075-8083.
12. OpenMP.org. [Электронный ресурс] URL: http://openmp.org/wp/(дата обращения: 21.02.2016).
13. Sinha S. Stability of small-world networks of neural populations // Physica A. 2005. - Vol. 346. - P. 147-153.
14. SNAP: Stanford Network Analysis Project. [Электронный ресурс] URL: http://snap.stanford.edu/index.html(дата обращения: 10.05.2015).
15. Watts D.J., Strogatz S.H. Collective dynamics of ’small-world’ networks. // Nature, 1998. - No. 393. - P. 440-442.
16. Беллман Р., Кук К.Л. Дифференциально-разностные уравнения. -М.: Мир, 1967.-548 с.
17. Иванов С.А. Устойчивость дискретных моделей стандартных конфигураций нейронных сетей с запаздывающими взаимодействиями: дис. ... канд. физ.-мат. наук: 05.13.18 / Иванов Сергей Александрович. - Челябинск, 2014. - 135 с.
18. Иванов С.А. Программный комплекс «Устойчивость матричных разностных уравнений с двумя запаздываниями»: свидетельство о регистрации электронного ресурса в ИНИМ РАО № 19417 / С.А. Иванов // ОФЭРНиО, 29.07.2013.
19. Седжвик Р. Фундаментальные алгоритмы на С++. Алгоритмы на графах: Пер. с англ./Роберт Седжвик. - СПб.: ООО «ДиаСофтЮП», 2002. - 496 с.
20. Технические характеристики вычислительного кластера «Торнадо ЮУрГУ». [Электронный ресурс] URL:
http://supercomputer.susu.ru/computers/tornado/(дата обращения: 25.04.2016).
21. Хохлова Т.Н. Устойчивость моделей нейронных сетей кольцевой и линейной конфигураций с запаздывающими взаимодействиями: дис. ... канд. физ.-мат. наук: 05.13.18 / Хохлова Татьяна Наилевна. - Челябинск, 2013. - 132 с.