金融工程筆記1:隨機分析基礎
目錄
- 布朗運動
- 伊藤積分
- 伊藤過程與隨機模擬
- 伊藤引理
1.布朗運動
是一個布朗運動,如果
a)
b)對於所有的
獨立
c)對於任意
2.伊藤積分
概念: 是一個布朗運動, 是一個隨機過程,在區間 上做任意劃分:
,令 。
若 存在,則將其記為:
性質:
a)
b)
註:當 為函數時性質也成立
c)
根據伊藤積分定義性質c)顯然,下面對性質a)、b)當 為函數時給出不嚴格說明:
1.當 (常數)
所以
2.當
可以直觀考慮 的期望和二階矩
3.當 為 上的函數時
根據定義
注意:這個並不嚴格
3.伊藤過程與隨機模擬
概念: 是伊藤過程,如果 為隨機微分方程:
或 的解。其中 為函數, 是布朗運動。 叫漂移率, 叫擴散率或標準差, 叫方差率。
注意: 區間可以推廣為
伊藤過程所滿足的隨機微分方程不一定都可以解出 的顯示錶達,不過一般而言可以通過隨機模擬的方法來模擬一個具體伊藤過程的隨機遊動,以及計算 。
伊藤過程隨機模擬:
根據伊藤過程定義,可以得到以下近似公式:
將時間 等分為 份:
則模擬伊藤過程的序列 滿足:
For i = 0 : N-1
當 時,用python伊藤過程進行模擬如下
代碼
import numpy as np
import matplotlib.pyplot as plt
plt.rcParams[font.sans-serif]=[SimHei]
plt.rcParams[axes.unicode_minus] = False
n = 5 #模擬伊藤過程個數
X = 2 #伊藤過程起點
N = 100 #各個伊藤過程模擬次數
T = 100 #伊藤過程從0時刻開始的結束時間
mu = 1
sigma = 4
e = np.random.randn(n,N) #5個伊藤過程,各100個點
Y = np.ones((n,N))
h = T/(N-1)
t = np.arange(N)*h
Y[:,0] = X*np.ones((n))
for i in np.arange(N-1):
Y[:,i+1]= Y[:,i] + mu*h + sigma*np.sqrt(h)*e[:,i]
for i in np.arange(Y.shape[0]):
plt.plot(t,Y[i,:],label= 伊藤過程+str(i+1))
plt.legend(loc=upper right)
plt.xlabel(t)
plt.ylabel(X)
plt.title(伊藤過程模擬)
plt.show()
運行結果