Loading [MathJax]/jax/output/CommonHTML/jax.js

3/01/2014

[隨機分析] Ito Integral 淺談 (I) - Ito Integral 的建構與 Ito Isometry property

這次要介紹隨機分析中的 Ito  integral 的建構:

目標:建立下面的(隨機)積分 or  Ito integral 
I(f)(ω)=T0f(ω,t)dBt    ()
其中 T 為固定時間, Bt 是標準布朗運動(Standard Brownian motion)。 f(ω,t)是一個隨機過程。
I(f)(ω) 表示積分為一個 mapping (之後會定義該從哪邊mapping到哪邊),且積分完畢之後會是一個隨機變數 (function of ω)

Comment:
你可能會問上面的積分跟一般積分有何不同!?
第一 是積分變數 f(ω,t) 不再是定數。此時的積分變數為一個隨機過程。
第二是 後方積分對象 dBt 亦為一個隨機過程 (標準布朗運動)。此時會使原本的Riemann -Stieltjes 積分無法定義(因為寫成sum之後左端點與右端點的答案不同)。 ( Ito Integral 選擇左端點因為之後會有較好的性質 ( Ito integral is (Local) martingale. )不過這是後話。)

建構 Ito 積分的想法如下:
先透過一類簡單的函數定義出上面的積分。再將其定義域擴展到更廣的函數類別。


積分變數(Integrand)需要那些條件?
為了要讓上述的積分可以make sense, () 式子中的積分變數 f(ω,t) 必須先滿足一些基本的可測性 (measurability) 與可積分 (integrability) 的條件。

首先考慮可測性(measurability):
B:= the smallest σ-algebra that contains all of the open subsets of [0,T]
{Ft}:= be standard Brownian filtration
且對所有的 t0Ft×B:= the smallest σ-algebra that contains all of the product sets A×B;其中 AFtBB

最後我們說 f(,) 是 可測的 (measurable) 如果 f(,)FT×B
f(,) 是 adapted 如果 對所有的 t[0,T], f(,t)Ft


再者考慮積分性(integrability):
為了建構前述 Ito 積分 (),我們會先將注意力放在 積分變數 f(w,t) 是落在 H2[0,T] 空間 ( H2 is a closed linear subspace of L2(dP×dt) )。其中 H2[0,T] 表示是由所有 measurable + adapted 函數 f組成的集合,且滿足下列積分限制
E[T0f(ω,t)2dt]<

Comments:
1. 上式積分限制為 L2(dP×dt)-norm ;
2. 上式積分限制為 一個雙重積分 ( 期望值由定義為 Lebesgue 積分;e.g.,  E[X]:=ΩXdP  ) 注意到 dP 指的是 Lebesgue 積分 對應於上式積分限制中的期望值部分, dt 指的上式積分限制中  Rieman 積分的部分,亦即對應於上式的  T0f(ω,t)2dt


現在我們可以開始定義積分。First step
回憶我們的idea:先建構 對基本函數的積分 再進一步拓展。所以在此我們先定義某一類基本函數所在的集合如下:

定義集合: H20H2的子集合,且此集合 H20 由下列基本 step 函數構成
f(ω,t):=n1i=0ai(ω)1(ti<tti+1)    (1)其中 aiFti, E[a2i]< , 1() 為 indicator function,且定義 partition
0=t0<t1<...<tn1<tn=T NOTE: 關於 Indicator function 1()有興趣的讀者請參閱 [機率論] Indicator function 

那麼,對於上式 (1) ,我們定義在 H20 上的 Ito Integral 為
I(f)(ω):=n1i=0ai(ω){Bti+1Bti}    (2)
Comment:
1. 注意到上述積分 I(f)(ω), 並未寫成一般常用的積分符號 是因為要強調此積分為某種mapping,且此積分目前只建構在基本step函數所構成的空間 H20,尚未值得引入標準的積分符號。
2. 上述積分的 ai(ω)是取左端點 i 的值 (Ito integral)

接著,我們要開始延伸積分的定義域從 H20 空間 到整個 H2 空間。
再延伸定義域之前,我們必須先確認 I:H20L2(dP×dt)L2(dP) 是一個 continuous mapping。下列引理建立了這個continuous mapping的關係

----------------------------
LEMMA 1 (Ito Isometry on H20)
對所有的 fH20,我們有如下關係
||I(f)||L2(dP)=||f||L2(dP×dt)-----------------------
Proof
idea: 分別計算左右兩邊的 L2 norm;為了方便起見我們計算 L2-norm的平方,;

給定 fH20; i.e., f(ω,t):=n1i=0ai(ω)1(ti<tti+1)

