UA MATH565C 随机微分方程II Wiener过程简介
UA MATH565C 隨機微分方程II Wiener過程簡介
- Wiener過程的簡單性質
- Wiener過程的定義
在上一講我們定義了WtW_tWt?:
dWt=ηtdt?Wt=∫0tηsdsdW_t = \eta_t dt \Leftrightarrow W_t = \int_{0}^t \eta_s dsdWt?=ηt?dt?Wt?=∫0t?ηs?ds
其中ηt\eta_tηt?是高斯白噪聲過程,WtW_tWt?被稱為Wiener過程或者Brown運動。這一part給出Wiener過程的構造理論,以及一些常用的性質與結論。
Wiener過程的簡單性質
因為Wiener過程根據高斯白噪聲構造,所以性質也主要由高斯白噪聲的性質推導得到。高斯白噪聲的性質為:
根據這三條,可以得到Wiener過程的三條簡單性質:
這一條是根據高斯白噪聲的性質2得到的,因為Wt?Ws=∫stηrdr,Wu?Wv=∫vuηrdrW_t-W_s=\int_s^t \eta_rdr,W_u-W_v=\int_v^u \eta_rdrWt??Ws?=∫st?ηr?dr,Wu??Wv?=∫vu?ηr?dr。
其實根據WtW_tWt?的定義還可以看出它的初值為0:
W0=∫00ηtdt=0W_0 = \int_0^0 \eta_t dt = 0W0?=∫00?ηt?dt=0
從Wiener過程的性質3得到,假設s<ts<ts<t,則
E[Ws(Wt?Ws)]=0?E[WsWt]=E[Ws2]=sE[W_s(W_t-W_s)]=0 \Rightarrow E[W_sW_t] = E[W_s^2] =sE[Ws?(Wt??Ws?)]=0?E[Ws?Wt?]=E[Ws2?]=s
這個結論可以寫成更一般的:
E[WsWt]=min(s,t),?s,t∈TE[W_sW_t]=min(s,t),\forall s,t \in \mathcal{T}E[Ws?Wt?]=min(s,t),?s,t∈T
其中WtW_tWt?的矩的計算可以用特征函數來做,這里給兩個結論:?k∈N\forall k \in \mathbb{N}?k∈N
因為Wt~N(0,t)W_t \sim N(0,t)Wt?~N(0,t),所以Wiener過程的矩可以直接套用正態分布的nnn階矩公式,公式的推導懶得打字,直接貼我的作業:
里面那個復合函數的高階導數公式可以參考柯召的組合論。
Wiener過程的定義
根據Wiener過程的簡單性質來定義它,如果隨機過程WtW_tWt?滿足下面三條性質:?s,t∈T\forall s,t \in \mathcal{T}?s,t∈T
比較正式的構造需要用Kolmogorov定理,為此WtW_tWt?的路徑函數依概率1連續,Wiener過程的構造放在下一講。
應用Wiener過程之前,需要根據這三條性質來驗證一個隨機過程是不是Wiener過程,下面舉幾個例子來說明。
如果WtW_tWt?是Wiener過程,則
例1
1、?Wt-W_t?Wt?是Wiener過程;
2、Xt=Wt0+t?WtX_t=W_{t_0+t}-W_tXt?=Wt0?+t??Wt?是Wiener過程,?t0∈T\forall t_0 \in \mathcal{T}?t0?∈T
第一條比較直觀,驗證一下第二條的性質3:
E[XsXt]=E[(Wt0+s?Ws)(Wt0+t?Wt)]=E[Wt0+sWt0+t]+E[WsWt]?E[Wt0+sWt]?E[Wt0+tWs]=min(t+t0,s+t0)+min(t,s)?min(t0+s,t)?min(t0+t,s)=min(t,s)E[X_sX_t] = E[(W_{t_0+s}-W_s)(W_{t_0+t}-W_t)] \\ = E[W_{t_0+s}W_{t_0+t}]+E[W_sW_t] - E[W_{t_0+s}W_t]-E[W_{t_0+t}W_s]\\=min(t+t_0,s+t_0)+min(t,s)-min(t_0+s,t)-min(t_0+t,s) =min(t,s) E[Xs?Xt?]=E[(Wt0?+s??Ws?)(Wt0?+t??Wt?)]=E[Wt0?+s?Wt0?+t?]+E[Ws?Wt?]?E[Wt0?+s?Wt?]?E[Wt0?+t?Ws?]=min(t+t0?,s+t0?)+min(t,s)?min(t0?+s,t)?min(t0?+t,s)=min(t,s)
例2 Yt=1cWc2tY_t=\frac{1}{c}W_{c^2t}Yt?=c1?Wc2t?是Wiener過程,這種構造叫Wiener Scaling
性質1、2都很顯然,驗證一下性質3:
E[YsYt]=E[1c2Wc2sWc2t]=1c2min(c2s,c2t)=min(s,t)E[Y_sY_t]=E[\frac{1}{c^2} W_{c^2s}W_{c^2 t}] = \frac{1}{c^2} min(c^2s,c^2t)=min(s,t)E[Ys?Yt?]=E[c21?Wc2s?Wc2t?]=c21?min(c2s,c2t)=min(s,t)
例3 Zt=tW1tZ_t=tW_{\frac{1}{t}}Zt?=tWt1??是Wiener過程
性質1、2同樣很顯然,先驗證一下性質3:
E[ZsZt]=E[stW1sW1t]=stmin(1s,1t)=min(s,t)E[Z_sZ_t] = E[stW_{\frac{1}{s}}W_{\frac{1}{t}}] = st min(\frac{1}{s},\frac{1}{t}) = min(s,t)E[Zs?Zt?]=E[stWs1??Wt1??]=stmin(s1?,t1?)=min(s,t)
需要注意的是這個隨機過程的時間在WtW_tWt?的基礎上做了一個反比例變換,導致t=0t=0t=0成為了一個奇點。因此需要驗證一下這個隨機過程的路徑函數是否以概率1連續,更精確一點,要驗證ZtZ_tZt?在t=0t=0t=0處能以概率1連續。
定義E={w:Zt→0,ast→0}E = \{w:Z_t \to 0, as\ t \to 0\}E={w:Zt?→0,as?t→0},用一個比較神奇的構造:
E=?n?m?q∈Q∩[0,1m]{w:∣Zq∣≤1n}E=\bigcap_n \bigcup_m \bigcap_{q \in \mathbb{Q}\cap [0,\frac{1}{m}]}\{w:|Z_q| \le \frac{1}{n}\}E=n??m??q∈Q∩[0,m1?]??{w:∣Zq?∣≤n1?}
要證明以概率1連續,需要
P(E)=1P(E)=1P(E)=1
因為Zq=W1q1/qZ_q = \frac{W_{\frac{1}{q}}}{1/q}Zq?=1/qWq1???,當q→0q \to 0q→0,1q→∞\frac{1}{q} \to \inftyq1?→∞,根據大數法則Zq→0a.s.Z_q \to 0\ a.s.Zq?→0?a.s.,用這個來證明P(E)=1P(E)=1P(E)=1就很直接了。再來看一下這個神奇的構造,要證明ZtZ_tZt?在t=0t=0t=0處連續,需要在T\mathcal{T}T上找一個收斂到0的子列,滿足按這個子列對ZtZ_tZt?的采樣依概率1收斂到0。因為ZtZ_tZt?中的時間在已知的Wiener過程中是1/t1/t1/t,所以這里用q=Q∩[0,1m],m=1,?,∞q=\mathbb{Q} \cap [0,\frac{1}{m}],m=1,\cdots,\inftyq=Q∩[0,m1?],m=1,?,∞是在確定一個收斂到0有理子列。∣Zq∣≤1n|Z_q| \le \frac{1}{n}∣Zq?∣≤n1?并對n=1,?,∞n=1,\cdots,\inftyn=1,?,∞取交集,ZqZ_qZq?其實就是按那個子列對ZtZ_tZt?做的一列采樣,如果這個東西的概率是1,就可以說明這個子列對ZtZ_tZt?的采樣依概率1收斂到0。
總結
以上是生活随笔為你收集整理的UA MATH565C 随机微分方程II Wiener过程简介的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH571B 试验设计V 析因
- 下一篇: UA MATH565C 随机微分方程II