Исследование модели фрактального броуновского движения
Исследование
модели фрактального броуновского движения
Реферат
Дипломная работа содержит 96 страниц, 24 рисунка, 5 таблиц, 11
источников, 1 приложение.
Фрактальное броуновское движение, фрактальный гауссовский
шум, автомодельность, оценка ковариационной функции, оценка параметра Харста,
интерполяция, экстраполяция, прогнозирование, фильтрация Калмана-Бьюси.
В данной работе рассматривается теоретические основы фрактального
броуновского движения (ФБД), вопросы статистического моделирования ФБД на
компьютере, а также применение теории ФБД при статистическом моделировании
процессов стохастической системы, описываемых линейным дифференциальным
уравнением с возмущениями в виде ФБД.
Оглавление
Реферат
Введение
. Теоретические основы фрактального броуновского движения
(ФБД)
.1 Свойство автомодельности
.2 Фрактальное броуновское движение
.3 Фрактальный гауссовский шум
. Моделирование ФБД
. Интерполяция, экстраполяция и прогнозирование процесса ФБД
по наблюдениям в двух точках
. Моделирование дифференциальной системы с возмущениями в
виде ФБД и оценка состояний дифференциальной системы с ФБД на основе фильтрации
Калмана-Бьюси
. Экономическая часть
. Охрана труда и окружающей среды
Заключение
Список использованной литературы
Приложение
Введение
При
статистическом анализе финансовых временных рядов давно было замечено,
что многие из них обладают свойствами (статистического) самоподобия,
проявляющимися в том, что, образно говоря, их «части устроены так же, как и
целое». Например, если 
- дневные значения финансового индекса, то
эмпирические плотности 
и 
, 
, найденные по большому ряду величин
и 
, 
соответственно,
оказываются такими, что
,
где
H - некоторая константа, которая (в отличие от
ожидаемой, согласно центральной предельной теореме, величины 
) значимым образом больше .
Эти наблюдения привели к развитию общей концепции (статистической)
автомодельности, приведшей к появлению понятий фрактального броуновоского
движения и фрактального гауссовского шума. Константа H, упомянутая выше, получила название
параметра Харста, в честь британского климатолога Гарольда Харста открывшего
эффект сильного последействия последовательности зависимых случайных величин
при анализе поведения флуктуаций годичной водности реки Нила. Позднее теория
фрактального броуновского движения получила широкое практическое применение при
анализе финансовых показателей (цен акций, обменных курсов валют), а также ряда
физических явлений, таких как турбулентность.
1.
Теоретические основы фрактального броуновского движения (ФБД)
1.1
Свойство автомодельности
Случайный
процесс 
со значениями в 
называется
автомодельным (самоподобным), или удовлетворяющим свойству (статистической)
автомодельности, если для каждого 
можно
найти такое 
, что
. (1)
Это
означает, что изменение временной шкалы 
приводит
к тому же самому результату, что и изменение фазовой шкалы 
.
Для
ненулевых строго устойчивых процессов существует константа H
такая, что 
. Таким образом, если в данном ранее определении
автомодельного процесса заменить , то
случайный процесс будет называться автомодельным с показателем Харста H,
или процессом, удовлетворяющим свойству статистической автомодельности с
показателем Харста H.
1.2
Фрактальное броуновское движение
Рассмотрим функцию:
. (2)
Эта
функция при 
является неотрицательно определенной, и,
следовательно, на некотором вероятностном пространстве существует гауссовский
процесс с нулевым средним и автоковариационной функцией
,
то
есть, с
. (3)
Отсюда видим, что
.
и,
значит,
.
Таким
образом, можно заключить, что рассматриваемый процесс 
является автомодельным с показателем Харста H.
Непрерывный
гауссовский процесс с нулевым средним и ковариационной функцией (3)
называется (стандартным) фрактальным броуновским движением с показателем
автомодельности Харста (в дальнейшем для такого процесса будет
использоваться обозначение 
).
Из
данного определения следует, что (стандартное) фрактальное броуновское движение
удовлетворяет следующим свойствам, которые можно было
бы также принять в качестве определения этого процесса:
)
, 
для всех
;
)
имеет стационарные приращения:
;
)
является гауссовским процессом,
,
где
.
)
имеет непрерывные траектории.
Из
этих свойств снова следует, что фрактальное броуновское движение обладает
свойством автомодельности.
В
случае 
(стандартное) фрактальное броуновское движение есть
ни что иное, как (стандартное) броуновское движение, или винеровский процесс.
1.3
Фрактальный гауссовский шум
В
прикладной теории вероятностей броуновское движение 
используется в качестве модели, дающей простой способ
получения белого шума. Если положить
, (4)
то
получаемая последовательность 
будет
гауссовской последовательностью независимых одинаково распределенных случайных
величин с 
, 
. Такая
последовательность называется белым (гауссовским) шумом.
По
аналогии с (4) положим
, (5)
и
будем называть последовательность фрактальным
(гауссовским) шумом с параметром Харста H, .
Из
формулы (3) для ковариационной функции (стандартного) процесса 
следует, что ковариационная функция 
имеет следующий вид:
. (6)
Отсюда
видно, что при 
. (7)
Тем
самым, в случае ковариация 
для 
, и 
образует
гауссовскую последовательность независимых случайных величин. Если же 
, то из (7) видно, что ковариация убывает с ростом 
достаточно медленно, что обычно интерпретируется, как
наличие «долгой памяти», или «сильного последействия».
Отметим
принципиальную разницу в случаях 
и 
.
Если
, то ковариация отрицательна (
, 
), при
этом 
.
Если
, то ковариация положительна (
, 
), при
этом 
.
Положительная
ковариация означает, что вслед за положительными (отрицательными) значениями 
следует ожидать также положительные (отрицательные)
значения. Отрицательность же ковариации означает, что вслед за положительными
(отрицательными) значениями следует ожидать отрицательные (положительные)
значения.
фрактальный броуновский моделирование стохастический
2. Моделирование ФБД
. (8)
Для того чтобы смоделировать ФБД по формуле (8), возьмем следующую оценку
фрактального гауссовского шума:
, (9)
где
- комплексно-значные независимые центрированные
гауссовские случайные величины с дисперсией 
.
Дисперсии
случайных величин определяются следующим образом:
, (10)
.
Спектральная
плотность имеет вид:
, (11)
где
ковариационная функция 
определяется по формуле (6).
Распишем
спектральную плотность в формуле (11):
(в
силу четности ковариационной функции)
. (12)
При моделировании ФБД будем использовать оценку спектральной плотности с
конечной суммой слагаемых в формуле (12):
, (13)
где
.
Оценим,
какую ковариационную функцию 
дает
такая оценка спектральной плотности:
Распишем
интеграл I:
Интеграл
II:
Таким
образом:
Поскольку
дальше при расчетах используется 
, можно
говорить о том, что оценка спектральной функции (13) является допустимой в
рамках данной работы.
В
формуле (10) выберем такие 
, 
чтобы 
. Для
этого построим функцию, равную интегралу от спектральной плотности, с
аргументом, равным верхнему пределу интеграла:
.
Поскольку
ряд в выражении (13) сходится, его можно интегрировать почленно:
Далее
для каждого 
, 
вычисляем
.
Поскольку
очевидно, что выбранное разбиение отрезка 
симметрично
относительно нуля, перейдем от , к 
, 
следующим образом:
,
, 
, (14)
.
Определим
случайные величины в формуле (9) следующим образом:
, (15)
где
- независимые одинаково распределенные случайные
величины.
Проверим,
что в случае представления (15) будут
действительно независимыми центрированными гауссовскими случайными величинами с
дисперсией :
)
независимость следует из независимости 
и 
,
)
центрированность следует из центрированности и ,
)
гауссовость следует из гауусовости и ,
)
, что и требовалось доказать.
С
учетом (14) и (15) формулу (9) для оценки фрактального гауссовского шума можно
преобразовать следующим образом:
. (16)
Далее с помощью специально разработанной на языке C++ программы, исходный код которой
представлен в Приложении, по формуле (16) производилось моделирование
реализаций фрактального гауссовского шума, а затем по формуле (8) вычислялись
значения ФБД.
Ниже
представлены графики ФБД для количества шагов 
и
различных значений параметра Харста.
Рис.
1. Реализация ФБД для 
Рис.
2. Реализация ФБД для
Рис.
3. Реализация ФБД для 
Рис.
4. Реализация ФБД для
Графики
ФБД, представленные на рис. 1-4, наглядно показывают различные типы
последействия для различных значений параметра Харста H.
Построим
теперь оценку ковариационной функции и параметра Харста смоделированного
процесса, чтобы убедиться, что смоделированный процесс действительно является
фрактальным броуновским движением.
Оценку
ковариационной функции можно построить по одной траектории, поскольку процесс
является стационарным:
, (17)
где
- наблюдения фрактального гауссовского шума.
Ниже представлены графики оценок ковариационных функций для
смоделированных ФБД, траектории которых представлены на рис. 1-4. На графиках
красным цветом обозначена оценка ковариационной функции, черным - точное
значение ковариационной функции, вычисленное по формуле (6). Графики построены
для 20 шагов.
Рис.
5. Оценка ковариационной функции ФБД для
(траектория процесса представлена на рис. 1)
Рис.
6. Оценка ковариационной функции ФБД для
(траектория процесса представлена на рис. 2)
Рис.
7. Оценка ковариационной функции ФБД для
(траектория процесса представлена на рис. 3)
Оценим
параметр Харста для смоделированных процессов по методу моментов. Подставим в
формулу (6) для точного значения ковариационной функции оценку параметра Харста
и приравняем полученное выражение оценке
ковариационной функции, рассчитанной по формуле (17), при 
.
,
,
.
Для смоделированных ФБД, траектории которых представлены
на рис. 1-4, получаем следующие оценки параметра Харста:
Таблица 1. Оценка параметра Харста методом моментов.
|
Процесс
|
Точное значение
параметра Харста H
|
Оценка параметра Харста
|
|
Рис. 1
|
0.8
|
0.7914
|
|
Рис. 2
|
0.8
|
0.7336
|
|
Рис. 3
|
0.2
|
0.2467
|
|
Рис. 4
|
0.2
|
0.2739
|
3. Интерполяция, экстраполяция и прогнозирование процесса ФБД
по наблюдениям в двух точках
Важным классом задач в теории случайных процессов является построение
оценки неизвестных значений процесса по нескольким известным значениям
(наблюдениям).
В
данной работе рассматривается задача построения с.к.-оптимальной оценки 
по наблюдениям 
. В
зависимости от значений 
можно выделить три подзадачи:
Подзадача
1 (интерполяция)
Рассмотрим
случай, когда 
:
Задачу
построения оценки по наблюдениям при будем называть задачей интерполяции.
Подзадача
2 (экстраполяция)
Рассмотрим
случай, когда 
:
Задачу
построения оценки по наблюдениям при будем называть задачей экстраполяции.
Подзадача
3 (прогнозирование)
Рассмотрим
случай, когда 
:
Задачу
построения оценки по наблюдениям при будем называть задачей прогнозирования.
Для
построения оценки воспользуемся теоремой о нормальной корреляции.
Пусть
-
оцениваемая случайная величина,
- вектор
наблюдений.
Тогда
по теореме о нормальной корреляции:
,
где
, (18)
, (19)
. (20)
Поскольку
, 
,
получаем
. (21)
Подставив в формулу (21) выражения (18), (19) и (20) и упростив
получившееся выражение, получаем:
,
где
.
Для
случая 
(обычное броуновское движение) в задаче интерполяции
аналитически получаем:
, (22)
то есть получаем линейную оценку неизвестного значения по двум
наблюдениям:
Рис.
9. Оценка по двум наблюдениям в задаче интерполяции в случае .
Для
задач экстраполяции и прогнозирования в случае также
получаем, что с.к.-оптимальной оценкой является линейная оценка.
Проанализируем
поведение с.к.-оптимальной оценки в случае 
.
Вычисление
оценки для случая
(график
процесса представлен на рис. 1)
Интерполяция:
|
   Точное
