Zagadnienia

14. Liniowe szeregi czasowe

Stacjonarność rzędu 2. Funkcje tworzące. Klasyczne modele liniowe MA, AR, ARMA i ARIMA. Ułamkowe ruchy Browna. (1 wykład)

14.1. Szeregi czasowe stacjonarne rzędu 2

Na początek ustalmy przestrzeń probabilistyczną

Ω,F,P.

Przez L2⁢Ω będziemy oznaczać zbiór zmiennych losowych X takich, że E⁢X2<∞. Jak łatwo zauważyć jest to przestrzeń liniowa nad R. Ponadto, gdy utożsamimy zmienne losowe, które są prawie wszędzie równe to funkcja

XL2=E⁢X2

jest normą, a L2⁢Ω/∼⁢ jest przestrzenią Banacha (por. [9] dodatek C.5 lub [3] §2.10).

Dla ustalenia uwagi przyjmiemy następującą definicję szeregu czasowego.

Definicja 14.1

Szeregiem czasowym nazywamy ciąg zmiennych losowych Xt o wartościach rzeczywistych, t∈Z

…,X-n,…,X-1,X0,X1,…,Xn,….
Definicja 14.2

Szereg czasowy Xt nazywamy stacjonarnym rzędu 2 gdy

a.∀t⁢⁢⁢⁢Xt∈L2⁢Ω;
b.∀t1,t2⁢⁢⁢⁢E⁢Xt1=E⁢Xt2;
c.∀t1,t2,n⁢⁢⁢⁢c⁢o⁢v⁢Xt1+n,Xt2+n=c⁢o⁢v⁢Xt1,Xt2.

Zauważmy, że czasami stacjonarność rzędu 2 implikuje silną stacjonarność.

Lemat 14.1

Gaussowski szereg czasowy stacjonarny rzędu 2 jest silnie stacjonarny w sensie definicji 8.6.

Dowód.
Jeśli Xt jest procesem Gaussowskim to dla kazdego s>0 łączny rozkład X1,…,Xs jest rozkładem normalnym N⁢μ,Σ, gdzie

μ1=…=μs=E⁢X1,⁢⁢⁢Σi,j=c⁢o⁢v⁢Xi,Xj.

Podobnie dla dowolnego n łączny rozkład Xn+1,…,Xn+s jest rozkładem normalnym N⁢μ′,Σ′, gdzie

μ1′=…=μs′=E⁢Xn+1,⁢⁢⁢Σi,j′=c⁢o⁢v⁢Xn+i,Xn+j.

Ponieważ

E⁢Xn+1=E⁢X1,⁢⁢⁢ i ⁢⁢⁢⁢c⁢o⁢v⁢Xn+i,Xn+j=c⁢o⁢v⁢Xi,Xj,

to oba rozkłady są identyczne.

□

Podstawowymi narzędziami służącymi do opisu stacjonarnych szeregów czasowych są funkcje autokowariancji. Dla każdego stacjonarnego rzędu 2 szeregu czasowego Xtt=-∞+∞ definiujemy funkcję γ

γ:Z⟶R,

przyporządkowującą liczbie całkowitej k - k-ty współczynnik autokowariancji (patrz definicja 11.1)

γ⁢k=c⁢o⁢v⁢X0,Xk.

Zauważmy, że ze stacjonarności wynika, że dla każdego n

c⁢o⁢v⁢Xn,Xn+k=γ⁢k,

a w szczególności

D2⁢Xn=γ⁢0.
Twierdzenie 14.1

Niech k≥0, wówczas funkcja autokowariancji stacjonarnego szeregu czasowego spełnia następujące warunki

1.γ⁢-k=γ⁢k;
2.γ⁢0≥0;
3.γ⁢k≤γ⁢0;
4.∀a-k,…,ak∈R⁢⁢⁢⁢∑i,j=-kkai⁢aj⁢γ⁢i-j≥0.
Uwaga 14.1

Funkcje spełniające warunek 4 dla każdego k nazywa się dodatnio określonymi.

Dowód.
Punkt 1 wynika z symetrii kowariancji.

