有限差分オプション価格評価：符号、安定性、障害条件、収束

時間の符号、離散作用素、境界・イベント条件、アメリカン型の相補性、ダンピング、価格とグリークスの収束を整合させ、監査可能な有限差分ソルバーを構築します。

教育目的のみであり、投資助言ではありません。投資により損失が生じる可能性があります。

直接の答え

有限差分オプション価格評価は、連続的な評価方程式を価格と時間のノードに置き換えます。有効な実装には細かいグリッドだけでなく、暦時間と残存期間の符号、空間作用素、終端・境界条件、イベント時のジャンプ、時間離散化、アメリカン型の障害問題ソルバー、残差、収束管理のすべてに整合性が必要です。

グリッドの細分化が減らすのは数値誤差だけです。誤ったペイオフ、確率モデル、配当処理、決済請求権、市場入力、執行仮定は修正できません。ソルバーは誤答へ高精度に収束することもあります。

コーディングの前にグリッドを導出する

請求権、ペイオフ Phi(S)、行使方式、決済、金利、借株費用、配当、ボラティリティ・モデル、座標系を確定します。暦時間 t では Black-Scholes PDE は V_t + L V = 0、残存期間 tau = T - t では V_tau = L V となり、V(S, 0) = Phi(S) から開始します。
価格または対数価格グリッド、切断領域、時間グリッドを選び、スポット、権利行使価格、バリア、配当、観測時点を整列させます。N が区間数かノード数かを明記し、S_max は Delta S と分けて検証します。
一様グリッド S_i = i Delta S 上で a_i = 0.5[sigma^2 i^2 - (r - q)i]、b_i = -(sigma^2 i^2 + r)、c_i = 0.5[sigma^2 i^2 + (r - q)i] を導出し、(L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1) とします。
時間ステップを厳密に定義します。陽解法は V^(n+1) = (I + Delta tau L_h)V^n、陰解法は (I - Delta tau L_h)V^(n+1) = V^n、Crank-Nicolson 法は (I - 0.5 Delta tau L_h)V^(n+1) = (I + 0.5 Delta tau L_h)V^n です。境界項を正しい右辺と時間層に組み込みます。
係数の符号を領域全体で確認し、CFL型制約、正値性、単調性、行列特性、振動を検査します。陽解法と後退 Euler 法はいずれも通常は時間一次精度です。陰解法は散逸性が強いものの、常に精度が低いわけではありません。Crank-Nicolson 法の線形安定性は、単調性や折れ点周辺の無振動性を保証しません。不連続点ごとに Rannacher ダンピングを使うか、その方法も含めて記録します。
アメリカン型の行使には、PSOR、ペナルティ法、仕様を明記した作用素分割法など、名称を明示した離散線形相補性解法を使用します。V - Phi >= 0、A V - b >= 0、(V - Phi) elementwise_times (A V - b) = 0 を満たし、緩和係数、反復上限、相補性残差を報告します。
スポットで価格とグリークスを補間し、領域、空間、時間、イベント整列、ソルバー許容誤差を細分化します。ヨーロピアン型は公式、アメリカン型は深い二項ツリーまたは独立ソルバーと比較し、価格、Delta、Gamma、自由境界、無裁定検査、残差を記録します。

境界条件は契約に固有です。標準的な仮定の下で、ヨーロピアン・プットには V(0, tau) = K exp(-r tau) を用いることができ、S が大きくなるとゼロへ近づきます。アメリカン・プットには通常 V(0, tau) = K を用います。配当付きヨーロピアン・コールの高価格側漸近式はこれとは異なります。離散現金配当には、V(S, t-) = V(max(S - D, 0), t+) のようなジャンプと補間が必要で、その後に正しい順序で行使判断を行います。

時間ステップが小さくても、中央差分で移流係数が負になる場合があります。その場合は風上差分、フィット型、対数価格離散化が必要になり得ます。拡散 PDE に価格ジャンプが自動的に加わることもありません。モデル、パラメータ、較正、数値、執行可能価格の誤差は分けて扱います。

4つの計算例

