Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
13
Добавлен:
16.01.2025
Размер:
7.83 Mб
Скачать

ряд идентичен ряду , то мы получим скейлограмму для ряда .

Представляет особый интерес вычисление кросс-вейвлет

преобразования в случае, когда используется комплексный

вейвлет (например, вейвлет Морле). Тогда

( )| ( )| ( )

( )

(

)

(

)

(

|

(

)|

 

|

(

)||

 

 

)| .

( )

( )/

Таким образом, вычисляя кросс-вейвлет преобразование

можно извлечь значение разности фаз между коэффициента-

ми вейвлет преобразования для двух врменных рядов.

На Рис. 32 показаны значения

( ) для радов Т и К.

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

а

б

Рис. 32 – Кросс-вейвлет преобразование с использованием

вейвлета мексиканская шляпа для рядов Т и К (а) и рядов Т и Х (б)

Если проинтегрировать коэффициенты вейвлет преобразования по времени, то получим вейвлет-кросс-корреляционную

121

меру (Wavelet Cross-Correlation Measure), которая зависит от

( )

|∫

( )

( ) |

 

|∫

( ) |

масштаба

∫| ( )|

∫| ( )|

 

∫| ( )|

∫| ( )|

 

 

мера помогает обнаружить корреляцию между сиг-

Такая

 

 

 

 

 

налами, которые содержат колебания с разной амплитудой или фазой, но, тем не менее, коррелированны между собой

(Рис. 33).

Также можно расширить определение вейвлет кросс-

корреляционной

меры, если ввести

зависимость от сдвига

между рядами

 

|∫

( )

(

)

|

 

(

) √∫| ( )|

∫| ( )|

 

Рис. 33 – Зависимоть вейвлет-кросс-корреляционной меры от

масштаба для временных рядов Т и К

Для исследования связи между компонентами временных рядов локально в плоскости преобразования, мы можем определить квадратичную оценку когерентности вейвлетов следующим образом:

122

где

( )

 

|

(

)

(

) |

 

обозначает локальную операцию сглаживания, как во

 

 

|

(

)|

|

(

)|

 

временной шкале, так и в масштабной, при этом сглаживание выполняется на компонентах преобразования.

Методы кросс-вейвлет анализа используют при исследова-

нии свойств нескольких временных рядов, зависимых между собой нетривиальным образом. Например, в геофизике возникает задача выявления причинно-следственных связей или

корреляции между метеорологическими или другими явлениями окружающей среды, которые происходят на большом расстоянии друг от друга. В [Maraun, 2004] анализируются особенности применения кросс-вейвлет преобразования и

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

[Adamowski, 2008] для изучения метеорологических временных рядов и данных об уровне потока реках. Из временных рядов двух типов выделяли компоненты, которые далее использовались в модели прогнозирования наводнений. В работе показано, что использование кросс-вейвлет анализа полез-

но в том случае, когда существует относительно стабильный сдвиг фазы между потоковым и метеорологическим временными рядами. С помощью кросс-вейвлет преобразования

определялась разность фаз между потоковыми и метеорологическими данными, что улучшило качество модели прогнозирования наводнений.

Другим примером является [Labat, 2010], где кросс-вейвлет

анализ проводился для климатических индексов и показателей сброса пресной воды в Африке. В этом случае, кросс-

вейвлет преобразование и оценка когерентности использовались для визуализации и анализа периодических колебаний в данных длиной в 2-8 лет. В [Kelly, 2003] подтверждается, что методы на основе кросс-вейвлет анализа могут быть эффек-

тивным инструментом в поиске квазипериодичности временного ряда.

Методы кросс-вейвлет анализа также нашли применение в медицине. В [Li, 2007] описывается использование кросс-

123

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

В [Aguiar-Conraria, 2008] используются инструменты кросс-

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

Данные, полученные с помощью кросс-вейвлет преобразо-

вания, также могут использоваться как исходные данные для алгоритмов классификации. В [Dey, 2010] коэффициенты кросс-вейвлет преобразования подавались на вход искусственной нейронной сети и классификатора Fuzzy.