γ⁢-k=c⁢o⁢v⁢X0,X-k=c⁢o⁢v⁢X-k,X0=c⁢o⁢v⁢X0,Xk=γ⁢k.

Punkt 2 jest oczywisty.

γ⁢0=c⁢o⁢v⁢X0,X0=D2⁢X0≥0.

Punkt 3 wynika ze związków między korelacją i kowariancją.

γ⁢k2=c⁢o⁢v⁢X0,Xk2≤D2⁢X0⁢D2⁢Xk=γ⁢02.

Punkt 4 wynika z nieujemności wariancji zmiennej losowej. Rozważmy zmienną losową

Xa=∑i=-kkai⁢Xi.

Korzystając z dwuliniowości kowariancji otrzymujemy

0≤D2⁢∑i=-kkai⁢Xi=∑i,j=-kkai⁢aj⁢c⁢o⁢v⁢Xi,Xj=∑i,j=-kkai⁢aj⁢γ⁢i-j.
□

Prawdziwe jest również twierdzenie odwrotne.

Twierdzenie 14.2

Jeśli funkcja γ:Z→R spełnia warunki 1,2,3 i 4 z powyższego twierdzenia to jest ona funkcją autokowariancji pewnego stacjonarnego szeregu czasowego.

Dowód
Patrz – [3] Theorem 1.5.1.

Zachowanie się funkcji autokowariancji dla dużych n ma istotne znaczenie w praktycznych zastosowaniach. Wyróżnia się dwie klasy szeregów czasowych stacjonarnych rzędu 2.

Definicja 14.3

Niech X będzie szeregiem czasowym stacjonarnym rzędu 2, a γX jego funkcją autokowariancji.
Gdy γX⁢n zbiega do zera dla dużych n w sposób wykładniczy

∃C>0⁢⁢⁢⁢∃δ∈0,1⁢⁢⁢⁢γX⁢n<C⁢δn,

to mówimy, że X jest procesem o krótkiej pamięci.
W przeciwnym przypadku mówimy, że X jest procesem o długiej pamięci.

Gdy szereg czasowy nie jest prawie na pewno stały czyli gdy γ⁢0>0 to definiujemy dodatkowo funkcję autokorelacji

ρ:Z⟶-1,1,⁢⁢⁢ρ⁢k=γ⁢kγ⁢0.

Kluczowym obiektem w teorii szeregów czasowych stacjonarnych rzędu 2 jest biały szum. Jest to przykład procesu o (bardzo!) krótkiej pamięci.

Definicja 14.4

Stacjonarny rzędu 2 szereg czasowy ε nazywamy białym szumem gdy

E⁢X0=0⁢⁢⁢⁢γ⁢0>0⁢⁢⁢⁢ i⁢⁢⁢γ⁢t=0⁢⁢⁢⁢ dla ⁢t≠0.

Gdy dodatkowo γ⁢0=1 to ε nazywamy unormowanym białym szumem.

Zbiór szeregów czasowych będących białym szumem o wariancji σ2=γ⁢0 będziemy oznaczać W⁢N⁢σ2.

Uwaga 14.2

Kolejne wyrazy białego szumu są nieskorelowane ale nie muszą być niezależne, chyba, że jest to gaussowski biały szum.

14.2. Sploty vel filtry

Na początek przypomnimy, kiedy ciąg liczb rzeczywistych jest klasy lp.

Definicja 14.5

Ciąg ann=0∞ należy do klasy lp gdy jest sumowalny w p-tej potędze.

ann=0∞∈lp⟺∑n=0∞anp<∞.
Uwaga 14.3

Zauważmy, że l2 zawiera l1. Ponadto dla p≥1 klasy lp są unormowanymi przestrzeniami liniowymi

anlp=∑n=0∞anp1/p.
Definicja 14.6

Niech X=Xtt=-∞+∞ będzie szeregiem czasowym, a a=ann=0∞ ciągiem o wyrazach rzeczywistych. Wówczas szereg czasowy Y=Ytt=-∞+∞, którego wyrazy dają się przedstawić jako sumy nieskończone w L2⁢Ω

Yt=∑i=0∞ai⁢Xt-i,

