(57)【特許請求の範囲】
【請求項１】
株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、Galerkin有限要素法とBTD法を用いて、対数株価確率分布の時間発展を数値解析する株価変動数値解析装置において、直前の過去の一定期間の株価データを入力するデータ入力手段と、前記データ入力手段から入力された直前の過去の一定期間の株価データを記憶するデータ記憶手段と、前記データ入力手段から入力された直前の過去の一定期間の株価データからトレンドμとボラティリティーσを求めて初期値とするパラメータ初期値計算手段と、前記パラメータ初期値計算手段で求めたトレンドμとボラティリティーσを初期値としてGalerkin有限要素法とBTD法を用いて株価確率分布の時間発展に関する数値解析を行い、さらに、確率の規格化条件により株価確率分布の補正を行う有限要素法数値計算手段と、前記有限要素法数値計算手段により数値計算して得られた対数株価に関する確率分布を基にして数値計算途中においてトレンドμとボラティリティーσを再評価し更新するパラメータ更新手段と、前記有限要素法数値計算手段により数値計算して得られた所定日数分の予想株価とボラティリティーσを出力するデータ出力手段とを具備することを特徴とする株価変動数値解析装置。
【請求項２】
前記パラメータ初期値計算手段は、トレンドμの初期値を過去の一定期間の株価から最小二乗法で計算する手段を備えたことを特徴とする請求項１記載の株価変動数値解析装置。
【請求項３】
形状関数を利用してLevy分布を数値計算内で再現し、対数株価確率分布の時間発展に関する数値解析を行うLevy分布有限要素法数値計算手段を備えたことを特徴とする請求項１記載の株価変動数値解析装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、株価変動数値解析装置に関し、特に、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、Galerkin有限要素法とBTD法を用いて、対数株価確率分布の時間発展を数値解析する株価変動数値解析装置に関する技術分野に属する。
【０００２】
【従来の技術】
経済や金融における現象は、その挙動が非常に複雑であることから、本質的構造を探ることや、予測を行うことは不可能であると考えられてきた。さらに、これらの現象は、限られた人々の経験と勘によって支配され、大多数にとって決して容易に関与できる領域ではないと捉えられてきた。しかしながら、昨今、経済や金融分野で起こる現象に対し、科学的な手法を用いてアプローチする研究が盛んに行われている。金融派生商品の価格評価方程式であるBlack-Scholesモデル以後、金融工学や数理ファイナンスにおける理論は急速に市場に浸透し、現在も進化を続けている。
【０００３】
幾何Brown運動で表現した株価モデルでは、株価変動における不確定要素をBrown運動による確率過程で表し、株価が変動傾向に従って動く確定的要素との直和として確率微分方程式で記述する。また、このモデルにより、株価変動は移流拡散型の現象として認識される。不確定要素に対するパラメータはボラティリティーと呼ばれている。ボラティリティーは、株価の一定期間における変動を年率換算したもので計測され、一般的にリスク指標として用いられる。また、確定的要素に対するパラメータは、トレンドと呼ばれており、株価変動の傾向成分を示すパラメータである。このトレンドとボラティリティーが、株価変動を扱う際の重要なパラメータとなる。
【０００４】
Brown運動を基礎とする確率過程は、その性質として拡散過程であることが知られている。このような確率過程に対して、確率密度関数に関する時間発展を記述するのがFokker-Planck方程式である。このFokker-Planck方程式は、確率の２重構造という特徴を有している。つまり、遷移過程において確率的であるのは勿論のこと、初期条件においても確率的である。また、Fokker-Planck方程式は、移流拡散型の方程式である。Fokker-Planck方程式は、確率密度関数の時間発展方程式であることから、確率過程を扱いながらも、確率統計的処理を施す必要がない。
【０００５】
確率微分方程式とは、ある物理量が確率的に時間変動する確率過程に対し、確率微分等を含んだ等式で表現し、その解として、確率過程を一意に決定付けできるものを呼ぶ。しかしながら、数学的に整備され体系化された議論が可能なのは、伊藤の確率微分方程式が唯一である。それを実現しているのが伊藤の補題であり、本来は微分不可能であるBrown運動に対して微分演算を可能としている。伊藤の確率微分方程式は、Brown運動を基礎とし、任意の確率変数z＝z(t)に対し、
dz＝g(z,t)dt＋h(z,t)dB
で表される。ただし、g(z,t)とh(z,t)は、確率変数ｚと時間ｔに従う任意の関数であり、dBは、Brown運動の微小変動である。また、伊藤の確率微分方程式に従う確率過程を伊藤過程と呼ぶ。伊藤の補題の詳細については、非特許文献２などを参照されたい。
【０００６】
確率変数zと時間ｔに関する任意の２変数関数ｆ＝f(z,t)をTaylor展開して、２次以上の微小項を無視すると、
df＝{(∂f/∂t)+(∂f/∂z)g(z,t)+1/2(∂2f/∂z2)h(z,t)2}dt
+(∂f/∂z)h(z,t)dB
＝(∂f/∂t)dt+(∂f/∂z)dz+(1/2)(∂2f/∂z2)h(z,t)2dt
となる。これは、伊藤の確率過程に従う確率変数zと時間ｔに関する任意の２変数関数ｆ＝f(z,t)が微分可能であるという伊藤の補題である。
【０００７】
株価をＳ＝S(t)としたとき、
dS＝μSdt+σSdB
が、幾何Brown運動で表現された株価モデルである。ただし、μはトレンドを表し、σはボラティリティーを表している。幾何Brown運動で表現された株価モデルは、伊藤の確率微分方程式に従う伊藤過程である。全ての株式銘柄に対して一般性を保つために、時刻ｔにおける株価S(t)を、初期時刻t＝0における初期株価S0＝S(0)を用いて、
x(t)≡S(t)/S0
と、規格化する。この規格化された株価ｘを用いて、株価モデルの式を書き直すと、
dx＝μxdt+σxdB
となる。
【０００８】
伊藤の補題を利用することで、規格化された株価ｘと時間ｔに関する任意の２変数関数ｆ＝f(x,t)に対して、
df＝{(∂f/∂t)+μx(∂f/∂x)+1/2(σx)2(∂2f/∂x2)}dt+σx(∂f/∂x)dB
の微分関係式が成り立つ。任意の２変数関数ｆを、
f(x,t)＝ln(x)
とすると、
d(ln(x))＝(μ-(1/2)σ2)dt+σdB
となる。さらに、
y(t)≡ln(x(t))  (x(t)≧0)
を、対数株価と呼ぶ。これらの式から、関係式
dy＝(μ-(1/2)σ2)dt+σdB
が得られる。この式から、対数株価ｙが、Brown運動に従うことがわかる。
【０００９】
本発明者らは、Black-Scholes方程式による金融派生商品の評価の正確さを高めるために、「金融派生商品評価システム」を提案した。特許文献１に開示されているように、金融派生商品の原資産のトレンドμをパラメータとして含む拡張Black-Scholes方程式を立て、その厳密解を求めて、記憶装置に格納しておく。入力手段から、金融市場における時系列価格データから求めた各統計データを入力し、記憶手段に格納する。価格演算手段で、各統計データに基づいて、厳密解から金融派生商品の価格を求め、出力手段から出力する。トレンド項を表現することによって、金融市場の短期トレンドを評価することが可能となり、従来のBlack-Scholes方程式よりも正確に経済物理の現象を捉えることができる。市場が現時刻において資産価格をどのように評価しているかを定量化できる。
【００１０】
【特許文献１】
特開2003-67581号公報
【非特許文献１】
棚橋隆彦著「流れの有限要素法解析I,II」（朝倉書店1997）
【非特許文献２】
保江邦夫著「数値確率解析入門」（朝倉書店2000）
【００１１】
【発明が解決しようとする課題】
しかし、従来の株価変動数値解析方法では、株価変動が移流拡散型の現象として認識されながら、移流を考慮しない拡散現象として解析されていたので、解析精度が不十分であるという問題があった。
【００１２】
本発明は、上記従来の問題を解決して、少量の株価データの入力のみで、株価の変動特性を自動的に推定し、株価変動を的確に解析して、株式投資における指針を示すことを目的とする。
【００１３】
【課題を解決するための手段】
上記の課題を解決するために、本発明では、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、Galerkin有限要素法とBTD法を用いて、対数株価確率分布の時間発展を数値解析する株価変動数値解析装置に、トレンドとボラティリティーの初期値を、過去の価格列から設定する手段と、数値計算内部において、株価変動の重要なパラメータであるトレンドとボラティリティーを更新する手段とを具備する構成とした。このように構成したことにより、長期間における株価変動の高精度な解析ができる。
【００１４】
【発明の実施の形態】
以下、本発明の実施の形態について、図１〜図４を参照しながら詳細に説明する。
【００１５】
（第１の実施の形態）
本発明の第１の実施の形態は、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、過去10日間の株価からトレンドとボラティリティーを計算して初期値とし、トレンドとボラティリティーを数値計算の途中で随時更新しながら、Galerkin有限要素法とBTD法を用いて株価確率分布の時間発展に関する数値解析を行う株価変動数値解析装置である。
【００１６】
図１（ａ）は、本発明の第１の実施の形態における株価変動数値解析装置の概念図である。図１（ｂ）は、株価変動数値解析装置で用いる解析格子の概念図である。図１（ａ）において、データ入力手段１は、直前の過去10日間の株価データを入力する手段である。データ記憶手段２は、直前の過去10日間の株価データを記憶する手段である。パラメータ初期値計算手段３は、直前の過去10日間の株価データから、パラメータμとσを求めて初期値とする手段である。有限要素法数値計算手段４は、Galerkin有限要素法とBTD法を用いて、株価確率分布の時間発展に関する数値解析を行う手段である。パラメータ更新手段５は、数値計算して得られた対数株価ｙに関する確率分布を基にし、トレンドμとボラティリティーσを再評価する手段である。データ出力手段６は、20日分の予想株価とボラティリティーσを出力する手段である。図２は、株価変動数値解析装置の動作手順を示す流れ図である。
【００１７】
上記のように構成された本発明の実施の形態における株価変動数値解析装置の機能と動作を説明する。最初に、図１を参照しながら、株価変動数値解析装置の動作の概略を説明する。株価変動を移流拡散型の現象として捉え、幾何Brown運動による株価モデルからFokker-Planck方程式を導く。直前の過去10日間の株価データを、データ入力手段１から入力して、データ記憶手段２に格納する。パラメータ初期値計算手段３で、直前の過去10日間の株価データから、トレンドμとボラティリティーσを求めて初期値とする。初期トレンド決定方法に最小二乗法を利用する。有限要素法数値計算手段４で、Fokker-Planck微分方程式を、Galerkin有限要素法とBTD法を用いて解き、対数株価確率分布の時間発展を数値解析する。さらに、確率の規格化条件により補正を行う。その結果に基いて、パラメータ更新手段５で、株価変動の重要なパラメータであるトレンドμとボラティリティーσを求めて更新する。これを20日分だけ繰り返す。データ出力手段６で、20日分の予想株価とボラティリティーσを出力する。
【００１８】
このように、株価変動を移流拡散型の現象として捉え、移流を考慮し、移流拡散方程式により数値解析して、数値計算途中でパラメータを再評価し更新することにより、少ない参照期間で長期間の推定が可能となる。過去10日分の株価データだけの少量の入力のみで、20日分の予想株価確率分布が簡便に得られるので、時々刻々と変化する株式市場を、的確かつ長期に亘り捉えることができ、予想株価のリスクも、ボラティリティーσから分かる。
【００１９】
次に、対数株価モデルからFokker-Planck方程式を導く方法を説明する。対数株価ｙが直接Brown運動に従うことに着目し、対数株価ｙで表した対数株価モデルからFokker-Planck方程式を導く。Fokker-Planck方程式を、Master方程式を利用して導出する。任意の確率過程ｚ＝z(t)の時間変化について考える。時刻ｔ0で値ｚ0の状態から、時刻ｔに値ｚとなる条件つき確率をP(z,t|z0,t0)で表すと、確率の規格化条件より、
∫-∞∞P(z,t|z0,t0)dz＝1
が成立する。この確率過程zが、Markov過程である場合を考える。単位時間あたりに、状態ｚから状態ｚ0へ遷移する確率として、遷移確率W(ｚ→ｚ0)を導入する。この遷移確率を用いて、Markov過程における条件つき確率の時間発展を記述した方程式が、Master方程式
(∂/∂t)P(z,t)＝-∫-∞∞W(z→z')P(z,t)dz'+∫-∞∞W(z'→z)P(z',t)dz'
である。ただし、今後は簡単のため、初期時刻ｔ0における初期状態ｚ0による条件つき確率P(z,t|z0,t0)は、単にP(z,t)と書くこととする。
【００２０】
このMaster方程式は、時刻ｔにおいて状態ｚとなる確率分布が、他の状態への遷移によって減少し、他の状態からの遷移によって増加することを意味している。さらに、このMaster方程式は、Kramers-Moyal展開により、偏微分方程式で書き表すことができる。遷移確率のｋ次モーメントCk(z)を、
Ck(z)≡∫-∞∞rkW(z→z+r)dr
で定義する。Master方程式は、
(∂/∂t)P(z,t)＝Σk=1∞{(1/k!)(-∂/∂z)k×Ck(z)P(z,t)}
で表すことができる。これがKramers-Moyal展開である。このKramers-Moyal展開を２次までで留めた形が、Fokker-Planck方程式
(∂/∂t)P(z,t)＝-(∂/∂z)C1(z)P(z,t)+(1/2)(∂2/∂z2)C2(z)P(z,t)
である。
【００２１】
株価モデルに、Fokker-Planck方程式を適用する。
dy＝(μ-(1/2)σ2)dt+σdB
に対する遷移確率の１次モーメントを計算すると、
C1(y)＝limΔt→0（1/Δt）E[Δy]＝μ-(1/2)σ2
が得られる。遷移確率の２次モーメントは、
C2(y)＝limΔt→0(1/Δt)E[(Δy)2]＝σ2
となる。３次以上のモーメントは零となるため、対数株価モデル式に対するFokker-Planck方程式は、
(∂/∂t)P(y,t)＝-(∂/∂y)(μ-(1/2)σ2)P(y,t)+(1/2)(∂2/∂y2)σ2P(y,t)
となる。
【００２２】
第３に、株価変動パラメータについて説明する。対数株価モデル
dy＝(μ-(1/2)σ2)dt+σdB
におけるボラティリティーσは、金融市場において値動きの標準偏差で計測され、一般的にリスクの概念となっており、非常に重要なパラメータとなっている。簡単のため、トレンドμとボラティリティーσが、空間領域ｙ全体で一定の値とする仮定を導入する。この仮定は、実際の市場における動きを想定することで、十分に妥当な仮定である。
【００２３】
トレンドμとボラティリティーσの初期値の算出方法について説明する。ヒストリカル・ボラティリティーによるボラティリティーσの算出方法を採用する。株価の変動傾向を測る方法として最小二乗法を援用し、トレンドμを算出する。過去の価格列から得られる情報は、新しいものほどその重要性は高い。したがって、過去の価格列を一律に扱うべきではない。直近の情報を重視するために、時間間隔が等しい価格列に対して重みを使用し、トレンドμとボラティリティーσの初期値を決定する。
【００２４】
トレンドμとボラティリティーσは、時間に関して変化するパラメータである。数値計算して得られた対数株価ｙに関する確率分布を基にし、トレンドμとボラティリティーσを再評価する。さらに、パラメータ更新時間ステップｍとパラメータ更新時間刻み幅Δｔmに従い時間進行することで、トレンドμとボラティリティーσを更新する。Fokker-Planck方程式の数値計算結果を受けて、トレンドμを随時更新していく。式
μm+1＝(1/2)(σm+1)2＋(1/Δtm)(Ｅ[ｙ]m+1−Ｅ[ｙ]m)
により、トレンドμを更新する。パラメータ更新時間ステップm+1におけるボラティリティーσm+1は、ボラティリティーσ更新式を用いて算出された値である。ボラティリティーσの更新式は
σm+1＝√(Πm/(2Δtm))
である。ただし、係数Πは、
Πm≡(Ｅ[ｙ]m)2＋(Ｅ[ｙ]m+1)2−2Ｅ[ｙ]mＥ[ｙ]m+1−2Var[ｙ]m＋2Var[ｙ]m+1
である。これらの式における肩付のｍなどは、時間ステップを示す添え字であり、べき乗の意味ではない。以下の同様な式においても同じである。
【００２５】
第４に、Fokker-Planck方程式を有限要素法で解く方法を説明する。有限要素法の詳細については、非特許文献１などを参照されたい。対数株価ｙにより表現された対数株価モデル式
dy＝(μ−(1/2)σ2)dt＋σdB
を基礎方程式とし、Fokker-Planck方程式
(∂/∂t)Ｐ(y,t)＝−(∂/∂y)(μ−(1/2)σ2)Ｐ(y,t)＋(1/2)(∂2/∂y2)σ2Ｐ(y,t)
を数値計算することで、対数株価確率分布の時間発展を追う。ただし、時間ｔに対数株価ｙとなる確率密度関数がＰ(y,t)であり、以後簡単のためＰと表示する。この式において、トレンドμとボラティリティーσは全空間領域で一定の仮定を用いることから、
(∂/∂t)Ｐ＝−ｕ(∂/∂y)Ｐ＋ν(∂2/∂y2)Ｐ
のように書き直せる。また、移流速度ｕと拡散係数νを、
ｕ≡μ−(1/2)σ2
ν≡(1/2)σ2
のように定義する。このFokker-Planck方程式を、数値解析における支配方程式として数値計算を行う。
【００２６】
また、確率過程を扱うことにより、
∫-∞∞Ｐ(ｙ,ｔ)dy＝１
に示す確率の規格化条件を満たす必要も生ずる。さらに、数値計算することで得られる対数株価確率分布から計算される期待値Ｅ[y]や分散Var[y]を用い、パラメータ更新時間ステップｍとパラメータ更新時間刻み幅Δtmに従い、トレンドμとボラティリティーσを、それぞれ更新式
μm+1＝(1/2)(σm+1)2＋(1/Δtm)(Ｅ[ｙ]m+1−Ｅ[ｙ]m)
σm+1＝√(Πm/(2Δｔｍ))
により時間進行し、更新を行う。前進Euler法を用いて、支配方程式の時間進行を行う。対数株価モデルから導出したFokker-Planck方程式を前進Euler法により時間進行することで、
(Ｐn+1−Ｐn)/Δｔ＝−ｕn(∂/∂y)Ｐn＋νn(∂2/∂y2)Ｐn
を得る。これらの式における肩付のｎなどは、時間ステップを示す添え字であり、べき乗の意味ではない。以下の同様な式においても同じである。
【００２７】
空間の離散化に関しては、Galerkin有限要素法を採用する。解析格子は、例として図１（ｂ）に示すように、指数関数を用いた不等分割格子である。支配方程式を前進Euler法により時間進行した式に対し、重みδωを乗じ、空間領域全体Ωで積分する。さらに全領域Ωを多数の要素に分割し、１つの要素Ωeに着目すると、
∫Ωeδω((Ｐn+1−Ｐn)/Δｔ)dΩ
＝−ｕn∫Ωeδω(∂/∂ｙ)ＰndΩ＋νn∫Ωeδω(∂2/∂ｙ2)ＰndΩ
となる。この式において、右辺最終項に対して部分積分を行い、Greenの定理を適用すると、
∫Ωeδω((Ｐn+1−Ｐn)/Δｔ)dΩ
＝−ｕn∫Ωeδω(∂/∂ｙ)(Ｐn)dΩ
＋νn∫Γeδω(∂/∂ｙ)(Ｐn)dΓ
−νn∫Ωe(∂/∂ｙ)(δω)(∂/∂ｙ)(Ｐn)dΩ
となる。要素内において確率密度関数Ｐを双１次関数で補間すると、
Ｐ＝ＮαＰα
の関係を得る。ここで、αは要素内での局所節点番号、Ｎαは形状関数であり、αについて総和の規約が適用される。
【００２８】
この形状関数を用いて、重み関数δωを、
δω＝Ｎαδωα
のように補間する。δωαの各成分は任意の定数であり、基本境界条件Γ1ではδωαの成分は零である。δωαの各成分の任意性を用いると、
∫ΩeＮαＮβdΩ(Ｐβn+1−Ｐβn)/(Δｔ)
＝−ｕn∫ΩeＮα(∂/∂ｙ)(Ｎβ)dΩＰβn
＋νn∫ΓeＮα(∂/∂ｙ)(Ｐn)dΓ
−νn∫Ωe(∂/∂ｙ)(Ｎα)(∂/∂ｙ)(Ｎβ)dΩＰβn
を導くことができる。さらに、係数行列を用い、質量の集中化を行うと、離散化式
Ｍ-αβ(Ｐβn+1−Ｐβn)/(Δｔ)
＝−ｕnＡαβＰβn−νnＤαβＰβn＋Λα
を得る。ただし、各係数行列は、
Ｍ-αβ＝∫ΩeＮαＮβdΩ（質量行列）
Ａαβ＝∫ΩeＮα(∂/∂ｙ)(Ｎβ)dΩ（移流行列）
Ｄαβ＝∫Ωe(∂/∂ｙ)(Ｎα)(∂/∂ｙ)(Ｎβ)dΩ（拡散行列）
で定義される。Ｍ-αβは集中対角化された質量行列を、Λαは境界積分項を表す。
【００２９】
２次の時間精度とするために、２階の時間微分項を、移流方程式で評価する。これは、CFDにおけるBTD法に対応する。まず、支配方程式において、移流項のみを考慮した移流方程式が、
(∂/∂ｔ)(Ｐ)＝−ｕ(∂/∂ｙ)(Ｐ)
で表される。確率密度関数Ｐに対して、着目時間ステップまわりにTaylor展開を行うことで、
(∂/∂ｔ)(Ｐ)|n＝(Ｐn+1−Ｐn)/Δｔ
−(Δt/2)(∂2/∂t2)(Ｐ)|n＋Ｏ(Δt2)
を得ることができる。この式を用いて、前進Euler法による時間進行に加えて、２階の時間微分項を考慮することで、時間精度を２次とすることができる。ここで、２階の時間微分項を、移流方程式で評価することで、
(∂2/∂t2)(Ｐ)＝−(dｕ/dｔ)(∂/∂y)(Ｐ)＋ｕ2(∂2/∂y2)(Ｐ)
を得る。ただし、トレンドμとボラティリティーσを時間変数としているため、移流速度ｕも、定義式より時間変数と見なせる。
【００３０】
よって、移流速度ｕの時間微分を、
(du/dt)≒Δu/Δt＝(un−un-1)/Δt
により近似する。したがって、時間ステップｎでのBTD項は、
BTDn＝(Δt/2){−(Δu/Δt)(∂/∂y)(Pn)＋(ｕn)2(∂2/∂ｙ2)(Ｐn)}
のように求めることができる。BTD項に対して、Galerkin有限要素法を適用し、係数行列を用いて表現すると、
BTDn＝(Δｔ/2){−(Δｕ/Δｔ)ＡαβＰβn−(ｕn)2ＤαβＰβn}
＋(Δｔ/2)νn∫ΓeＮα(∂/∂ｙ)(Ｐn)dΓ
となる。離散化式にBTD項を加えると、離散化式
Ｍ-αβ(Ｐβn+1−Ｐβn)/(Δｔ)
＝−(ｕn＋(1/2)Δｕ)ＡαβＰβn
−{νn＋(1/2)Δｔ(ｕn)2}ＤαβＰβn＋Λα
となる。このBTD項も含めた離散化式に従い、数値計算を行う。ただし、境界積分項はまとめてΛαとして表示している。
【００３１】
確率過程を扱う上で、確率の規格化条件を満たすことは非常に重要である。そのため、数値計算を行う際にも、計算領域内で総確率量を保存させ、確率の規格化条件を満足させなければならない。境界条件は、計算領域の両端で初期状態のまま固定させる。つまり、基本境界条件Γ1を採用していることから、Galerkin有限要素法の性質から境界積分項を計算する必要は生じない。また、計算領域をある程度広く取ることで、境界を固定する方法が、数値計算に与える影響を最小限度に抑えることができる。確率変数を数値解析する際には、確率の規格化条件を満足させながら数値計算を進める必要がある。以下の方法により総確率量を保存させている。
【００３２】
任意の時間ステップ後における状態を考える。各要素での確率量をＰeとすると、計算領域全体での総確率量ＰΩは、
ＰΩ＝∫ΩＰedΩ
により求まる。したがって、確率の規格化条件より、確率の誤差量Pεは、
Ｐε＝１−ＰΩ
と求まる。この確率の誤差量Pεを、各要素における確率量に応じて、計算領域全体で補正する。したがって、確率の規格化条件を適用した後の各要素における確率量Ｐe*は、
Ｐe*＝Ｐe＋Ｐε(Ｐe/ＰΩ)
として計算できる。以上の一連の計算を毎ステップ計算後に行い、確率の規格化条件を満足させている。
【００３３】
陽解法で離散化された非定常移流拡散方程式に対して、安定的に数値計算を進めるためには、Courant数と拡散数の制約を受ける。移流速度ｕ、拡散係数νに対して、Courant数をｃ、拡散数をｄと書くと、各条件はそれぞれ、
ｃ＝|ｕΔｔ/Δｙ|≦１
ｄ＝νΔｔ/(Δｙ)2≦1/2
で与えられる。したがって、離散化式を数値計算する際には、この条件を満たす必要がある。
【００３４】
第５に、図２を参照しながら、数値解析の手順を説明する。解析対象は、日次変動である。初期パラメータを決定する際に利用する価格列{Ｓi}等は、日次データとなっている。したがって、数値計算内におけるパラメータ更新は、１日毎としている。まず、ステップ１において、計算開始時刻以前の過去の価格列を入力する。ステップ２において、計算開始時刻以前の過去の価格列から、数値解析に用いる初期パラメータを計算する。その後、初期条件を設定して、時間ステップｎの繰り返し計算を開始する。繰り返し計算による時間進行の過程は、はじめに、ステップ３において、Galerkin有限要素法とBTD法により、離散化した支配方程式を計算する。得られた確率密度関数を、ステップ４において、確率の規格化条件により補正を行う。ステップ５において、設定したパラメータ更新時間刻み幅Δtm毎の計算が終了したか検査する。終了していれば、ステップ６において、ボラティリティーσを、更新式に基づき再計算する。ステップ７において、トレンドμを、更新式に基づき再計算する。この一連の計算を繰り返すことで時間進行し、指定の時刻まで計算を行う。
【００３５】
上記のように、本発明の第１の実施の形態では、株価変動数値解析装置を、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、過去10日間の株価からトレンドとボラティリティーを計算して初期値とし、トレンドとボラティリティーを数値計算の途中で随時更新しながら、Galerkin有限要素法とBTD法を用いて株価確率分布の時間発展に関する数値解析を行う構成としたので、少量の株価データの入力のみで、株価の変動特性を自動的に推定し、株価変動を的確に解析して、株式投資における指針を示すことができる。
【００３６】
この株価変動数値解析装置により、株価の変動特性を推定し、株式投資における指針を示すことができる。確率分布の時間発展とトレンドμとボラティリティーσの推移を得ることで、具体的な投資活動に対して各種の情報を示すことができる。株価変動を推定することにより、将来における株価の確率分布を得ることができ、ある一定の範囲内にどの程度の確率で株価が収まるかを示すことができる。株価変動におけるリスクの推定や管理に役立つ。
【００３７】
（第２の実施の形態）
本発明の第２の実施の形態は、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、過去10日間の株価からトレンドとボラティリティーを計算して初期値とし、トレンドとボラティリティーを数値計算の途中で随時更新しながら、Galerkin有限要素法とBTD法を用いて、形状関数を利用してLevy分布を数値計算内で再現し、株価確率分布の時間発展に関する数値解析を行う株価変動数値解析装置である。
【００３８】
図３（ａ）は、本発明の第２の実施の形態における株価変動数値解析装置の概念図である。図３（ｂ）は、株価変動数値解析装置で用いる形状関数の概念図である。図３（ａ）において、Levy分布有限要素法数値計算手段７は、Galerkin有限要素法とBTD法を用いて、形状関数を利用してLevy分布を数値計算内で再現し、株価確率分布の時間発展に関する数値解析を行う手段である。その他は、第１の実施の形態と同じである。図４は、株価変動数値解析装置の動作手順を示す流れ図である。
【００３９】
上記のように構成された本発明の第２の実施の形態における株価変動数値解析装置の機能と動作を説明する。最初に、図３を参照しながら、株価変動数値解析装置の動作の概略を説明する。株価変動を移流拡散型の現象として捉え、幾何Brown運動による株価モデルからFokker-Planck方程式を導く。直前の過去10日間の株価データを、データ入力手段１から入力して、データ記憶手段２に格納する。パラメータ初期値計算手段３で、直前の過去10日間の株価データから、トレンドμとボラティリティーσを求めて初期値とする。初期トレンド決定方法に最小二乗法を利用する。ここまでは、第１の実施の形態と同じである。Levy分布有限要素法数値計算手段７で、Fokker-Planck微分方程式を、形状関数を利用してLevy分布を数値計算内で再現し、Galerkin有限要素法とBTD法を用いて解き、対数株価確率分布の時間発展を数値解析する。さらに、確率の規格化条件により補正を行う。その結果に基いて、パラメータ更新手段５で、株価変動の重要なパラメータであるトレンドμとボラティリティーσを求めて更新する。これを20日分だけ繰り返す。データ出力手段６で、20日分の予想株価とボラティリティーσを出力する。
【００４０】
図３（ｂ）に示す形状関数を利用することで、Levy分布を数値計算内で再現する。幾何Brown運動によるモデルでは、株価リターンの分布が正規分布により与えられることを前提としている。しかしながら、実際の市場においては、正規分布が確認されることは稀であり、その分布は、正規分布と比較して「ハイピーク・ファットテール」形状に分布することが知られている。またその分布は、Levy分布に従う。Levy分布では、分布の中央部分で正規分布よりも尖度が高く、分布の裾野部分が厚くなるといった特徴を持っている。これはつまり、変動幅が小さい株式においては、正規分布と比較して変動が抑制され、逆に変動幅が大きい株式については、正規分布と比較して変動が促進される現象と捉えることができる。
【００４１】
この抑制と促進の相反する力が働く境界にあたる点が、期待値を中心とした±１σの値である。なぜなら，ボラティリティーσは、株価変動の標準偏差で与えられ、正規分布において期待値を中心とした±１標準偏差の値は、分布関数の変曲点となっており、その内側（外側）では、分布を上（下）に凸とする方向に力が働いている。すなわち、Levy分布では、これらの力が各点でより一層強まることで、「ハイピーク・ファットテール」形状の分布が形成される。これらのLevy分布を形成する要因を株価モデルへ導入することで、より高精度に株価変動を捉えることができる。
【００４２】
数値解析手法として、Galerkin有限要素法を採用している。このGalerkin有限要素法における特徴のひとつが、要素内補間と重み関数に使用する形状関数である。この形状関数を利用し、Levy分布を数値計算内で再現する。Levy分布では、期待値を中心とした±１σ内外において、それぞれ別方向の力が働き、その分布を形作っている。期待値を中心として±１σの内側においては、中心方向へ重みを強くする形状関数を選択し、±１σの外側では、逆に外部方向への重みを強くする形状関数を用いて、重み関数を構成することで、Levy分布を再現する。
【００４３】
Ｎα＝(1/2)(１＋ξαξ)＋γξα
で表される形状関数を用いる。Levy分布を再現するための形状関数の形状を、図３（ｂ）に示す。この形状関数は、有限要素法における上流化で用いられる重み関数のひとつと同一の形状である。γが正（負）であれば、図３（ｂ）において、左側（右側）の重みが大きくなる性質を持つ。この形状関数を用いて、期待値を中心とした±１σとの位置関係により、適宜γの値を調節し、有限要素法における重み関数を構成することで、数値計算内におけるLevy分布の再現ができる。
【００４４】
この形状関数を重み関数として使用し、数値解析の支配方程式を、有限要素法とBTD法により離散化した場合、離散化式はすでに示した形と同等であるが、係数行列は、
Ｍαβ＝∫ΩeＮαＮβdΩ＝Ｊe((1/6)ξαξβ＋(1/2)＋γξα)（質量行列）
Ａαβ＝∫ΩeＮα(∂/∂ｙ)ＮβdΩ＝((1/2)＋γξα)ξβ（移流行列）
Ｄαβ＝∫Ωe(∂/∂ｙ)(Ｎα)(∂/∂ｙ)(Ｎβ)dΩ
＝(1/Ｊe)((1/2)ξαξβ)（拡散行列）
となる。Jeは、要素平均で求めるJacobianである。
【００４５】
図４を参照しながら、数値解析の手順を説明する。解析対象は、日次変動である。初期パラメータを決定する際に利用する価格列{Ｓi}等は、日次データとなっている。したがって、数値計算内におけるパラメータ更新は、１日毎としている。まず、ステップ11において、計算開始時刻以前の過去の価格列を入力する。ステップ12において、計算開始時刻以前の過去の価格列から、数値解析に用いる初期パラメータを計算する。その後、初期条件を設定して、時間ステップｎの繰り返し計算を開始する。繰り返し計算による時間進行の過程は、はじめに、ステップ13において、形状関数を利用してLevy分布を数値計算内で再現し、Galerkin有限要素法とBTD法により、離散化した支配方程式を計算する。得られた確率密度関数を、ステップ14において、確率の規格化条件により補正を行う。ステップ15において、設定したパラメータ更新時間刻み幅Δtm毎の計算が終了したか検査する。終了していれば、ステップ16において、ボラティリティーσを、更新式に基づき再計算する。ステップ17において、トレンドμを、更新式に基づき再計算する。この一連の計算を繰り返すことで時間進行し、指定の時刻まで計算を行う。
【００４６】
上記のように、本発明の第２の実施の形態では、株価変動数値解析装置を、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、過去10日間の株価からトレンドとボラティリティーを計算して初期値とし、トレンドとボラティリティーを数値計算の途中で随時更新しながら、Galerkin有限要素法とBTD法を用いて、形状関数を利用してLevy分布を数値計算内で再現し、株価確率分布の時間発展に関する数値解析を行う構成としたので、より高精度に株価変動を捉えることができ、少量の株価データの入力のみで、株価の変動特性を自動的に推定し、株価変動を的確に解析して、株式投資における指針を示すことができる。
【００４７】
【発明の効果】
以上の説明から明らかなように、本発明では、株価変動を移流拡散型の現象として捉え、幾何Brown運動に基づく株価モデルからFokker-Planck方程式を導き、Galerkin有限要素法とBTD法を用いて、対数株価確率分布の時間発展を数値解析する株価変動数値解析装置に、トレンドとボラティリティーの初期値を、過去の価格列から設定する手段と、数値計算内部において、株価変動の重要なパラメータであるトレンドとボラティリティーを更新する手段とを具備する構成としたので、長期間における株価変動の高精度な解析ができる。
【図面の簡単な説明】
【図１】本発明の第１の実施の形態における株価変動数値解析装置の概念図、
【図２】本発明の第１の実施の形態における株価変動数値解析装置の動作手順を示す流れ図、
【図３】本発明の第２の実施の形態における株価変動数値解析装置の概念図、
【図４】本発明の第２の実施の形態における株価変動数値解析装置の動作手順を示す流れ図である。
【符号の説明】
１  データ入力手段
２  データ記憶手段
３  パラメータ初期値計算手段
４  有限要素法数値計算手段
５  パラメータ更新手段
６  データ出力手段
７  Levy分布有限要素法数値計算手段
【図１】
0004461245.tif 000002
【図２】
0004461245.tif 000003
【図３】
0004461245.tif 000004
【図４】
0004461245.tif 000005