Full Text

Введение Одним из наиболее перспективных подходов к изучению поведения финансовых рынков является фрактальный подход. Он основан на гипотезе фрактального рынка Б. Мандельброта [1]. Эта гипотеза гласит, что «процесс ценообразования на рынках глобально детерминирован, зависим от “начальных условий” [2], прошлые цены влияют на будущие, а графики ценовых рядов финансовых инструментов обладают фрактальными свойствами: самоподобием, масштабированием по степенному закону и фрактальной размерностью» [2, 3]. По мнению ученых - приверженцев данной теории, исследование финансовых рядов методами фрактального анализа позволяет с достаточно высокой точностью судить об их будущем поведении: предсказывать критические точки (кризисы) и строить точечные и/или интервальные прогнозы. Изначально фрактальный подход, включающий в себя методы физики, использовался при исследовании физических объектов и пространственных структур. Затем его начали применять при исследовании экономических процессов, где в качестве пространственных структур стали рассматривать графики временных рядов [4-7]. Сегодня фрактальный подход является одним из методов особого направления в экономической науке - эконофизике, возникшей на рубеже XX-XXI вв. [8]. Строго говоря, эконофизика - это наука, которая применяет методы физики для изучения экономических систем. В настоящий момент продолжаются исследования по эконофизике на кафедре информационных систем и математических методов в экономике экономического факультета Пермского государственного научного исследовательского университета (ПГНИУ). На кафедре ПГНИУ с помощью методов эконофизики изучаются модели финансовых пузырей, модель взаимодействия агентов на рынке [9]. Также на кафедре в рамках эконофизики изучаются фрактальный подход, представленный в работах [3, 10, 11], и р-адический метод, представленный в работах [12, 13]. Далее дадим характеристику фрактальному анализу. Сущность фрактального анализа состоит «в обработке масштабированием … фрактальной структуры и описании распределения какой-либо структурной характеристики - или меры - при этом масштабировании» [14]. Для этого определяется фрактальная размерность структуры D - численный показатель, характеризующий то, как объект заполняет пространство, его сложность. Применительно к временным рядам говорят о фрактальной размерности их графиков как пространственных структур. В отличие от других структур, для временных рядов существует собственная мера фрактальности - показатель Херста H, характеризующий свойство ряда, названное персистентностью (способностью ряда сохранять свою тенденцию) [2, 3]. Персистентность тесно связана с наличием в ряду так называемой длинной памяти - зависимости между далеко отстоящими друг от друга во времени значениями исследуемой переменной. Более подробно связь самоподобия (другими словами, фрактальности) временного ряда, наличия у него длинной памяти и персистентности показана в работах Б. Мандельброта [15]. Если исследуемый временной ряд является реализацией гауссовского случайного процесса, то показатель H связан с фрактальной размерностью D простым линейным соотношением H = 2 - D [16]. Ввиду этого фрактальный анализ временных рядов заключается в определении их фрактальных свойств как с помощью оценки размерности, так и с помощью оценки показателя Херста. Оценки фрактальных показателей (фрактальной размерности и показателя Херста) исследуемого временного ряда или его фрагмента (в этом случае речь идет о локальной размерности - понятии, введенном Н.В. Старченко [4]) позволяют судить о поведении этого ряда в будущем. И, конечно, при вычислении этих оценок важно получать значения, близкие к истинным. В настоящее время существует большое количество методов и способов оценивания фрактальных показателей, например метод R/S-анализа, метод Фурье, метод максимумов модулей коэффициентов вейвлет-преобразования, методы минимального и клеточного покрытия и другие [17]. На фоне большого количества известных методов работ, посвященных их подробному описанию и сравнению друг с другом, немного. В данной статье предпринята попытка детальной проработки одного из самых популярных и в то же время простых, интуитивно понятных методов оценки фрактальной размерности - метода клеточного покрытия. Новизна работы заключается в том, что русскоязычные авторы используют этот метод, но подробно его не рассматривают, хотя его применение имеет определенные сложности. Например, в работе [16] сказано о том, что клеточный метод «реализуется при графическом выводе фрактала на экран PC как массива пикселей» (под PC автором подразумевается персональный компьютер). На наш взгляд, анализ временного ряда сложно свести к анализу пикселей его графика, так как пиксели - всего лишь техническая особенность визуализации ряда на экране компьютера. Характеристики массива пикселей зависят от характеристик экрана, в частности от его разрешения, в то время как временной ряд существует сам по себе и никаким образом не должен быть связан с аппаратной конфигурацией компьютера исследователя. Более оправданным выглядит подход к определению размерности, основанный на аналитических вычислениях. Таким образом, для достижения цели исследования дано описание метода клеточного покрытия; предложен механизм его адаптации к временным рядам; разработана программная реализация метода на языке Python; проведена апробация разработанной программы на сгенерированных временных рядах и выполнено сравнение с реализациями клеточного метода и R/S-анализа из библиотек языка R. 1. Теоретический обзор Рассмотрим основные результаты и выводы по методу клеточного покрытия тех ученых, которые занимались проблемой применения этого метода для фрактального анализа временных рядов. Среди отечественных авторов популярность клеточного метода как средства исследования не пространственных, а временных структур невелика. К числу наиболее значимых работ можно отнести диссертацию Н.В. Старченко, в которой им разработана модификация клеточного метода, получившая название метода минимального покрытия. Автором утверждается, что минимальное покрытие точнее клеточного и что использование минимальных покрытий вместо клеточных позволяет существенно уменьшить количество данных, необходимое для определения фрактальной размерности [4]. Предложенный им метод пользуется популярностью у других исследователей [18-20], ссылающихся на выводы Старченко о преимуществах его подхода, изложенные в диссертации и сопутствующих работах. Другим исследователем - Е.К. Кривоносовой - в работе [16] применяются методы как клеточного, так и минимального покрытия для анализа и прогнозирования котировок акций ПАО «Лукойл» и выручки предприятий Приволжского федерального округа. Автором сделан вывод о том, что определенная обоими методами фрактальная размерность может быть использована в качестве индикатора стабильности работы предприятий. Среди англоязычных авторов клеточному методу как подходу к фрактальному анализу временных рядов уделяется больше внимания. Так, в работе [7] выполнен обзор различных алгоритмов оценки фрактальной размерности как меры гладкости (шероховатости) временных рядов и пространственных данных, в том числе рассмотрены различные вариации клеточного метода (box-counting). Авторами также разработана библиотека fractaldim для языка программирования R, содержащая функции, реализующие рассмотренные в указанной работе методы. Подробный анализ подходов к оценке фрактальной размерности временных рядов и возникающих при этом сложностей выполнен I. Pilgrim и R. Taylor. В их публикации [21] озвучены проблемы, характерные для клеточного метода, которые, как правило, не затрагиваются в отечественных источниках. (Следует отметить, что в процессе программной реализации метода мы столкнулись с теми же трудностями, что и Pilgrim и Taylor в своей работе.) Авторами сделан вывод, что клеточный и другие методы фрактального анализа временных рядов должны применяться с осторожностью к данным, «ограниченным по длине или спектральному содержанию», так как подобные ограничения могут препятствовать выявлению их фрактальных свойств. В работе [22] N.C. Kenkel подробно рассмотрены требования, предъявляемые к величине конечной выборки данных, по которой оценивается размерность клеточным методом. Автором сделан вывод, что для надежной оценки показателя требуются очень большие размеры выборки: обычно на порядки больше, чем используемые в большинстве эмпирических исследований, однако отмечено, что, если в качестве данных выступает временной ряд, объем выборки определяется просто как длина этого ряда. Также автором предложена вероятностная модель для вычисления минимально допустимого размера клетки, в которой общее число покрывающих фигуру клеток является функцией от размера клетки и величины клеточной размерности геометрического объекта. Согласно этой модели, минимальный размер клетки выбирается таким образом, чтобы отклонение размерности от величины ее оценки составляло не более e = 0,001. E. Rosenberg в работе [23] предложенную Kenkel модель обобщает на случай произвольного значения e. Проведенный анализ литературы дает основания полагать, что в русскоязычных источниках клеточный метод освещен далеко не так подробно, как в зарубежных, при этом авторами он сводится к анализу растрового изображения графика временного ряда. В настоящей работе произведена попытка вспомнить этот пробел, логичным образом адаптировать метод к временным рядам, устранив необходимость формирования компьютерного изображения графика ряда, и изложить возникшие при этом проблемы. Но, прежде чем приступить к описанию разработанной методики, следует кратко изложить теоретические основы метода клеточного покрытия как способа определения фрактальной размерности. Термин «размерность» означает число измерений некоторого геометрического объекта. Топологической размерностью называют минимальное количество координат, необходимое для описания объекта [3]. Эта величина может принимать только целые значения, поэтому подходит лишь для описания «идеальных», гладких объектов (прямая, куб, сфера), однако в природе такие объекты почти не встречаются. Ф. Хаусдорф обратил внимание, что если геометрический объект покрывать шарами радиусом δ, то минимальное число шаров N(δ) с уменьшением их радиуса возрастает по степенному закону Показатель DH получил название размерности Хаусдорфа. Его формула имеет вид [4] (1) Размерность Хаусдорфа обобщила понятие топологической размерности. Показатель DH может принимать как целые, так и дробные значения, что оказалось весьма полезным свойством для описания объектов реального мира. DH показывает, насколько плотно и равномерно объект заполняет пространство, например при малой изрезанности кривой значение ее DH приближается к 1, при большой - к 2 (рис. 1). а б в Рис. 1. Геометрические фигуры различной размерности: а - прямая, б - кривая Пеано в - полностью заполненное пространство Позднее Б. Мандельбротом был введен термин «фрактальная размерность» для обозначения частного случая размерности Хаусдорфа, в котором она принимает дробные значения, превышающие значение топологической размерности. Объекты, обладающие такой размерностью, Мандельброт назвал фракталами [1]. Метод клеточного покрытия весьма прост и основан на идее Хаусдорфа о степенной зависимости между количеством фигур элементарной формы, полностью покрывающих геометрический объект, и размером фигуры. Чем меньше фигура, тем больше таких фигур необходимо для полного покрытия объекта. Как следует из названия метода, в качестве фигуры в двумерном пространстве используется клетка, т.е. прямоугольник. В трехмерном пространстве фигурой является «коробка» (box, отчего англоязычное название метода - box-counting method). Далее будет рассмотрена разновидность метода клеточного покрытия, используемая для анализа объектов двумерного пространства. а б Рис. 2. Графики фрактального объекта (кривой Коха), где серым цветом выделены клетки, содержащие хотя бы одну точку анализируемого объекта. Переходя от а к б, можно заметить, что при уменьшении размера клетки со стороной δ происходит увеличение числа клеток и точности (аккуратности) покрытия объекта Алгоритм клеточного метода состоит в следующем. Плоскость, содержащая анализируемое изображение объекта, покрывается сеткой, которая делит ее на квадратные клетки стороной δ. Подсчитывается количество клеток N(δ), содержащих в себе хотя бы одну точку изображения. Затем поступательно уменьшается δ и подсчитываются соответствующие значения N(δ) (рис. 2). После в двойном логарифмическом масштабе строится график функции N(δ), который аппроксимируется прямой с помощью метода наименьших квадратов. Тогда величина клеточной размерности D, являющаяся тангенсом угла наклона данной прямой относительно оси Ox, вычисляется как взятый с противоположным знаком коэффициент при независимой переменной уравнения линейной регрессии где b - свободный член. В работе [21] предлагается модификация данного подхода, названная авторами вариационным методом клеточного покрытия (variational box-counting method), суть которой состоит в том, что клетки каждого столбца располагаются не строго по сетке, а таким образом, чтобы минимизировать их количество, необходимое для покрытия фигуры (рис. 3). Рис. 3. График фрактального объекта (кривой Коха), где серым цветом выделены клетки, содержащие хотя бы одну точку анализируемого объекта. Количество клеток, по сравнению с рис. 2 б, меньше В результате модификации количество клеток, соответствующих определенной величине δ, уменьшается. Однако, по утверждению авторов, внесение данных изменений в метод не оказывает значимого влияния на вычисленное значение размерности. 2. Адаптация метода клеточного покрытия для его применения к временным рядам Как было отмечено во введении, обычно в публикациях предлагается использовать клеточный метод для анализа графика временного ряда как изображения. Нами предпринята попытка отойти от замены анализа временного ряда анализом компьютерного изображения его графика и применить аналитические расчеты. Опишем ниже механизм адаптации к временным рядам метода клеточного покрытия. Пусть дан дискретный временной ряд V, состоящий из элементов vt, где Очевидно, что для реализации клеточного метода от понятия «график ряда» избавиться не удастся, однако под графиком следует понимать не растровое изображение, а совокупность точек ряда, соединенных между собой. Поскольку временной ряд является численным дискретным описанием непрерывного динамического процесса, значение показателя которого зафиксировано в моменты времени t, видится логичным принять в качестве формы траектории движения показателя наиболее простую из возможных - прямую (линейный тренд), т.е. представить график ряда как некоторую кусочно-линейную функцию, непрерывную и определенную на интервале Не менее очевидно, что в подсчете должны участвовать клетки, покрывающие не только точки, но и отрезки между ними. Прежде всего, необходимо ввести абстрактную плоскость (холст), вмещающую в себя график ряда. Пусть W - ширина холста, H - высота (безразмерные величины). Левый нижний угол холста лежит в точке с координатами (0; 0), правый верхний угол имеет координаты (W; H). Плоскость следует покрыть сеткой, т.е. разбить на квадратные клетки со стороной δ. CW - ряд, содержащий границы клеток по горизонтали, состоит из элементов CH - аналогичный ряд границ клеток по вертикали, состоящий из элементов Правые границы рядов вычисляются таким образом, потому что длина и высота холста могут не быть кратными стороне клетки, в то же время клетки, покрывающие график, несомненно, должны быть учтены даже в том случае, если они выходят за границы холста. Следующим шагом (или, скорее, комплексом шагов) является отображение на холсте графика временного ряда. Каждое значение ряда vt является точкой с координатами Возникает необходимость соотнести координатную сетку данных с введенной ранее абстрактной координатной сеткой холста, разделенного на клетки. Для этого следует «подогнать» первую ко второй, т.е. отмасштабировать ее координаты таким образом, чтобы они соответствовали координатам холста. Известно, что «плоскость данных» ограничена диапазонами по горизонтали и по вертикали, где - соответственно минимальное и максимальное значения ряда, а холст с клетками ограничен диапазонами и Чтобы привести «шкалу данных» к «шкале холста», имеет смысл воспользоваться следующими формулами: (2) (3) где - приведенные к единицам измерения холста координаты точки, соответствующей значению ряда Теперь известны границы всех клеток (ячеек сетки), на которые разделен холст, и координаты всех точек ряда, по которым можно определить, входит ли та или иная точка в любую из клеток. Однако, как было сказано выше, недостаточно выявить клетки, покрывающие точки, необходимо также найти и учесть те клетки, внутри которых находятся отрезки (целиком или фрагментом), соединяющие соседние точки. С одной стороны, координатами концов этих отрезков являются координаты самих точек, поэтому клетки, покрывающие точки, не имеет смысла учитывать отдельно: они войдут в набор клеток, пересекаемых отрезками. С другой стороны, необходимо выявить клетки, содержащие отрезки, соединяющие соседние точки. Следует рассматривать не каждую очередную точку в отдельности, а пары соседних точек. Пусть A - точка с координатами соответствующая значению ряда B - точка с координатами соответствующая значению ряда Точки A и B, помещенные на сетку холста, попадают в клетки, которые являются углами прямоугольника P. Среди всех клеток, составляющих этот прямоугольник, следует выделить те, которые покрывают отрезок AB (рис. 4). Рис. 4. Прямоугольник P и составляющие его клетки. Серым цветом выделены клетки, покрывающие отрезок AB При этом нужно учесть, что точка A может находиться как в левой верхней клетке, так и в левой нижней, а B - соответственно в правой нижней или верхней, что зависит от того, какая из точек на холсте расположена выше (какое значение из vt, vt+1 больше). Таким образом, нужно перебрать все клетки прямоугольника и для каждой определить, содержит ли она фрагмент отрезка AB либо весь отрезок целиком. Определить границы клетки, содержащей точку с координатами можно по следующим формулам. Пусть (ряд, содержащий границы клеток по горизонтали) состоит из элементов где - общее число элементов ряда (ряд, содержащий границы клеток по вертикали) состоит из элементов где - общее число элементов ряда Тогда (4) где - номер элемента ряда ближайшего, не меньшего p; - последний элемент ряда (5) где - номер элемента ряда ближайшего, не меньшего - последний элемент ряда Номера элементов ряда являющихся горизонтальными границами клетки, покрывающей точку, можно рассчитать по формулам (6) (7) где - номера левой и правой границ соответственно. Аналогичным образом можно вычислить номера элементов ряда являющихся вертикальными границами клетки: (8) (9) где - номера нижней и верхней границ соответственно. Подставив в формулы (4)-(9) координаты точек A и B, можно определить номера границ клеток, покрывающих эти точки: - номера левой и правой, нижней и верхней границ клетки, покрывающей точку A, и соответственно - номера границ клетки, покрывающей точку B. Тогда границы прямоугольника P можно вычислить следующим образом: где - номера левой, правой, нижней и верхней границ соответственно. Теперь, ограничив ряд диапазоном а ряд диапазоном можно перебрать все клетки прямоугольника P, координатами углов которых будут являться те или иные элементы этих рядов. Количество клеток прямоугольника определяется как произведение числа элементов первого ряда, уменьшенного на единицу, и числа элементов второго ряда, также уменьшенного на единицу. Далее, для определения того, является ли клетка закрашенной, т.е. содержащей внутри себя фрагмент отрезка AB, следует выявить факт пересечения отрезка AB и хотя бы одной из сторон этой клетки. Стоит отметить, что, несмотря на возможность ситуации, когда отрезок AB находится целиком внутри клетки и не касается ни одной из ее сторон, такой подход к выявлению закрашенных клеток в конечном счете сработает всегда: даже если на одной из итераций, при определенной паре точек A, B, обе из них попадут внутрь клетки, в результате чего она не будет учтена как закрашенная, на следующей (или предыдущей) итерации данная клетка все равно окажется учтена, потому что хотя бы одна ее грань будет пересечена отрезком, соединяющим предыдущую точку и A либо B и следующую точку. В случае же если клетка единственная и покрывает весь холст целиком, ее левая и верхняя границы совпадают с первым значением ряда V и его максимальным значением а следовательно, как минимум два из возможных отрезков AB пересекутся с одной из сторон этой клетки и она окажется учтена. Чтобы определить, пересекаются ли отрезки, можно воспользоваться методом, основанным на косом произведении векторов [24]. Пусть - отрезок, соответствующий стороне клетки, его концы имеют координаты Известно, что концы отрезка AB имеют координаты Можно ввести следующие обозначения максимумов и минимумов координат точек: Тогда условием пересечения отрезков AB и CD будет являться выполнение следующих неравенств (в левой части неравенств 2 и 3 приведены косые произведения векторов): 1) 2) 3) В случае, если отрезок AB пересекается хотя бы с одной из сторон клетки, эта клетка считается закрашенной и учитывается в подсчете величины N(δ). Разумеется, одна и та же клетка может быть учтена лишь один раз. Результатом перебора всех значений ряда V и выполнения описанных шагов является пара значений: N(δ) - подсчитанное количество закрашенных клеток и δ - размер стороны клетки. Дальнейшая последовательность действий полностью соответствует алгоритму метода клеточного покрытия, применяемому к анализу пространственных объектов: многократно уменьшается величина δ, вычисляются соответствующие значения N(δ), затем логарифмируются полученные ряды и определяется величина фрактальной размерности D из уравнения линейной регрессии В том случае, если размах ряда V равен нулю, т.е. каждый элемент ряда а график представляет собой отрезок прямой, параллельной оси Ox, для вычисления отмасштабированной координаты любой точки ряда согласно формуле (3), необходимо разделить на ноль, что приведет к невозможности реализации приведенной методики. Это частный случай входных данных, который нами предлагается учитывать особо и для такого ряда принять D равной величине теоретической размерности линии - единице. Несмотря на то, что изложенная методика вычисления клеточной фрактальной размерности временного ряда вполне соответствует сути метода и может быть использована, при реализации ее на практике возникает ряд проблем, требующих разрешения. Помимо самого ряда V, входными данными алгоритма являются диапазон изменения размера клетки δ и ширина и высота холста W, H. Вопросы определения этих параметров в работах отечественных авторов обычно опускаются. В то же время данные величины влияют на получаемую оценку размерности D, поэтому их произвольный выбор может привести к искажению этой оценки. Как отмечает E. Rosenberg в работе [23], проблема определения диапазона значений δ широко обсуждается в работах иностранных авторов. Очевидно, что необходимо выбрать такой диапазон, при котором становится возможным исследовать тенденцию масштабирования. Согласно работам [22, 23], не следует использовать величину δ, при которой где N - общее количество точек изображения. Поскольку величины δ и имеют обратную зависимость, минимальный размер клетки достигается при выполнении равенства Другими словами, при реализации метода можно уменьшать δ до тех пор, пока не будет выполняться упомянутое равенство. Однако в случае временного ряда остается неясным, какой величине соответствует показатель N. Выше было введено понятие холста абстрактного размера на который «помещался» график ряда V. Визуализация метода клеточного покрытия на компьютере помогает придать величинам W и H вполне конкретный смысл: их можно рассматривать как ширину и высоту окна с графиком и клетками в пикселях, в этом случае Но точку, из множества которых состоит холст, нельзя приравнивать к пикселю изображения на экране монитора, она не является шагом дискретизации ряда V ни по горизонтали, ни по вертикали. Во-первых, это следует из того факта, что количество значений ряда T может превышать величину W, а значит, «расстояние» между соседними наблюдениями vt и будет меньше единицы. Во-вторых, вопрос о дискретизации размаха значений ряда V, который будет рассмотрен ниже, и вовсе является трудноразрешимым. Исходя из этого, на наш взгляд, предложенный способ определения применим для анализа пространственных объектов, но не временных рядов. Авторы работы [7] предлагают определять δ путем деления длины ряда V на степени двойки (в нашем случае делимым выступает ширина холста W; обозначим количество фрагментов холста, на которые он делится по ширине, как Тогда где K таково, что При этом на диапазон значений γ накладываются ограничения: отбрасываются те γ, для которых а также и что улучшает статистическую и вычислительную эффективность алгоритма. Данный подход к определению размеров клеток и предлагается использовать в настоящей статье. Таким образом, размер клетки для анализа графиков временных рядов вычисляется как где таково, что Выполнение указанного неравенства и будет являться условием прекращения итераций. Помимо сложностей с определением диапазона размеров клетки, существует еще один вопрос, остро встающий при адаптации метода клеточного покрытия для его применения к временным рядам: как выбрать значения W и H - ширину и высоту абстрактного холста, на который помещается график ряда? Поскольку размер клетки δ зависит от W, правильнее говорить не об абсолютных значениях W и H (сами по себе они могут быть произвольными), а об их соотношении, т.е. соотношении сторон холста. Строго говоря, величина W может принимать любое значение (логичным выглядит «привязать» его к ширине окна, которое отображает график временного ряда, измеряемой в пикселях). В свою очередь, величина H зависит от величины W и поэтому должна вычисляться как произведение W и некоторого коэффициента причем (чтобы холст не оказался в высоту больше, чем в ширину). От выбора значения этого коэффициента зависит, чему будет равна величина D - оценка фрактальной размерности ряда. Проблема выбора значения коэффициента r рассмотрена в работе [21]. В ней также приведена иллюстрация, демонстрирующая, что при различных соотношениях сторон холста величина различается даже при неизменном размере клетки δ (рис. 5). Рис. 5. Графики биржевой цены закрытия, выраженной в долларах США. Количество клеток, покрывающих анализируемый временной ряд, при разных соотношениях сторон холста отличается [21] Вопрос, каково истинное соотношение W и H, т.е. каково истинное значение r, ответа не имеет. Очевидно, что ширина графика W связана с длиной ряда T, высота H - с размахом Однако единицы измерения графика ряда по горизонтали и по вертикали не сопоставимы друг с другом: мерой ширины графика ряда служит время единицей измерения W являются, например, торговые дни на бирже; мерой высоты графика ряда служат биржевые показатели (цены, доходности и другие), единицей измерения H являются, например, рубли, если речь идет о ценах, или безразмерная величина, если речь идет о доходностях. Величину r можно зафиксировать. Например, если то холст будет квадратным, т.е. иметь соотношение сторон 1:1. Согласие об использовании определенного фиксированного значения r обеспечит сопоставимость результатов при оценке размерности различных временных рядов, однако «абсолютные значения размерности, полученные в результате такого анализа, будут оставаться по существу произвольными» [21]. Ввиду этого Pilgrim, Taylor рекомендуют для анализа временных рядов использовать методы оценки размерности, специально адаптированные для них, в частности метод Херста (R/S-анализ). В то же время, как показал анализ литературы, классический клеточный метод продолжает использоваться для оценки размерности временных рядов, а следовательно, отказаться от его применения (во всяком случае, от сравнения с другими методами), на наш взгляд, неправильно. Мы предлагаем три способа определения значения коэффициента r: 1. Первый способ состоит в использовании но не единицы, а, например, 0,7, так как обычно графики временных рядов растянуты по ширине сильнее, чем по высоте. 2. Второй способ заключается в том, чтобы попытаться соотнести «расстояние» между соседними значениями ряда по горизонтали с «расстоянием» между ними же по вертикали Шаг по горизонтали фиксированный и равен единице. В качестве шага по вертикали может быть взята минимальная абсолютная величина изменения между соседними значениями ряда V, а именно: при этом количество шагов вычисляется как Коэффициент r рассчитывается как Однако при коэффициент что приводит к вытянутости графика по вертикали, а при коэффициент r вычислить не представляется возможным. На наш взгляд, для более приближенной к реальности ситуации в качестве шага по вертикали необходимо брать не минимальное, а среднее изменение исследуемой переменной по модулю, т.е. В этом случае 3. Третий способ состоит в том, что коэффициент r определяется как отношение абсолютного минимального значения ряда к размаху: т.е. предполагается, что высота холста достигает его ширины в том случае, если рассматриваемый показатель на протяжении рассматриваемого временного интервала увеличивается вдвое относительно своего минимума. Наконец, стоит отметить, что величина r может принимать значения как меньше, так и больше единицы. Ввиду этого, чтобы высота холста не превысила его ширину, предлагается вычислять высоту следующим образом: Далее перейдем к апробации описанной методики нахождения оценки фрактальной размерности с помощью адаптированного к временным рядам метода клеточного покрытия. 3. Апробация разработанной методики Было решено для программной реализации методики написать функцию на языке программирования Python 3.7.5. Выбор этого высокоуровневого языка обусловлен его простотой, популярностью и существованием большого количества подключаемых библиотек для научных и инженерных расчетов. На рис. 6 приведены скриншоты окна разработанной программы, содержащего холст с графиком некоторого временного ряда длиной 300 при различных размерах клетки. На рис. 7 представлен график линейной регрессии для полученный в результате анализа того же ряда. а б Рис. 6. Визуальная реализация метода клеточного покрытия, где количество фрагментов холста по ширине: а - γ = 5; б - γ = 26 Рис. 7. График линейной регрессии логарифма количества закрашенных клеток от логарифма длины стороны клетки Несмотря на то, что для написания программы был выбран язык Python, при ее тестировании нами было решено не ограничиваться одним этим инструментом и воспользоваться также возможностями языка R, что связано с наличием в нем некоторых библиотек с подходящими для наших целей готовыми алгоритмами. Разумным видится в качестве данных взять сгенерированные временные ряды, которые характеризуются различными состояниями и для которых заранее известна величина фрактальной размерности. В зависимости от характера развития случайного процесса временные ряды могут находиться в одном из трех состояний, каждое из которых определяется величиной фрактальной размерности D: 1) тренд 2) случайное блуждание 3) флэт [3]. Также известно, что размерность является мерой изрезанности или гладкости графика. С нашей точки зрения, имеет смысл проанализировать ряды со следующими характеристиками: 1) ряд, порожденный нестационарным процессом с тенденцией, низкой волатильностью и гладким графиком (его размерность должна стремиться к единице, т.е. 2) ряд, порожденный нестационарным дробно-интегрированным процессом с длинной памятью 3) ряд, порожденный нестационарным процессом «случайное блуждание» 4) ряд, порожденный стационарным дробно-интегрированным процессом 5) ряд, порожденный стационарным процессом с очень высокой волатильностью и изрезанным графиком Опишем методику апробации. Чтобы сделать вывод о состоятельности реализованного алгоритма, следует с его помощью оценить размерность сгенерированных с определенными параметрами рядов и сравнить оценку с величиной априорной размерности. Ряды 1 и 5, графики которых представляют собой с некоторыми оговорками прямую и заштрихованную плоскость, можно сформировать с помощью простых формул в единственных вариантах, в то время как ряды 2-4 генерируются случайным образом и при каждой генерации отличаются, поэтому имеет смысл сформировать множество вариантов каждого ряда, оценить размерность каждого из них и вычислить ее математическое ожидание, которое и сравнивать с априорной величиной D. Анализ каждого ряда следует проводить с применением различных способов определения коэффициента r, чтобы посмотреть, как данный параметр влияет на результат. Также было решено сравнить результаты нашего алгоритма адаптированного метода клеточного покрытия с результатами других методов вычисления фрактальных показателей, а именно R/S-анализа (данный метод позволяет вычислить показатель Херста H, и размерность определяется по формуле и box-counting метода, описанного в работе [7]. Для этого были использованы функции compute_Hc() из библиотеки hurst (Python) и fd.estim.boxcount() из библиотеки fractaldim (R). Длина генерируемого ряда - 4096, количество генераций - 5000. Ряд 1 (обозначим его V1), представляющий собой гладкую прямую, сформируем с помощью следующей модели: Ряд - заштрихованная область: Ряды нами решено генерировать в виде фрактального броуновского движения (ФБД), которое, с одной стороны, может быть смоделировано с заданной величиной размерности, а с другой стороны, способно правдоподобно описывать реальные процессы на финансовых рынках (например, динамику цен акций) [25]. Для симуляции таких рядов воспользуемся функцией fbm() из библиотеки somebm (язык R). Входным параметром данной функции является величина H, которая связана с D через уже упомянутое соотношение Ряд сгенерирован с В последнем случае ряд порожден винеровским процессом, являющимся моделью случайного блуждания с непрерывным временем. На рис. 8 приведены графики тестовых рядов всех видов (DT означает теоретическую, ожидаемую величину фрактальной размерности D). а б в г д Рис. 8. Графики сгенерированных временных рядов с заданной величиной фрактальной размерности: а - ряд б - один из вариантов ряда в - один из вариантов ряда г - один из вариантов ряда д - ряд Также для рядов 2-4 был проведен ADF-тест (расширенный тест Дики - Фуллера на наличие единичного корня), т.е. была выполнена проверка стационарности. Для этого использовалась функция adf.test() из библиотеки tseries (R) со следующими параметрами: количество лагов - 15, альтернативная гипотеза - ряд стационарен. Тест проводился для каждого из 5000 вариантов рядов, после чего было вычислено математическое ожидание p-значения статистики с 99%-ным доверительным интервалом. Результаты: В первых двух случаях для отклонения гипотезы о нестационарности оснований недостаточно, поэтому ряд можно считать нестационарным; в последнем случае ряд стационарен. В следующем подразделе представлены полученные результаты и их анализ. 4. Результаты и их обсуждение В таблице ниже приведены полученные по итогам проведенных вычислений значения фрактальной размерности. Оценка фрактальной размерности тестовых рядов различными методами Ряд Эмпирически оцененная величина размерности DE Наша реализация (r = 0,7) Наша реализация (r = 1) Наша реализация (r вычислен способом 2) Наша реализация (r вычислен способом 3) R/S-анализ (готовая функция) Box-counting метод (готовая функция) V1 1,024 1,028 1,028 -** 0,920 1,036 V2 1,187 ± 0,002* 1,211 ± 0,002* 1,108 ± 0,001* -** 1,243 ± 0,001* 1,211 ± 0,002* V3 1,382 ± 0,002* 1,414 ± 0,002* 1,096 ± 0,000* -** 1,450 ± 0,001* 1,415 ± 0,002* V4 1,590 ± 0,002* 1,623 ± 0,002* 1,062 ± 0,000* 1,546 ± 0,003* 1,658 ± 0,001* 1,624 ± 0,002* V5 1,982 2,000 -** 2,000 1,911 2,000 Примечание: *Для рядов, сгенерированных случайным образом по 5000 раз, приведено математическое ожидание размерности с 99%-ным доверительным интервалом. **При некоторых вариантах рядов вычисление DE оказалось невозможным. Можно констатировать, что факт чувствительности метода клеточного покрытия к входным параметрам в ходе экспериментов подтвердился. Коэффициент r, величиной которого определяется соотношение сторон холста, достаточно сильно влияет на полученную оценку размерности: чем его значение выше (т.е. холст более вытянут по высоте), тем больше оказывается DE. Очевидно, это связано с увеличением количества клеток на плоскости и, как следствие, ростом числа закрашенных клеток N(δ) при неизменном размере клетки δ, что, в свою очередь, приводит к увеличению угла наклона прямой регрессии в логарифмических координатах. Это обусловлено спецификой методики: клетки квадратные и величина δ зависит только от ширины холста W, но не от его высоты. Также можно утверждать, что на сгенерированных данных предложенные нами способы автоматического расчета r (способы определения r под номерами 2 и 3) оказались неприменимы. Оба из них базируются на попытке соотнести между собой несопоставимые величины: длину ряда и размах его значений. Способ 2, основанный на сопоставлении длины ряда и числа вхождений в размах среднего абсолютного прироста показателя, в результате дал очень низкий коэффициент высоты холста, так как при большой выборке (4096 в нашем случае) величина оказывается чрезвычайно мала, а низкое значение r приводит к вне зависимости от характера анализируемого ряда. В некоторых же случаях, как можно заметить на примере ряда вычисление D с таким подходом к определению коэффициента высоты холста и вовсе невозможно: в каком-то из вариантов ряда значение r получилось настолько малым, что H оказалась меньше единицы, а так как в нашей программной реализации безразмерная величина H соответствует высоте холста в пикселях, т.е. целому числу, то высота составила 0, что и привело к невозможности продолжения работы алгоритма. Что касается способа 3, согласно которому r определяется по абсолютному минимальному значению ряда, отнесенному к размаху, он в нашем случае также оказался несостоятельным по той же причине, что и предыдущий способ: в каких-то из вариантов рядов принимало значения, близкие к нулю, что приводило к и невозможности вычисления Кроме того, при данном подходе r зависит не только от размаха значений ряда, но от минимального значения: если к каждому члену ряда прибавить некоторую константу, это не изменит форму графика и вероятностные характеристики ряда, но повлияет на оценку D, что является некорректным. Использование рассмотренных подходов, несостоятельных в случае анализа сгенерированных рядов, может иметь смысл при исследовании рядов, отражающих реальные экономические процессы (динамику курсов акций, валютных пар и так далее), так как при разработке этих подходов предполагалось, что характер изменения значений анализируемых переменных (процессов и явлений) подчиняется природным закономерностям. Например, сложно себе представить, чтобы цена акции или биржевой индекс многократно увеличивались относительно своего минимума на рассматриваемом временном интервале при условии, что процесс на этом интервале существенно не меняет законов своего развития, т.е. исследователем производится локальный анализ процесса (рассмотрен в работе [4]). Вопросы применимости методов фрактального анализа к рядам различной длины, отражающим динамику процесса на разных временных интервалах, требуют отдельного исследования, а пока можно сделать вывод о неприменимости предложенных способов определения соотношения сторон холста при анализе сгенерированных рядов и необходимости их тестирования на реальных экономических рядах. Из двух вариантов r меньшую разницу между и в четырех случаях из пяти показало значение уступив лишь при анализе ряда где графиком является прямая линия. Поскольку ряды 2-4, несомненно, являются более приближенными к реальным процессам, можно заключить, что квадратная форма холста позволяет получать более точные оценки размерности. Наиболее точные результаты были получены с помощью R/S-анализа - метода, специально разработанного для оценки фрактальных характеристик хаотических рядов с длинной памятью. Значения вычисленные на основе этого метода, оказались наиболее близки к хотя и не достигли их с высокой точностью: абсолютная разница между математическим ожиданием и составила от 0,007 до 0,092. В случае ряда величина оказалась некорректной (плоская фигура не может иметь размерность меньше единицы, как и значение показателя Херста H не должно выходить за границы диапазона но это можно объяснить весьма специфическим характером ряда: на практике такие временные ряды, где графиком является идеальная прямая, вряд ли существуют. Следует заметить, что для генерации фрактального броуновского движения использовалась функция в R, а для оценки H была применена функция из библиотеки Python. Эта же библиотека имеет функцию для генерации ФБД, и возможно, что оценки H сгенерированных с ее помощью рядов оказались бы ближе к величинам, указанным при генерации, однако нами сознательно было решено для большей объективности исследования использовать программные библиотеки разных авторов. Другой вывод, который вытекает из того, что данные для исследования сформированы случайным образом, но с явным указанием вида распределения и параметра H, состоит в следующем. Даже в «идеальных» условиях (т.е. в условиях искусственного характера исследуемых временных рядов) оценивание фрактального показателя с высокой точностью является неординарной задачей. При анализе рядов, отражающих реальные процессы (т.е. в «естественных» условиях), сложность этой задачи, несомненно, повышается в связи с неизвестностью истинного значения H этих рядов. Также стоит отметить, что методика R/S-анализа требует для достаточно надежного определения H нескольких тысяч наблюдений [4]. В нашем случае использовались ряды из 4096 значений. При исследовании более коротких рядов точность R/S-анализа относительно точности клеточного метода может снизиться, однако этот момент требует провести сравнение нескольких методов в разных условиях. Наконец, реализация метода клеточного покрытия из библиотеки fractaldim (R) показала значения отличающиеся от вычисленных нашим алгоритмом с на уровне погрешности. Особенностью реализации метода клеточного покрытия в R является то, что по умолчанию для анализа используют только первые значений входного ряда, остальные отбрасывают. В эксперименте мы использовали ряды длиной а в действительности сгенерировано было по 4097 значений и ряды именно такой длины подавались на вход функции в R. Если бы количество значений в ряду было меньше хотя бы на единицу, тогда величина оценивалась бы по 2049 наблюдениям, т.е. практически по половине всего ряда. По нашему мнению, отбрасывание такого объема данных является некорректным, может привести к искаженным оценкам фрактальных показателей, поэтому в предложенной нами методике используются все значения ряда вне зависимости от того, является длина ряда степенью двойки или нет. Таким образом, разработанный алгоритм метода клеточного покрытия, адаптированный к временным рядам, корректно воплощает в себе идею метода оценивания фрактальной размерности; его программная реализация на языке Python может быть использована при анализе финансовых временных рядов. Проведенные эксперименты показали, что получаемые с помощью нашей реализации оценки D не уступают в точности оценкам, полученным с помощью другой существующей в R реализации метода. Однако наша программа уступает программе в R по вычислительной эффективности, что выражено в большей длительности ее выполнения, поэтому требуется ее оптимизация. Заключение В настоящее время моделирование и прогнозирование финансовых временных рядов является одним из наиболее актуальных и необходимых направлений исследований, так как позволяет участникам, которых непосредственно затрагивают финансовые процессы, определять благоприятные исходы и вырабатывать перспективные решения. В качестве инструментов, помогающих выработать данные решения, можно использовать методы фрактального анализа временных рядов. Один из таких методов - метод клеточного покрытия - нами был изучен и детально проработан. В результате была разработана методика, адаптирующая метод клеточного покрытия для оценки размерности временных рядов без формирования компьютерного изображения графика ряда, рассмотрены возникающие при такой адаптации трудности, предложены способы их разрешения. Также была разработана программа для использования данной методики на практике и произведена ее апробация на сгенерированных данных. По итогам проведенной работы можно сделать следующие основные выводы: 1. Метод клеточного покрытия для фрактального анализа финансовых временных рядов достаточно часто используется в исследованиях, проводимых отечественными авторами. Однако последними не уделяется должного внимания особенностям и проблемам применения метода к временным рядам. 2. В то же время эти проблемы являются достаточно важными. Во-первых, метод обладает высокой чувствительностью к регулирующим параметрам, а именно: соотношению сторон плоскости с графиком, покрываемой клетками; используемому диапазону значений длины сторон клетки; количеству разбиений на клетки плоскости. Во-вторых, выбор наиболее подходящих значений этих параметров неочевиден. Наибольшие трудности вызывает определение корректного соотношения сторон, потому что оно требует сопоставления несоизмеримых величин. Результаты проведенного эксперимента свидетельствуют, что наиболее точными являются оценки размерности при соотношении сторон 1:1, т.е. квадратной форме плоскости. При этом предложенные нами способы автоматического определения соотношения сторон плоскости применительно к сгенерированным рядам показали себя несостоятельными. 3. R/S-анализ лишен вышеупомянутых проблем, так как разработан непосредственно для анализа временных рядов. В нашем эксперименте он позволил получить более точные оценки размерности. Однако у этого метода есть ограничение, которое заключается в необходимости исследования выборок большого объема. Таким образом, утверждения о его более высокой эффективности во всех случаях преждевременны, необходимо выполнить сравнение точности различных методов в тех или иных условиях, а именно - при анализе временных рядов различной длины. Итак, в ходе изучения метода клеточного покрытия для оценивания фрактальной размерности пространственных структур была разработана методика, адаптирующая этот метод к временным рядам. По результатам проведенного исследования видно, что разработанный метод на сгенерированных рядах уступает по точности оценивания фрактальных характеристик широко известному методу в области фрактального анализа данных - R/S-анализу, однако не уступает реализации метода клеточного покрытия из библиотеки fractaldim языка R. Вполне закономерными дальнейшими шагами в области исследований фрактального анализа временных рядов могут являться следующие: выполнить сравнение точности различных методов оценивания фрактальной размерности на сгенерированных и реальных временных рядах различной длины; осуществить сравнение методов фрактального анализа и p-адического анализа при моделировании и прогнозировании финансовых временных рядов; определить условия применимости различных методов; провести апробацию разработанных способов определения соотношения сторон плоскости с графиком на реальных экономических данных.

About the authors

R. V Garafutdinov

Реrm State University

S. A Akhunyanova

Реrm State University

References

Statistics

Views

Abstract - 73

PDF (Russian) - 26

Refbacks

• There are currently no refbacks.