nazywamy splotem ciągu a i szeregu X.

Operację splotu będziemy oznaczać przez ”⋆”

Y=a⋆X.
Uwaga 14.4

W niektórych źródłach operację splotu z ciągiem an nazywa się filtrem o współczynnikach an.

Podamy teraz dwa warunki gwarantujące istnienie splotu Yt ciągu a=an i szeregu Xt.

Twierdzenie 14.3

Jeżeli zachodzi jeden z poniższych warunków

1.an∈l1⁢⁢⁢⁢ i ⁢⁢⁢⁢∀t∈Z⁢⁢⁢⁢E⁢Xt≤μ,⁢⁢⁢D2⁢Xt≤σ2.
2.an∈l2,⁢⁢⁢∀t∈Z⁢⁢⁢⁢E⁢Xt=0,⁢⁢⁢D2⁢Xt≤σ2,⁢⁢⁢ oraz ⁢⁢⁢⁢Xt⁢ nieskorelowane ze sobą.

to istnieje splot ciągu a=an i szeregu Xt.

Dowód.
Skorzystamy z zasady majoryzacji.
Ad 1.

YtL2≤∑i=0∞ai⁢Xt-iL2=∑i=0∞ai⁢Xt-iL2≤∑i=0∞ai⁢σ2+μ2=σ2+μ2⁢al1.

Ad 2.

YtL22=D2⁢∑i=0∞ai⁢Xt-i=∑i=0∞ai2⁢D2⁢Xt-i≤σ2⁢∑i=0∞ai2=σ2⁢al22.
□
Twierdzenie 14.4

Gdy szereg czasowy Xt jest stacjonarny rzędu 2, a ciąg a=an należy do l1 to ich splot Yt jest stacjonarny rzędu 2.

Dowód.
Pokażemy, że wartość oczekiwana i autokowariancje szeregu Yt nie zmieniają się przy przesunięciu. Oznaczmy przez μ wartość oczekiwaną Xt a przez γX funkcję autokowariancji.

E⁢Yk=E⁢∑i=0∞ai⁢Xk-i=∑i=0∞ai⁢E⁢Xk-i=∑i=0∞μ⁢ai=∑i=0∞ai⁢E⁢X-i=E⁢Y0.
c⁢o⁢v⁢Yn,Yk+n=∑i,j=0∞ai⁢aj⁢c⁢o⁢v⁢Xn-i,Xn+k-j=∑i,j=0∞ai⁢aj⁢γX⁢k+i-j=
=∑i,j=0∞aiajcov(X-i,Xk-j)=cov(Y0,Yk).
□
Twierdzenie 14.5

Gdy szereg czasowy εt jest białym szumem, a ciąg a=an należy do l2 to ich splot Yt jest stacjonarny rzędu 2. Ponadto

E⁢Yt=0,⁢⁢⁢γY⁢k=σ2⁢∑i=0∞ai⁢ak+i,

gdzie σ2 wariancja ε0.

Dowód.
Powtarzamy rozumowanie z poprzedniego dowodu i wstawiamy E⁢εt=0 oraz γε⁢0=σ2 i γε⁢t=0 dla t≠0.

□

Okazuje się, że z dokładnością do pewnego ”nieistotnego” składnika, wszystkie szeregi czasowe stacjonarne rzędu 2 można przedstawić w postaci splotu z białym szumem.

Twierdzenie 14.6

Niech X będzie szeregiem czasowym stacjonarnym rzędu 2 o zerowej wartości oczekiwanej, E⁢X0=0. Wówczas istnieją ciąg a klasy l2 i unormowany biały szum ε takie, że

X=a⋆ϵ+V,

gdzie szereg czasowy V jest nieskorelowany z ε

∀s,t⁢⁢⁢⁢E⁢Vs⁢εt=0

i ponadto jest zawarty w przecięciu domknięć (w L2) podprzestrzeni generowanych przez początkowe Xs

∀tVt∈⋂n=-∞+∞ClL2(lin{Xs;s≤n}).

Dowód.
Powyższe twierdzenie jest wnioskiem z ”rozkładu Wolda” - [3] Theorem 5.7.1.