значение  Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
70
|
120
|
75
|
33,02
|
31,24
|
31,69
|
|
70
|
100
|
36,85
|
35,4
|
36,06
|
|
70
|
120
|
110
|
39,82
|
37,07
|
37,42
|
Экстраполяция:
|
Точное
значение Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
100
|
120
|
70
|
30,41
|
34,03
|
26,95
|
|
100
|
120
|
80
|
35,03
|
34,97
|
30,65
|
|
100
|
120
|
90
|
35,88
|
35,91
|
34,05
|
Прогнозирование:
|
Точное
значение Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
70
|
100
|
110
|
39,82
|
39,0
|
38,75
|
|
70
|
100
|
120
|
38,73
|
41,14
|
40,58
|
|
70
|
100
|
130
|
39,47
|
43,29
|
42,35
|
Вычисление
оценки для случая
(график процесса представлен на рис. 3)
Интерполяция:
|
Точное
значение Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
70
|
110
|
80
|
2,382
|
-0,288
|
-0.356
|
|
70
|
110
|
90
|
0,047
|
-0,614
|
-0.563
|
|
70
|
110
|
100
|
0,257
|
-0,941
|
Экстраполяция:
|
Точное
значение Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
100
|
120
|
70
|
0,039
|
1,295
|
0,016
|
|
100
|
120
|
80
|
2,382
|
0,949
|
0,031
|
|
100
|
120
|
90
|
0,047
|
0,603
|
0,062
|
Прогнозирование:
|
Точное
значение Линейная оценкаС.к.-оптимальная оценка
|
|
|
|
|
|
|
70
|
100
|
110
|
-1,267
|
0,33
|
0,187
|
|
70
|
100
|
120
|
-0,435
|
0,402
|
0,172
|
|
70
|
100
|
130
|
1,701
|
0,475
|
0,163
|
4.
Моделирование дифференциальной системы с возмущениями в виде ФБД и оценка
состояний дифференциальной системы на основе фильтрации Калмана-Бьюси
Рассмотрим следующую дифференциальную систему:
(23)
в
которой 
, 
, 
, 
(такая
система описывает процесс Орнштейна-Уленбека). 
-
стандартные процессы фрактального броуновского движения. 
, являются
независимыми.
Для
решения системы (23) перейдем от дифференциалов к конечным разностям:
(24)
Здесь
, 
-
фрактальный гауссовский шум.
, 
.
, , -
независимые.
Рассчитаем
получившуюся рекуррентную формулу (24) при следующих данных: 
, 
, 
, 
, 
.
Для
получаем:
Рис.
10. Решение разностного уравнения (24) . Синим
цветом нарисован график наблюдаемого процесса 
, черным
- оцениваемого процесса .
Для
получаем:
Рис.
11. Решение разностного уравнения (24) для . Синим
цветом нарисован график наблюдаемого процесса , черным
- оцениваемого процесса .
Задача
фильтрации состоит в определении с.к.-оптимальной оценки процесса 
по наблюдениям процесса 
. С.к.-оптимальной оценкой является условное
математическое ожидание:
,
где
есть
-алгебра, порожденная значениями процесса 
до текущего момента времени .
Классическая
теорема Калмана-Бьюси дает уравнение, описывающее эволюцию математического
ожидания 
в случае, когда в системе (23) вместо фрактального
броуновского движения присутствуют обычные винеровские процессы.
Для
рассмотренного случая процесса Орнштейна-Уленбека дифференциальное уравнение
для математического ожидания имеет
вид:
, (25)
где
.
Исследуем
поведение фильтра в случае системы (23) с возмущениями в виде ФБД.
В
выражении (25) перейдем от дифференциалов к конечным разностям:
. (26)
Рассчитаем значения фильтра Калмана-Бьюси вместе со значениями состояния
системы (24):
, :
Ошибка
оценки 
:
, 
:
Ошибка
оценки :
, 
:
Ошибка
оценки :
, 
:
Ошибка
оценки :
, 
:
Ошибка
оценки :
, :
Ошибка
оценки :
Как
видно из представленных графиков, в случае при
уменьшении 
фильтр Калмана дает более точные значения. В случае
же этого не происходит - фильтр Калмана не дает оценки
процесса с приемлемой точностью.
5.
Экономическая часть
В
современном мире информационные технологии играют большую
(а иногда и решающую) роль в успешности бизнеса. Информационные технологии
позволяют автоматизировать производственные процессы, управлять в реальном
времени сложными динамическими системами и организовывать почти мгновенный
обмен информацией, как между людьми, так и между различными информационными
системами.
Интернет
в современном мире является одним из самых мощных коммуникационных инструментов
(уступая пока только телефону). Именно поэтому в настоящее время коммерческие
структуры обращают все больше внимания на использование «Всемирной Паутины» как
эффективного канала коммуникации с потребителями.
Веб-сайты
являются виртуальными представительствами компаний в Интернете. Для многих
компаний веб-сайт является реальным маркетинговым или PR-
инструментом, без которого их бизнес был бы неэффективным или вообще не смог бы
существовать. Веб-сайт позволяет доносить до потенциальных потребителей товаров
или услуг компании нужную информацию именно тогда, когда эта информация нужна,
и с минимальными затратами.
Чем
больше развиваются информационные технологии (в частности Интернет-технологии)
и чем больше проникают эти технологии в жизнь обычных людей, тем более
комплексные задачи ставятся перед веб-сайтами компаний, тем сложнее с точки
зрения производства становятся веб-сайты. Современный веб-сайт - это высоко
технологичный продукт, в котором сочетаются сложные программные решения и
креативные дизайнерские идеи. Если в начале 90-х годов разработка сайта могла
осуществляться от начала и до конца всего лишь одним человеком, то сейчас
производство обычного сайта - это сложный процесс, в который могут быть
вовлечены более 10-15 разработчиков.
В
данной работе я рассмотрю пример разработки сайта для компании «ИксСтрахование»
- одного из крупнейших страховщиков в России на данный момент.
Задача
Разработать корпоративный сайт для компании «ИксСтрахование».
Сайт должен решать следующие задачи:
· Информирование потенциальных и существующих клиентов компании
(как физических, так и юридических лиц) о страховых услугах, которые оказывает
«ИксСтрахование»;
· Предоставление удобных сервисов по расчету стоимости
страховых полисов для посетителей сайта;
· Построение коммуникации с представителями СМИ и
бизнес-аналитиками;
· Поддержка рекламной активности компании в Интернете.
Сетевой график
В рамках разработки сайта выполняются следующие работы:
|
№
|
Название работы
|
Исполнитель
|
|
Проектная аналитика
|
|
1
|
Сбор требований к
разрабатываемому сайту у Заказчика
|
Аналитик
|
|
2
|
Разработка технического
задания
|
Технический писатель
|
|
3
|
Разработка прототипов
страниц
|
Технический писатель
|
|
4
|
Разработка технической
спецификации программных модулей
|
Аналитик
|
|
Дизайн
|
|
5
|
Разработка дизайн-концепции
|
Арт-директор
|
|
6
|
Разработка дизайн-шаблонов
|
Дизайнер
|
|
Программирование
|
|
7
|
Разработка программных
модулей
|
Программист
|
|
8
|
Сборка сайта
|
Программист
|
|
Верстка
|
|
9
|
Верстальщик
|
|
Публикация
|
|
10
|
Тестирование сайта на
соответствие ТЗ и ТС
|
Тестировщик
|
|
11
|
Наполнение сайта
информацией
|
Контент-менеджер
|
|
12
|
|
|
|
|
|
|