Stochastyka ściąga

Mowimy,że {Wt : t>=0} jest procesem Wienera, gdy1) W0 = 0,2) proces ma przyrosty niezależne,3) Wt.Ws ma rozkład normalny z E(Wt-Ws)=0 i Var(Wt-Ws)=t-s dla t>s>=0,4) proces ma ciągłe realizacje. Gdy własności1-4są spełnione dla t∈[0,a],mowimy opr Wienera na[0,a]. Pr Wienera jako pr gaussowski. Mowimy,że wektor los(Y1, . . . , Yn) ma rozkład normalny, gdy zmienna suma od k=1 do n z aYk ma rozkład nor lub jest stałą dla dow liczb rzeczyw a1, . ,an. Mowimy,że proces (Yt) jest gaussowski, gdy wektor(Yt1, . . . , Ytn) ma rozkład normalny dla każdego rosnącegociąguchwil.DOWÓD:udowodnimy że pr Wien jest Gauss.Istotnie,przyrosty pr Wienera są NZL i mają rozkład nor.Suma NZL o rozkładzie nor ma rozkład nor, zatem wektor(Wt1 ,Wt2 -Wt1 , .,Wtn -Wtn.1) ma rozkład nor dla t1< . <tn.wektor (Wt1 ,Wt2 ,..,Wtn) ma rownież rozkład nor,gdyż Wtk = Wt1+(Wt2 -Wt1)+. +(Wtk -Wtk.1)dla dow k,zatem pr Wienerajest gaussowski. Def1’Pr Wienera naz prgaussowski (Wt)t>=0 o ciągłych realizacjach z W0=0,EWt=0 i Cov(Ws,Wt)= min(s,t) Tw(Paley-Wiener-Zygmund)P(A) = 0. Ruch BrownaWroćmy do opisu cząstki poruszającej się pod wpływem zderzeń z cząsteczkami cieczy,w której jest zanurzona.Gdy ciecz jest w ruchu, warto.oczekiwana jej położenia będzie funkcją czasu.Gdy temp cieczy zmienia się w czasie, wariancja przebytej drogi rownież będzie się zmieniała. Aby opisać tego typu procesy wprowadzimy pojęcie ruchu Browna z dryfem µ(t) i dyfuzją σ(t): Bt:=µ(t)+σ(t)Wt, t>=0,gdzie µ(t) i σ(t)> 0 są danymi funkcjami. Geometryczny ruch Browna Ct =C0exp{(µ-σ2/2)t+σWt, t>=0. gdzie Ct oznacza cenę w chwili t oraz liczby C0,µ,σ są dodatnie. logarytmiczne stopy zwrotu ln(Ct/Cs) są NZL i mają rozkład normalny. Mowimy,że pr (Yt), określony na (Ω,A,P), jest całkowalny,gdy E|Yt|<∞ dla dow t. Filtracją naz dow rosnącą rodzinę σ-ciał (Ft)⊂A,tzn. Fs⊂Ft dla s<t.Zbior Ft składa się z tych zdarzeń,o ktorych w chwili t wiemy czy zaszły,czy nie. Mowimy,że pr(Yt) jest zgodny z filtracją (Ft)(inaczej (Ft)zgodny, adaptowalny do (Ft)),gdy zm Yt jest Ft-mierzalna dla każdego t. Mowimy,że pr(Yt) jest martyngałem względem filtracji (Ft),gdy jest zgodny,całkowalny oraz dla t>s i A∈Fs E((Yt -Ys)1A)=0. Gdy pr (Yt) jest ciągiem,wystarczy przyjąć że dla dow n i A∈Fn,E((Yn+1-Yn)1A)=0.Istotnie,dla k>n E((Yk -Yn)1A) =suma od i=n do k-1 z E((Yi+1-Yi)1A)=0,Gdy A∈F ⊂Fi dla i>=n. Charakteryzacja martyngałów Tw Całkowalny i zgodny pr Yt jest martyngałem wtwg dla t>s E((Yt -Ys)Zs)=0, gdzie Zs jest dowolną Fs-mierzalną zm los taką,że (Yt-Ys)Zs ∈ L1. Włas o przyros nieskorelowanych Każdy martyngał (Yt)⊂L2 ma przyrosty nieskorelowane,tzn. dla t1< t2<=t3< t4 Cov(Yt2 -Yt1 ,Yt4 -Yt3)=0.Istotnie, wart oczekiwana martyngału jest funkcją stałą,zatem Cov(Yt2 -Yt1 ,Yt4 -Yt3) = E((Yt4 -Yt3)Z),gdzie Z = Yt2 -Yt1 oraz zm (Yt4 -Yt3 )Z jest całkowalna i Ft2-mierzalna(z nierowności Schwarza).

FunkcjęT :Ω→ {0,1,..}∪{∞} naz momentem Markowa, względem filtracji (Fn), gdy dla n<∞ {T=n}∈Fn. Mowimy,że moment Markowa jest mom stopu,gdy T<∞.Gdy T<=m dla stałej m<∞, funkcję T naz ograniczonym momentem stopu.

Tw(O losowym wyborze, Doob)

Załóżmy,że (Yn)n>=0

jestmartyngałem,T zaś jest mom stopu.Wówczas EYT=EY0, jeśli jeden z poniższych warunków jest spełniony:a) T jest ogr momstopu, b)supn>=0 |Yn|<=c, gdzie c jest liczbą, c)suma od n=1 do niesk z E(|Yn- Yn-1|1{T>=n} )< ∞,d) suma od n=1 do T z E|Yn-Yn-1| < ∞ i zm Yn-Yn-1 jest NZL od Fn-1 dla każ n. Pierwsza tożsamość Walda Niech T będzie mom stopu wzg naturalnej filtracji ciągu i.i.d. Xn.Jeśli TЄL1 to EST = ET EX. Druga tożsamość Walda Niech X,X1,X2,..będzie ciągiem i.i.d. z EX=0 i EX2<∞. Jeśli T∈L1 jest mom stopu wzgl naturalnej filtracji ciągu Xn, to EST2=ET EX2. Całkowalny i zgodny pr (Yt) naz submartyngałem wzg (Ft), gdy E((Yt-Ys)1A)>= 0 dla t>s i A∈Fs. Dekompozycja Dooba Mowimy,że ciąg (Kn)n>=0 jest prognozowalny,gdy zm Kn jest Fn-1mierzalna dla n>=1. Je.eli filtracja (Fn) jest generowana przez jakiś ciąg (Xn), to (Kn) jest prognozowany wtwg istnieją takie funkcje borelowskie gn,że

