**布朗运动**,又称**维纳过程**(Wiener Process),是连续时间、连续状态的最基本随机过程之一

布朗运动及其应用概述

布朗运动(Brownian Motion)

布朗运动,又称维纳过程(Wiener Process),是连续时间、连续状态的最基本随机过程之一。它最初由植物学家罗伯特·布朗在1827年观察花粉微粒在液体中的不规则运动时发现,后经爱因斯坦(1905)和诺伯特·维纳(1923)等人建立数学模型,成为现代概率论与随机分析的核心工具。


一、定义

一个随机过程 $ {B(t)}_{t \geq 0} $ 称为 标准布朗运动(Standard Brownian Motion),如果满足以下四个条件:

  1. 初始值为零:$ B(0) = 0 $ 几乎必然
  2. 独立增量:对任意 $ 0 \leq s < t \leq u < v $,增量 $ B(t) - B(s) $ 与 $ B(v) - B(u) $ 独立
  3. 平稳正态增量:对任意 $ 0 \leq s < t $,
    B ( t ) − B ( s ) ∼ N ( 0 , t − s ) B(t) - B(s) \sim N(0, t - s) B(t)B(s)N(0,ts)
  4. 路径连续:样本路径 $ t \mapsto B(t) $ 是连续函数(几乎处处成立)

📌 若去掉 $ B(0)=0 $ 的限制,则称其为从某点 $ x $ 出发的布朗运动,记作 $ B_x(t) = x + B(t) $


二、推广形式

1. 带漂移的布朗运动

X ( t ) = μ t + σ B ( t ) X(t) = \mu t + \sigma B(t) X(t)=μt+σB(t)

  • $ \mu $:漂移系数(趋势项)
  • $ \sigma > 0 $:扩散系数(波动率)
  • 应用于金融资产价格建模(如股票)

2. 几何布朗运动(GBM)

S ( t ) = S 0 exp ⁡ ( ( μ − σ 2 2 ) t + σ B ( t ) ) S(t) = S_0 \exp\left( \left(\mu - \frac{\sigma^2}{2}\right)t + \sigma B(t) \right) S(t)=S0exp((μ2σ2)t+σB(t))

  • 满足随机微分方程:
    d S ( t ) = μ S ( t ) d t + σ S ( t ) d B ( t ) dS(t) = \mu S(t) dt + \sigma S(t) dB(t) dS(t)=μS(t)dt+σS(t)dB(t)
  • Black-Scholes 期权定价模型的基础

三、关键性质

性质描述
马尔可夫性给定当前状态 $ B(t) $,未来演化与过去无关 ⇒ 是马尔可夫过程
鞅性$ E[B(t) \mid \mathcal{F}_s] = B(s),\ s < t $ ⇒ 是连续时间鞅
无记忆性(来自正态分布)增量仅依赖于时间差,具有时齐性
不可微性路径处处连续但几乎处处不可导 ⇒ 极端震荡
二次变差非零$ B,B = t $,这是伊藤积分理论的基础
尺度不变性(自相似性)对任意 $ c > 0 $,过程 $ {c^{-1/2} B(ct)} $ 也是标准布朗运动
反射原理$ P\left(\max_{0 \leq s \leq t} B(s) \geq a\right) = 2P(B(t) \geq a) $,用于首达时间分析

四、首达时间与最大值分布

定义 首达时间(First Passage Time):
τ a = inf ⁡ { t > 0 : B ( t ) = a } , a > 0 \tau_a = \inf\{t > 0 : B(t) = a\},\quad a > 0 τa=inf{t>0:B(t)=a},a>0
则其分布为:
P ( τ a ≤ t ) = 2 P ( B ( t ) ≥ a ) = 2 ( 1 − Φ ( a t ) ) P(\tau_a \leq t) = 2P(B(t) \geq a) = 2\left(1 - \Phi\left(\frac{a}{\sqrt{t}}\right)\right) P(τat)=2P(B(t)a)=2(1Φ(t a))
密度函数为:
f τ a ( t ) = a 2 π t 3 e − a 2 / ( 2 t ) , t > 0 f_{\tau_a}(t) = \frac{a}{\sqrt{2\pi t^3}} e^{-a^2/(2t)},\quad t > 0 fτa(t)=2πt3 aea2/(2t),t>0
这称为逆高斯分布(Inverse Gaussian Distribution)。

📌 特别地:$ E[\tau_a] = \infty $,即平均首次到达时间为无穷!


五、多维布朗运动

设 $ B_1(t), B_2(t), \ldots, B_d(t) $ 是 $ d $ 个相互独立的标准布朗运动,则向量过程:
B ( t ) = ( B 1 ( t ) , B 2 ( t ) , … , B d ( t ) ) \mathbf{B}(t) = (B_1(t), B_2(t), \ldots, B_d(t)) B(t)=(B1(t),B2(t),,Bd(t))
称为 $ d $-维标准布朗运动。

  • 每一分量独立且服从一维布朗运动
  • 常用于物理粒子在空间中的随机扩散建模

🔁 Polya 定理指出:当 $ d=1,2 $ 时,d维布朗运动是常返的;当 $ d \geq 3 $ 时为暂态


六、与其它随机过程的关系

关系说明
随机游走的极限根据 Donsker 定理,适当缩放后的对称随机游走弱收敛到布朗运动
泊松过程对比
  • 布朗运动:连续状态、路径连续、无限变差
  • 泊松过程:离散跳跃、有限变差、不连续路径 |
    | 伊藤过程的一般化 | 布朗运动是所有扩散过程的基本“驱动噪声” |
    | 作为高斯过程 | 布朗运动是均值为0、协方差函数为 $ \min(s,t) $ 的高斯过程 |

七、应用领域

领域应用方式
金融工程股价建模(几何布朗运动)、利率模型、期权定价
物理学微观粒子的热运动、扩散方程解的概率表示
生物学动物迁徙路径、细胞内分子运动
信号处理噪声建模(白噪声积分)
机器学习扩散模型(Diffusion Models)中前向过程基于布朗运动设计
偏微分方程利用费曼-卡茨公式将PDE解表示为布朗运动的期望

八、模拟方法(Python 示例)

import numpy as np
import matplotlib.pyplot as plt

def simulate_brownian_motion(T, N):
    """模拟 [0, T] 上的布朗运动,N 为步数"""
    dt = T / N
    dW = np.random.normal(0, np.sqrt(dt), N)  # 增量 ~ N(0, dt)
    W = np.cumsum(dW)
    W = np.insert(W, 0, 0)  # 加上 B(0) = 0
    t = np.linspace(0, T, N+1)
    return t, W

# 参数设置
T = 1
N = 500

t, B = simulate_brownian_motion(T, N)

plt.plot(t, B, label="Brownian Motion")
plt.xlabel("Time $t$")
plt.ylabel("B(t)")
plt.title("Simulation of Standard Brownian Motion")
plt.grid(True)
plt.legend()
plt.show()

在这里插入图片描述

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Bol5261

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值