Operacja splotu przeprowadza szeregi czasowe stacjonarne rzędu 2 na szeregi czasowe stacjonarne rzędu 2. Zatem sploty można iterować.

Twierdzenie 14.7

Niech ciągi a i b należą do przestrzeni l1, a X będzie szeregiem czasowym stacjonarnym rzędu 2. Wówczas

a⋆b⋆X=c⋆X,

gdzie c jest ciągiem o wyrazach

cn=∑i=0nai⁢bn-i.
Uwaga 14.5

Ciąg c z powyższego twierdzenia nazywa się iloczynem Cauchy'ego ciągów a i b. Jak łatwo sprawdzić c też należy do l1. Zatem splot definiuje działanie algebry l1 na zbiorze szeregów czasowych, które są stacjonarne rzędu 2. Ponadto mnożenie Cauchy'ego jest przemienne gdyż

∑i=0nai⁢bn-i=∑j=0nan-j⁢bj,

zatem dla a,b∈l1

b⋆a⋆X=a⋆b⋆X.

Dowód twierdzenia.
Oznaczmy przez Y splot b i X, a przez Z splot a i Y. Wówczas

Yt=∑i=0∞bi⁢Xt-i,⁢⁢⁢Zt=∑j=0∞aj⁢Yt-j.

Przedstawiamy szereg Z w zależności od szeregu X

Zt=∑j=0∞aj⁢∑i=0∞bi⁢Xt-j-i=∑i,j=0∞aj⁢bi⁢Xt-j-i.

Podstawiamy k=i+j i porządkujemy powyższą sumę względem Xt-k.

Zt=∑k=0∞∑j=0kaj⁢bk-j⁢Xt-k=∑k=0∞ck⁢Xt-k.
□

Gdy szereg X z twierdzenia 14.7 jest białym szumem to możemy osłabić założenia dotyczące ciągu b. W analogiczny sposób jak twierdzenie 14.7 dowodzi się następujące twierdzenie.

Twierdzenie 14.8

Niech ciąg a należy do przestrzeni l1, ciąg b do l2, a ε będzie białym szumem. Wówczas

a⋆b⋆ε=c⋆ε,

gdzie c jest ciągiem o wyrazach

cn=∑i=0nai⁢bn-i.

14.3. Funkcje tworzące

Mnożenie Cauchy'ego ciągów jest ściśle związane z mnożeniem szeregów potęgowych. Otóż

∑n=0∞an⁢zn⁢∑n=0∞bn⁢zn=∑n=0∞cn⁢zn,

gdzie ciąg c jest iloczynem Cauchy'ego ciągów a i b

cn=∑i=0nai⁢bn-i.

Dlatego przyporządkujemy ciągom funkcje tworzące czyli sumy szeregów potęgowych. Pozwala to wykorzystać aparat analizy zespolonej do badania szeregów czasowych.

Definicja 14.7

Funkcję holomorficzną

A⁢z=∑n=0∞an⁢zn

nazywamy funkcją tworzącą ciągu liczbowego a=ann=0∞.

Zauważmy, że istnieje zależność między promieniem zbieżności RA szeregu potęgowego A⁢z a klasą ciągu a.

Lemat 14.2
RA>1⟹a∈l1,
RA<1⟹a∉l2.

Dowód.
Gdy promień zbieżności szeregu potęgowego A jest większy od 1 to jest on zbieżny bezwzględnie w punkcie z∗=1. Zatem

∑n=0∞an=∑n=0∞an⋅z∗n<∞.

Natomiast gdy ciąg a jest klasy l2 to musi on zbiegać do 0

limn→∞⁡an=0.

Zatem

lim supn→∞⁡ann≤1.

Z czego wynika, że promień zbieżności szeregu potęgowego A jest nie mniejszy niż 1

RA=lim supn→∞⁡ann-1≥1.

Więc jeśli promień zbieżności RA jest mniejszy od 1, to ciąg a nie należy do l2.

□

Funkcje tworzące dla funkcji autokowariancji określamy jako sumy szeregów Laurenta (por. [15] Cz. I, §6).