Kn=gn(X1,..,Xn-1) dla każd n. Dow submartyngał (Yn)n>=0 można jednoznacznie przedstawić w postaci Yn=Mn+Kn, gdzie Mn jest martyngałem, Kn zaś jest rosnącym i prognozowalnym ciągiem losowym z K0(w)=0 dla każd w.Rozkład Yn=Mn+Kn naz dekompozycją Dooba.O ciągu Kn mowimy,że jest kompensatorem, a ciąg Mn naz pr innowacji. MPWL tw Jeśli ciąg (Sn)n>=0⊂L2 z S0=0 jest martyngałem i suma od k=1 do niesk z ak2E(Sk-Sk-1)2<∞,gdzie an↓0, to anSn→0, gdy n→∞.

Tw Doob Jeżeli supn>=1 E(Yn)+<∞, to submartyngału Yn jest zbieżny prawie na pewno do całkowalnej zm los. Mowimy,że ciąg (Yn)n>=0 jest łańcuchem Markowa, gdy

(Yn+1|Y0,.., Yn)=(Yn+1|Yn) dla dow n, tzn. P(Yn+1=s|Y0=s0,.., Yn=sn)=P(Yn+1=s|Yn=sn)

dla takich s0,..,sn, s ∈ S i n>=1,że P(Y0=s0,..,Yn=sn)>0. Gdy zbior S jest skończony, mowimy,że łańcuch jest skończony. Łańcuch (Yn)n>=0 naz jednorodnym (wzg czasu), gdy rozkład Yn+1|Yn nie zależy od n.

