確率論および関連分野で使用される数学的手法
確率論および関連分野において、 マリアヴァン計算は、 変分法の数学 的分野を決定論的関数から 確率過程 に拡張する数学的手法とアイデアの集合である 。特に、 ランダム変数 の 導関数の計算を可能にする。マリアヴァン計算は 確率 変分 法とも呼ばれる。P. マリアヴァンは、最初に無限次元空間での計算を開始した。その後、S. 楠岡、D. ストロック、 JM. ビスムット 、 渡辺真三 、I. 重川などの重要な貢献者によって、 最終的に基礎が完成された。
マリアヴァン計算は、ヘルマンダーの条件が確率微分方程式の解の密度の存在と滑らかさを意味することを証明したポール・マリアヴァンにちなんで名付けられました 。 ヘルマン ダー の 最初 の 証明 は 偏 微分方程式 の理論に基づいていました 。この計算は 確率偏微分方程式 にも
適用されています。
この計算は ランダム変数 の 部分積分を可能にし、 数理ファイナンスでは 金融派生商品 の感度を計算するのに 使われます 。この計算は、例えば 確率的フィルタリング に応用されています。
概要と歴史
マリアヴァンは、ヘルマンダーの条件が 確率微分方程式 の解の 密度 の存在を意味することを 確率 的に証明するためにマリアヴァン計算を導入した。 ヘルマンダーの最初の証明は 偏微分方程式 の理論に基づいていた 。彼の計算により、マリアヴァンは解の密度の正則性境界を証明することができた。この計算は 確率偏微分方程式 に適用されている 。
ガウス確率空間
ウィーナー関数( 古典的なウィーナー空間 からの関数) を考え 、その微分を求めることを考えてみましょう。自然なアイデアは ガトー微分を使うことです。
ふ
{\displaystyle F}
だ
グ
ふ
:=
d
d
τ
ふ
[
ふ
+
τ
グ
]
|
τ
=
0
、
{\displaystyle D_{g}F:=\left.{\frac {d}{d\tau }}F[f+\tau g]\right|_{\tau =0},}
しかし、これは常に存在するわけではありません。したがって、方向を制限することによって、そのような空間の新しい微分積分を見つけることは理にかなっています。
マリアバン計算のおもちゃのモデルは、既約な ガウス確率空間 です。これは、すべてが 平均 0 のガウス変数であり、 であるような閉じた 部分空間 を伴う(完全な) 確率空間です 。 の基底を選択すると、 数値モデル と 呼ばれます 。一方、任意の可分ヒルベルト空間に対して、のガウス部分空間を持つ、 シーガルモデル ( アーヴィング・シーガル にちなんで名付けられました)と 呼ばれる標準的 な既約なガウス確率空間が存在します 。この場合、 に対して、 の 関連するランダム変数を と表記し ます 。
バツ
=
(
Ω
、
ふ
、
ポ
、
H
)
{\displaystyle X=(\Omega ,{\mathcal {F}},P,{\mathcal {H}})}
(
Ω
、
ふ
、
ポ
)
{\displaystyle (\Omega ,{\mathcal {F}},P)}
H
⊂
ら
2
(
Ω
、
ふ
、
ポ
)
{\displaystyle {\mathcal {H}}\subset L^{2}(\Omega ,{\mathcal {F}},P)}
H
∈
H
{\displaystyle H\in {\mathcal {H}}}
ふ
=
σ
(
H
:
H
∈
H
)
{\displaystyle {\mathcal {F}}=\sigma (H:H\in {\mathcal {H}})}
H
{\displaystyle {\mathcal {H}}}
バツ
{\displaystyle X}
グ
{\displaystyle {\mathcal {G}}}
セグ
(
グ
)
{\displaystyle \operatorname {Seg} ({\mathcal {G}})}
グ
{\displaystyle {\mathcal {G}}}
グ
∈
グ
{\displaystyle g\in {\mathcal {G}}}
セグ
(
グ
)
{\displaystyle \operatorname {Seg} ({\mathcal {G}})}
わ
(
グ
)
{\displaystyle W(g)}
ガウス確率空間の特性のうち、基底の特定の選択に依存しないものは、 固有 と呼ばれ、 は 外部 の 選択に依存するものと呼ばれます 。 [1] 実空間の可算無限積を と表します 。
R
いいえ
=
∏
私
=
1
∞
R
{\displaystyle \mathbb {R} ^{\mathbb {N} }=\prod \limits _{i=1}^{\infty }\mathbb {R} }
キャメロン・マーティン定理 の現代版を思い出してください
円筒ガウス測度を持つ 局所凸ベクトル空間を 考える 。位相双対の要素について、 平均までの距離を定義する。
え
{\displaystyle E}
γ
{\displaystyle \gamma}
ふ
∈
え
′
{\displaystyle f\in E'}
t
γ
(
ふ
)
:=
ふ
−
∫
え
ふ
(
x
)
γ
(
d
x
)
、
{\displaystyle t_{\gamma }(f):=f-\int _{E}f(x)\gamma (\mathrm {d} x),}
これは写像であり 、閉包を 次のように表す。
t
γ
:
え
′
→
ら
2
(
え
、
γ
)
{\displaystyle t_{\gamma }\colon E'\to L^{2}(E,\gamma )}
ら
2
(
え
、
γ
)
{\displaystyle L^{2}(E,\gamma )}
え
γ
′
:=
クローズ
{
t
γ
(
ふ
)
:
ふ
∈
え
′
}
{\displaystyle E_{\gamma }':=\operatorname {clos} \left\{t_{\gamma }(f)\colon \ f\in E'\right\}}
平行移動を と表記する 。 すると、 その上の 共変性演算子は それぞれ、キャメロン・マーティン空間と呼ばれる 再生核ヒルベルト空間を 誘導し、任意の に対して が 同値と なる 。
γ
メートル
:=
γ
(
⋅
−
メートル
)
{\displaystyle \gamma _{m}:=\gamma (\cdot -m)}
メートル
∈
え
{\displaystyle m\in E}
え
γ
′
{\displaystyle E_{\gamma }'}
R
γ
:
え
γ
′
→
(
え
γ
′
)
∗
{\displaystyle R_{\gamma }:E_{\gamma }'\to (E_{\gamma }')^{*}}
R
{\displaystyle R}
m
∈
R
{\displaystyle m\in R}
γ
m
∼
γ
{\displaystyle \gamma _{m}\sim \gamma }
[2]
実際、 フェルドマン・ハイェックの定理 を使用すると、他の そのような測度では が特異になることが分かります。
h
∉
R
{\displaystyle h\not \in R}
を標準ガウス測度とすると、 からのキャメロン・マーティン定理を 数値 モデル に移すことで 、 の加法群は 上の準自己同型群を定義します 。構築は次のように行うことができます。 の正規直交基底を選択し 、 上 の変換を で表し 、 キャメロン・マーティン空間への写像を で表し 、 と表します。
γ
{\displaystyle \gamma }
(
R
N
,
B
(
R
N
)
,
γ
N
=
⊗
n
∈
N
γ
)
{\displaystyle (\mathbb {R} ^{\mathbb {N} },{\mathcal {B}}(\mathbb {R} ^{\mathbb {N} }),\gamma ^{\mathbb {N} }=\otimes _{n\in \mathbb {N} }\gamma )}
X
{\displaystyle X}
H
{\displaystyle {\mathcal {H}}}
Ω
{\displaystyle \Omega }
H
{\displaystyle {\mathcal {H}}}
τ
α
(
x
)
=
x
+
α
{\displaystyle \tau _{\alpha }(x)=x+\alpha }
R
N
{\displaystyle \mathbb {R} ^{\mathbb {N} }}
α
{\displaystyle \alpha }
j
:
H
→
ℓ
2
{\displaystyle j:{\mathcal {H}}\to \ell ^{2}}
L
∞
−
0
(
Ω
,
F
,
P
)
=
⋂
p
<
∞
L
p
(
Ω
,
F
,
P
)
{\displaystyle L^{\infty -0}(\Omega ,{\mathcal {F}},P)=\bigcap \limits _{p<\infty }L^{p}(\Omega ,{\mathcal {F}},P)\quad }
そして
q
:
L
∞
−
0
(
R
N
,
B
(
R
N
)
,
γ
N
)
→
L
∞
−
0
(
Ω
,
F
,
P
)
,
{\displaystyle \quad q:L^{\infty -0}(\mathbb {R} ^{\mathbb {N} },{\mathcal {B}}(\mathbb {R} ^{\mathbb {N} }),\gamma ^{\mathbb {N} })\to L^{\infty -0}(\Omega ,{\mathcal {F}},P),}
自己準同型 に作用する 加法群の標準的な表現は 次のように定義
される。
ρ
:
H
→
End
(
L
∞
−
0
(
Ω
,
F
,
P
)
)
{\displaystyle \rho :{\mathcal {H}}\to \operatorname {End} (L^{\infty -0}(\Omega ,{\mathcal {F}},P))}
ρ
(
h
)
=
q
∘
τ
j
(
h
)
∘
q
−
1
.
{\displaystyle \rho (h)=q\circ \tau _{j(h)}\circ q^{-1}.}
の作用は 外在的であり、つまり の基底の選択に依存しないことを示すことができる 。
さらに
の 、そして の 無限小生成子 の
ρ
{\displaystyle \rho }
H
{\displaystyle {\mathcal {H}}}
ρ
(
h
+
h
′
)
=
ρ
(
h
)
ρ
(
h
′
)
{\displaystyle \rho (h+h')=\rho (h)\rho (h')}
h
,
h
′
∈
H
{\displaystyle h,h'\in {\mathcal {H}}}
(
ρ
(
h
)
)
h
{\displaystyle (\rho (h))_{h}}
lim
ε
→
0
ρ
(
ε
h
)
−
I
ε
=
M
h
{\displaystyle \lim \limits _{\varepsilon \to 0}{\frac {\rho (\varepsilon h)-I}{\varepsilon }}=M_{h}}
ここで は 恒等演算子であり、 は ランダム変数による乗算演算子(自己準同型に作用する)を表す 。任意のヒルベルト空間 とシーガル模型 の 場合、 (したがって)。すると、上の極限は に関連付けられた ランダム変数による乗算演算子になる 。 [3]
I
{\displaystyle I}
M
h
{\displaystyle M_{h}}
h
∈
H
{\displaystyle h\in {\mathcal {H}}}
G
{\displaystyle {\mathcal {G}}}
Seg
(
G
)
{\displaystyle \operatorname {Seg} ({\mathcal {G}})}
j
:
G
→
ℓ
2
{\displaystyle j:{\mathcal {G}}\to \ell ^{2}}
ρ
:
G
→
End
(
L
∞
−
0
(
Ω
,
F
,
P
)
)
{\displaystyle \rho :{\mathcal {G}}\to \operatorname {End} (L^{\infty -0}(\Omega ,{\mathcal {F}},P))}
W
(
g
)
{\displaystyle W(g)}
g
∈
G
{\displaystyle g\in {\mathcal {G}}}
と に対して、 方向微分を定義する。
F
∈
L
∞
−
0
(
Ω
,
F
,
P
)
{\displaystyle F\in L^{\infty -0}(\Omega ,{\mathcal {F}},P)}
h
∈
H
{\displaystyle h\in {\mathcal {H}}}
⟨
D
F
,
h
⟩
=
D
h
F
=
lim
ε
→
0
(
ρ
(
ε
h
)
−
I
)
F
ε
.
{\displaystyle \langle DF,h\rangle =D_{h}F=\lim \limits _{\varepsilon \to 0}{\frac {\left(\rho (\varepsilon h)-I\right)F}{\varepsilon }}.}
ヒルベルト空間と ガウス空間を持つ シーガルモデルが与えられている。 部分積分公式
を導くことができる。
H
{\displaystyle H}
Seg
(
H
)
{\displaystyle \operatorname {Seg} (H)}
H
=
{
W
(
h
)
:
h
∈
H
}
{\displaystyle {\mathcal {H}}=\{W(h):h\in H\}}
F
∈
L
∞
−
0
(
Ω
,
F
,
P
)
{\displaystyle F\in L^{\infty -0}(\Omega ,{\mathcal {F}},P)}
E
[
D
h
F
]
=
E
[
M
W
(
h
)
F
]
=
E
[
W
(
h
)
F
]
{\displaystyle \mathbb {E} [D_{h}F]=\mathbb {E} [M_{W(h)}F]=\mathbb {E} [W(h)F]}
[4 ]
不変性原理
実数直線全体にわたるルベーグ積分 に対する通常の不変性原理は 、任意の実数 ε と積分可能な関数 f に対して、次の式が成り立つという
ものである。
∫
−
∞
∞
f
(
x
)
d
λ
(
x
)
=
∫
−
∞
∞
f
(
x
+
ε
)
d
λ
(
x
)
{\displaystyle \int _{-\infty }^{\infty }f(x)\,d\lambda (x)=\int _{-\infty }^{\infty }f(x+\varepsilon )\,d\lambda (x)}
そしてそれゆえ
∫
−
∞
∞
f
′
(
x
)
d
λ
(
x
)
=
0.
{\displaystyle \int _{-\infty }^{\infty }f'(x)\,d\lambda (x)=0.}
これは部分積分公式を 導くために使用できる 。なぜなら、 f = gh と設定すると、
0
=
∫
−
∞
∞
f
′
d
λ
=
∫
−
∞
∞
(
g
h
)
′
d
λ
=
∫
−
∞
∞
g
h
′
d
λ
+
∫
−
∞
∞
g
′
h
d
λ
.
{\displaystyle 0=\int _{-\infty }^{\infty }f'\,d\lambda =\int _{-\infty }^{\infty }(gh)'\,d\lambda =\int _{-\infty }^{\infty }gh'\,d\lambda +\int _{-\infty }^{\infty }g'h\,d\lambda .}
同様の考え方は、キャメロン・マーティン・ギルサノフ方向の微分に対する確率解析にも適用できる。実際、 を平方積分可能な 予測可能な過程 とし、
h
s
{\displaystyle h_{s}}
φ
(
t
)
=
∫
0
t
h
s
d
s
.
{\displaystyle \varphi (t)=\int _{0}^{t}h_{s}\,ds.}
がウィーナー過程 である 場合 、 ギルサノフの定理は 不変性原理の次のような類似物を導きます。
X
{\displaystyle X}
E
(
F
(
X
+
ε
φ
)
)
=
E
[
F
(
X
)
exp
(
ε
∫
0
1
h
s
d
X
s
−
1
2
ε
2
∫
0
1
h
s
2
d
s
)
]
.
{\displaystyle E(F(X+\varepsilon \varphi ))=E\left[F(X)\exp \left(\varepsilon \int _{0}^{1}h_{s}\,dX_{s}-{\frac {1}{2}}\varepsilon ^{2}\int _{0}^{1}h_{s}^{2}\,ds\right)\right].}
両辺をεについて微分し、ε=0で評価すると、次の部分積分公式が得られます。
E
(
⟨
D
F
(
X
)
,
φ
⟩
)
=
E
[
F
(
X
)
∫
0
1
h
s
d
X
s
]
.
{\displaystyle E(\langle DF(X),\varphi \rangle )=E{\Bigl [}F(X)\int _{0}^{1}h_{s}\,dX_{s}{\Bigr ]}.}
ここで、左辺は 方向の ランダム変数の マリアバン微分 であり、右辺に現れる積分は 伊藤積分 として解釈されるべきである。
F
{\displaystyle F}
φ
{\displaystyle \varphi }
マリアヴァン計算の最も有用な結果の 1 つは、 マルチンゲール表現定理 のプロセスを 明示的に識別すること
を可能にする クラーク - オコーネ定理です。この定理の簡略版は次のとおりです。
正準空間 上の標準ウィーナー測度 とその正準フィルタリングについて考える。を 満たすに対しては リプシッツであり、 F は C [0,1]
において に対して となる意味で強い微分核を持つ。
C
[
0
,
1
]
{\displaystyle C[0,1]}
F
:
C
[
0
,
1
]
→
R
{\displaystyle F:C[0,1]\to \mathbb {R} }
E
(
F
(
X
)
2
)
<
∞
{\displaystyle E(F(X)^{2})<\infty }
φ
{\displaystyle \varphi }
lim
ε
→
0
1
ε
(
F
(
X
+
ε
φ
)
−
F
(
X
)
)
=
∫
0
1
F
′
(
X
,
d
t
)
φ
(
t
)
a
.
e
.
X
{\displaystyle \lim _{\varepsilon \to 0}{\frac {1}{\varepsilon }}(F(X+\varepsilon \varphi )-F(X))=\int _{0}^{1}F'(X,dt)\varphi (t)\ \mathrm {a.e.} \ X}
それから
F
(
X
)
=
E
(
F
(
X
)
)
+
∫
0
1
H
t
d
X
t
,
{\displaystyle F(X)=E(F(X))+\int _{0}^{1}H_{t}\,dX_{t},}
ここで、 Hは F '( x , ( t ,1])の可視射影であり 、これは関数 F をその定義域の一部( t ,1])にわたってプロセス X を適切に平行移動させたときの微分として見ることができる。
これをもっと簡潔に表現すると、
F
(
X
)
=
E
(
F
(
X
)
)
+
∫
0
1
E
(
D
t
F
∣
F
t
)
d
X
t
.
{\displaystyle F(X)=E(F(X))+\int _{0}^{1}E(D_{t}F\mid {\mathcal {F}}_{t})\,dX_{t}.}
マリアヴァン計算の正式な開発における作業の多くは、上で使用した微分カーネルを、 結果の上記の記述で示された「 マリアヴァン微分 」に置き換えることによって 、この結果を関数 Fの可能な限り最大のクラスに拡張することである。 [ 引用が必要 ]
D
t
{\displaystyle D_{t}}
スコロホッド積分
慣例的にδと表記されるスコロホッド積分作用素は、ヒルベルト空間が空間であるときの白色ノイズの場合のマリアヴァン微分の随伴作用素として定義される 。 したがって 、 の 部分集合である作用素の領域内のuに対して 、 マリアヴァン微分の領域内の
Fに対して、次式が成立する必要がある。
L
2
{\displaystyle L^{2}}
L
2
(
[
0
,
∞
)
×
Ω
)
{\displaystyle L^{2}([0,\infty )\times \Omega )}
E
(
⟨
D
F
,
u
⟩
)
=
E
(
F
δ
(
u
)
)
,
{\displaystyle E(\langle DF,u\rangle )=E(F\delta (u)),}
ここで内積は 、
L
2
[
0
,
∞
)
{\displaystyle L^{2}[0,\infty )}
⟨
f
,
g
⟩
=
∫
0
∞
f
(
s
)
g
(
s
)
d
s
.
{\displaystyle \langle f,g\rangle =\int _{0}^{\infty }f(s)g(s)\,ds.}
この随伴作用素の存在は、 ヒルベルト空間 上の線型作用素に対する リースの表現定理 から導かれます。
u が適応される
と、
δ
(
u
)
=
∫
0
∞
u
t
d
W
t
,
{\displaystyle \delta (u)=\int _{0}^{\infty }u_{t}\,dW_{t},}
ここで積分は伊藤の意味で理解される。したがって、これは伊藤積分を非適応積分関数に拡張する方法を提供する。
アプリケーション
この計算は、 確率変数 の 部分積分を可能にし、 数理ファイナンスにおいて 金融派生商品 の感度を計算するために 使用されます 。この計算は、例えば 確率的フィルタリング に応用されています。
参考文献
^ ポール・マリアビン (1997)。 確率的分析 。 Grundlehren der mathematischen Wissenschaften。ベルリン、ハイデルベルク:シュプリンガー。 4~15ページ。 ISBN 3-540-57024-1 。
^ Bogachev, Vladimir (1998). ガウス測度 ロードアイランド州: アメリカ数学会 .
^ ポール・マリアビン (1997)。 確率的分析 。 Grundlehren der mathematischen Wissenschaften。ベルリン、ハイデルベルク:シュプリンガー。 20~22ページ。 ISBN 3-540-57024-1 。
^ ポール・マリアビン (1997)。 確率的分析 。 Grundlehren der mathematischen Wissenschaften。ベルリン、ハイデルベルク:シュプリンガー。 p. 36.ISBN 3-540-57024-1 。
楠岡 誠、ストロック 大澤 剛 (1981)「マリアビン微積分の応用 I」、 確率解析、谷口国際シンポジウム堅田・京都 1982 講演論文集、271~306 ページ
楠岡 誠、ストロック 大 (1985)「マリアビン微積分の応用 II」、 東京大学理学部誌、1A 数学編 、32 pp 1–76
楠岡 誠、ストロック 大 (1987)「マリアビン微積分の応用 III」、 東京大学理学部誌、第 1A 部、数学 、34 pp 391–442
ポール・マリアヴィン、アントン・タルマイアー共著『 数理ファイナンスにおける変分法の確率計算』 、Springer 2005、 ISBN 3-540-43431-3
ヌアルト、デイヴィッド (2006)。 マリアヴァン計算と関連トピック (第 2 版)。Springer- Verlag。ISBN 978-3-540-28328-7 。
ベル、デニス。(2007) マリアヴァン微積分 、ドーバー。ISBN 0-486-44994-7 ; 電子 書籍
Sanz-Solé, Marta (2005) Malliavin Calculus, with applications to stochastic partial different equations . EPFL Press、CRC Press、Taylor & Francis Group が配布。
シラー、アレックス (2009) 金融アプリケーションにおけるモンテカルロシミュレーションのためのマリアバン計算。プリンストン大学数学部論文
Øksendal, Bernt K. (1997) マリアビン微積分入門と経済学への応用。講義ノート、オスロ大学数学部 (論文と補遺を含む Zip ファイル)
Di Nunno, Giulia 、Øksendal, Bernt、 Proske , Frank (2009)「Lévy Processes のマリアヴァン計算と金融への応用」、Universitext、Springer。ISBN 978-3-540-78571-2
外部リンク
ウィキクォートにおけるマリアヴァン計算に関する引用
Friz, Peter K. (2005-04-10). 「マリアヴァン微積分入門」 (PDF) 。2007-04-17 に オリジナル (PDF)からアーカイブ 。2007-07-23 に取得 。 講義ノート、43ページ
Zhang, H. (2004-11-11). 「マリアヴァン計算」 (PDF) . 2004-11-11 に閲覧。 論文、100ページ