МГУ им М В Ломоносова http chaos phys

МГУ им. М. В. Ломоносова http: //chaos. phys. msu. ru

Содержание 1. 2. 3. 4. 5. Введение Энтропия и размерность Показатели Ляпунова Анализ временного ряда Анализ электрокардиограмм

1. Введение Два основных подхода к исследованию природных процессов: 1. Моделирование 2. Анализ данных Идентификация изучаемого процесса

x. Исследуемый процесс: n 10 20 30 40 50 60 70 80 80 90 100 110 120 130 n 140 . . . Конечная реализация: временной ряд Задача идентификации: определить свойства системы. Задача прогноза: предсказать или l >1. При этом порождающий механизм неизвестен

Формализация понятия «временной ряд» Измерения: отображение множества натуральных чисел в некоторое множество. Основной пример: эволюция динамической системы с компактным пространством состояний X, определяемое отображением Временные реализации имеют вид.

Задачи: • Разделение случайного и детерминированного сигнала • Определение эффективного числа степеней свободы • Выявление скрытой динамики системы (периодических и квазипериодических составляющих) • Разладка. Сегментация • Возможность управление ситуацией • Прогноз

Разладка Резкое изменение свойств наблюдаемого ряда (размерности, энтропии, показателей Ляпунова), происходящее в неизвестный заранее момент времени.

Априорная ? Апостериорная

2. Энтропия и размерность Основные вопросы • Прогноз развития ситуации. • Как отличить случайность от хаоса? • Каково число степеней свободы? • Какие характеристики нужны для этого? ЭНТРОПИЯ ПРОЦЕССА

Энтропия S – мера беспорядка в системе. Пример, когда S растет: Pi – вероятности оказаться системе в состоянии i.

Эта задача прояснит смысл энтропии На раздаче в баре совершено убийство; яд находится водном из 1200 стаканов. Сколько надо сделать проб для определения, в каком стакане яд? Если положение шара заранее неизвестно и мы его узнаём, то тем самым мы увеличиваем суммарную информацию на 1 бит. Необходимо два вопроса, т. е. максимальное количество информации 2 бит. Или от числа возможных состояний, которое равно 4. Рост беспорядка связан с ростом незнания о состоянии системы:

Поскольку беспорядок есть понятие теории информации, энтропию, показывающую, насколько система хаотична (неупорядочена), можно определить соотношением Шеннона: это средняя скорость потери информации о состоянии динамической системы с течением времени Как рассчитать энтропию и другие характеристики реального наблюдаемого процесса, т. е. временного ряда?

РАЗМЕРНОСТЬ Размерность D – геометрическая характеристика. Она показывает, сколько переменных необходимо взять, чтобы однозначно описать систему. D=1 x

D = 1, 2 или 3

Каковы энтропия и размерность случайного, регулярного и хаотического процессов? Эти величины характеризуют сложность. Для случайных систем Для хаотических процессов 0

– независимые величины преобразование Бернулли 1 Это детерминированное преобразование порождает случайные последовательности для почти любого н/у.

Теорема (F. Takens). Для любой дифференцируемой динамической системы с компактным пространством состояний энтропия либо , либо Система называется хаотической, если хаотической в пространстве ее состояний существует такое подмножество положительной меры, что для почти любой траектории с выполняется

Как рассчитать энтропию H и размерность D реального процесса, т. е. временного ряда? xn n 10 20 30 40 50 60 70 Временной ряд 80 80 90 100 110 120 130 140 . . .

Назовем (n, ε)-разделенными два отрезка, и , траектории длины n, если существует такое, что . ,

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

Мера того, как растет при увеличении n, называется энтропией процесса: Энтропия говорит о ! горизонте предсказуемости

«Проклятие» размерности Неожиданный пример: 1 r

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

Это означает, что минутные данные после преобразования масштаба не должны отличаться от часовых. a f ( f ( … (x))) +h= Xn+1 = f (Xn)

Шум Каждая точка окружена «облаком» :

Здесь наклон с возрастанием размерности не меняется Здесь наклон с ростом увеличивается

3. Показатели Ляпунова Основное свойство траекторий: W W s u

Показатели Ляпунова x 1(t) ξ(0) ξ(t) x 2(t) Теорема Я. Песина: Для однородных подмножеств

4. Анализ временного ряда Как найти энтропию наблюдаемого процесса? • К 1 – энтропия Колмогорова-Синая, • K 2 – корреляционная энтропия (нижняя граница К 1).

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

Предельную зависимость K 2 можно представить в виде K 2 -энтропия

Как найти D? • Функциональный метод • Корреляционный интеграл • Метод ложных соседей Необходимо ввести спектр размерностей Dq

Функциональный метод Зафиксируем одну из последовательностей, начинающуюся с номера Определим Наличие функциональной зависимости: При выполняется

При выполняется Система Лоренца

Метод ложных соседей

5. Анализ электрокардиограмм В данной части приведены результаты, полученные в Рентгенохирургическом центре интервенционной кардиологии Главного военного клинического госпиталя им. акад. Н. Н. Бурденко совместно с группой кардиохирургов во главе с проф. А. В. Ардашевым. Все больные перед проведением исследования в течение 12 ч. не курили, не употребляли напитков, содержащих кофеин; в течение 36 часов не употребляли алкоголь; за 1 ч. до исследования рекомендовался легкий прием пищи. Исследование проводилось после 15 -минутного отдыха в горизонтальном положении в условиях покоя в тихой комнате при спокойном освещении от 11. 00 до 16. 00.