Uogólniony schemat BernoulliegoNiech zm Yn przyjmuje tylko dwie wartości: 0 (”porażka”) i 1 (”sukces”) dla dow n. Jeżli Yn=i, to Yn+1=1 albo Yn+1=0 odpowiednio z pstwem pi i 1-pi dla i=0, 1. Liczby p0, p1 mogą zależeć od n.Ciąg (Yn) jest łańcuchem Markowa o przestrzeni stanów S={0,1}. Macierz przejścia ma postać P=(1-p0 p01-p1 p1)

Gdy p0=p1 i p1 nie zależy od n,otrzymujemy schematBernoulliego.

Rozkład łańcucha Markowa Wykażemy,że opis probabilistyczny jednorodnego łańcucha Markowa jest pełny,gdy podamy rozkład początkowy i macierz przejścia. Ze wzoru na pstwo warunkowe obliczamy rozkład wektora (Y0,Y1): P(Y0=i,Y1=j)=P(Y0=i)P(Y1=j|Y0=i)=P(Y0=i)pij .Rozkład wektora (Y0,Y1,Y2) ma postać:

P(Y0=i,Y1=j,Y2=k)=P(Y0=i)P(Y1=j|Y0=i)P(Y2=k|Y0=i,Y1=j)

=P(Y0=i)P(Y1=j|Y0=i)P(Y2=k|Y1=j)=P(Y0=i)·pij·pjk.W ostatnim przejściu skorzystaliśmy z jednorodności łańcucha. Postępując podobnie wyznaczamy rozkład wektora (Y0,Y1,.., Yn): P(Y0=i0,Y1=i1,..,Yn=in) =P(Y0=i0)·pi0i1 pi1i2·..·pin.1in,

gdy P(Y0=i0,Y1= i1,.., Yn=in)>0.

Rozkład zmiennej Yn Ozn rozkład zm Yn przez П(n): П(n)=(П0(n), П1(n),..),gdzie Пi(n)=P(Yn=i).Ze wzoru na pstwo całkowite (zu-pełne):

Пj(1)=P(Y1=j)=suma po i∈S z

P(Y1=j|Y0=i)P(Y0=i)=suma i∈S

pijПi(0), j∈S.Stąd wynika,że

П(1)=П(0)P.Powtarzając powyższe wyprowadzenie dla Yn otrzymamy wzor rekurencyjny:П(n)=П(n-1)P, n=2,3,..Stąd П(n)=П(0)Pn, n=1,2,..Stany chwilowe i powracająceNiech fij(n) oznacza pstwo przejścia ze stanu i do j po raz pierwszy po n kolejnych krokach: fij(n)=P(Y1≠j,..,Yn-1≠j,Yn=j|Y0=i) Pstwo fij ,że łańcuch wychodząc ze stanu i trafi kiedykolwiek do stanu j rowna się fij =suma od n=1 do niesk z fij(n).Stan i∈ S naz powracającym,gdy fii=1.Gdy

fii<1, stan i nazyw chwilowym.

Ozn przez pij(n) pstwo przejścia ze stanu i do j po n kolejnych krokach: pij(n)=P(Yn=j|Y0=i).Tw fii = 1 wtwg Pi=∞.Jeśli Pi<∞,to fii=Pi/(1+Pi).

