UA MATH565C 随机微分方程III Ito积分简介
UA MATH565C 隨機微分方程III Ito積分簡介
- Wiener過程的分割
- 與Riemann積分的對比
- Ito積分的構造
在隨機微分方程解的構造中,積分
∫0tσ(Xs)dWs\int_{0}^t \sigma(X_s)dW_s∫0t?σ(Xs?)dWs?
被叫做Ito積分,這一part的目標是構造Ito積分并給出一些Ito積分常用的性質。
Riemann積分和Lebesgue積分的構造是大家耳熟能詳的,Riemann積分的構造中最重要的思想是積分就是基于分割定義的部分和的極限;Lebesgue積分的構造中主要的思路是先定義簡單可測函數的積分,再用簡單可測函數序列逼近一般可測函數。把這些思想用在Ito積分的構造中,首先我們研究分割和部分和的性質,然后定義簡單隨機過程的積分,再用簡單隨機過程的序列逼近一般隨機過程。
Wiener過程的分割
給定t∈Tt \in \mathcal{T}t∈T,定義tj=jnt,j=0,?,n?1t_j = \frac{j}{n}t,j=0,\cdots,n-1tj?=nj?t,j=0,?,n?1,則
[0,t]=[t0,t1]?(t1,t2]???(tn?1,tn][0,t]=[t_0,t_1] \sqcup (t_1,t_2] \sqcup \cdots \sqcup (t_{n-1},t_n][0,t]=[t0?,t1?]?(t1?,t2?]???(tn?1?,tn?]
這是時域上的一個equal-partition。從而
Wt=∑j=0n?1(Wtj+1?Wtj)W_t = \sum_{j=0}^{n-1} (W_{t_{j+1}}-W_{t_j})Wt?=j=0∑n?1?(Wtj+1???Wtj??)
定理 ∑j=0n?1(Wtj+1?Wtj)2→L2t,asn→∞\sum_{j=0}^{n-1} (W_{t_{j+1}}-W_{t_j})^2 \to_{L^2}t,as\ n\to \inftyj=0∑n?1?(Wtj+1???Wtj??)2→L2?t,as?n→∞
證明
等價敘述是
E[∑j=0n?1(Wtj+1?Wtj)2?t]2→0E\left[ \sum_{j=0}^{n-1} (W_{t_{j+1}}-W_{t_j})^2 - t \right]^2 \to 0E[j=0∑n?1?(Wtj+1???Wtj??)2?t]2→0
根據Wiener過程的性質,不同的增量ΔjWt\Delta_j W_tΔj?Wt?是獨立的。定義兩個符號
ΔjWt=Wtj+1?Wtj~N(0,Δjt)Δjt=tj+1?tj=tn\Delta_j W_t = W_{t_{j+1} }-W_{t_j} \sim N(0,\Delta_j t) \\ \Delta_j t = t_{j+1}-t_j = \frac{t}{n}Δj?Wt?=Wtj+1???Wtj??~N(0,Δj?t)Δj?t=tj+1??tj?=nt?
在構造積分取部分和的極限時,求的是部分和在分割的模趨于0的時候部分和的極限。Riemann積分中分割的模一般用所有分割的Lebesgue測度的上確界。但構造Ito積分的時候分割ΔjWt\Delta_j W_tΔj?Wt?是隨機過程,當nnn趨近于無窮時,ΔjWt→L20\Delta_j W_t \to_{L^2} 0Δj?Wt?→L2?0,感覺上符合分割的模趨于零這個想法的。
∑j=0n?1(Wtj+1?Wtj)2?t=∑j=0n?1[(ΔjWt)2?Δjt]\sum_{j=0}^{n-1} (W_{t_{j+1}}-W_{t_j})^2 - t = \sum_{j=0}^{n-1} [(\Delta_j W_t)^2 - \Delta_j t]j=0∑n?1?(Wtj+1???Wtj??)2?t=j=0∑n?1?[(Δj?Wt?)2?Δj?t]
定義Zj=(ΔjWt)2?ΔjtZ_j = (\Delta_j W_t)^2 - \Delta_j tZj?=(Δj?Wt?)2?Δj?t,顯然EZj=0EZ_j=0EZj?=0。則
(∑j=0n?1Zj)2=∑j=0n?1Zj2+∑j≠kZjZk(\sum_{j=0}^{n-1} Z_j)^2 = \sum_{j=0}^{n-1}Z_j^2 + \sum_{j \ne k} Z_jZ_k(j=0∑n?1?Zj?)2=j=0∑n?1?Zj2?+j?=k∑?Zj?Zk?
注意到j≠kj \ne kj?=k時,ZjZ_jZj?與ZkZ_kZk?是獨立的,因此
E(∑j=0n?1Zj)2=E∑j=0n?1Zj2+∑j≠kEZjEZk=∑j=0n?1EZj2E(\sum_{j=0}^{n-1} Z_j)^2 = E\sum_{j=0}^{n-1}Z_j^2 + \sum_{j \ne k} EZ_jEZ_k = \sum_{j=0}^{n-1}EZ_j^2E(j=0∑n?1?Zj?)2=Ej=0∑n?1?Zj2?+j?=k∑?EZj?EZk?=j=0∑n?1?EZj2?
計算
Zj2=[(ΔjWt)2?Δjt]2=(ΔjWt)4?2(ΔjWt)2Δjt+(Δjt)2Z_j^2 = [(\Delta_j W_t)^2 - \Delta_j t]^2 = (\Delta_j W_t)^4 - 2(\Delta_j W_t)^2\Delta_j t + (\Delta_j t)^2Zj2?=[(Δj?Wt?)2?Δj?t]2=(Δj?Wt?)4?2(Δj?Wt?)2Δj?t+(Δj?t)2
上一講給出了求Wiener過程nnn階矩的公式,直接套用可以得到
E(ΔjWt)4=3(Δjt)2,E(ΔjWt)2=(Δjt)2E(\Delta_j W_t)^4 = 3 (\Delta_j t)^2,E(\Delta_j W_t)^2 = (\Delta_j t)^2E(Δj?Wt?)4=3(Δj?t)2,E(Δj?Wt?)2=(Δj?t)2
帶入計算
E(∑j=0n?1Zj)2=∑j=0n?12(Δjt)2=2t2n→0,asn→∞E(\sum_{j=0}^{n-1} Z_j)^2 =\sum_{j=0}^{n-1} 2(\Delta_j t)^2 = \frac{2t^2}{n} \to 0, as\ n \to \inftyE(j=0∑n?1?Zj?)2=j=0∑n?1?2(Δj?t)2=n2t2?→0,as?n→∞
之所以需要這個定理是因為在[0,t][0,t][0,t]上WtW_tWt?的方差為ttt,作了分割∑j=0n?1(Wtj+1?Wtj)\sum_{j=0}^{n-1} (W_{t_{j+1}}-W_{t_j})∑j=0n?1?(Wtj+1???Wtj??)之后,要保證當nnn趨近于無窮的時候,分割的方差要趨近于ttt。大于ttt意味著做了分割之后隨機過程包含的信息反而增多了,顯然是不合理的;小于ttt說明這個分割會造成信息損失。
與Riemann積分的對比
我們可以考慮一下與∫0tσ(Xs)dWs\int_{0}^t \sigma(X_s)dW_s∫0t?σ(Xs?)dWs?相似的Riemann積分:假設g∈C1(R)g \in C^1(\mathbb{R})g∈C1(R)
∫0tf(s)dg(s)\int_0^t f(s)dg(s)∫0t?f(s)dg(s)
同樣用以下分割
[0,t]=[t0,t1]?(t1,t2]???(tn?1,tn][0,t]=[t_0,t_1] \sqcup (t_1,t_2] \sqcup \cdots \sqcup (t_{n-1},t_n][0,t]=[t0?,t1?]?(t1?,t2?]???(tn?1?,tn?]
定義Δjg(t)=g(tj+1)?g(tj)\Delta_j g(t) = g(t_{j+1}) - g(t_j)Δj?g(t)=g(tj+1?)?g(tj?),要與Ito積分對比的話需要考慮分割的平方和:
∑j=0n?1(Δjg)2=∑j=0n?1(g(tj+1)?g(tj))2\sum_{j=0}^{n-1} (\Delta_j g)^2 = \sum_{j=0}^{n-1} (g(t_{j+1}) - g(t_j))^2j=0∑n?1?(Δj?g)2=j=0∑n?1?(g(tj+1?)?g(tj?))2
用Lagrange中值定理,?tj?∈(tj,tj+1)\exists t_j^* \in (t_j,t_{j+1})?tj??∈(tj?,tj+1?)
∑j=0n?1(g(tj+1)?g(tj))2=∑j=0n?1(g′(tj?)Δjt)2≤[sup?jg′(tj?)]2∑j=0n?1(Δjt)2=[sup?jg′(tj?)]2t2n→0\sum_{j=0}^{n-1} (g(t_{j+1}) - g(t_j))^2 = \sum_{j=0}^{n-1} (g'(t_{j}^*) \Delta_j t)^2 \\ \le [\sup_j g'(t_j^*)]^2 \sum_{j=0}^{n-1} (\Delta_j t)^2 = \frac{[\sup_j g'(t_j^*)]^2 t^2}{n} \to 0j=0∑n?1?(g(tj+1?)?g(tj?))2=j=0∑n?1?(g′(tj??)Δj?t)2≤[jsup?g′(tj??)]2j=0∑n?1?(Δj?t)2=n[supj?g′(tj??)]2t2?→0
Wiener過程的分割結論是:
∑j=0n?1[(Wtj+1?Wtj)2?Δjt]→L20\sum_{j=0}^{n-1} [(W_{t_{j+1}}-W_{t_j})^2-\Delta_j t] \to_{L^2}0j=0∑n?1?[(Wtj+1???Wtj??)2?Δj?t]→L2?0
和Riemann積分的分割相比主要有兩點區別:
如果只有第二個區別還可以直接把Riemann積分的構造套用在Ito積分上,但第一個區別導致這種生搬硬套注定要發散,所以需要另外的構造方式。
Ito積分的構造
用一個簡單的例子,∫0tWsdWs\int_{0}^t W_sdW_s∫0t?Ws?dWs?來說明Ito積分的部分和與Riemann積分和Lebesgue積分的區別。
Riemann積分的部分和是對每個分割的模與分割內選取的任一點函數值的乘積求和;簡單函數的Lebesgue積分在每個分割內的函數值是常數,不存在選取的問題。Ito積分的部分和應該在分割內選取哪一點是一個值得討論一下的問題,這里給出兩種可能的方案:
計算
In+Jn=∑j=0n?1(Wtj+Wtj+1)ΔjWt=∑j=1n?1(Wtj+i2?Wtj2)=Wt2Jn?In=∑j=0n?1(ΔjWt)2→L2tI_n+J_n = \sum_{j=0}^{n-1} (W_{t_j}+W_{t_{j+1}})\Delta_j W_t = \sum_{j=1}^{n-1}(W_{t_{j+i}}^2 - W_{t_j}^2) = W_t^2 \\ J_n - I_n = \sum_{j=0}^{n-1} (\Delta_j W_t)^2 \to_{L^2} tIn?+Jn?=j=0∑n?1?(Wtj??+Wtj+1??)Δj?Wt?=j=1∑n?1?(Wtj+i?2??Wtj?2?)=Wt2?Jn??In?=j=0∑n?1?(Δj?Wt?)2→L2?t
所以
In→L212Wt2?t2Jn→L212Wt2+t2I_n \to_{L^2} \frac{1}{2}W_t^2 - \frac{t}{2}\\ \ J_n\to_{L^2} \frac{1}{2}W_t^2 + \frac{t}{2}In?→L2?21?Wt2??2t??Jn?→L2?21?Wt2?+2t?
這兩種不同的選取導致部分和的極限完全不同,說明Ito積分的部分和不能選取分割內的任意點來構造。一個更關鍵的問題是應該選取哪個點才是合理的。
我們再來審視一下Ito積分與Riemann積分的區別。Riemann積分中用來作分割的函數ggg原象空間RRR就表示是一個參數空間,其中的點是沒有先后順序的。Ito積分中用來作分割的Wiener過程的原象空間是(Ω×T,F?B(T),P?λ)(\Omega \times \mathcal{T},\mathcal{F} \otimes \mathcal{B}(\mathcal{T}),P\otimes\lambda)(Ω×T,F?B(T),P?λ),T=[0,∞)\mathcal{T}=[0,\infty)T=[0,∞)。事實上這并不是一個合理的定義,因為在這個乘積可測空間中時間集合T\mathcal{T}T的σ\sigmaσ代數用的是Borel σ\sigmaσ代數,這意味著T\mathcal{T}T的子集在映射中會像ggg的參數那樣被同等地對待。但是關于隨機過程WtW_tWt?,一個常識是ttt比較小的WtW_tWt?值一定比ttt比較大的先出現。因此T\mathcal{T}T中的元素與[0,∞)[0,\infty)[0,∞)有本質區別。
考慮到這個區別,Ito積分做了一個假設:Ito積分的被積函數不考慮當前時點以后的情況。也就是說在考慮分割ΔjWt\Delta_j W_tΔj?Wt?的時候,我們永遠假設只知道WtjW_{t_j}Wtj??的行為,而tjt_jtj?以后的行為是隨機的。并且還沒發生,所以部分和按InI_nIn?的方式來定義:
假設f:(Ω×T,F?B(T),P?λ)→(Rm,B(Rm))f:(\Omega \times \mathcal{T},\mathcal{F} \otimes \mathcal{B}(\mathcal{T}),P\otimes\lambda) \to (\mathbb{R}^m,\mathcal{B}(\mathbb{R}^m))f:(Ω×T,F?B(T),P?λ)→(Rm,B(Rm)),
∫0tfsdWs=lim?n→∞∑j=0n?1ftjΔjWt\int_0^t f_s dW_s = \lim_{n \to \infty} \sum_{j=0}^{n-1} f_{t_j}\Delta_j W_t∫0t?fs?dWs?=n→∞lim?j=0∑n?1?ftj??Δj?Wt?
顯然這個積分就和Lebesgue積分比較相似了,下一講用和Lebesgue積分的構造類似的思路,先考慮step processes的Ito積分,再用step processes序列逼近一般隨機過程。但做這個之前還需要考慮一下對時間集T\mathcal{T}T的σ\sigmaσ代數的重新定義。上面討論了用Borel σ\sigmaσ代數是不行的,所以下一講還要給T\mathcal{T}T重新定義一個σ\sigmaσ代數。
總結
以上是生活随笔為你收集整理的UA MATH565C 随机微分方程III Ito积分简介的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH565C 随机微分方程II
- 下一篇: UA MATH564 概率论VI 数理统