日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

UA MATH565C 随机微分方程II Wiener过程简介

發布時間:2025/4/14 编程问答 46 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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過程根據高斯白噪聲構造,所以性質也主要由高斯白噪聲的性質推導得到。高斯白噪聲的性質為:

  • Eηt=0,?t∈TE\eta_t=0,\forall t \in \mathcal{T}Eηt?=0,?tT
  • Eηtηs=δ(t?s),?t,s∈T,σ2<+∞E\eta_t\eta_s=\delta(t-s),\forall t,s \in \mathcal{T},\sigma^2 < +\inftyEηt?ηs?=δ(t?s),?t,sT,σ2<+
  • 高斯過程
  • 根據這三條,可以得到Wiener過程的三條簡單性質:

  • 高斯過程
  • EWt=E∫0tηsds=∫0tEηsds=0EW_t = E\int_{0}^t \eta_s ds = \int_{0}^t E\eta_s ds=0EWt?=E0t?ηs?ds=0t?Eηs?ds=0
  • ?0≤s<t<v<u,E[(Wt?Ws)(Wu?Wv)]=0\forall 0\le s<t<v<u,E[(W_t-W_s)(W_u-W_v)]=0?0s<t<v<u,E[(Wt??Ws?)(Wu??Wv?)]=0
    這一條是根據高斯白噪聲的性質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,tT
    其中WtW_tWt?的矩的計算可以用特征函數來做,這里給兩個結論:?k∈N\forall k \in \mathbb{N}?kN

  • E[Wt]2k=(2k)!2kk!tkE[W_t]^{2k} = \frac{(2k)!}{2^k k!} t^kE[Wt?]2k=2kk!(2k)!?tk
  • E[Wt]2k?1=0E[W_t]^{2k-1} = 0E[Wt?]2k?1=0
  • 因為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,tT

  • 高斯過程;
  • EWt=0EW_t=0EWt?=0;
  • E[WtWs]=min(s,t)E[W_tW_s]=min(s,t)E[Wt?Ws?]=min(s,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?t0},用一個比較神奇的構造:
    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??qQ[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 0q01q→∞\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过程简介的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。