【FinE】测度变换及Radon-Nicodym导数
導航
- Probability Space
- Random Variable
- 測度變換
- Radon-Nicodym
- 例題
- 解析
- 參考資料
Probability Space
使用三元組(Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P})(Ω,F,P)表示概率空間.
- Ω\OmegaΩ表示樣本空間(sample space),Ω={w1,w2,…,wn}\Omega=\{w_1,w_2, \dots, w_n\}Ω={w1?,w2?,…,wn?}
- F\mathcal{F}F是一個σ\sigmaσ代數,表示域流,包含所有可測的事件,有以下性質
(1) ?∈F\empty \in \mathcal{F}?∈F.
(2) A∈F→Ac∈FA\in \mathcal{F}\to A^c\in \mathcal{F}A∈F→Ac∈F.
(3) A1,A2,?∈F→?i=1∞Ai∈FA_1, A_2, \dots \in \mathcal{F}\to \bigcup\limits_{i=1}^\infty A_i\in \mathcal{F}A1?,A2?,?∈F→i=1?∞?Ai?∈F
(4)A1c,A2c,…,∈F→(?i=1∞Aic)∈F→(?i=1∞Aic)c∈F→((?i=1nAi)c)c=?i=1nAi∈FA_1^c, A_2^c, \dots, \in \mathcal{F}\to(\bigcup\limits_{i=1}^\infty A_i^c)\in \mathcal{F}\to (\bigcup\limits_{i=1}^\infty A_i^c)^c\in \mathcal{F}\to((\bigcap\limits_{i=1}^n A_i)^c)^c=\bigcap\limits_{i=1}^n A_i\in \mathcal{F}A1c?,A2c?,…,∈F→(i=1?∞?Aic?)∈F→(i=1?∞?Aic?)c∈F→((i=1?n?Ai?)c)c=i=1?n?Ai?∈F(德摩根律) - P\mathbb{P}P是概率測度(Probability Measure),所有測度空間上的事件被賦予了[0,1][0, 1][0,1]之間的值. 滿足以下性質(Measure →\to→ Probability Measure)
(1) P(Ω)=1\mathbb{P}(\Omega)=1P(Ω)=1
(2) 如果A1,A2,…A_1, A_2, \dotsA1?,A2?,…是disjoint sets, 那么P(?i=1∞Ai)=∑i=1∞P(Ai)\mathbb{P}(\bigcup\limits_{i=1}^\infty A_i)=\sum\limits_{i=1}^\infty \mathbb{P}(A_i)P(i=1?∞?Ai?)=i=1∑∞?P(Ai?)(可數可加性)
有以下推論
(1) P(?)=0\mathbb{P}(\empty)=0P(?)=0
(2) P(?i=1NAi)=∑i=1NP(Ai)\mathbb{P}(\bigcup\limits_{i=1}^NA_i)=\sum\limits_{i=1}^N\mathbb{P}(A_i)P(i=1?N?Ai?)=i=1∑N?P(Ai?)
(3) P(A)+P(Ac)\mathbb{P}(A)+\mathbb{P}(A^c)P(A)+P(Ac)=1
Random Variable
X:Ω→RX:\Omega\to \mathbb{R}X:Ω→R,如果是高維隨機變量,則為X:Ω→RkX:\Omega\to \mathbb{R}^kX:Ω→Rk.
滿足約束:X?1(B)={w∈Ω:X(w)∈B}∈FX^{-1}(B)=\{w\in \Omega: X(w)\in B\}\in \mathcal{F}X?1(B)={w∈Ω:X(w)∈B}∈F,其中B∈B(R)B\in\mathcal{B(\mathbb{R})}B∈B(R)(Borel Sets).
數學期望
可數無窮
EX=∑w∈ΩX(w)P(w)\mathbb{E}X=\sum\limits_{w\in\Omega}X(w)\mathbb{P}(w) EX=w∈Ω∑?X(w)P(w)
不可數無窮,令w→w+dww\to w+dww→w+dw,劃分為區間積分
E∫ΩX(w)dP(w)\mathbb{E}\int_\Omega X(w)d\mathbb{P}(w) E∫Ω?X(w)dP(w)
在樣本空間的子集上A?Ω,A?FA\subset \Omega, A\subset \mathcal{F}A?Ω,A?F積分
∫AX(w)dP(w)=∫ΩIA(w)X(w)dP(w)\int_AX(w)d\mathbb{P}(w)=\int_\Omega \mathbb{I}_A(w)X(w)d\mathbb{P}(w) ∫A?X(w)dP(w)=∫Ω?IA?(w)X(w)dP(w)
應用:歐式看漲期權定價
設當前時刻為000時刻,期權開始時刻為ttt,行權時刻為TTT,期權價格CtC_tCt?
Ct=EQ[e?r(T?t)(ST?K)+∣Ft]C_t=\mathbb{E}^\mathbb{Q}[e^{-r(T-t)}(S_T-K)^+\mid \mathcal{F}_t] Ct?=EQ[e?r(T?t)(ST??K)+∣Ft?]
性質
如果A?Ω,B?Ω,A,B∈F,A∩B=?A\subset \Omega, B\subset \Omega, A, B\in \mathcal{F}, A\cap B=\emptyA?Ω,B?Ω,A,B∈F,A∩B=?
可以推出
∫A∪BX(w)dP(w)=∫ΩI{A∪B}X(w)dP(w)=∫Ω(IA+IB)X(w)dP(w)=∫ΩIAX(w)dP(w)+∫ΩIBX(w)dP(w)=∫AX(w)dP(w)\begin{aligned} \int_{A\cup B}X(w)d\mathbb{P}(w)&=\int_\Omega\mathbb{I}_{\{A\cup B\}}X(w)d\mathbb{P}(w)\\ &=\int_\Omega(\mathbb{I}_A+\mathbb{I}_B)X(w)d\mathbb{P}(w)\\ &=\int_\Omega\mathbb{I}_AX(w)d\mathbb{P}(w)+\int_\Omega\mathbb{I}_BX(w)d\mathbb{P}(w)\\ &=\int_AX(w)d\mathbb{P}(w) \end{aligned} ∫A∪B?X(w)dP(w)?=∫Ω?I{A∪B}?X(w)dP(w)=∫Ω?(IA?+IB?)X(w)dP(w)=∫Ω?IA?X(w)dP(w)+∫Ω?IB?X(w)dP(w)=∫A?X(w)dP(w)?
期望收斂
要求
E∣X∣=∫Ω∣X(w)∣dP(w)<∞\mathbb{E}|X|=\int_\Omega|X(w)|d\mathbb{P}(w)<\infty E∣X∣=∫Ω?∣X(w)∣dP(w)<∞
一些分布函數的期望,如Laplace Distribution,計算發現期望并不收斂.
測度變換
測度變換的本質是對樣本點的權重進行轉換
| w1w_1w1? | 1/31/31/3 | 1/61/61/6 |
| w1w_1w1? | 1/31/31/3 | 1/21/21/2 |
| w1w_1w1? | 1/31/31/3 | 1/31/31/3 |
考慮不同測度P\mathbb{P}P和Q\mathbb{Q}Q之間的聯系,不難發現存在如下轉換
P×(1/23/21)=Q\mathbb{P}\times\left( \begin{matrix} 1/2\\ 3/2\\ 1 \end{matrix} \right) =\mathbb{Q} P×???1/23/21????=Q
Radon-Nicodym導數可以建立不同測度之間關系
Z(w)=Q(w)P(w)Z(w)=\frac{\mathbb{Q}(w)}{\mathbb{P}(w)} Z(w)=P(w)Q(w)?
具有如下性質
EP[Z(w)]=1EQ[X]=EP[XZ]EP[X]=EQ[XZ],Z>0\begin{aligned} &\mathbb{E}^\mathbb{P}[Z(w)]=1\\ &\mathbb{E}^\mathbb{Q}[X]=\mathbb{E}^\mathbb{P}[XZ]\\ &\mathbb{E}^\mathbb{P}[X]=\mathbb{E}^\mathbb{Q}[\frac{X}{Z}], Z>0 \end{aligned} ?EP[Z(w)]=1EQ[X]=EP[XZ]EP[X]=EQ[ZX?],Z>0?
定理:令(Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P})(Ω,F,P)為一個概率空間,隨機變量Z≥0Z\geq 0Z≥0, 且EP[Z]=1\mathbb{E}^\mathbb{P}[Z]=1EP[Z]=1,對于A?FA\subset \mathcal{F}A?F,定義
Q(A)=∫AZ(w)dP(w)\mathbb{Q}(A)=\int_AZ(w)d\mathbb{P}(w) Q(A)=∫A?Z(w)dP(w)
那么Q\mathbb{Q}Q是一個概率測度. 并且,如果隨機變量X≥0X\geq 0X≥0,那么
EQ[X]=EP[XZ]\mathbb{E}^\mathbb{Q}[X]=\mathbb{E}^\mathbb{P}[XZ] EQ[X]=EP[XZ]
如果隨機變量Z>0Z\gt 0Z>0,那么
EP[X]=EQ[XZ]\mathbb{E}^\mathbb{P}[X]=\mathbb{E}^\mathbb{Q}[\frac{X}{Z}] EP[X]=EQ[ZX?]
Radon-Nicodym
定義:概率測度(Ω,F,P)(\Omega, \mathcal{F}, \mathbb{P})(Ω,F,P),令Q\mathbb{Q}Q為域流F\mathcal{F}F上的另一個概率測度,如果測度Q\mathbb{Q}Q和P\mathbb{P}P相等,需要滿足這樣的條件,對于隨機變量Z(w)>0,a.s.Z(w)>0, a.s.Z(w)>0,a.s.和事件A∈FA\in \mathcal{F}A∈F
Q(A)=∫AZ(w)dP(w)\mathbb{Q}(A)=\int_AZ(w)d\mathbb{P}(w) Q(A)=∫A?Z(w)dP(w)
那么ZZZ就是Radon-Nicodym導數,且
Z(w)=dQ(w)dP(w)Z(w)=\frac{d\mathbb{Q}(w)}{d\mathbb{P}(w)} Z(w)=dP(w)dQ(w)?
在P\mathbb{P}P測度下計算隨機變量Z(w)Z(w)Z(w)的期望得到
EP[Z(w)]=∫ΩZ(w)dP(w)=∫ΩZ(w)dQZ(w)=∫ΩdQ=1\begin{aligned} \mathbb{E}^\mathbb{P}[Z(w)]&=\int_\Omega Z(w)d\mathbb{P}(w) \\ &=\int_\Omega Z(w)\frac{d\mathbb{Q}}{Z(w)}\\ &=\int_\Omega d\mathbb{Q}\\ &=1 \end{aligned} EP[Z(w)]?=∫Ω?Z(w)dP(w)=∫Ω?Z(w)Z(w)dQ?=∫Ω?dQ=1?
Radon-Nicodym導數計算
使用極限的思想計算,設概率測度P\mathbb{P}P為均勻測度,且0≤a≤b≤10\leq a\leq b\leq 10≤a≤b≤1,
P[a,b]=b?a\mathbb{P}[a, b]=b-a P[a,b]=b?a
概率測度Q\mathbb{Q}Q滿足
Q[a,b]=b2?a2\mathbb{Q}[a,b]=b^2-a^2 Q[a,b]=b2?a2
對于小區間[w,w+Δw][w, w+\Delta w][w,w+Δw]計算ΔQ\Delta \mathbb{Q}ΔQ
ΔQ=2wΔw+(Δw)2\Delta\mathbb{Q}=2w\Delta w+(\Delta w)^2 ΔQ=2wΔw+(Δw)2
ΔP\Delta \mathbb{P}ΔP
ΔP=Δw\Delta\mathbb{P}=\Delta w ΔP=Δw
計算R-N導數Z(w)Z(w)Z(w)
Z(w)=dQdP=2w+Δw=lim?Δw→02wZ(w)=\frac{d\mathbb{Q}}{d\mathbb{P}}=2w+\Delta w\xlongequal{\lim \Delta w\to 0}2w Z(w)=dPdQ?=2w+ΔwlimΔw→02w
在新測度下,每個樣本點的權重均發生了變化(發生了拉伸或者壓縮)
例題
隨機變量X(w)=wX(w)=wX(w)=w,在概率測度P\mathbb{P}P下為標準正態分布NP(0,1)\mathcal{N}^\mathbb{P}(0, 1)NP(0,1),隨機變量Y(w)=w+θ,(θ>0)Y(w)=w+\theta, (\theta>0)Y(w)=w+θ,(θ>0)在概率測度P\mathbb{P}P下為正態分布N(θ,1)\mathcal{N}(\theta, 1)N(θ,1). 求一個概率測度Q\mathbb{Q}Q使得隨機變量Y(w)Y(w)Y(w)滿足NQ(0,1)\mathcal{N}^\mathbb{Q}(0, 1)NQ(0,1).
解析
根據定義
Z(w)=dQdP(w)Z(w)=\frac{d\mathbb{Q}}{d\mathbb{P}}(w) Z(w)=dPdQ?(w)
標準正態分布的pdf為
φ(x)=12πe?x22,x∈R\varphi(x)=\frac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2}}, x\in\mathbb{R} φ(x)=2π?1?e?2x2?,x∈R
cdf為
Φ(x)=∫∞x12πe?s22ds\Phi(x)=\int_\infty^x\frac{1}{\sqrt{2\pi}}e^{-\frac{s^2}{2}}ds Φ(x)=∫∞x?2π?1?e?2s2?ds
ΔP(l)=P(l≤w≤l+Δl)=Φ(l+Δl)?Φ(l)\Delta \mathbb{P}(l)=\mathbb{P}(l\leq w\leq l+\Delta l)=\Phi(l+\Delta l)-\Phi(l) ΔP(l)=P(l≤w≤l+Δl)=Φ(l+Δl)?Φ(l)
隨機變量YYY滿足Y~NQ(0,1)Y\sim\mathcal{N}^\mathbb{Q}(0, 1)Y~NQ(0,1)
ΔQ(l)=Q(l≤w≤l+Δl)\Delta\mathbb{Q}(l)=\mathbb{Q}(l\leq w\leq l+\Delta l)ΔQ(l)=Q(l≤w≤l+Δl)
因為Y(w)=X(w)+θ=w+θY(w)=X(w)+\theta=w+\thetaY(w)=X(w)+θ=w+θ
ΔQ?(l)=Q(l≤Y(w)≤l+Δl)=Φ(l+Δl)?Φ(l)=Q(l?θ≤w≤l+Δl?θ)\Delta\mathbb{Q}^*(l)=\mathbb{Q}(l\leq Y(w)\leq l+\Delta l)=\Phi(l+\Delta l)-\Phi(l)=\mathbb{Q}(l-\theta\leq w\leq l+\Delta l-\theta)ΔQ?(l)=Q(l≤Y(w)≤l+Δl)=Φ(l+Δl)?Φ(l)=Q(l?θ≤w≤l+Δl?θ).
由于lll值任意,令l=l+θl=l+\thetal=l+θ,可以得到
Q(l≤w≤l+Δl)=Φ(l+θ+Δl)+Φ(l+θ)=ΔQ(l)\mathbb{Q}(l\leq w\leq l+\Delta l)=\Phi(l+\theta+\Delta l)+\Phi(l+\theta)=\Delta \mathbb{Q}(l) Q(l≤w≤l+Δl)=Φ(l+θ+Δl)+Φ(l+θ)=ΔQ(l)
所以
Z(l)=dQdP=lim?Δl→0Φ(l+θ+Δl)?Φ(l+θ)Φ(l+Δl)?Φ(l)=lim?Δl→0φ(l+θ+Δl)φ(l+Δl)=φ(l+θ)φ(l)\begin{aligned} Z(l)&=\frac{d\mathbb{Q}}{d\mathbb{P}}=\lim_{\Delta l \to 0}\frac{\Phi(l+\theta+\Delta l)-\Phi(l+\theta)}{\Phi(l+\Delta l)-\Phi(l)}\\ &=\lim_{\Delta l\to 0}\frac{\varphi(l+\theta+\Delta l)}{\varphi(l+\Delta l)}\\ &=\frac{\varphi(l+\theta)}{\varphi(l)} \end{aligned} Z(l)?=dPdQ?=Δl→0lim?Φ(l+Δl)?Φ(l)Φ(l+θ+Δl)?Φ(l+θ)?=Δl→0lim?φ(l+Δl)φ(l+θ+Δl)?=φ(l)φ(l+θ)??
求得Radon-Nicodym導數為
Z(w)=φ(w+θ)φ(w)=exp?(?(w+θ)22+w22)=exp?(?θ22?wθ)=exp?[?θ22?X(w)θ]\begin{aligned} Z(w)&=\frac{\varphi(w+\theta)}{\varphi(w)}\\ &=\exp(-\frac{(w+\theta)^2}{2}+\frac{w^2}{2})\\ &=\exp(-\frac{\theta^2}{2}-w\theta)\\ &=\exp[-\frac{\theta^2}{2}-X(w)\theta] \end{aligned} Z(w)?=φ(w)φ(w+θ)?=exp(?2(w+θ)2?+2w2?)=exp(?2θ2??wθ)=exp[?2θ2??X(w)θ]?
參考資料
Radon-Nicodym Derivative Jerry Xu
總結
以上是生活随笔為你收集整理的【FinE】测度变换及Radon-Nicodym导数的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 免流服务器系统怎么选,免流云服务器选购
- 下一篇: Android自定义View的多点触控