тает сравнительно медленно-но она работает. На рис. 12.2 показаны корреляционные интегралы для аттрактора Хенона, полученные с использованием восстановленного фазового пространства по одной переменной способом, рассмотренным выше. График представляет собой зависимость Сщ от Я в двойных логарифмических координатах для размерности вложения, равной 3. Регрессия представляет собой прямую линию. Оценка D составляет 1.25 против 1.26, полученной при измерении методом оконного скейлинга. Метод Грассберге-ра и Прокаччи предлагает надежный, относительно простой способ оценки фрактальной размерности в том случае, когда наблюдается только одна динамическая переменная. Этот метод требует обработки большого количества данных и соответственно много машинного времени, однако его результаты отличаются надежностью.

ПОКАЗАТЕЛИ ЛЯПУНОВА

Мы не можем рассчитать полный спектр показателей Ляпунова, поскольку нам неизвестны уравнения движения. Однако Уолфом и др. (Wolf et all., 1985) был разработан метод для расчета наибольшего показателя Ляпунова Li по экспериментальным данным. Показатель Li, больший О, свидетельствует о том, что имеется чувствительная зависимость от начальных условий и что в системе существует странный аттрактор. С помощью этого метода измеряется разбегание близлежащих точек в восстановленном фазовом пространстве и указывается, как изменяется скорость разбегания за пределами фиксированных иитррва TTOR вррмртги

Сначала выбираются две точки, отстоящие одна от другой по крайней мере на один орбитальный период. После установленного интервала времени («эволюционного периода») измеряется расстояние между этими двумя точками. Если расстояние становится слишком большим, ищется точка замены с той же угловой ориентацией движения, что и у исходной точки. Относительная угловая ориентация новой пары точек должна быть по возможности ближе к ориентации пары исходной.

Точка замены является необходимой, так как мы измеряем только растяжение, или разбегание в фазовом пространстве. Если эти точки отстоят слишком далеко, они сложатся одна в другую, что будет мерой конвергенции. Конвергенция не

Рис. 12.3. Схематическое представление алгоритма Уолфа для оценки наибольшего показателя Ляпунова по временному ряду (воспроизведено с разрешения Financial Analysts Journal).

есть часть Li. На рис. 12.3 схематически представлен этот алгоритм. Формально он записывается следующим образом:

(12.4)

Теоретически при бесконечном количестве незашумленных данных уравнение (12.4) эквивалентно уравнению (11.3). В реальности мы должны иметь дело с ограниченным количеством зашумленных данных - это означает, что размерность вложения т, временной лаг t и допустимые максимальное и минимальные дистанции между точками должны выбираться с осмотрительностью.

К счастью, Уолф дал несколько «правил большого пальца» для обработки экспериментальных данных. Во-первых, размсрпогтт. Б.тожгпи" должна бьггт. болтлпе, чг?.т у ф-зкл-го пространства аттрактора - это не является неожиданным, учитывая вьшюизложенное относительно фрактальных размерностей. Что ново, однако, так это то, что размерность вложения должна быть больше, чем следующее ближайшее целое число, поскольку грубая поверхность часто выглядит более гладкой в более высокой размерности. Эта размерность, однако, не должна быть слишком высокой, ибо данные становятся слишком редкими, когда мы восстанавливаем фазовое пространство. В результате оказывается слишком мало кандидатур на точки замены.

Временной лаг может быть рассчитан по соотношению (12.1). Уолф с соавторами установили, что эволюция длинЫ не должна превышать 10% длины аттрактора в фазовом пр

странстве. В сущности максимальная длина не должна превышать 10% разности между максимальной и минимальной величинами временного ряда. Уолф с соавторами пришли к это-jiy числу экспериментально, так как логика, стоящая за этими рассуждениями, по-видимому, не поддается точному количественному анализу. В целом же, на мой взгляд, это вполне работоспособное правило. Временная эволюция должна быть достаточно долгой для того, чтобы измерить растяжение, но не складывание. И здесь нет правила - просто чем короче, тем лучше. Здесь должен быть соблюден компромисс: хотя короткие периоды эволюции требуют больше вычислений, при этом имеет место меньше точечных замен, и в результате сходимость более устойчива.

После завершения вычислений по длинному временному ряду получается устойчивая величина Ly. Если пет устойчивой сходимости, это означает, что параметры выбраны неверно, или недостает циклов данных для анализа, или система в действительности не является нелинейной.

Требования к данным, обусловленные использованием алгоритма Уолфа, варьируют в зависимости от сложности системы. Как минимум, нам нужно 10 точек и 10" орбитальных периодов. Поэтому, если размерность аттрактора равна 2, необходимо 100 точек данных; если же размерность равна 6, то требуется один миллион точек. Определение размерности является решающим условием, если предпринимается попытка расчета показателя Ляпунова.

В Приложении 5 описана бейсик-программа для расчета показателя Li по данным временного ряда.

Ня рис 12.4 показана сходимость показателя 1-, к зпаче нию 0.45 для аттрактора Хенона, который при использовании уравнений движения проявил спектр показателей Ляпунова 2 диапазоне 0.4 - 1.6.

В заключение следует сказать о природе «экспериментальных данных» в теории инвестиций, экономической теории и физических науках. В физических науках экспериментальные Данные получаются из управляемого эксперимента. В задаче •конвекции жидкости, например, данные собирают только тог-Да, когда температура достаточно высока для индуцирования Турбулентности. Эти данные анализируют с целью установления того, является ли турбулентность действительно хаотиче-Кой, со странным аттрактором или просто случайной.