時計を変えると時間の符号も変わります。 V_t + L V = 0 と tau = T - t から、連鎖律により V_tau = L V を得ます。したがって、tau = 0 の満期ペイオフから大きな tau へ進むときは正の L 作用素を使います。暦時間で後ろ向きに進む際の負号をそのまま移すと、別の問題を解くことになります。
一つのノードだけなら局所的に安全に見えることがあります。 sigma = 25%、r = 4%、q = 0、i = 95、Delta tau = 0.0005 とすると、a_i = 280.13125、b_i = -564.1025、c_i = 283.93125 です。陽解法の重みは 0.140065625、0.71794875、0.141965625 で、合計は 0.99998 です。隣接値 6, 5, 4 は 4.998 に更新されます。障害値が 5 のアメリカン・プットでは陽解法の射影により 5 を残しますが、陰解法や Crank-Nicolson 法では結合した LCP 残差の確認がなお必要です。
制約は領域全体で満たす必要があります。 同じ入力で i = 200 とすると、陽解法の重みは 0.623、-0.25002、0.627 です。スポット近傍のノードが妥当に見えても、中央の負の重みにより単調性が失われます。ステップを Delta tau = 0.0002 に縮めると、重みは 0.2492、0.499992、0.2508 となり、すべて非負で合計は 0.999992 です。
価格とグリークスには独立した管理が必要です。 ヨーロピアン・プットを S = $95、K = $100、T = 0.5、sigma = 25%、r = 4%、q = 0 とすると、Black-Scholes は d_1 = -0.088633257689、d_2 = -0.265409952985、価格 $8.413258404601、Delta -0.535313311861、Gamma 0.023662224614 を与えます。グリッドが3つの出力へ個別に収束してからでなければ、アメリカン・プレミアムや自由境界を信頼できません。

7段階の管理で起こり得る不備

暦時間 t と残存期間 tau を整合しない符号で進めることがあります。
日・年、金利、ボラティリティ、配当の単位を混同することがあります。
ペイオフ、行使方式、決済請求権、終端層を誤って実装することがあります。
連続配当利回りと離散現金配当を二重計上したり、処理順序を誤ったりすることがあります。
S_max がスポットに近すぎることがあります。これは Delta S が小さく見える場合も同じです。
低価格・高価格の境界式が別の請求権やキャリー前提のものになっていることがあります。
配当、バリア、フィキシング、観測が時間ノード間に入り、補間を誤ることがあります。
ノード数と区間数、または非一様グリッドの微分式にオフバイワンが生じることがあります。
ドリフトにより中央空間係数が負になり、単調性が損なわれることがあります。
陽解法の重みがスポットから離れた位置で全域安定性または正値性を満たさないことがあります。
Crank-Nicolson 法がペイオフの折れ点、バリア、イベント時のジャンプ周辺で振動することがあります。
Rannacher その他のダンピングを省略したり、誤ったイベントに適用したり、記録しなかったりすることがあります。
価値曲面が滑らかに見えても、線形系の残差が大きいことがあります。
PSOR の緩和係数、停止許容誤差、初期値、反復上限が不安定なことがあります。
一度の無制約解にノードごとの最大値を適用するだけで、厳密な陰解法または Crank-Nicolson LCP を解いたと誤認することがあります。
スポットでの補間が価値、Delta、Gamma、自由境界を偏らせることがあります。
価格が収束して見えても、Gamma と行使境界にはノイズが残ることがあります。
滑らかな漸近領域に入る前に、観測された収束次数や Richardson 外挿を使うことがあります。
拡散、ローカル・ボラティリティ補間、ジャンプ、パラメータ、較正が構造的に誤っていることがあります。
理論価値を執行可能な買い気配、売り気配、ヘッジ、裁定機会として示すことがあります。

よくある誤解

「グリッドを細かくすれば価格が正しいと証明できる。」契約もモデルも証明できません。
「陰解法は本質的に陽解法より精度が低い。」どちらの Euler 法も通常は時間一次精度で、安定性と誤差特性が異なります。
「Crank-Nicolson 法は無条件安定だから、単調で振動しない。」これらは別々の性質です。
「アメリカン型の価値は、無制約解を一度求めて各点で最大値を取れば得られる。」結合した陰解法系には LCP 解法または定量化された近似が必要です。
「一つの気配値に一致すればソルバーとグリークスを検証できる。」較正と市場ノイズは数値誤差や構造誤差を隠すことがあります。