犀角(Diceros Horn)

とくながの「書き散らかし」です

ここは研究・調査・開発などの記録がメインのページです。 日常の雑事、読書記録は はてなダイアリー の方に書いています。よろしければそちらもどうぞ。

Sun, 19 Feb 2006

ベータ分布

ベータ分布乱数生成プログラム(Java Web Start)

ベルヌイ分布に関するベイズ推測を行うと、 事後分布としてベータ分布が現れる。そこでベータ分布の 勉強をかねて、乱数生成プログラムを作成。

履歴

2003/11/03
乱数生成アルゴリズムおよび簡単なインタフェイス作成

実行形式

a,b の値を入力し、[SET PARAMETER] ボタンを押すと、 確率分布が更新される。START ボタンを押せば乱数が生成されます。結果はヒストグラム表示され ます。試行回数によってヒストグラムの幅は自動的に調整されます。

これだけで終わっては仕方ないので、KL距離を計算したり、AIC を計算したりする プログラムに発展させていきたいな。

posted at 15:04 | category: /Program/Math | 固定リンク(ベータ分布 )

正規分布

正規分布の数値計算(Java Web Start)

正規分布表が手元にないときに代わりに数値計算するプログラムです。

posted at 14:53 | category: /Program/Math | 固定リンク(正規分布 )

平方剰余を求める

平方剰余を計算するプログラム(Java Web Start)

いわゆる Legendre 記号を計算する。

posted at 14:42 | category: /Program/Math | 固定リンク(平方剰余を求める )

Sun, 12 Feb 2006

ジュリア集合

ジュリア集合を描く(Java Web Start)

マンデルブロー集合とならんで有名なフラクタル図形であるジュリア集合。 描画する複素平面の左上の点とその大きさを指定して マンデルブロー集合を描画する。マウスによって拡大する領域を選択すること も可能。左上の座標が x + yi でパラメータλ= a + bi です。

posted at 21:16 | category: /Program/Math | 固定リンク(ジュリア集合 )

マンデルブロー集合

マンデルブロー集合を描く(Java Web Start)

有名なフラクタル図形であるマンデルブロー集合。 描画する複素平面の左上の点とその大きさを指定して マンデルブロー集合を描画する。マウスによって拡大する領域を選択すること も可能。

Save Image は Java Web Start では動きません。

posted at 21:07 | category: /Program/Math | 固定リンク(マンデルブロー集合 )

Fri, 10 Feb 2006

Pell 方程式を解く

Pell 方程式の最小解を求めるプログラム(Java Web Start)

x2 - N y2 = ±1
を満たす最小の整数解 (x,y) を求める。

posted at 03:00 | category: /Program/Math | 固定リンク(Pell 方程式を解く )

ブラックショールズの公式

Black Scholes の公式(Java Web Start)

ヨーロッパ型コールオプションの価格をブラック・ショールズの 公式を使って計算するプログラム。ブラック・ショールズの公式とは、 現在の価格が Y の証券を満期期間 T 後に行使価格 K で売るオプションの 現在での合理的価格 C(Y) を求める公式である。ただし、非危険利子率を r 、 ボラティリティを σ とする。

  C(Y) = YΦ(g)-Ke-rTΦ(h)

   g = ((ln(Y/K) + (r + σ2/2) T) / σ T1/2

   h = g - σ T1/2
Φは標準正規分布の累積分布関数で、数値計算には不完全Γ関数を用いる方法を使っ た。

参考文献

  1. 「ファイナンスのための確率微分方程式」 トーマス・ミコシュ 著 (東京電機大学出版局)
  2. 「C言語による最新アルゴリズム辞典」 奥村晴彦 著 (技術評論社)
ブラック=ショールズの公式については[1]が詳しい。 [2]は標準正規分布の分布関数の計算に用いた。

posted at 02:47 | category: /Program/Math | 固定リンク(ブラックショールズの公式 )