ItoProcess

ItoProcess[{a,b},x,t]

伊藤過程 を表す.ただし,である.

ItoProcess[{a,b,c},x,t]

伊藤過程 を表す.ただし,である.

ItoProcess[,{x,x0},{t,t0}]

初期条件 を使う.

ItoProcess[,,,Σ]

共分散 Σ のウィナー(Wiener)過程 を使用する.

ItoProcess[proc]

可能な場合は常に proc を標準伊藤過程に変換する.

ItoProcess[sdeqns,expr,x,t,wdproc]

確率微分方程式 sdeqns,出力式 expr,状態 x,時間 t で指定され,過程 dproc に従う w によって駆動される伊藤過程を表す.

詳細とオプション

  • ItoProcessは伊藤拡散あるいは確率微分方程式としても知られている.
  • ItoProcessは連続時間・連続状態のランダム過程である.
  • ドリフト a 次元ベクトルで拡散 b× 次元行列の場合,過程は 次元で 次元WienerProcessによって駆動される.
  • 係数 a および b によく使われる指定
  • a スカラー,b スカラー
    a スカラー,b ベクトル
    a ベクトル,b ベクトル
    a ベクトル,b 行列
  • 確率微分方程式 は積分方程式 として書かれることもある.
  • デフォルトの初期時間 t0 は0であるとみなされる.デフォルトの初期状態 x0 は0である.
  • デフォルトの共分散 Σ は恒等行列である.
  • 一般共分散 Σ については,ItoProcessは拡散行列 bb.Σ1/2に変換することで過程を正規化する.Σ1/2は可能な場合は Σ の下コレスキー因子である. »
  • 標準伊藤過程は,微分状態 の部分集合からなる出力 を持つ.
  • 標準ItoProcess形式に変換できる過程 proc には,OrnsteinUhlenbeckProcessGeometricBrownianMotionProcessStratonovichProcessItoProcessが含まれる.
  • ItoProcessを標準形式に変換するためには,伊藤の補助定理が自動的に使われる.
  • sdeqns における確率微分方程式は の形式を取ることができる,ここで,\[DifferentialD]であり,ddを使って入力することができる.微分 は伊藤微分であると解釈される.
  • 出力式 exprx[t]t を含む任意の式でよい.
  • 駆動する過程 dproc は標準伊藤過程に変換できる任意の過程でよい.
  • 次は,ItoProcessの関連特性である.
  • "Drift"ドリフト項
    "Diffusion"拡散行列
    "Output"出力状態
    "TimeVariable"時間変数
    "TimeOrigin"時間変数の原点
    "StateVariables"状態変数
    "InitialState"初期状態値
    "KolmogorovForwardEquation"Kolmogorov前進方程式(Fokker-Planckの方程式)
    "KolmogorovBackwardEquation"Kolmogorov後退方程式
    "Derivative"伊藤導関数
    "FeynmanKacFormula"Feynman-Kacの公式から得られた偏微分方程式
  • ItoProcessに特有のRandomFunctionにおけるMethod設定  »
  • "EulerMaruyama"オイラー・丸山(次数1/2,デフォルト)
    "KloedenPlatenSchurz"KloedenPlatenSchurz(次数3/2)
    "Milstein"Milstein(次数1)
    "StochasticRungeKutta"3段階Rossler SRKスキーム(次数1)
    "StochasticRungeKuttaScalarNoise"スカラーノイズのための3段階Rossler SRKスキーム(次数3/2)
  • ItoProcessは,RandomFunctionCovarianceFunctionPDFExpectation等の関数で使うことができる.

例題

すべて開くすべて閉じる

  (1)

過程をその確率微分方程式で定義する:

過程のシミュレーションを行う:

平均値関数を計算する:

共分散関数を計算する:

スコープ  (19)

基本的な用法  (10)

確率微分方程式 からのドリフト と拡散 でウィナー過程を定義する:

パラメトリック過程から直接変換する:

過程 を定義する.ただし,である:

微分表記を使って同じ過程を定義する:

ベクトル過程 を出力 で定義する:

微分表記を使う:

ベクトル過程 を定義する.ただし,である:

微分表記を使う:

ベクトル過程 を定義する.ただし,である:

微分表記を使う:

相関する2つのウィナー過程で駆動される過程を定義する:

正規化された過程はと等しい拡散行列を持つ. は正規化前の拡散行列である:

確率微分方程式 に対応するスカラー過程 を定義する:

確率微分方程式 および に対応するベクトル過程 および を定義する:

2D相関ウィナー過程に相当する過程を定義する:

2D相関ウィナー過程に駆動されるベクトル過程を定義する:

さまざまなメソッドでItoProcessの経路のシミュレーションを行う:

シミュレーションメソッドと対応する次数:

シミュレーションメソッドをRandomFunctionのオプションとして指定する:

過程特性の抽出  (2)

伊藤過程をその確率微分方程式で定義する:

伊藤過程の使用可能な特性:

過程のドリフトと拡散:

Kolmogorov前進方程式:

ここではInactiveを使って偏導関数の展開を防いでいる.Activateを使って式を展開する:

Kolmogorov後退方程式:

関数 の伊藤導関数を計算する.出力はドリフト項と拡散項からなるリストである:

特性"FeynmanKacFormula"は,解 が条件付き期待値 と終端条件 を満足する偏微分方程式を与える:

一般化された状態には追加的な引数を与えることができる.特性"FeynmanKacFormula"は,追加引数 を与えられると解が条件付き期待値 と同じ終端条件を満足する偏微分方程式を与える:

特性"FeynmanKacFormula"は,第3引数 があると解が条件付き期待値 と同じ終端条件を満足する偏微分方程式を与える:

ItoProcessでHestonモデルを定義する:

ドリフト:

(正規化後の)拡散行列:

Kolmogorov前進方程式:

Kolmogorov後退方程式:

伊藤導関数の公式:

特殊伊藤過程  (5)

WienerProcessに対応する伊藤過程:

GeometricBrownianMotionProcessに対応する伊藤過程:

BrownianBridgeProcessに対応する伊藤過程:

OrnsteinUhlenbeckProcessに対応する伊藤過程:

CoxIngersollRossProcessに対応する伊藤過程:

過程スライス特性  (2)

ヤコビの拡散過程を定義する:

時間スライス分布の低次数キュムラントを計算する:

無限時間枠の極限を求める:

一様分布のキュムラントと比較する:

線形確率微分方程式の系で与えられたベクトル過程を定義する:

時間スライス分布の確率密度関数を求める:

の相互共分散を計算する:

無限時間範囲極限は のときにのみ存在する:

アプリケーション  (11)

計算特性  (3)

OrnsteinUhlenbeck過程 ともとになっているウィナー過程 の相互共分散を計算する:

過程 のモーメントを計算する.ただし, は標準ウィナー過程である:

スカラーノイズに駆動されたベクトル伊藤過程(ホワイトノイズに駆動された1D振動子):

過程経路のシミュレーションを行う:

平均値関数と分散関数を計算する:

平均値関数と標準偏差帯を生成された経路とともにをプロットする:

マルチンゲール  (3)

過程 がマルチンゲールである の値を決定する.ただし, は標準ウィナー過程である:

標準形式に変換する:

標準形式のゼロドリフト係数は,がマルチンゲールになるための必要条件である:

ベクトルWiener過程に駆動されたスカラー伊藤過程:

同じ過程を確率方程式によって定義する:

2つのWiener過程によって駆動されたスカラー過程を構築する:

の二次変動:

Lévy条件によると,はブラウン運動である.過程の平均は初期状態に等しい:

モデリング  (2)

熱変動の影響下での自由粒子のダイナミクスは,Langevin運動方程式 によってモデル化できる.ここで は標準のWienerProcessであり, は熱ノイズの強度である.ここでは,にのみ依存し,速度の方程式に焦点を当てることができると想定されている.運動方程式の統合には,伊藤定式化とStratonovich定式化という一般的な2つの方法があるが,それらは次の方法で定義できる:

が定数のとき,2つの定式化は同一で,のときに同じ定常分布に至る:

が速度に依存する場合は,WienerProcessの性質のために,は非零の二次変動を持ち,2つの定式化は異なる結果をもたらす.Stratonovichの公式を伊藤の公式相当に変換する:

Stratonovichの定式化でのドリフトは伊藤の定式化でのドリフトとは異なる:

Gompertz曲線は,腫瘍の成長のような成長過程のモデリングによく使われる.成長過程の対数にガウスノイズを仮定すると,モデルを確率微分方程式として書くことができる:

過程の平均は通常のGompertz曲線である:

時点 における過程のスライス分布はLogNormalDistributionに従う:

および からまでで過程のシミュレーションを行う:

サンプルされた経路を可視化する:

同じ条件で1000個のサンプルを生成し,次に における経路とスライスデータを可視化する:

伊藤過程表現  (3)

ItoProcessを使って標準的なWienerProcessを表す:

はマルチンゲールである.伊藤の補題を使って の導関数を計算する:

CoxIngersollRossProcessを作成し,ItoProcessでこれを表す:

Kolmogorov前進方程式を得る:

における局所化された初期条件,ディリクレ境界条件で方程式を数値的に解く:

Kolmogorov前進方程式の における解をプロットし,これを閉じた形の密度関数と比較する:

Animateで解のダイナミクスを可視化する:

GeometricBrownianMotionProcessItoProcessで表す.ここで, は無リスク金利を, は変動性を表す:

割引過程 についての伊藤導関数を計算する:

古典的なブラック・ショールズ(BlackScholes)方程式は,若干の単純化を加えてドリフトを0にすることで得られる:

この方程式はFeynmanKacの公式から直接得ることもできる:

DSolveでブラック・ショールズ方程式を記号的に解く:

特性と関係  (2)

StratonovichProcessItoProcessに変換する:

変換してもとに戻す:

変換ウィナーはItoProcessに関連している:

平均値関数と分散関数は一致する:

考えられる問題  (2)

ItoProcessはランダムな初期条件をサポートしないので,表現することができない:

しかし,固定された初期条件を持つ過程はサポートする:

駆動された過程の初期時間はItoProcessにマッチする必要がある:

初期時間がマッチすると,これを表現することができる:

Wolfram Research (2012), ItoProcess, Wolfram言語関数, https://reference.wolfram.com/language/ref/ItoProcess.html (2016年に更新).

テキスト

Wolfram Research (2012), ItoProcess, Wolfram言語関数, https://reference.wolfram.com/language/ref/ItoProcess.html (2016年に更新).

CMS

Wolfram Language. 2012. "ItoProcess." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/ItoProcess.html.

APA

Wolfram Language. (2012). ItoProcess. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/ItoProcess.html

BibTeX

@misc{reference.wolfram_2024_itoprocess, author="Wolfram Research", title="{ItoProcess}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/ItoProcess.html}", note=[Accessed: 21-November-2024 ]}

BibLaTeX

@online{reference.wolfram_2024_itoprocess, organization={Wolfram Research}, title={ItoProcess}, year={2016}, url={https://reference.wolfram.com/language/ref/ItoProcess.html}, note=[Accessed: 21-November-2024 ]}