Definicja 14.8

Funkcję holomorficzną

Γ⁢z=∑-∞+∞γ⁢n⁢zn

określoną na pierścieniu

1r<z<r,⁢⁢⁢r>1,

nazywamy funkcją tworzącą ciągu liczbowego γ=γ⁢n-∞+∞.

Uwaga 14.6

Współczynniki rozwinięcia funkcji Γ w szereg Laurenta na pierścieniu

1r<z<r,⁢⁢⁢r>1,

są wyznaczone jednoznacznie (por. [15] Cz. I, §6, Twierdzenie 2).

Lemat 14.3

Niech γ będzie funkcją autokowariancji szeregu czasowego X stacjonarnego rzędu 2. Wówczas następujące warunki są równoważne:
1. Szereg Laurenta Γ⁢z jest zbiezny w pewnym pierścieniu

1r<z<r,⁢⁢⁢r>1.

2. Szereg czasowy X jest procesem o krótkiej pamięci.

Dowód.
1⇒2
.
Z nierówności Cauchy'ego dla współczynników szeregu Laurenta ([15] s.120) otrzymujemy, że dla każdego ρ∈1,r

∀n⁢⁢⁢⁢γ⁢n≤Mρn,

gdzie M to maksimum modułu funkcji Γ na okręgu z=ρ.

2⇒1.
Gdy szereg czasowy X jest procesem o krótkiej pamięci to istnieją stałe C i δ, C>0, δ∈0,1, takie, że

∀n⁢⁢⁢⁢γ⁢n<C⁢δn.

Zatem

lim supn→∞⁡γ⁢-nn=lim supn→∞⁡γ⁢nn≤δ<1.

Zatem szereg Laurenta Γ⁢z jest zbieżny na pierścieniu

δ<z<1δ

(por. [15] s.117).

□

14.4. Operator przesunięcia

Oznaczymy przez L operator przesunięcia szeregu czasowego o 1 w prawo

Y=L⁢X⁢⁢⁢⟺⁢⁢⁢∀t∈Z⁢⁢⁢⁢Yt=Xt-1.

Pozwoli to nam zapisać splot Y=a⋆X jako szereg potęgowy iterowanych operatorów przesunięcia

∀t∈Z⁢⁢⁢⁢Yt=∑i=0∞ai⁢Xt-i⁢⁢⁢⟹⁢⁢⁢Y=∑i=0∞ai⁢Li⁢X=A⁢L⁢X.

Z twierdzeń 14.7 i 14.8 wynika następujący wniosek.

Wniosek 14.1

Funkcja tworząca złożenia operatorów A1⁢L i A2⁢L jest iloczynem funkcji A1⁢z i A2⁢z

A1⁢L∘A2⁢L=A1⋅A2⁢L.

Funkcja tworząca funkcji autokowariancji szeregu czasowego X będącego obrazem unormowanego białego szumu ε, X=A⁢L⁢ε wynosi

Γ⁢z=A⁢z⁢A⁢z-1.
Uwaga 14.7

Gdy promień zbieżności szeregu potęgowego A⁢z jest większy od 1 to szereg Γ⁢z=A⁢z⁢A⁢z-1 jest zbieżny na pierścieniu

z:RA-1<z<RA.

Z powyższej uwagi i lematu 14.3 otrzymujemy:

Wniosek 14.2

Gdy promień zbieżności szeregu potęgowego A⁢z jest większy od 1 to szereg czasowy X=A⁢L⁢ε, ⁢ε∈W⁢N⁢σ2, jest procesem o krótkiej pamięci.

14.5. Przykłady

Przedstawimy teraz kilka najpopularniejszych szeregów czasowych stacjonarnych rzędu 2. Niech ε=εt-∞+∞ będzie unormowanym białym szumem.

1. Wielomianowa funkcja tworząca.
Szereg czasowy X postaci

X=A⁢L⁢ε,

gdzie A⁢z wielomian stopnia p nazywa się szeregiem średnich ruchomych rzędu p – M⁢A⁢p.

∀t⁢⁢⁢⁢Xt=a0⁢εt+a1⁢εt-1+…+ap⁢εt-p.

