2. Brownian Motion

Tags:

https://www.math.uchicago.edu/~lawler/inprogress

2.1 Limits of sums of independent variables


2024.02.03

Theorem 2.1.1. Central Limit Theorem

mean of μ\mu, variance of σ2<\sigma^2 < \infty 의 i.i.d.인 r.v. X1,X2,...X_1, X_2, ... 에 대해
Zn=(X1+...+Xn)nμσ(n)Z_n = \frac{(X_1 + ... + X_n) - n\mu}{\sigma\sqrt(n)} 라 하면 limnP{aZnb}=Φ(b)Φ(a)\lim_{n\rightarrow\infty}\mathbb{P}\{a \le Z_n \le b\}=\Phi(b)-\Phi(a) 이다. (35p)
이 책에서는 Φ\Phi 가 standard normal distribution의 pdf

  • finite variance라면 XnX_n 이 어떤 분포를 가지든 sample mean은 잘 scale하면 normal로 근사 가능

Theorem 2.1.2.

Convergence to the Poisson
limnP{Yn=k}=eλλkk!\lim_{n\rightarrow\infty}\mathbb{P}\{Y_n=k\}=e^{-\lambda}\frac{\lambda^k}{k!}
Poisson distribution의 경우가 nonnormal limit의 예시. 가정에 따라 극한을 취한 결과는 다를 수 있음

2.2 Multivariate normal distribution


2024.02.03

sequence of r.v. (X1,...,Xn)(X_1, ... ,X_n) 이 joint/multivariate normal distribution을 가진다는 것은, 각각이 여러 independent standard normal r.v.들의 linear combination으로 나타내질 수 있다는 것.

  • 즉 i.i.d.인 ZnZ_n 과 상수 mj,ajkm_j, a_{jk} 에 대해 Xj=mj+aj1Z1+aj2Z2+...+ajmZmX_j=m_j+a_{j1}Z_1+a_{j2}Z_2+...+a_{jm}Z_m 인 것.
  • ajka_{jk}행렬로 나타내면 그것의 곱이 covariance matrix가 될 수 있음
    • mean-zero joint normal is determined by covariance matrix
  • joint normal r.v.s의 special property는 orthogonal하면 independent하다는 것.
    그래서 E(X1,X2)=0\mathbb{E}(X_1,X_2)=0 이면 둘은 independent. (42p)

2.3 Limits of random walks


2024.02.03

12\frac{1}{2} 확률로 1 또는 -1이 나오는 r.v. XjX_j 가 있다고 하자. 이 때 E[Xj]=112+(1)12=0,Var[Xj]=(10)212+(10)212=1\mathbb{E}[X_j]=1 \cdot \frac{1}{2} + (-1) \cdot \frac{1}{2} = 0, Var[X_j]=(1 - 0)^2 \cdot \frac{1}{2} + (-1 - 0)^2 \cdot \frac{1}{2} = 1 이다.

기본적으로는 time increment Δt=1\Delta t = 1, space increment Δx=1\Delta x = 1 이다.
이 때 Δt=1N\Delta t = \frac{1}{N} 이라 하면, space increment ±Δx\pm\Delta x 1N\frac{1}{N} 만큼의 시간 변화에 대응된다.

시간이 1=NΔt1=N\Delta t 만큼 변화할 때 이 process의 값은 W1N=Δx(X1+...+XN)W_1^{N}=\Delta x (X_1 + ... + X_N) 이다.
이 때 우리가 원하는 건, Δx\Delta x 를 적당한 값으로 정해서 Var[W1N]=1Var[W_1^{N}]=1 이 되게 만드는 것이다.

Var[W1N]=Var[Δx(X1+...+XN)]=(Δx)2[Var(X1)+...+Var(XN)]=(Δx)2N=1 \begin{aligned} Var[W_1^{N}] & = Var[\Delta x (X_1 + ... + X_N)] \\ & = (\Delta x)^2[Var(X_1) + ... + Var(X_N)] \\ & = (\Delta x)^2 N \\ & = 1 \\ \end{aligned}

따라서 1=(Δx)2N1 = (\Delta x)^2 N 이고, 1N=(Δx)2\frac{1}{N}=(\Delta x)^2 , Δt=(Δx)2\Delta t=(\Delta x)^2 이다.

2.4 Brownian motion


2024.02.04

BM is a model of random continuous motion.

  • random continuous motion이 되기 위한 세 가지 가정(assumption):

    1. Stationary increments
    2. Independent increments
    3. Continuous paths
  • 위 가정들을 만족하는 BtB_t 가 i.i.d.라면, s<ts < t 에 대해

    • E[Bt]=E[Bs]+E[BtBs]=E[Bs]+E[Bts]\mathbb{E}[B_t] = \mathbb{E}[B_s] + \mathbb{E}[B_t - B_s] = \mathbb{E}[B_s] + \mathbb{E}[B_{t-s}] ,
    • Var[Bt]=Var[Bs]+Var[BtBs]=Var[Bs]+Var[Bts]Var[B_t] = Var[B_s] + Var[B_t - B_s] = Var[B_s] + Var[B_{t-s}] .

