资本市场三体动力学 SDE 方程组推导
2026-02-07
基于《资本市场的三体动力学假说》一文,本文推导出完整的随机微分方程(SDE)系统,并逐一验证其是否满足原文的所有定性约束。
一、状态变量定义
系统有两个时间尺度:
快变量(秒到天):
- 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(Position)>0 |
动量资本定义 |
| C2 |
V 对价格负反馈:dSd(Position)<0 |
价值资本定义 |
| C3 |
L 对价格无方向反馈:dSd(Position)≈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 |
三体对波动率的影响系数 |
| ξ |
vol-of-vol 系数 |
| η |
动量噪声强度(新信息冲击) |
| ai,bi,ci |
收益/风险/成本系数 (i=M,V,L) |
| ρi |
拥挤惩罚系数 |
| ε |
时间尺度分离参数(≪1) |
五、约束检验
5.1 约束 C1:M 对价格正反馈
原文:动量资本对价格变化产生正反馈,dSd(Position)>0(涨时买入,跌时卖出)。
方程对应:方程 (II) 中的 αM⋅m⋅p 项。
验证:
- 当 p>0(价格上涨)时,该项为正,使 dp/dt 增大,即加速上涨
- 当 p<0(价格下跌)时,该项为负,使 dp/dt 更负,即加速下跌
- 这正是"追涨杀跌"的数学表达:动量的变化与动量本身同向
结论:✓ 通过
5.2 约束 C2:V 对价格负反馈
原文:价值资本对价格变化产生负反馈,dSd(Position)<0(价格上涨时减仓,价格下跌时加仓)。
方程对应:方程 (II) 中的 −αV⋅v⋅x 项。
验证:
- 当 x>0(价格高于内在价值)时,该项为负,产生向下的力,抑制上涨
- 当 x<0(价格低于内在价值)时,该项为正,产生向上的力,抑制下跌
- 这正是均值回归的力学机制
结论:✓ 通过
5.3 约束 C3:L 对价格无方向反馈
原文:流动性资本对价格变化无方向性反应,dSd(Position)≈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,这是一个二阶系统。当 α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 的正 Lyapunov 指数方向,系统对初值敏感
- 同时 −αVvx/l 提供非线性耦合
- 线性不稳定 + 非线性耦合 + 噪声 = 混沌行为的经典配方
结论:✓ 通过
十、统计特性验证
波动率聚集
原文:波动率聚集是市场的 stylized fact。
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 比值为分岔参数,绘制分岔图
- 参数校准:用真实市场数据估计参数
- 平均化方法:利用时间尺度分离,分析慢变量的有效动力学