2. Funkcja tworząca jest odwrotnością funkcji wielomianowej.
Szereg czasowy X postaci

X=C⁢L⁢ε,⁢⁢⁢C⁢z=1B⁢z,

gdzie B⁢z wielomian stopnia q, który nie zeruje się na kole jednostkowym

z≤1⁢⁢⁢⟹⁢⁢⁢B⁢z≠0,

nazywa się szeregiem autoregresyjnym rzędu q – A⁢R⁢q.
Szereg X spełnia równanie

B⁢L⁢X=ε,

zatem

∀t⁢⁢⁢⁢b0⁢Xt+b1⁢Xt-1+…+bq⁢Xt-q=εt.

3. Funkcja tworząca jest ilorazem funkcji wielomianowych.
Szereg czasowy X postaci

X=D⁢L⁢ε,⁢⁢⁢D⁢z=A⁢zB⁢z,

gdzie A⁢z wielomian stopnia p a B⁢z wielomian stopnia q, który nie zeruje się na kole jednostkowym

z≤1⁢⁢⁢⟹⁢⁢⁢B⁢z≠0,

nazywa się autoregresyjnym szeregiem średnich ruchomych rzędu q,p – A⁢R⁢M⁢A⁢q,p.
Szereg X spełnia równanie

B⁢L⁢X=A⁢L⁢ε,

zatem

∀t⁢⁢⁢⁢b0⁢Xt+b1⁢Xt-1+…+bq⁢Xt-q=a0⁢εt+a1⁢εt-1+…+ap⁢εt-p.

Biorąc pod uwagę, że funkcje tworzące w trzech powyższych przykładach mają promień zbieżności większy od 1, to na mocy wniosku 14.2 otrzymujemy:

Lemat 14.4

Szeregi czasowe M⁢A, A⁢R i A⁢R⁢M⁢A są procesami o krótkiej pamięci.

Uwaga 14.8

Niech γ będzie funkcją autokowariancji szeregu czasowego X=A⁢L⁢ε, ε∈W⁢N⁢1. Gdy X jest klasy A⁢R to operator A jest wyznaczony przez funkcje autokowariancji γ z dokładnością do znaku. Natomiast dla szeregów czasowych M⁢A i A⁢R⁢M⁢A taka jednoznaczność zachodzi tylko przy dodatkowym warunku, że funkcje tworzące nie zerują się w kole jednostkowym.

Przykład A⁢R⁢M⁢A⁢1,1.
Rozważmy szereg czasowy

X=D⁢L⁢ε,⁢⁢⁢D⁢z=3⁢z+1z+3,⁢⁢⁢ε∈W⁢N⁢1.

Funkcja D⁢z jest holomorficzna na całej płaszczyźnie zespolonej poza punktem z=-3 gdzie ma biegun. Zatem promień zbieżności jej szeregu Taylora w 0 wynosi 3.
Okazuje się, że funkcja tworząca funkcji autokowariancji jest stała

Γ⁢z=D⁢z⁢D⁢z-1=3⁢z+1z+3⁢3⁢z-1+1z-1+3=3⁢z+1z+3⁢3+1⁢z1+3⁢z=1.

Czyli szereg X jest białym szumem. Zatem reprezentacja A⁢R⁢M⁢A szeregu czasowego o zadanej funkcji autokowariancji nie jest jednoznaczna.

14.6. Procesy o przyrostach stacjonarnych

Niech Y będzie dowolnym szeregiem czasowym. Szereg X taki, że

∀t⁢⁢⁢⁢Xt=Yt-Yt-1,

czyli

X=1-L⁢Y,

nazywamy szeregiem przyrostów szeregu Y.

Operację brania przyrostów można iterować. Szereg

X=1-Lk⁢Y=P⁢L⁢Y,⁢⁢⁢P⁢z=1-zk,

nazywamy szeregiem k-tych przyrostów szeregu Y.

Gdy k-te przyrosty szeregu czasowego Y są stacjonarne rzędu 2 i należą do klasy A⁢R⁢M⁢A⁢q,p to mówimy, że szereg Y jest klasy A⁢R⁢I⁢M⁢A⁢q,k,p.