ZBIEŻNOŚĆ DO ROZKŁADU STACJONARNEGO Jeśli (Xn) jest łańcuchem Mark, to rozkład П(n) zm Yn jest dany Wzor П(n)=П(0)Pn, ale wyznaczenie potęgi Pn jest zwykle trudnym zadaniem. Gdy istnieje granica limn→∞ Pn, istnieje rownież granica П=limn→∞ П(n); granica macierzy jest macierzą złożoną z granic. Ponadto,П=ПP, ponieważ П=limn→∞П(n+1)=limn→∞П(n)P=ПP.Jeżli rozw П=(П0,П1,..) rownania П=ПPjest rozkładem pstwa na przestrzeni stanow S, to można je uznać za przybliżenie dla П(n),gdy n jest duże. Rozkład П naz stacjonarnym (lub niezmien-

niczym), gdyż jeśli Y0 ma rozkład П,to Y1 ma rozkład ПP=П,a zatem Yn ma rozkład П dla każd n. Tw Markow Załóżmy,że istnieje stała c>0 oraz k-elementowy zbiór Sk⊂S taki,że pij>=c dla i∈S,j∈Sk. Wówczas istnieje tylko jeden rozkład stacjonarny П oraz dla dow rozkładu początkowego П(0) i n>=1: ||П(n)-П||<= 2/kc*(1-kc)n, gdzie П(n) oznacza rozkład zm Yn.Łańcuch Mar naz nieprzywiedlnym,jeśli wszystkie jego stany komunikują się tzn dla dow i,jЄS istnieje liczba n taka że pij(n)>0.Okresem stanu iЄS naz najwiekszy wspólny dzielnik zbioru tych n>=1 dla których pii(n)>0.Jeśli okres równa się 1to stan i nazywamy nieokresowym. Tw:Załóżmy, że łańcuch Yn o zbiorze stanów S jest nieprzywiedlny i nieokresowy.Jeśli zbiór S jest skończony, to istnieje n takie,

że mini,j∈S pij(n) >0,a zatem istnieje tylko jeden rozkład stacjonarny П i ||П(n)-П||→ 0,gdy n→∞, gdzie П(n) oznacza rozkład wektora Yn.

Mowimy że F-mierzalna zm los Y z daszkiem jest warunkową

wartością oczekiwaną całkowalnej zm Y wzg σ-ciała F, gdy dla każd A∈F E((Y-Yz daszk)1A)=0.Własności WWO

W1Liniowość. Jeśli a,b∈R,to

E(aX+bY|F)= aE(X|F)+bE(Y|F).W2Jeśli zm X jest F-mierzalna i XY ∈L1,to

E(XY|F)=XE(Y|F).W szczeg, E(X|F)=X, gdy X jest Fmierzalna. W3Wł iteracyjna E(E(Y|F))=EY. W4Jeżeli Y jest NZL od F (tzn zm Y i 1A są NZLdla dow A∈F),to E(Y|F)=EY . W5Monotoniczność Jeśli X<=Y to E(X|F)<=E(Y|F).W6Uogólniona wł iteracyjna Jeśli F1⊂F2⊂A,to

E(E(Y|F2)|F1)= E(E(Y|F1)|F2)=E(Y|F1). W7Nierówność Jensena Jeżeli g:R→R jest funkcją wypukłą oraz g(Y)∈L1,to E(g(Y)|F)>=g(E(Y|F)).Zastos wł WWO Całkowalny i zgodny z (Ft) pr (Yt) jest martyngałem, gdy dla t>s E(Yt|Fs)=Ys.

PROCESEM LOSOWYM (STOCHASTYCZNYM) nazyw dow rodzinę zm los {Yt:t Є T}.Wartością procesu los w chwili t nie jest liczba,ale zm los Yt czyli funk z Ω w R. Dla każ ustalonego w∈Ω , funkcję t → Yt(w),gdzie t ∈ T,nazyw realizacją (lub: trajektorią) procesu {Yt; t ∈ T}. Zbior S := {Yt(w); t ∈ T, w∈Ω } jest przestrzenią stanów. Ciąg (Xn) złożony z nzl zm los o tym samym rozkładzie (iid) nazyw szumem,gdy EX1=0 i VarX1<∞ i ozn εn.Błądzeniem losowym nazyw ciąg Yn=Yn-1+Xn, n=1,2,… gdzie Y0 jest dane.Jeśli dodatkowo Y0 jest liczba całkowitą i zmienne Xn przyjmują wartości -1 lub 1, to ciąg Yn nazyw błądzeniem los po zb liczb całkowitych.O procesie (Yt)t∈T mowimy,że ma przyrosty niezależne,gdy zmienne losowe Yt2 -Yt1 , . . . , Ytn –Ytn-1 są niezależne dla dowolnego rosnącego ciągu chwil (ti).Mowimy, że proces (Yt)t∈T ma przyrosty stacjonarne, gdynie zmieni się rozkład wektora (Yt2 . Yt1 , . . . , Ytn . Ytn.1 ) po przesunięciu dow rosnącego ciągu chwil (ti) o dow liczbę h tak że ti + h ∈ T dla każdego i, tzn(Yt2 -Yt1 , . . . , Ytn –Ytn-1 ) =(Yh+t2 -Yh+t1, . . . , Yh+tn . Yh+tn.1).O procesach (Xt)t∈T i (Yt)t∈T′ mowimy,że są niezależne, gdy wektory los (Xt1 , . . . ,Xtk ) i (Ys1, . . . , Ysn) są niezależne dla dow chwil (ti) ⊂ T, (sj) ⊂ T′.

