UA MATH565C 随机微分方程III Ito积分简介
UA MATH565C 隨機(jī)微分方程III Ito積分簡(jiǎn)介
- Wiener過(guò)程的分割
- 與Riemann積分的對(duì)比
- Ito積分的構(gòu)造
在隨機(jī)微分方程解的構(gòu)造中,積分
∫0tσ(Xs)dWs\int_{0}^t \sigma(X_s)dW_s∫0t?σ(Xs?)dWs?
被叫做Ito積分,這一part的目標(biāo)是構(gòu)造Ito積分并給出一些Ito積分常用的性質(zhì)。
Riemann積分和Lebesgue積分的構(gòu)造是大家耳熟能詳?shù)?#xff0c;Riemann積分的構(gòu)造中最重要的思想是積分就是基于分割定義的部分和的極限;Lebesgue積分的構(gòu)造中主要的思路是先定義簡(jiǎn)單可測(cè)函數(shù)的積分,再用簡(jiǎn)單可測(cè)函數(shù)序列逼近一般可測(cè)函數(shù)。把這些思想用在Ito積分的構(gòu)造中,首先我們研究分割和部分和的性質(zhì),然后定義簡(jiǎn)單隨機(jī)過(guò)程的積分,再用簡(jiǎn)單隨機(jī)過(guò)程的序列逼近一般隨機(jī)過(guò)程。
Wiener過(guò)程的分割
給定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?]
這是時(shí)域上的一個(gè)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→∞
證明
等價(jià)敘述是
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
根據(jù)Wiener過(guò)程的性質(zhì),不同的增量ΔjWt\Delta_j W_tΔj?Wt?是獨(dú)立的。定義兩個(gè)符號(hào)
Δ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?
在構(gòu)造積分取部分和的極限時(shí),求的是部分和在分割的模趨于0的時(shí)候部分和的極限。Riemann積分中分割的模一般用所有分割的Lebesgue測(cè)度的上確界。但構(gòu)造Ito積分的時(shí)候分割ΔjWt\Delta_j W_tΔj?Wt?是隨機(jī)過(guò)程,當(dāng)nnn趨近于無(wú)窮時(shí),ΔjWt→L20\Delta_j W_t \to_{L^2} 0Δj?Wt?→L2?0,感覺(jué)上符合分割的模趨于零這個(gè)想法的。
∑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時(shí),ZjZ_jZj?與ZkZ_kZk?是獨(dú)立的,因此
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?
計(jì)算
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過(guò)程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
帶入計(jì)算
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→∞
之所以需要這個(gè)定理是因?yàn)樵?span id="ozvdkddzhkzd" class="katex--inline">[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??)之后,要保證當(dāng)nnn趨近于無(wú)窮的時(shí)候,分割的方差要趨近于ttt。大于ttt意味著做了分割之后隨機(jī)過(guò)程包含的信息反而增多了,顯然是不合理的;小于ttt說(shuō)明這個(gè)分割會(huì)造成信息損失。
與Riemann積分的對(duì)比
我們可以考慮一下與∫0tσ(Xs)dWs\int_{0}^t \sigma(X_s)dW_s∫0t?σ(Xs?)dWs?相似的Riemann積分:假設(shè)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積分對(duì)比的話需要考慮分割的平方和:
∑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過(guò)程的分割結(jié)論是:
∑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積分的分割相比主要有兩點(diǎn)區(qū)別:
如果只有第二個(gè)區(qū)別還可以直接把Riemann積分的構(gòu)造套用在Ito積分上,但第一個(gè)區(qū)別導(dǎo)致這種生搬硬套注定要發(fā)散,所以需要另外的構(gòu)造方式。
Ito積分的構(gòu)造
用一個(gè)簡(jiǎn)單的例子,∫0tWsdWs\int_{0}^t W_sdW_s∫0t?Ws?dWs?來(lái)說(shuō)明Ito積分的部分和與Riemann積分和Lebesgue積分的區(qū)別。
Riemann積分的部分和是對(duì)每個(gè)分割的模與分割內(nèi)選取的任一點(diǎn)函數(shù)值的乘積求和;簡(jiǎn)單函數(shù)的Lebesgue積分在每個(gè)分割內(nèi)的函數(shù)值是常數(shù),不存在選取的問(wèn)題。Ito積分的部分和應(yīng)該在分割內(nèi)選取哪一點(diǎn)是一個(gè)值得討論一下的問(wèn)題,這里給出兩種可能的方案:
計(jì)算
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?
這兩種不同的選取導(dǎo)致部分和的極限完全不同,說(shuō)明Ito積分的部分和不能選取分割內(nèi)的任意點(diǎn)來(lái)構(gòu)造。一個(gè)更關(guān)鍵的問(wèn)題是應(yīng)該選取哪個(gè)點(diǎn)才是合理的。
我們?cè)賮?lái)審視一下Ito積分與Riemann積分的區(qū)別。Riemann積分中用來(lái)作分割的函數(shù)ggg原象空間RRR就表示是一個(gè)參數(shù)空間,其中的點(diǎn)是沒(méi)有先后順序的。Ito積分中用來(lái)作分割的Wiener過(guò)程的原象空間是(Ω×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,∞)。事實(shí)上這并不是一個(gè)合理的定義,因?yàn)樵谶@個(gè)乘積可測(cè)空間中時(shí)間集合T\mathcal{T}T的σ\sigmaσ代數(shù)用的是Borel σ\sigmaσ代數(shù),這意味著T\mathcal{T}T的子集在映射中會(huì)像ggg的參數(shù)那樣被同等地對(duì)待。但是關(guān)于隨機(jī)過(guò)程WtW_tWt?,一個(gè)常識(shí)是ttt比較小的WtW_tWt?值一定比ttt比較大的先出現(xiàn)。因此T\mathcal{T}T中的元素與[0,∞)[0,\infty)[0,∞)有本質(zhì)區(qū)別。
考慮到這個(gè)區(qū)別,Ito積分做了一個(gè)假設(shè):Ito積分的被積函數(shù)不考慮當(dāng)前時(shí)點(diǎn)以后的情況。也就是說(shuō)在考慮分割ΔjWt\Delta_j W_tΔj?Wt?的時(shí)候,我們永遠(yuǎn)假設(shè)只知道WtjW_{t_j}Wtj??的行為,而tjt_jtj?以后的行為是隨機(jī)的。并且還沒(méi)發(fā)生,所以部分和按InI_nIn?的方式來(lái)定義:
假設(shè)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?
顯然這個(gè)積分就和Lebesgue積分比較相似了,下一講用和Lebesgue積分的構(gòu)造類似的思路,先考慮step processes的Ito積分,再用step processes序列逼近一般隨機(jī)過(guò)程。但做這個(gè)之前還需要考慮一下對(duì)時(shí)間集T\mathcal{T}T的σ\sigmaσ代數(shù)的重新定義。上面討論了用Borel σ\sigmaσ代數(shù)是不行的,所以下一講還要給T\mathcal{T}T重新定義一個(gè)σ\sigmaσ代數(shù)。
總結(jié)
以上是生活随笔為你收集整理的UA MATH565C 随机微分方程III Ito积分简介的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: UA MATH565C 随机微分方程II
- 下一篇: UA MATH564 概率论VI 数理统