Przykłady
∙
Błądzenie przypadkowe bez dryfu

Xt=Xt-1+σ⁢εt

jest procesem klasy A⁢R⁢I⁢M⁢A⁢0,1,0.

∙ Błądzenie przypadkowe z dryfem

Xt=Xt-1+μ+σ⁢εt

jest procesem klasy A⁢R⁢I⁢M⁢A⁢0,2,1.

∙ Trend liniowy

Xt=a⁢t+b+σ⁢εt

jest procesem klasy A⁢R⁢I⁢M⁢A⁢0,2,2.

∙ Trend wielomianowy stopnia k

Xt=ak⁢tk+ak-1⁢tk-1+…⁢a0+σ⁢εt

jest procesem klasy A⁢R⁢I⁢M⁢A⁢0,k+1,k+1.

14.7. Ułamkowy ruch Browna

Definicja 14.9

Gaussowski proces BtH, ⁢t∈Z, ⁢H∈0,1, taki, że

∀t⁢⁢⁢⁢E⁢Xt=0;
∀t,s⁢⁢⁢⁢c⁢o⁢v⁢BtH,BsH=σ22⁢t2⁢H+s2⁢H-t-s2⁢H,⁢⁢⁢σ>0,

nazywamy ułamkowym ruchem Browna.

Pokażemy, że szereg BH jest procesem o przyrostach stacjonarnych. Niech

Xt=BtH-Bt-1H.
Lemat 14.5

Szereg czasowy XtH jest stacjonarny o zerowej wartości oczekiwanej, E⁢XtH=0, oraz funkcji autokowariancji

γ⁢n=σ22⁢n+12⁢H+n-12⁢H-2⁢n2⁢H.

Dowód.
Ponieważ BH jest szeregiem o zerowej wartości oczekiwanej to to samo zachodzi dla szeregu przyrostów. Natomiast

c⁢o⁢v⁢XtH,Xt+nH=c⁢o⁢v⁢BtH-Bt-1H,Bt+nH-Bt+n-1H=
=σ22((|t|2⁢H+|t+n|2⁢H-|n|2⁢H)+(|t-1|2⁢H+|t+n-1|2⁢H-|n|2⁢H)-
-(|t+n|2⁢H+|t-1|2⁢H-|n+1|2⁢H)-(|t+n-1|2⁢H+|t|2⁢H-|n-1|2⁢H))=
=σ22⁢n+12⁢H+n-12⁢H-2⁢n2⁢H.

Zatem XH jest stacjonarny z funkcją autokowariancji

γ⁢n=σ22⁢n+12⁢H+n-12⁢H-2⁢n2⁢H.
□
Uwaga 14.9

Dla wszystkich XH ⁢⁢γ⁢0=σ2. Dodatkowo dla H=12 i n≠0 ⁢⁢γ⁢n=0. Zatem X1/2 jest gaussowskim białym szumem, a B1/2 jest gaussowskim błądzeniem przypadkowym.

Lemat 14.6

Dla H≠12

limn→∞⁡γ⁢n⁢n2-2⁢H=H⁢2⁢H-1⁢σ2.

Dowód.
Dla n>1

γ⁢nσ2=12⁢n+12⁢H+n-12⁢H-2⁢n2⁢H=n2⁢H2⁢1+1n2⁢H+1-1n2⁢H-2=
=n2⁢H2(1+2⁢Hn+2⁢H⁢2⁢H-12⁢n2+1-2⁢Hn+2⁢H⁢2⁢H-12⁢n2-2+O(n-3))=
=H⁢2⁢H-1⁢n2⁢H-2+O⁢n2⁢H-3.
□
Uwaga 14.10

Gdy H≠12, to szereg czasowy XH jest procesem stacjonarnym o długiej pamięci.

Treść automatycznie generowana z plików źródłowych LaTeXa za pomocą oprogramowania wykorzystującego LaTeXML.

Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego.

Projekt współfinansowany przez Ministerstwo Nauki i Szkolnictwa Wyższego i przez Uniwersytet Warszawski.