2025-11-10T02:38:09.804207

Stochastic Simulation and Monte Carlo Method

Mirzaei

These lecture notes are intended to cover some introductory topics in stochastic simulation for scientific computing courses offered by the IT department at Uppsala University, as taught by the author. Basic concepts in probability theory are provided in the Appendix A, which you may review before starting the upcoming sections or refer to as needed throughout the text.

academic

確率的シミュレーションとモンテカルロ法

基本情報

論文ID: 2501.00997
タイトル: 確率的シミュレーションとモンテカルロ法
著者: Davoud Mirzaei (ウプサラ大学)
分類: math.NA cs.NA stat.CO stat.OT
発表日時: 2024年11月1日 (第2版)
論文リンク: https://arxiv.org/abs/2501.00997

要約

本資料はウプサラ大学IT学部の科学計算コース向けの確率的シミュレーション教授用講義ノートである。内容は確率的シミュレーションの入門的主題を網羅しており、モンテカルロ法、確率変数生成、確率過程、マルコフ連鎖モンテカルロ法(MCMC)などの中核内容を含む。付録Aは確率論の基礎概念を提供し、後続章の理論的基盤として機能する。

研究背景と動機

核心的問題

本講義は科学計算における決定論的モデルと確率的モデルの選択と応用に関する問題の解決を目指している:

モデリング方法の選択: 決定論的モデリングと確率的モデリングの間で適切な方法をいかに選択するか
確率的シミュレーション技術: 複雑なシステムシミュレーションのための確率変数をいかに効果的に生成・利用するか
パラメータ推定: 不確実性環境下においてベイズ的パラメータ推定をいかに実施するか

重要性分析

確率的シミュレーションは現代科学計算において重要な意義を有する:

実システムのモデリング: 多くの実システムは内在的な確率性と不確実性を含む
高次元積分計算: モンテカルロ法は高次元積分計算において独特の利点を有する
複雑システムの分析: 生物システム、金融市場、疾病伝播などの複雑システムは確率的モデリングを必要とする

既存方法の限界

決定論的方法: システムの確率的変動と不確実性を捉えることができない
従来の数値方法: 高次元問題において計算複雑度が急速に増加する
パラメータ推定の困難性: 複雑な事後分布は直接計算が困難である

核心的貢献

体系的教授フレームワーク: 基礎確率論から高度なMCMC法までの完全な教授体系を提供
実用的アルゴリズム実装: 確率変数生成、Gillespieアルゴリズム、MCMCなどの主要アルゴリズムのPython実装を提供
多領域応用例: 放射性崩壊、疾病伝播、金融オプション価格設定、生化学反応などの応用シナリオを網羅
理論と実践の統合: 数学理論と実際のプログラミング実装を有機的に結合

方法の詳細説明

タスク定義

本講義は以下の核心的タスクの解決に主に取り組む:

確率変数生成: 与えられた分布から効率的に確率標本を生成
モンテカルロ積分: 確率的サンプリング法を用いて複雑な積分を計算
確率過程のシミュレーション: マルコフ連鎖、ブラウン運動などの確率過程をモデル化
パラメータ推定: MCMC法によるベイズ的パラメータ推論

核心的方法アーキテクチャ

1. 確率変数生成法

逆変換法 (Inverse Transform Method)

# 基本的考え方: U ~ U(0,1)ならば、X = F^(-1)(U) ~ f
def inverse_transform_sampling(cdf_inverse, n):
    U = np.random.uniform(0, 1, n)
    return cdf_inverse(U)

受理-棄却法 (Acceptance-Rejection Method)

提案分布g(x)と定数Cを用いて目標分布f(x)を界定
受理確率: α = min{f(X)/(Cg(X)), 1}
効率: P(受理) = 1/C

2. モンテカルロ積分

基本的モンテカルロ推定 積分 I = ∫g(x)f(x)dx に対して:

I ≈ (1/N) Σ g(xi), xi ~ f

重要度サンプリング (Importance Sampling)

I = ∫g(x)f(x)dx = ∫g(x)[f(x)/ℓ(x)]ℓ(x)dx
I ≈ (1/N) Σ g(xi)w(xi), xi ~ ℓ, w(xi) = f(xi)/ℓ(xi)

3. 確率過程の生成

マルコフ連鎖の生成

状態遷移行列P = (pij)
定常分布πは πP = π を満たす

ブラウン運動の生成

# Wt+1 = Wt + √(Δt) * Z, Z ~ N(0,1)
def brownian_motion(t_vec, dim):
    W = np.zeros([dim, len(t_vec)])
    for k in range(len(t_vec)-1):
        Z = np.random.normal(0, 1, dim)
        dt = t_vec[k+1] - t_vec[k]
        W[:, k+1] = W[:, k] + np.sqrt(dt) * Z
    return W

4. Gillespieアルゴリズム (SSA)

核心的ステップ:

総倾向関数を計算: a(y) = Σ wj(y)
待機時間を生成: τ ~ Exp(a(y))
反応を選択: k ~ DD(1,...,m, p1,...,pm)
状態を更新: y ← y + vk

技術的革新点

次元無関の収束: モンテカルロ法の収束率O(N^(-1/2))は次元に無関
適応的サンプリング: 重要度サンプリングは適切な提案分布の選択により効率を向上
確率微分方程式の求解: オイラー・丸山法による拡散過程の求解
MCMC収束性: 詳細釣り合い方程式により目標分布への収束を保証

実験設定

応用シナリオ

放射性崩壊モデル
- 決定論的モデル: dy/dt = -λy(t)
- 確率的モデル: y →^λ z (Gillespieアルゴリズム)
SIR疾病伝播モデル
- 状態: S(感受性)、I(感染)、R(回復)
- パラメータ: μ(出生死亡率)、β(感染率)、γ(回復率)
金融オプション価格設定
- 幾何ブラウン運動: dSt = μStdt + σStdWt
- ヨーロッパ型コールオプション: C0 = e^(-rT)Emax(ST-K, 0)