Black Scholes Equation ブラック ショールズ方程式 知的財産仲裁センター知的財産価値評価人候補者 序 弁理士堀城之 これが ブラック ショールズ方程式であり フィッシャー ブラックとマイロン ショールズにより発明された偏微分方程式である オプション取引 ( ヨーロピアンオプショ

Similar documents
パソコンシミュレータの現状

Microsoft Word - NumericalComputation.docx

Microsoft Word - å“Ÿåłžå¸°173.docx

DVIOUT-SS_Ma

受信機時計誤差項の が残ったままであるが これをも消去するのが 重位相差である. 重位相差ある時刻に 衛星 から送られてくる搬送波位相データを 台の受信機 でそれぞれ測定する このとき各受信機で測定された衛星 からの搬送波位相データを Φ Φ とし 同様に衛星 からの搬送波位相データを Φ Φ とす

微分方程式による現象記述と解きかた

線積分.indd

( 慣性抵抗 ) 速度の 2 乗に比例流体中を進む物体は前面にある流体を押しのけて進む. 物 aaa 体の後面には流体が付き従う ( 渦を巻いて ). 前面にある速度 0 の流体が後面に移動して速度 vとなったと考えてよい. この流体の質量は単位時間内に物体が押しのける体積に比例するので,v に比例

1.民営化

多変量解析 ~ 重回帰分析 ~ 2006 年 4 月 21 日 ( 金 ) 南慶典

<4D F736F F D2094F795AA95FB92F68EAE82CC89F082AB95FB E646F63>

Microsoft PowerPoint - e-stat(OLS).pptx

Probit , Mixed logit

Microsoft Word - 1B2011.doc

以下 変数の上のドットは時間に関する微分を表わしている (ex. 2 dx d x x, x 2 dt dt ) 付録 E 非線形微分方程式の平衡点の安定性解析 E-1) 非線形方程式の線形近似特に言及してこなかったが これまでは線形微分方程式 ( x や x, x などがすべて 1 次で なおかつ

NumericalProg09

解析力学B - 第11回: 正準変換

(Microsoft Word - 10ta320a_\220U\223\256\212w\223\301\230__6\217\315\221O\224\274\203\214\203W\203\201.docx)

スライド 1

データ解析

<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63>

不偏推定量

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

Microsoft Word - 高森-伊藤_優先株JAROS_2006.doc

<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63>

様々なミクロ計量モデル†

スライド 1

行列、ベクトル

19年度一次基礎科目計算問題略解

Microsoft Word - thesis.doc

Microsoft PowerPoint - 第7章(自然対流熱伝達 )_H27.ppt [互換モード]

0. はじめに ここでは 金融工学の基礎であるブラックショールズの公式を導くまでの過程を説明する そのためには ランダムウォークから派生したブラウン運動と確率積分の概念の理解は必要不可欠である そしてそこから求まる伊藤の公式を用いて確率微分方程式を解き ブラックショールズ過程について紹介する 1.

Microsoft PowerPoint - aep_1.ppt [互換モード]

応用数学Ⅱ 偏微分方程式(2) 波動方程式(12/13)

OpenFOAM(R) ソースコード入門 pt1 熱伝導方程式の解法から有限体積法の実装について考える 前編 : 有限体積法の基礎確認 2013/11/17 オープンCAE 富山富山県立大学中川慎二

<4D F736F F D20824F B834E835882CC92E8979D814690FC90CF95AA82C696CA90CF95AA2E646F63>

FEM原理講座 (サンプルテキスト)

memo

3 数値解の特性 3.1 CFL 条件 を 前の章では 波動方程式 f x= x0 = f x= x0 t f c x f =0 [1] c f 0 x= x 0 x 0 f x= x0 x 2 x 2 t [2] のように差分化して数値解を求めた ここでは このようにして得られた数値解の性質を 考

講義「○○○○」

微分方程式補足.moc

PowerPoint Presentation

統計学 - 社会統計の基礎 - 正規分布 標準正規分布累積分布関数の逆関数 t 分布正規分布に従うサンプルの平均の信頼区間 担当 : 岸 康人 資料ページ :

ナビエ・ストークス方程式

構造力学Ⅰ第12回

ヤコビ楕円関数とはなにか

2015-2017年度 2次数学セレクション(複素数)解答解説

4STEP 数学 Ⅲ( 新課程 ) を解いてみた関数 1 微分法 1 微分係数と導関数微分法 2 導関数の計算 272 ポイント微分法の公式を利用 (1) ( )( )( ) { } ( ) ( )( ) ( )( ) ( ) ( )( )

s とは何か 2011 年 2 月 5 日目次へ戻る 1 正弦波の微分 y=v m sin ωt を時間 t で微分します V m は正弦波の最大値です 合成関数の微分法を用い y=v m sin u u=ωt と置きますと dy dt dy du du dt d du V m sin u d dt

PowerPoint Presentation

OCW-iダランベールの原理

ギリシャ文字の読み方を教えてください

Microsoft PowerPoint - H17-5時限(パターン認識).ppt

2014年度 千葉大・医系数学

Microsoft PowerPoint - 1章 [互換モード]

2018年度 東京大・理系数学

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (

ミクロ経済学・基本講義 第2回

PowerPoint プレゼンテーション

講義資料サンプル

2016年度 筑波大・理系数学

14 化学実験法 II( 吉村 ( 洋 mmol/l の半分だったから さんの測定値は くんの測定値の 4 倍の重みがあり 推定値 としては 0.68 mmol/l その標準偏差は mmol/l 程度ということになる 測定値を 特徴づけるパラメータ t を推定するこの手

PowerPoint プレゼンテーション

ÿþŸb8bn0irt

Microsoft PowerPoint - H22制御工学I-2回.ppt

微分方程式 モデリングとシミュレーション

Laplace2.rtf

耳桁の剛性の考慮分配係数の計算条件は 主桁本数 n 格子剛度 zです 通常の並列鋼桁橋では 主桁はすべて同じ断面を使います しかし 分配の効率を上げる場合 耳桁 ( 幅員端側の桁 ) の断面を大きくすることがあります 最近の桁橋では 上下線を別橋梁とすることがあり また 防音壁などの敷設が片側に有る

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように

ニュートン重力理論.pptx

多次元レーザー分光で探る凝縮分子系の超高速動力学

Microsoft Word - 補論3.2

Microsoft Word - t30_西_修正__ doc

Microsoft PowerPoint - 熱力学Ⅱ2FreeEnergy2012HP.ppt [互換モード]

2016年度 京都大・文系数学

Transcription:

Black Scholes Equation ブラック ショールズ方程 知的財産仲裁センター知的財産価値評価人候補者 序 弁理士堀城之 これが ブラック ショールズ方程であり フィッシャー ブラックとマイロン ショールズにより発明された偏微分方程である オプション取引 ( ヨーロピアンオプション ) における理論価格の決定に用いられるものとして広く普及した この功績により 1997 年のノーベル経済学賞の受賞対象になった 現在では ブラック ショールズ方程が 知的財産価値評価 (pl-x 社の TRRU(TM) 等 ) パテントポートフォリオ 知的財産価値デューデリジェンスなどにも適用されている また 知財過度な節税防止海外移転後 高収益で再課税 という記事も発表され (2017 年 3 月 10 日日本経済新聞電子版 ) 当該知財評価にも用いられる これらの手法においてはブラック ショールズ方程は基本の基であり且つ大変重要なであり 知的財産価値評価にも用いられる しかし ブラック ショールズ方程に関する数学的アプローチが書かれた書籍を見たことが無く 確率 統計などになじみが無い方には数学的に難しい また経済用語などは我々弁理士にはとっつきにくい そこで 可能な限り分かり易くブラック ショールズ方程を纏めてみた 以下 ブラック ショールズ方程の導出 解法を行う 1.Winner 過程 1

確率統計でよく用いられるものであるが 難しいことはなくブラウン運動で あり よっぱらいが歩いた過程や株価である Winner 過程では 確率過程 W(ti) の変化量 W(ti)- W(t(i-1)) は 正規分布 N(0, σ =(tk- t(k-1))) に従う ΔX=aΔt + bδw (1) を Winner 過程の一般化という なぜ一般化と言うかというと 筆者が思うに マクロ的には株価等はある直線に従って増減すると共に ミクロ的には Winner 過程で増減しているからだと解釈している また Δとしているのは 微分不可能だから 以上が ブラック ショールズ方程の前提の前提 2. 伊藤のレンマ元京都大学教授の伊藤清氏が考案したレンマであり ブラック ショールズ方程の前提である 伊藤のレンマを導くための前提が伊藤過程である まずは一般的な解法を記載する (1) 伊藤過程とは 上記一般化した Winner 過程である (1) を以下のように変形したものをいう 以下 文系の読者にも分かり易いように丁寧に記載する 2

dxt=a(t, Xt)dt+b (t, Xt)dBt (2) 伊藤過程では 前提として 微分可能であり 2 階のテーラー展開が可能な ので として多変数テーラー展開すると 三次以降及び二乗微小項を排除すると ここで は 上記のごとく確率過程であるから その期待値は Winner 過程は 上記の如く正規分布に従うから ( 期待値 0) 分散の公から 正規分布 から分散は dt なので (4) 及び (5) から は正の期待値であり 大数の法則 を鑑みると と書いても構わない (7) を (3) にいれると 3

これが伊藤のレンマである 大数の法則とは Winner 過程を前提に記載すると 試行を繰り返した結果が正規分布 N(0, σ =(tk- t(k-1))) に従うこととなるということ 例えば 試行結果が正規分布になるように算数の同一テストを沢山行ったということではなく 日本人全員でテストしたら結果が正規分布に従ったということである 換言すれば良いテストでしたということ 実際には 全世界で正規分布となるテストを日本で行うと 識字率に教養の高さが比例することが真であるとして日本人は右正規分布に偏ったえぼし岩のような結果になるかもしれない (2) 流体力学的アプローチ検証する意味で 筆者の専門の一つである流体でも導出できた ( 筆者が調べた限りにおいて 世界で筆者だけの方法 ) 伊藤過程は 微分可能なので流線と等価であると仮定することができるとする 流線を 以下の如く時間微分 場所微分が可能な関数で表す 微小点 流線 4

流線における微小点を拡大すると以下になる この図は 微小時間経過すると X 方向に進み その結果 はだけ 増加することを示す この図をで表すと 2 次までテーラー展開すると ここで 下付き添え字は当該文字での微分を示す 微小項を無視すると 伊藤過程を示す (1), (2) から 微小項を無視すると 大数の法則から ( (9) は (8) と数学的等価であり 是すなわち伊藤のレンマである 流体においても砕波 滝等に水が落ちた場合の水滴の動きに適用できるかもしれない 3. ブラック ショールズモデル 5

株価の確率過程を S 時刻 t における株価をとすると 微小時間経過後の収益は 単位時間当たりの儲けは 株価を 1. で述べたように比例部分と変動部分 ( 確率過程 ) とで表すと (11) を (10) に代入すると ( ) これをブラック ショールズモデルという : ドリフト : ボラティリティボラティリティは 確率過程における株価変動率であり 大きければ値動きが激しく リスクも大きいということになる 4. ブラック ショールズ方程 今 オプション ( 売買権 ) を として伊藤のレンマ ( (8)) を適用すると 微分項で整理すると ( ) (12) から単位時間当たりの価値変化率は ( ) (14) 右辺に (13),(12) を代入すると 6

( ) だから 利子率を r とすると ( 儲かるためには r > 0) だから このの左辺に (16) を入れると 両辺を時間積分して整理すると これがブラック ショールズ方程である 5. ブラック ショールズ方程の解法 ブラック ショールズ方程は wikipedia 等にも記載されているように熱 伝導方程に変換できる 以下簡便に記載する と仮定すると T: ポートフォリオ満了日 t: 初期値 0 から経過時間 7

とすると (18) (19) を (17) に入れると V の解を と仮定すると (21)~(23) を (20) に代入すると 両辺を で除して 微分階数で整理すると 熱伝導方程にするには 第 3 項及び第 4 項の [ ] の中をゼロにする必要があるから連立方程をとして a,b を計算しそれを代入するれば良い すると (24) は となる 8

と置くと (25) は となる 因みに左辺を 2 階にすれば振動方程である いずれも 変数分離で解くことができる 熱伝導方程の解は 熱力学の書物に出ているのでそれらを参照されたい ブラック ショールズモデルは ヨーロピアン コールオプション ( 売買権 ) の価格を決定するものなので ヨーロピアンコールオプションの形で解を記載すると以下になる ここで N: 標準正規累積分布関数, : 株価 : 行使価格 : ボラティリティ : r: 利子率 T: ポートフォリオ満了日 t: 初期値 0からの経過時間さらに使いやすい汎用ブラック ショールズモデルと呼ばれている形も記載する q : 原資産利回り Tt : 期間 (1 年なら 1 半年なら 0.5 ) 計算するには 各係数を widows の関数電卓で計算して (28) に入れれば良いだけである エクセルで計算を作成し 各係数を入力するようにすれば簡単に計算できる N については エクセルの normsdist 関数で求めることができる 計算例 9

金融大学 (http://www.findai.com/) から引用許可を得た計算例を以下に記載する S: 原資産価格 ( 株価 ) ----- 100 円 K: 行使価格 ----- 99 円 σ: ボラティリティ ----- 10% ボラティリティが唯一の未知数です Tt: オプションの有効期間 ----- 0.5( 半年 ) r: 安全利子率 ----- 3% Ln( d = 1 1 2 ) ( r+ ) t 2 t 100 1 2 Ln( ) (0.03+ 0.1 ) 0.5 99 = 2 0.1 0.5 =0.3896 2 1 s k d =d - t =0.3896-0.1 0.5 =0.3189 N(d 1 )=N(0.3896)=0.6516 N(d 2 )=N(0.3189)=0.6251 x rt C SN( d1) e KN( d2) 100 0.6516-0.9851 99 0.6251 =4.1953 これがコール価格である 6. 纏めブラック ショールズ方程は ナビエストークス方程のようなエレガントさはない また 伊藤のレンマが無ければ絶対にできなかった 数学的には伊藤のレンマの方が格段に興味をそそられる 10

ただ 伊藤のレンマは上記の如く N(0, σ =(tk- t(k-1))) の正規分布に従うことを前提にしているので ( 当該正規分布に従わなければ数学的に解けない ) 大数の法則で記載したえぼし岩の如く試行結果は正規分布に従わない場合があり 故に リーマンショックが生まれたのであろう もしかしたら 正規分布の代わりに以下の確率密度関数で表される分布を基 礎とするとボラティリティの発散も解消でき第 2 リーマンショックの発生も防 止できるかもしれない 終わり 11