先計算右邊的 ||f||L2(dP×dt),由 L2-norm 定義:
||f||L2(dP×dt):=(E[T0f(ω,t)2dt])1/2也就是 ||f||2L2(dP×dt);照定義計算:
f2(ω,t)=(n1i=0ai(ω)1(ti<tti+1))2f2(ω,t)=n1i=0a2i(ω)1(ti<tti+1)

||f||2L2(dP×dt)=E[T0f(ω,t)2dt]=E[T0;n1i=0a2i(ω)1(ti<tti+1)dt]

f2(ω,t)=E[ti+1tin1i=0a2i(ω)dt]

f2(ω,t)=E[n1i=0a2i(ω)(ti+1ti)]

f2(ω,t)=n1i=0E[a2i(ω)](ti+1ti))    ()


接著再計算左邊的 ||I(f)||L2(dP),由定義可知:
||I(f)||L2(dP):=E[I2(f)(ω)]1/2

由於我們對基本step函數定義的積分為
I(f)(ω):=n1i=0ai(ω){Bti+1Bti}
I(f)2(ω)=(n1i=0ai(ω){Bti+1Bti})2

=n1i=0a2i(ω){Bti+1Bti}2+n1ijai(ω)aj(ω){Bti+1Bti}{Btj+1Btj}

兩邊同取期望值 + 利用 Brownian motion 定義 (independent increment 上式的交叉項取期望值 =0) + aiFti可得

E[I(f)2(ω)]=n1i=0E[a2i(ω){Bti+1Bti}2]
=n1i=0E[a2i(ω)]E[{Bti+1Bti}2]=n1i=0E[a2i(ω)](ti+1ti)    ()

現在比較 可知兩式相等,故得証。

------

現在我們得知 I(f)(ω) 確實從 H20 映射到 L2(dP),且因為 I(f)(ω) 在映射過程保持 了 norm  (isometry property),這暗示了 I(f)(ω) 亦將 Cauchy sequence in H20 映射到 a Cauchy sequence in L2(dP)

故下面的引理說明了我們可以用在 H20 中的元素來逼近 落在 H2中的元素。

--------
LEMMA 2 ( H20  is Dense in H2)
對任意 fH2L2(dP×dt),存在一個序列 {fn}fnH20  使得當  n
||ffn||L2(dP×dt)0---
Proof: omitted.

至此我們便完成Ito 積分的第一步,用 simple step function H20 逼近(取極限) 落在 H2的元素,並且將此極限稱為 Ito integral in H2

另外由於現在 Ito 積分已經被定義在 H2 空間,之前的 Ito isometry 也被拓展到 H2,此為Ito 積分建構中的重要結果,我們將其寫為下面的定理:

--------
Theorem 3 Ito Isometry in H2)
對任意 fH2L2(dP×dt),我們有下列的 Isometry 性質
||I(f)(ω)||L2(dP)=||f||L2(dP×dt)---
Proof:
因為定理陳述為對任意的 fH2,故我們令 fH2 ;則由 LEMMA2 可知,存在一個序列 {fn}fnH20  使得當  n
||ffn||L2(dP×dt)0再者,因為上述為 L2-norm,(反向)三角不等式 | ||X||||Y|| |||XY||告訴我們
| ||fn||L2(dP×dt)||f||L2(dP×dt) |||fnf||L2(dP×dt)0||fn||L2(dP×dt)||f||L2(dP×dt) 同樣的, I(fn)I(f) in L2(dP),由三角不等式我們可知
||I(fn)||L2(dP)||I(f)||L2(dP) 但由 LEMMA1,我們知道,||fn||L2(dP×dt)=||I(fn)||L2(dP)    (),且由前述討論可知左右兩式的極限都存在,亦即
||fn||L2(dP×dt)||f||L2(dP×dt)||I(fn)||L2(dP)||I(f)||L2(dP)故對式 () 兩邊取極限可得
||I(f)(ω)||L2(dP)=||f||L2(dP×dt)

上式 即為所求。


=====
延伸閱讀
[隨機分析] Ito Integral 淺談 (II) - 再論 Ito 積分的建構

ref:
J. M. Steele, Stochastic Calculus and Financial Applications, Springer

沒有留言:

張貼留言

[人工智慧] 本地端 DeepSeek R1 快速安裝:以 Macbook Pro M4 Chip為例

最近火熱的 DeepSeek R1 模型由於採用了 distill 技術,可以大幅降低計算成本,使得一般人有機會在自家筆電上跑性能逼近 Open AI ChatGPT o1的大語言模型。本文簡單介紹一步安裝在 Macbook Pro 的方法以及使用方法,以下測試採用 Macboo...