確率共鳴をExcel VBAで計算してみる。

この「微分方程式による計算科学入門」の確率微分方程式のところのっていた例題。

地球の百万年程度のスパンでの温度変動を表す簡単な数理モデルとして、

dX(t) = [ X (1-X^2) + a cos(ωt)] dt + σ dW(t)

が上げられている。ここでX(t)は地球の気温で、-1が氷期、1が間氷期を表す。

計算はオイラー・丸山スキームで。上からσ=0.1,0.2,0.3,0.4の時の計算結果。

なるほど。こういう風に確率共鳴が起きるのか。