資本市場三体力学 SDE 方程式系の導出
2026年2月7日
『資本市場の三体力学仮説』に基づき、本稿では完全な確率微分方程式(SDE)系を導出し、原文のすべての定性的制約を満たすかどうかを逐一検証します。
一、状態変数の定義
システムには2つの時間スケールがあります:
速い変数(秒から日):
- x(t)=ln(S/S∗) — 対数プレミアム(S は価格、S∗ は内在価値)
- p(t) — モメンタム(価格変化率の状態変数)
- σ(t) — 瞬間ボラティリティ
遅い変数(週から年):
- m(t) — モメンタム資本の規模
- v(t) — バリュー資本の規模
- l(t) — 流動性資本の規模
派生変数:
- δ=(S−S∗)/S∗=ex−1≈x(x が小さい場合)
- μ=dS/dt(SDE ではドリフト項として表現)
二、核心制約の抽出
原文から形式化可能な制約を抽出:
| 番号 |
制約 |
出典 |
| C1 |
M の価格への正帰還:dSd(ポジション)>0 |
モメンタム資本の定義 |
| C2 |
V の価格への負帰還:dSd(ポジション)<0 |
バリュー資本の定義 |
| C3 |
L の価格への無方向帰還:dSd(ポジション)≈0 |
流動性資本の定義 |
| C4 |
M↑→σ↑ |
正帰還ループ |
| C5 |
σ↑→L↓(流動性の退避) |
正帰還ループ |
| C6 |
L↑→ 価格ショック↓ |
L→M 制約によるショック |
| C7 |
V は S<S∗ で買い、S>S∗ で売る |
V のアンカー機構 |
| C8 |
$V\uparrow \to |
S - S^* |
| C9 |
M の収益 ∝μ、リスク ∝δ、コスト ∝σ |
収益行列 |
| C10 |
V の収益 ∝δ、リスク ∝σ、コスト ∝μ |
収益行列 |
| C11 |
L の収益 ∝σ、リスク ∝μ、コスト ∝δ |
収益行列 |
| C12 |
資本過剰 → 収益率低下 → 規模縮小 |
自然選択機構 |
三、SDE 方程式系
速い変数サブシステム
(I) 対数プレミアム:
dx=p⋅dt+lβσ⋅dW1
(II) モメンタム:
dp=l1[αM⋅m⋅p−αV⋅v⋅x−γ⋅p]dt+η⋅dW2
整理すると:
dp=lαMm−γ⋅p⋅dt−lαVv⋅x⋅dt+η⋅dW2
(III) ボラティリティ:
dσ=κσ(σˉ−σ)⋅dt+λMm∣p∣⋅dt−λVv⋅dt−λLl⋅dt+ξσ⋅dW3
遅い変数サブシステム(ε≪1)
(IV) モメンタム資本:
dm=ε⋅m⋅[aM∣p∣−bM∣x∣−cMσ−ρMm]dt+εσmm⋅dW4
(V) バリュー資本:
dv=ε⋅v⋅[aV∣x∣−bVσ−cV∣p∣−ρVv]dt+εσvv⋅dW5
(VI) 流動性資本:
dl=ε⋅l⋅[aLσ−bL∣p∣−cL∣x∣−ρLl]dt+εσll⋅dW6
ここで W1,…,W6 は独立な標準ブラウン運動です。
四、パラメータ表
| パラメータ |
意味 |
| αM |
M の正帰還強度(順張りの積極性) |
| αV |
V の負帰還強度(価格回帰の強さ) |
| γ |
自然減衰(摩擦、情報減衰) |
| β |
流動性による価格ショック緩衝指数 |
| κσ |
ボラティリティの平均回帰速度 |
| σˉ |
ボラティリティの長期平均 |
| λM,λV,λL |
三体がボラティリティに与える影響係数 |
| ξ |
ボラティリティのボラティリティ係数 |
| η |
モメンタムノイズ強度(新情報ショック) |
| ai,bi,ci |
収益/リスク/コスト係数 (i=M,V,L) |
| ρi |
混雑ペナルティ係数 |
| ε |
時間スケール分離パラメータ(≪1) |
五、制約検証
5.1 制約 C1:M の価格への正帰還
原文:モメンタム資本は価格変化に対して正帰還を生み出す、dSd(ポジション)>0(上昇時に買い、下落時に売る)。
方程式対応:方程式 (II) の αM⋅m⋅p 項。
検証:
- p>0(価格上昇)のとき、この項は正であり、dp/dt を増加させる、すなわち上昇を加速する
- p<0(価格下落)のとき、この項は負であり、dp/dt をより負にする、すなわち下落を加速する
- これはまさに「順張り」の数学的表現である:モメンタムの変化はモメンタム自体と同じ方向
結論:✓ 合格
5.2 制約 C2:V の価格への負帰還
原文:バリュー資本は価格変化に対して負帰還を生み出す、dSd(ポジション)<0(価格上昇時にポジションを減らし、価格下落時にポジションを増やす)。
方程式対応:方程式 (II) の −αV⋅v⋅x 項。
検証:
- x>0(価格が内在価値を上回る)のとき、この項は負であり、下方への力を生み出し、上昇を抑制する
- x<0(価格が内在価値を下回る)のとき、この項は正であり、上方への力を生み出し、下落を抑制する
- これはまさに平均回帰の力学機構である
結論:✓ 合格
5.3 制約 C3:L の価格への無方向帰還
原文:流動性資本は価格変化に対して方向性のある反応を示さない、dSd(ポジション)≈0。
方程式対応:方程式 (I) および (II) において、l は分母にのみ現れ、方向性のあるドリフト項を生み出さない。
検証:
- l が影響するのは価格ショックの大きさ(1/l および σ/lβ)であって、方向ではない
- マーケットメーカーは方向性を賭けず、流動性の緩衝を提供するのみ
結論:✓ 合格
5.4 制約 C4:M↑→σ↑
原文:モメンタム資本の増加はボラティリティの上昇をもたらす。
方程式対応:方程式 (III) の +λM⋅m⋅∣p∣ 項。
検証:
- m が増加すると、λM⋅m⋅∣p∣ が増加する
- これは dσ/dt のドリフト項を直接増加させる
- モメンタム資本が多く、モメンタムが強いほど、ボラティリティは高くなる
結論:✓ 合格
5.5 制約 C5:σ↑→L↓(崩壊スパイラル中)
原文:高ボラティリティ時にはマーケットメーカーが退避する。
方程式対応:方程式 (VI) の aLσ−bL∣p∣。
検証:
- 表面上は、aLσ は L の収益項であり、σ↑ は L に利益をもたらすはずである
- 重要な理解:原文の論理は、崩壊スパイラルにおいて、高 σ は高 ∣p∣(強いトレンド)を伴うというものである
- bL∣p∣>aLσ のとき、L の純収益は負となり、l の縮小をもたらす
- これはまさに「高ボラティリティ+強トレンド」環境下でのマーケットメーカー退避の機構である
パラメータ条件:崩壊スパイラルには、トレンドリスクがボラティリティ収益を上回るように bL/aL が十分に大きい必要がある。
結論:✓ 合格(適切なパラメータ条件下で)
5.6 制約 C6:L↑→ 価格ショック↓
原文:十分な流動性はモメンタム資本による価格へのショックを緩和する。
方程式対応:
- 方程式 (I) の σ/lβ:l が大きいほど、価格ノイズは小さい
- 方程式 (II) の 1/l:l が大きいほど、同じ力によるモメンタム変化は小さい
検証:
- 市場深度 l は価格ショック係数の分母である
- 深い市場では、大口注文でも価格を動かすことは難しい
- 流動性は「ショックアブソーバー」である
結論:✓ 合格
5.7 制約 C7:V は S<S∗ で買い、S>S∗ で売る
原文:バリュー資本は内在価値をアンカーとして逆張り操作を行う。
方程式対応:方程式 (II) の −αV⋅v⋅x 項。
検証:
- x=ln(S/S∗)
- S<S∗ のとき、x<0、この項 −αVvx>0 は上方への力(買い圧力)を生み出す
- S>S∗ のとき、x>0、この項 −αVvx<0 は下方への力(売り圧力)を生み出す
- 力の大きさは乖離度 ∣x∣ とバリュー資本規模 v に比例する
結論:✓ 合格
5.8 制約 C8:V↑→∣S−S∗∣↓
原文:バリュー資本の介入は価格を内在価値に回帰させる。
方程式対応:方程式 (II) の −αV⋅v⋅x 項。
検証:
- v が大きいほど、回帰力 ∣αVvx∣ は強い
- 強い回帰力は x をより速くゼロに近づける
- バリュー資本はシステムの「安定化装置」である
定量的分析:遅い変数を凍結した場合、方程式 (II) の x に関する部分は:
p˙≈−lαVvx−lγ−αMmp
x˙=p と組み合わせると、これは2次システムである。αVv>0 かつ γ>αMm のとき、システムは安定で、x は振動しながらゼロに収束する。
結論:✓ 合格
5.9 制約 C9:M の収益 ∝μ、リスク ∝δ、コスト ∝σ
原文:モメンタム資本はトレンドから利益を得、プレミアムはリスク、ボラティリティはコストである。
方程式対応:方程式 (IV) の収益率項 rM=aM∣p∣−bM∣x∣−cMσ。
検証:
- aM∣p∣:トレンド継続(∣p∣ 大)= M の利益 ✓
- −bM∣x∣:プレミアム大は反転の予兆であり、リスク ✓
- −cMσ:高ボラティリティ時にはストップロスが頻発し、コスト ✓
結論:✓ 合格
5.10 制約 C10:V の収益 ∝δ、リスク ∝σ、コスト ∝μ
原文:バリュー資本は価値乖離から利益を得、ボラティリティはリスク、トレンドはコストである。
方程式対応:方程式 (V) の収益率項 rV=aV∣x∣−bVσ−cV∣p∣。
検証:
- aV∣x∣:プレミアム大 = V の機会 ✓
- −bVσ:高ボラティリティ時には評価損が大きくなり、アンカーも不安定になりうる ✓
- −cV∣p∣:トレンド継続時には V はより長く待つ必要があり、資金効率が低下する ✓
結論:✓ 合格
5.11 制約 C11:L の収益 ∝σ、リスク ∝μ、コスト ∝δ
原文:流動性資本はボラティリティから利益を得、トレンドはリスク、プレミアムはコストである。
方程式対応:方程式 (VI) の収益率項 rL=aLσ−bL∣p∣−cL∣x∣。
検証:
- aLσ:ボラティリティ大 = 取引機会多、マーケットメイク収益高 ✓
- −bL∣p∣:トレンド強時には在庫が一方向に蓄積し、方向性損失に直面する ✓
- −cL∣x∣:プレミアム大時には自己防衛のためより大きなスプレッドが必要となり、効率が低下する ✓
結論:✓ 合格
5.12 制約 C12:資本過剰 → 収益率低下 → 規模縮小
原文:収益率駆動の自然選択により、三種類の資本は長期的に共存する。
方程式対応:方程式 (IV-VI) の混雑項 −ρi⋅i(i=m,v,l)。
検証:
- M を例にとる:dm/dt∝m⋅(rM−ρMm)
- m が過大になると、ρMm が支配的となり、実効収益率は負に転じ、m は縮小する
- これはロジスティック成長の混雑効果である
- 同じ機構が V と L にも適用され、いずれかの資本が無限に膨張するのを防ぐ
結論:✓ 合格
六、正帰還ループの完全な追跡
原文:M↑→σ↑→L↓→価格ショック↑→σ↑→M↑(または暴落)
SDE 追跡:
m↑:外部ショックまたは収益率による吸引
→σ↑:方程式 (III) の λMm∣p∣ が増加、ボラティリティ上昇
→l↓:方程式 (VI) において、bL∣p∣>aLσ のとき(強トレンドがボラティリティ収益を上回る)、l が縮小
→ 価格ショック↑:方程式 (II) の 1/l が増加、同じ力でより大きな p の変化が生じる
→σ↑:方程式 (III) の λMm∣p∣ がさらに増加
→m↑(または暴落):方程式 (IV) において:
- aM∣p∣ が支配的の場合:m は成長を継続
- bM∣x∣+cMσ が支配的の場合(プレミアム過大、ボラティリティ過高):m は縮小(暴落)
結論:✓ 正帰還ループは完全に実現されている
七、負帰還ループの完全な追跡
原文:∣S−S∗∣↑→σ↑→V↑→∣S−S∗∣↓→σ↓→L↑
SDE 追跡:
∣x∣↑:外部ショックにより価格が内在価値から乖離
→σ↑:価格乖離は通常ボラティリティを伴う(∣p∣ を介して伝播)
→v↑:方程式 (V) の aV∣x∣ が増加、V の収益が増加し、より多くのバリュー資本を吸引
→∣x∣↓:方程式 (II) の −αVvx が強化、回帰力が増大し、x はゼロに収束
→σ↓:∣x∣ と ∣p∣ が減少、方程式 (III) のボラティリティが低下
→l↑:方程式 (VI) において、低 ∣p∣ は L のリスクを低下させ、流動性の回帰を吸引
結論:✓ 負帰還ループは完全に実現されている
八、相状態分析の検証
冷状態(000):δ 低, μ 低, σ 低
対応:x≈0, p≈0, σ≈σˉlow
各資本収益率:
- rM=aM⋅0−bM⋅0−cMσˉ≈−cMσˉ<0
- rV=aV⋅0−bVσˉ−cV⋅0≈−bVσˉ<0
- rL=aLσˉ−bL⋅0−cL⋅0≈aLσˉ(微正または微負)
原文の記述:三者ともに利益を得られず、市場は縮小する。
検証:✓ 合格
熱状態(111):δ 高, μ 高, σ 高
対応:∣x∣ 大, ∣p∣ 大, σ 大
各資本収益率:三項とも大きく、符号はパラメータ比率に依存し、不確定。
原文の記述:三者ともに極端な環境に直面し、高収益高リスク、システムは臨界点にある。
検証:✓ 合格
九、三体類推の検証
原文:市場は三体が拮抗するためにカオス的振る舞いを示し、初期値に敏感である。
SDE 分析:
方程式 (II) の p のドリフト項:
lαMm−γ⋅p
- αMm>γ のとき、これは p の正のリアプノフ指数方向であり、システムは初期値に敏感
- 同時に −αVvx/l が非線形結合を提供
- 線形不安定 + 非線形結合 + ノイズ = カオス的振る舞いの古典的レシピ
結論:✓ 合格
十、統計的特性の検証
ボラティリティ・クラスタリング
原文:ボラティリティ・クラスタリングは市場の様式化された事実である。
SDE 実装:方程式 (III) の ξσdW3 の乗法的ノイズ。
機構:高 σ 時にはノイズが大きく、σ は高水準を維持しやすい。λMm∣p∣ の正帰還と相まって、クラスタリング効果は増幅される。
ファットテール分布
原文:収益率分布はファットテールを示す。
SDE 実装:価格ノイズ σ/lβ⋅dW1 において、σ と l はともに確率的である。
機構:確率的ボラティリティ自体がファットテールを生み出す;l が極端な状況で急落する(流動性退避)ことがテールをさらに厚くする。
十一、まとめ
| 制約 |
対応方程式項 |
検証結果 |
| C1 M 正帰還 |
αMmp |
✓ |
| C2 V 負帰還 |
−αVvx |
✓ |
| C3 L 無方向 |
l は分母のみ |
✓ |
| C4 M↑→σ↑ |
$\lambda_M m |
p |
| C5 σ↑→L↓ |
$a_L\sigma - b_L |
p |
| C6 L↑→ショック↓ |
1/l, σ/lβ |
✓ |
| C7 V アンカー機構 |
−αVvx |
✓ |
| C8 V↑→回帰 |
−αVvx |
✓ |
| C9 M 収益行列 |
$a_M |
p |
| C10 V 収益行列 |
$a_V |
x |
| C11 L 収益行列 |
$a_L\sigma - b_L |
p |
| C12 混雑効果 |
−ρi⋅i |
✓ |
12 の制約すべてに合格。
十二、今後の研究方向
- 数値シミュレーション:遅い変数を固定し、速い変数サブシステムをシミュレートし、アトラクター、リミットサイクル、カオス領域を探索する
- 分岐解析:m/v 比率を分岐パラメータとして、分岐図を描く
- パラメータ較正:実市場データを用いてパラメータを推定する
- 平均化法:時間スケール分離を利用し、遅い変数の有効力学を解析する