Для кардинального лечения разного рода аритмий был использован достаточно прогрессивный метод радиочастотной абляции (РЧА) сердечной ткани. Задача: изучить практическую значимость характеристик, известных из теории временных рядов. Фазовый портрет сердечного ритма практически здорового человека

Фазовый портрет ритма пациентов, страдающих атриовентрикулярной тахикардией. Исходные количественные характеристики, рассчитанные на основе анализа ритмов у больных с АВТ и лиц контрольной группы Параметры Больные с АВУРТ (n=25) Группа контроля (n=20) D 0 2, 45 0, 26 2, 48 0, 18 D 1 2, 55 0, 34 2, 52 0, 32 Н 0, 18 0, 08 0, 66 0, 29 D 2 2, 52 0, 21 2, 54 0, 26 Результаты свидетельствуют о большей упорядоченности ритма у больных с данными нарушениями сердечного ритма.

Фазовый портрет ритма пациентов, страдающих атриовентрикулярной тахикардией, в первые сутки после РЧА. Изменение характеристик хаотичности сердечного ритма у больных с атриовентрикулярной тахикардией после РЧА Параметры Исходно Через 6 ч. Через 24 ч. Через 2 мес. Через 6 мес. Через 12 мес. D 0 2, 45 0, 26 2, 56 0, 27 2, 64 0, 14 2, 70 0, 26 2, 56 0, 16 2, 49 0, 19 D 1 2, 55 0, 34 2, 48 0, 38 2, 52 0, 34 2, 55 0, 35 2, 37 0, 35 2, 41 0, 44 Н 0, 18 0, 08 0, 47 0, 40 0, 28 0, 16 0, 42 0, 14 0, 42 0, 22 0, 43 0, 26 D 2 2, 52 0, 21 2, 44 0, 26 2, 42 0, 27 2, 45 0, 24 2, 45 0, 29 2, 46 0, 15 Происходит увеличении степени хаотичности колебаний

Сравнение количественных мер хаотичности сердечного ритма у пациентов с рецидивами и без рецидивов после РЧА Параметр ы D 0 D 1 Н D 2 Подгруппа Исходно Через 6 ч. Через 24 ч. Через 2 мес. Через 6 мес. Через 1 год А 2, 47 0, 25 2, 60 0, 23 2, 51 0, 27 2, 49 0, 26 2, 52 0, 22 2, 42 0, 40 Б 2, 50 0, 30 2, 63 0, 21 2, 60 0, 22 2, 55 0, 22 2, 60 0, 15 2, 44 0, 18 А 2, 51 0, 34 2, 43 0, 32 2, 46 0, 33 2, 42 0, 36 2, 36 0, 38 2, 41 0, 33 Б 2, 44 0, 41 2, 43 0, 40 2, 50 0, 42 2, 38 0, 42 2, 37 0, 43 2, 45 А 0, 57 0, 14 0, 38 0, 14 0, 48 0, 11 0, 56 0, 15 0, 66 0, 25 0, 62 0, 28 Б 0, 37 0, 17 0, 35 0, 14 0, 21 0, 09 0, 50 0, 21 0, 61 0, 22 0, 70 0, 2 А 2, 55 0, 26 2, 47 0, 30 2, 56 0, 29 2, 52 0, 17 2, 45 0, 31 2, 41 0, 22 Б 2, 45 0, 32 2, 55 0, 35 2, 52 0, 13 2, 49 0, 23 2, 52 0, 26 2, 44 0, 2 А – подгруппа пациентов без рецидивов; Б – подгруппа пациентов с рецидивами.

Фазовый портрет ритма пациентов, страдающих трепетанием предсердий. Параметры Группа пациентов с ТП до РЧА (n=23) Группа клинического сравнения (n=20) D 0 2, 61 0, 30 2, 48± 0, 18 D 1 2, 28 0, 36 2, 52± 0, 32 H 0, 47 0, 25 0, 66± 0, 29 D 2 2, 53 0, 19 2, 54± 0, 26 Показатели свидетельствуют о большей упорядоченности сердечного ритма у больных с ТП.

Фазовый портрет ритма пациентов, страдающих трепетанием предсердий, в первые сутки после РЧА. Параметры Исходно Через 6 ч. Через 24 ч. Через 2 мес. Через 6 мес. Через 1 год DF 2, 61 0, 30 2, 47 0, 24 2, 55 0, 31 2, 67 0, 19 2, 75 0, 16 2, 62 0, 32 DI 2, 28 0, 36 2, 59 0, 31 2, 65 0, 35 2, 59 0, 26 2, 66 0, 34 2, 26 0, 33 H 0, 47 0, 15 0, 32 0, 09 0, 42 0, 09 0, 54 0, 23 0, 83 0, 25 0, 71 0, 40 DC 2, 53 0, 19 2, 52 0, 14 2, 54 0, 22 2, 66 0, 16 2, 48 0, 35 2, 53 0, 17