Definition

  • 다음을 만족하는 stochastic process BtB_tBrownian motion with drift mm and variance σ2\sigma^2 라고 한다.

    • B0=0B_0=0.
    • s<ts < t 에 대해 BtBsB_t - B_s 는 평균이 m(ts)m(t-s) 이고 분산이 σ2(ts)\sigma^2(t-s) 인 정규분포를 따른다.
    • s<ts < t 이면 r.v. BtBsB_t - B_srsr \le sBrB_r 의 값들과 독립이다.
    • 1의 확률로 함수 tBtt \mapsto B_ttt 에 대한 연속 함수이다.
  • Standard normal을 따르는 ZN(0,1)Z \sim N(0,1) 를 이용해 Y=σZ+mY = \sigma Z + m 이라 나타내면 YN(m,σ2)Y \sim N(m, \sigma^2) 이듯, m=0,σ2=1m=0, \sigma^2=1 인 standard Brownian motion BtB_t 를 이용해 Yt=σBt+mtY_t = \sigma B_t + mt 는 drift가 mm 이고 variance가 σ2\sigma^2 인 Brownian motion이 된다.

  • timestep이 uncountable한 경우에 BM의 존재성에 대한 이야기 (45-46p)

2.5 Construction of Brownian motion


2024.02.04

(존재성에 대한 내용이 궁금할 경우에 보면 됨)

2.6 Understanding Brownian motion


2024.02.04

  • 시뮬레이션을 위해 시간 tt 를 discretize 해서 Δt\Delta t 에 대해 생각해보면, NkN(0,1)N_k \sim N(0,1)에 대해
    B(k+1)ΔtBΔt=ΔtNkB_{(k+1)\Delta t}-B_{\Delta t} = \sqrt{\Delta t}N_k 이다.
  • 그럼 ΔBt=ΔBt+ΔtBtΔt\lvert \Delta B_t \rvert = \lvert \Delta B_{t+\Delta t} - B_t\rvert \approx \sqrt{\Delta t} 라는 것인데
  • 이렇게 되면 도함수의 정의를 생각해봤을 때
    • limΔt0ΔBt+ΔtBtΔt\lim_{\Delta t \rightarrow 0}\frac{\Delta B_{t+\Delta t} - B_t}{\Delta t} 의 분자는 Δt\sqrt{\Delta t} 로 수렴한다.
    • 이 때 작은 값 Δt\Delta t 에 대해 Δt\sqrt{\Delta t}Δt\Delta t 보다 훨씬 크다.
    • 따라서 도함수가 존재할 수 없다.

Theorem 2.6.1

With probability one, the function tBtt \mapsto B_t is nowhere differentiable.

Theorem 2.6.2

앞으로 사용하는 BM 모델의 Hölder exponent는 1/2

2.6.1

Brownian motion은 continuous martingale 이고

2.6.2

Brownian motion은 Markov process 이고

2.6.3

Brownian motion은 Gaussian process 이다.

2.7 Computations for Brownian motion


2024.02.04

skip

2.8 Quadratic variation


2024.02.04

  • Quadratic variation은 다음과 같이 정의된다.
    Xt=limnjtn[X(jn)X(j1n)]2\langle X \rangle _t = \lim_{n \rightarrow \infty}{\sum_{j \le tn}{\biggr[ X(\frac{j}{n}) - X(\frac{j-1}{n}) \biggr]^2}}
  • 즉 increment를 제곱해서 모두 더한 것인데, 이 때 increment의 제곱은 (Δx)2=Δt(\Delta x)^2 = \Delta t 이므로 Δt\sum{\Delta t} 와 같다.
  • Quadratic variation은 drift와 무관하다.
  • Quadratic variation은 tt 시점까지의 randomness의 총량 혹은 총 베팅 금액으로 볼 수 있다. (Section 3.2, 97p)

Theorem 2.8.1

drift가 mm 이고 variance가 σ2\sigma^2 인 Brownian motion BtB_t 에 대해 Xt=σ2t\langle X \rangle _t=\sigma^2t 이다.

Theorem 2.8.2

skip

2.9 Multidimensional Brownian motion


2024.02.04

많은 asset의 가치를 동시에 고려하기 위해 d-dimensioinal의 Multidimensional Brownian motion에 대해 다룬다.

  • d-dimensional process Bt=(Bt1,...,Btd)B_t = (B_t^1, ..., B_t^d) 이 따르는 조건들도 1차원 Brownian motion이 따르는 조건들을 multivariate하게 바꾼 것들이다. (67p)

  • 두 processes X,YX, Y 에 대해 CovariationXt=limnjtn[X(jn)X(j1n)][Y(jn)Y(j1n)]\langle X \rangle _t = \lim_{n \rightarrow \infty}{\sum_{j \le tn}{\biggr[ X(\frac{j}{n}) - X(\frac{j-1}{n}) \biggr]\biggr[ Y(\frac{j}{n}) - Y(\frac{j-1}{n}) \biggr]}} 와 같이 정의한다.

    • 자연스럽게 X,Xt=Xt\langle X, X \rangle _t = \langle X \rangle _t
  • Quadratic variation

2.10 Heat equation and generator


2024.02.04

  • random continuous motion을 하는 heat particle들을 상정하여 Heat flow에 대해 생각해볼 수 있다. heat particle들의 밀도(density)에 의해 1차원 막대 각 지점의 온도가 결정된다고 하자.
  • pt(x)p_t(x) 가 지점 xx 의 시점 tt 에서의 온도라 하고, Rpt(x)dx=1\int_R p_t(x) dx = 1 라 하면 pt(x)p_t(x) 를 Brownian motion의 probability density로 볼 수 있다. 즉 BtB_t 에 대응되는 확률분포를 pt(x)p_t(x) 로 알 수 있는 것

skip

2.11 Exercises


(Skip)

References