Mowimy, że {Nt : t >=0} jest procesem Poissona z intensyw λ>0,gdy 1) N0 = 0, 2) proces ma przyrosty niezależne,

3) Nt -Ns ma rozkład Poissona z wartością oczekiwaną λ(t-s) dla t >s>=0, 4) realizacje są rosnące, prawostronnie ciągłe i schodkowe o skokach równych 1. Czas czekania na skoki Oznaczmy przez Tn chwilę pojawienia się n-tego skoku.Wykażemy, że czasy czekania na kolejne skoki: X1:=T1,X2 :=T2-T1,X3 :=T3-T2, . . . , tworzą ciąg i.i.d. o rozkładzie wykładniczym G(1,λ),tzn.P(X1 > t) = exp{-λt}. Oczywiście,dla t >0 P(X1 > t) = P(Nt = 0) = exp(-λt),zatem X1 ∼ G(1, λ). Dalej,P(X1 > s,X2 > t) = całka od0 doniesk z P(X1 > s,X2 > t|X1 = x)λexp{-λx}dx=całka od s do niesk z

P(brak zgłoszeń w (x,x + t)|X1 = x)λexp{-λx}dx= całka od s do ∞ z P(brak zgłoszeń w (x,x + t))λexp{-λx}dx=całka od s do ∞ z

Exp{-λt}λexp{-λx}dx = exp{-λt}exp{-λs}więc X1,X2 są niezależne i mają rozkład G(1,λ).Metoda indukcji kończy dowód,że Xn jest ciągiem i.i.d o rozkładzie wykładniczym. Konstrukcja procesu Poissona Niech (Xn)n>=1 będzie ciągiem i.i.d. o rozkładzie wykładniczym G(1,λ). Położmy Tn =X1+. . + Xn dla n>=1.Niech Nt będzie liczbą wskaźnikow n>=1,dla ktorych Tn<=t, czyli Nt =suma od n=1 do niesk z1(Tn<=t).(1) Tw.Proces(1)jest pr Poissona z intensywnością λ.Nietrywialne własności procesu Poissona.prawo wielkich liczb (MPWL)limt→∞ Nt/t=λ, centralne twierdzenie graniczne (CTG):(Ntλt)/ Pier(λt)→N(0,1) przy t→∞.


Wyszukiwarka

Podobne podstrony:
1 sciaga ppt
metro sciaga id 296943 Nieznany
ŚCIĄGA HYDROLOGIA
AM2(sciaga) kolos1 id 58845 Nieznany
Narodziny nowożytnego świata ściąga
finanse sciaga
Jak ściągać na maturze
Ściaga Jackowski
Aparatura sciaga mini
OKB SCIAGA id 334551 Nieznany
Przedstaw dylematy moralne władcy i władzy w literaturze wybranych epok Sciaga pl
fizyczna sciąga(1)
Finanse mala sciaga
Podział węży tłocznych ze względu na średnicę ściąga
OLIMPIADA BHP ŚCIĄGAWKA
Opracowanie Sciaga MC OMEN
Finanse Sciaga3 (str 7) id 171404

więcej podobnych podstron