Фрактальный анализ временных рядов. Любушин А.А. - 12 стр.

UptoLike

Составители: 

Рубрика: 

12
(iv) Пусть .
0k
Тогда:
0.5 ( ) 0;
00.5 ()
0.5 1 ( ) 0
XX
XX
XX
Hk
Hk
Hk
0;
γ
γ
γ
=⇒
<< <
<< >
=
(v) Пусть . Тогда 0.5H
222
() (2 1)| | , | |
H
XX
kHHk k
γσ
⋅→
(vi) Пусть
1
() ()
2
ik
XX XX
k
Ske
ω
ωγ
π
+∞
=−
=
- спектральная плотность стационарной случайной
последовательности
k
X
; ( ) ( )
ik
XX XX
kSed
π
ω
π
γ
ω
+
=
ω
. Тогда , т.е. если
(2 1)
() , 0
H
XX
S
ωω ω
−−
0.5 1 ( ) , 0
XX
HS
ω
ω
<< и сигнал
k
X
носит низкочастотный характер.
Соответственно, спектральная оценка конечной выборки из фрактального броуновского
движения будет иметь поведение:
(2 1)
() , 0
HH
H
BB
S
ωω ω
−+
5.
Оценка постоянной Херста по экспериментальным данным. Наиболее очевидная
оценка непосредственно из определения показателя, данного самим автором. Пусть
( ), 1,...,
=
x
tt Nанализируемый временной ряд;
<
N
длина скользящего временного
окна;
τ
- номер отсчета правого конца скользящего окна, то есть мы рассматриваем
моменты времени , которые удовлетворяют условию
t 1
+≤Lt
τ
τ
. Пусть - длина
внутреннего временного окна, которое используется внутри текущего основного окна для
операций усреднения. Мы рассматриваем длины внутренних окон, удовлетворяющие
условию: . Пусть
s
/5sL
()
,
1
1
(
=
1)
=
+−
s
su
t
xxut
s
τ
(41)
выборочная оценка среднего значения на интервале длиной отсчетов, который лежит
внутри текущего основного окна и начинается в точке . Следующий шаг состоит в
вычислении отклонений от среднего значения (41), их накопленной суммы и размаха
накопленной суммы:
s
u
() () () () () () ()
,,,,,,
1
() () , () (), max () min ()
=
Δ= =Δ =
t
su su su su su su su
t
t
v
,
x
txtx t xvR t t
ττττττ
ξξ
τ
ξ
(42)
                                                                            12
(iv) Пусть k ≠ 0 .
                                    H = 0.5               ⇒ γ XX (k ) = 0;
Тогда:                              0 < H < 0.5 ⇒ γ XX (k ) < 0;
                                    0.5 < H < 1 ⇒ γ XX (k ) > 0


(v) Пусть H ≠ 0.5 . Тогда γ XX (k ) ∼ σ 2 H (2 H − 1)⋅ | k |2 H − 2 , | k |→ ∞
                                      +∞
                                1
(vi) Пусть S XX (ω ) =
                               2π
                                     ∑γ
                                     k =−∞
                                             XX   (k )e − ikω - спектральная плотность стационарной случайной

                                                     +π
последовательности X k ; γ XX (k ) =
                                                     −
                                                      ∫π S   XX   (ω )eikω dω . Тогда S XX (ω ) ∼ ω − (2 H −1) , ω → 0 , т.е. если

 0.5 < H < 1 ⇒ S XX (ω ) → ∞, ω → 0                           и           сигнал           Xk      носит       низкочастотный                характер.
Соответственно, спектральная оценка конечной выборки из фрактального броуновского
движения будет иметь поведение:
S BH BH (ω ) ∼ ω − (2 H +1) , ω → 0

    5. Оценка постоянной Херста по экспериментальным данным. Наиболее очевидная
оценка непосредственно из определения показателя, данного самим автором. Пусть
x(t ), t = 1,..., N – анализируемый временной ряд; L < N – длина скользящего временного
окна; τ - номер отсчета правого конца скользящего окна, то есть мы рассматриваем
моменты времени t , которые удовлетворяют условию τ − L + 1 ≤ t ≤ τ . Пусть s - длина
внутреннего временного окна, которое используется внутри текущего основного окна для
операций усреднения. Мы рассматриваем длины внутренних окон, удовлетворяющие
условию: s ≤ L / 5 . Пусть


                                                                               1 s
                                                                  xs(,τu) =      ∑ x(u + t − 1)
                                                                               s t =1
                                                                                                                                                  (41)



− выборочная оценка среднего значения на интервале длиной s отсчетов, который лежит
внутри текущего основного окна и начинается в точке u . Следующий шаг состоит в
вычислении отклонений от среднего значения (41), их накопленной суммы и размаха
накопленной суммы:


                                                                    t
           Δ xs(τ, u) (t ) = x(t ) − xs(,τu) , ξ s(,τu) (t ) =    ∑Δ x τ
                                                                   v =1
                                                                              ( )
                                                                              s, u   (v), Rs(τ, u) = max ξ s(,τu) (t ) − min ξ s(,τu) (t )
                                                                                                          t                  t
                                                                                                                                                  (42)