Дискретное вейвлет преобразование и приближение функций с помощью ряда

Дискретное вейвлет преобразование определяется та-

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

быстрому вычислению вейвлет преобразования и его обратно-

 

( )

 

 

 

 

 

 

го [Addison, 2017].

 

 

 

 

 

Пусть

 

принадлежит пространству

периодических

квадратично

интегрируемых функций.

Тогда

( )

можно

представить в виде ряда Фурье

 

 

 

 

 

( )

 

 

 

 

где коэффициенты имеют вид:

)

 

 

 

 

 

 

∫ (

 

 

 

 

 

 

124

 

 

 

 

Набор функций

 

( )

– это ортонормированный ба-

зис

в пространстве

, построенный с помощью мас-

( )

преобразования

( ) из базовой функции

штабного .

 

 

(

)( )

Пусть теперь

 

 

. Базисной функций в пространстве

 

должна быть

функция, которая достаточно быстро убы-

 

 

( )

 

вает к 0 на

. Поэтому, для построения базиса используют-

ся( вейвлеты)

– хорошо

локализованные солитоноподобные

функции. Для того чтобы покрыть вейвлетами всю действительную ось, используют перенос вдоль оси. Для простоты

можно использовать целые сдвиги и аналоги синусоидаль-

ной частоты, как степени двойки

 

. Вейвлет

 

 

семейство функ-

 

 

называется ортогональным, если(

)

 

ций

{

( ) образует ортонормированный базис в

 

.

 

2.6}.5. Корреляция с шаблоном

( )

С помощью непрерывного вейвлет-преобразования выяв-

ляются участки исследуемого ряда, которые по форме наиболее похожи на вейвлет (Рис. 34). Идея состоит в том, чтобы

сравнить части ряда с некоторым шаблоном на разных масштабах (Рис. 35). При этом вейвлет как функция должен об-

ладать определенными математическими свойствами, в частности быстро убывать к нулю на бесконечности. В некоторых случаях полезно использовать шаблон, который не соответствует требованиям к вейвлету. Для этого вместо вейвлет-

преобразования будем вычислять корреляцию между частью

 

(

)

 

 

(

̅)(

 

̅)

 

 

временного ряда и некоторым шаблоном

 

 

 

Полученный

 

√∑

 

(

 

̅)∑ (

 

̅)

 

 

 

коэффициент

 

зависит от значений

.

То есть параметр

отвечает сдвигу шаблона, а

( )

 

 

 

 

параметр

соответствует количеству точек в шаблоне и в

рассматриваемом отрезке ряда. Параметр

 

в данном случае

является аналогом масштаба s, который использовали при вейвлет-преобразовании.

125

Если при вычислении коэффициента( ) вейвлет-

преобразования всегда использовался весь временной ряд, то в данном случае для вычисления используются точек ряда и шаблон длины .

Рис. 34 – Отрезок временного ряда с наложенным шаблоном

Рис. 35 – Шаблон «змея» с разным количеством точек

126

Полученные корреляционные коэффициенты ( ) представим на графике, который похож на скейлограмму (Рис.

36).

а

б

 

в

 

 

Рис. 36 – Корреляционные коэффициенты

 

вычисленные

для рядов Т (а), К (б) и Х (в) с

использованием шаблона, пока-

 

( )

 

занного на Рис. 3.22

2.6.6. Фрактальный анализ

Термин фрактал ввел и популяризировал Бенуа Мандельброт. Чаще всего фракталами называют геометрические объекты, которые имеют сильно изрезанную форму и обладают свойством самоподобия.

Строго и общепринятого определения фрактала в данный момент не существует, хотя Бенуа Мальдеброт использовал несколько пробных определений. Одно из них, введенное в

[Mandelbrot, 1982], звучит так:

127

Фракталом называется множество, размерность Хаусдор- фа-Безиковича которого строго больше его топологической

размерности.

Строгое определение размерности Хаусдорфа-Безиковича

или фрактальной размерности будет введено позже. Суть такого определения сводится к тому, чтобы выделить класс сильно изрезанных объектов, для описания которых недостаточно топологической размерности. Например, существуют кривые, топологическая размерность которых всегда равна 1, но они изогнуты таким сложным образом, что заполнят плоскость или пространство. Так кривые Пеано, проходят через любую точку единичного квадрата. Другой пример – траектория броуновской частицы, которая не является гладкой ни в одной точке.

Первое определение, хотя и является строгим, но исключает многие физические фракталы, и поэтому не используется. Было предложено следующее определения фрактала:

Фракталом называется структура, состоящая из частей, которые в каком-то смысле подобны целому.

Второе определение подчеркивает, что отличительным признаком для фрактала является самоподобие. Приведем строгое определение самоподобного множества, которое используется в математике. Для этого потребуется несколько предварительных определений.

Пусть

множество

замкнуто.

|

Тогда

отображение

То есть

 

| ( )

( )|

|

 

 

 

называется

отображением подобия (similarity) на

,

если

 

 

 

 

 

 

 

 

 

отображение подобия

превращает множество

в

геометрически подобное множество.

 

 

 

 

Рассмотрим набор отображение подобия

. Множе-

ство

является инвариантным относительно преобразо-

ваний , если

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

Множество, которое является инвариантным относительно набора отображений подобия, называется самоподобным.

128

Определение самоподобного множества можно понять интуитивно. Действительно, по определению множество самоподобно, если его можно «собрать» из кусочков, которые по-

добны целому множеству. Тогда простейшим примером само-

и

. Тогда , -

(,

, -) -

(, -).

подобного множества будет отрезок

. Возьмем, например

Очевидно, что простого самоподобия не достаточно для того, чтобы назвать объект фракталом. В самом деле, не будем же мы считать фракталами отрезок прямой или листок бумаги в клеточку. Мы будем рассматривать фрактальные объекты, которые обладают свойством самоподобия, а также сложной структурой.

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

 

 

Удалим из

среднюю треть, останется множество

0,

 

-1 0

 

1. Множество

состоит из двух отрезков, из

 

 

каждого из которых теперь удалим среднюю треть, получим множество . Продолжая повторять* + эту процедуру, получим последовательность множеств . Множество Кантора –

это пересечение

Заметим, что множество

самоподобно. Возьмем отобра-

( )

 

 

 

 

 

 

 

 

 

 

жения подобия

 

 

и

 

 

 

 

,

тогда

 

. С другой

стороны известно, что мера Лебега множества

( )

 

 

( )

 

 

 

 

 

( )

Кантора равна 0, точно также как и для точки, или любого счетного множества. Но очевидно, что структура множества Кантора намного более сложная, и здесь уже возникает идея, что для описания такого множества нужна специальная мера, которой и станет фрактальная размерность множества.

129

Фрактальная размерность

Вернемся к определению размерности Хаусдорфа-

Безиковича. Для этого нам понадобится понятие покрытия множества.

Пусть

 

– непустое множество в

. Диаметр множества

по определению равен

| |

 

 

жеств

*

+

 

 

 

 

Если

 

 

и| |

*|

для|

любого , то набор мно-

 

 

+

 

 

 

называется

покрытием для множества .

Пусть – подмножество некоторого замкнутого множества

в

. Для произвольных

|

*

и

+

определим

}

 

( )

{∑|

 

 

где

инфинум

берется

по

всем

возможным

покрытиям

множества . По определению s-размерная мера Хаусдорфа

Такой предел

существует для любого множества

, с

( )

( )

той оговоркой, что часто он равен нулю или бесконечности. Размерность Хаусдорфа-Безиковича множества опреде-

ляется как

 

( )

+

(* )

( )

+

-же самое

Или, что то( )

*

Для примера вычислим размерность множества Кантора .

( )

{

( )

Выше была описана процедура построения, и согласно с ней на n-м шаге имеется отрезков длины каждый и далее

множество только уменьшается. Поэтому в качестве диаметра

покрытия можно взять величину

и использовать

множеств в покрытии. По определению

 

от предела по

, к пределу

и теперь можно перейти ( )

( )

 

130

 

Соседние файлы в папке !Экзамен зачет 25-26 год