Главная > ЗНАКОМСТВО С НЕЛИНЕЙНОЙ ДИНАМИКОЙ (В.С.Анищенко)
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
След.
Макеты страниц

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

Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике

ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO

Предположим, что в результате эксперимента получены значения некоторой физической величины a, то есть набор ai=a(iΔt), i=1,,N. Будем считать, что временной ряд ai порождается некоторой ДС с непрерывным или дискретным временем, представляя собой дискретизованную с шагом Δt одномерную проекцию фазовой траектории. Эта траектория принадлежит аттрактору системы, размерность которого равна d. Согласно Такенсу, задать вектор состояния можно по методу (9.2). Размерность пространства вложения n определяется на основании формулы Манэ (9.4), то есть требует предварительного вычисления величины d. На практике значение n, определяемое формулой (9.4), часто оказывается завышенным, и можно ограничиться пространством меньшей размерности (n<dd).

Вначале необходимо определить размерность аттрактора d. Для ее оценки обычно вычисляют так называемую корреляционную размерность, используя формулу
Dc=limϵlimNlgC(ϵ,N)lgϵ,

где C(ϵ,N)=N2ieqjv(ϵ|xixj|) — корреляционный интеграл, ϵ — размер ячейки разбиения фазового пространства, N — число точек, используемых для оценки размерности, v функция Хевисайда, xi=x(iΔt). Для определения Dc строят зависимость lgC(ϵ,N) от lgϵ и ищут на ней линейный участок, нак.тон которого и определяет искомое значение размерности. Кроме того, иногда анализируют зависимость Dc(n) и увеличивают n до тех пор, пока Dc не достигнет насыщения.

Известно, что при вычислении Dc существуют ограничения на величину ϵ. Если ϵ приближается к размерам аттрактора ϵmax, то линейная зависимость lgC(ϵ,N) от lgϵ пропадает, что объясняется влиянием границ аттрактора, где число соседей у каждой точки обычно меньше, чем в \»середине\». В пределе, если ϵ=ϵmax, то lgC(ϵ,N)=1. С другой стороны, при уменьшении значения ϵ существует некоторое ϵmin, такое, что для ϵ<ϵmin структура аттрактора остается неразрешенной. Как следствие, вновь нарушается линейная зависимость lgC(ϵ,N) от lgϵ.

Одной из проблем при расчете размерности является выбор величин N и Δt. Целесообразнее их подбирать не по отдельности, а с учетом размера временного окна NΔt. При этом следует принимать во внимание существование фундаментальных ограничений на величину Dc, определяемых формулой
Dmax=2lgNlg(1/r),r=ϵϵmax,

которая означает, что алгоритм расчета размерности не может дать значение больше, чем Dmax при заданном числе точек N. Иными словами, если r=0.1 и N=1000, то Dmax6; если N=100000, то Dmax10. В результате возникают определенные сложности в случае, когда требуется отличить детерминированный процесс в системе с достаточно большим числом степеней свободы от истинно случайного процесса.

Знание размерности позволит нам реконструировать аттрактор методом (9.2), как было описано выше.

Метод задержки Такенса является наиболее известным, но не единственным способом задания вектора состояния. Альтернативой ему служит так называемый метод последовательного дифференцирования, имеющий определенные преимущества при решении задачи реконструкции. Идея данного метода следующая. Пусть имеется временной ряд a(iΔt)=ai,i=1,,N. Задание вектора состояния в фазовом пространстве производится следующим образом:
x(t)=(a(t),da(t)/dt,,dn1a(t)/dtn1)=(x1,x2,,xn).

Поскольку известны значения ai только в дискретные моменты времени iΔt, координаты xj вектора x определяются путем численного дифференцирования исходного временного ряда по приближенным математическим формулам. Очевидно, что точность вычисления производных будет определяться малостью величины шага дискретизации Δt. Недостатком метода является повышенная чувствительность к шуму, что ограничивает его применимость для пространств вложения большой размерности (по крайней мере, без проведения предварительной процедуры фильтрации).

1
Оглавление
email@scask.ru