σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft Academy路线图学习路线
Ctrl K
AO
+1.2k 位学习者
登录加入 ↗
路线图学习路线
AO
+1.2k 位学习者
登录
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft Academy路线图学习路线
Ctrl K
AO
+1.2k 位学习者
登录加入 ↗
路线图学习路线
AO
+1.2k 位学习者
登录
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Particle filtering and non-stationary time series 学习路线›模块 1 · 会变动的隐含状态›主题 1.2·~50分钟高级

卡尔曼滤波回顾

路线进度模块 1/3 · 主题 2/2
1.2 卡尔曼滤波回顾
—实验得分— / 100总 XP—连续天数—

钟形曲线的滤波器

当状态做线性变动、噪声都是钟形时,会出现一个好结果:关于状态的信念永远是一条钟形曲线。钟形曲线只需要两个数,均值 mmm 和方差 PPP,所以滤波器只要带着这两个数就行。这就是卡尔曼滤波,在这种设定下它给出的是精确答案。它每天走两步。预测:均值不变,方差增加 qqq。更新:观测 yyy 把均值朝自己拉近一个比例 KKK(称为增益),方差随之减小。

P−=P+q,K=P−P−+r,m=m−+K (y−m−),P=(1−K) P−P^{-}=P+q,\quad K=\frac{P^{-}}{P^{-}+r},\quad m=m^{-}+K\,(y-m^{-}),\quad P=(1-K)\,P^{-}P−=P+q,K=P−+rP−​,m=m−+K(y−m−),P=(1−K)P−(1)

更新其实就是加权平均:新均值为 (1−K) m−+K y(1-K)\,m^{-}+K\,y(1−K)m−+Ky。增益 KKK 在 000 与 111 之间。当预测很不确定(P−P^{-}P− 大)或观测很精确(rrr 小)时,KKK 接近 111,你更相信数据;反之接近 000,你更相信预测。差值 y−m−y-m^{-}y−m− 是“意外”,称为新息。每次更新都会减小方差(因为 1−K<11-K<11−K<1),每次预测又使其增加 qqq。运行足够久之后,两者在一个稳定水平上达到平衡。

两天的滤波取信念 N(0,2)N(0,2)N(0,2)、q=0.5q=0.5q=0.5、r=1r=1r=1,第一个观测 y1=2y_1=2y1​=2。预测:P−=2.5P^{-}=2.5P−=2.5。增益:K=2.5/3.5=0.714K=2.5/3.5=0.714K=2.5/3.5=0.714。更新:m1=0+0.714×2=1.429m_1=0+0.714\times2=1.429m1​=0+0.714×2=1.429,P1=(1−0.714)×2.5=0.714P_1=(1-0.714)\times2.5=0.714P1​=(1−0.714)×2.5=0.714。第二天 y2=1y_2=1y2​=1:预测 P−=0.714+0.5=1.214P^{-}=0.714+0.5=1.214P−=0.714+0.5=1.214,K=1.214/2.214=0.548K=1.214/2.214=0.548K=1.214/2.214=0.548。意外为 1−1.429=−0.4291-1.429=-0.4291−1.429=−0.429,所以 m2=1.429−0.548×0.429=1.194m_2=1.429-0.548\times0.429=1.194m2​=1.429−0.548×0.429=1.194。长期来看,预测方差稳定在满足 P−=q+P−r/(P−+r)P^{-}=q+P^{-}r/(P^{-}+r)P−=q+P−r/(P−+r) 的位置,对这组数恰好是 111,稳态增益为 0.50.50.5。

同一个滤波器还能跟踪随时间变化的贝塔。状态是像随机游走一样漂移的贝塔 βt\beta_tβt​,观测是股票收益 yt=βt mt+ety_t=\beta_t\,m_t+e_tyt​=βt​mt​+et​,其中 mtm_tmt​ 是市场收益。唯一的改变是增益变成 K=P−mt/(P−mt2+r)K=P^{-}m_t/(P^{-}m_t^{2}+r)K=P−mt​/(P−mt2​+r):市场大幅波动的日子对贝塔很有信息,mt≈0m_t\approx0mt​≈0 的日子几乎没有信息。2005 年以来苹果对 SPY 的全样本贝塔是 1.071.071.07,而 q=10−3q=10^{-3}q=10−3 的卡尔曼滤波给出:2008 年底 0.770.770.77,2018 年底 1.421.421.42,最后一天 0.400.400.40,标准误差约 0.250.250.25。

  • 每天两个数:一切都是线性且钟形时,均值和方差就够了。
  • 增益是信任旋钮:KKK 接近 111 信任新数据,接近 000 信任预测。
  • 局限:如果波动率本身是隐藏的,或噪声有厚尾,信念就不再是钟形,需要用粒子。

实验:评分内容

  • 第一个增益15 分
  • 更新后的均值15 分
  • 更新后的方差10 分
  • 第二天15 分
  • 稳态水平10 分
  • 噪声更大的观测10 分
  • 配对术语10 分
  • 为什么更新会减小方差15 分
0/ 已答 8 题
任务
1数值作答15 分未作答

预测方差 P−=2.5P^{-}=2.5P−=2.5,测量噪声 r=1r=1r=1,增益 K=P−/(P−+r)K=P^{-}/(P^{-}+r)K=P−/(P−+r) 是多少?保留三位小数。

K=K =K=允许误差:± 0.002
2数值作答15 分未作答

预测均值为 000,观测 y1=2y_1=2y1​=2,增益 K=0.714K=0.714K=0.714。更新后的均值 m1=m−+K(y1−m−)m_1=m^{-}+K(y_1-m^{-})m1​=m−+K(y1​−m−) 是多少?保留三位小数。

m=m =m=允许误差:± 0.003
3数值作答10 分未作答

取 P−=2.5P^{-}=2.5P−=2.5、K=0.714K=0.714K=0.714,更新后的方差 P1=(1−K)P−P_1=(1-K)P^{-}P1​=(1−K)P− 是多少?保留三位小数。

P=P =P=允许误差:± 0.003
4数值作答15 分未作答

第二天预测均值为 1.4291.4291.429,增益 K=0.548K=0.548K=0.548,新观测 y2=1y_2=1y2​=1。更新后的均值 m2m_2m2​ 是多少?保留三位小数。

m=m =m=允许误差:± 0.003
5数值作答10 分未作答

长期来看,预测方差 aaa 满足 a=q+a r/(a+r)a=q+a\,r/(a+r)a=q+ar/(a+r)。对 q=0.5q=0.5q=0.5、r=1r=1r=1 检验哪个值是解:a=1a=1a=1 行吗?请输入 a=1a=1a=1 时 q+a r/(a+r)q+a\,r/(a+r)q+ar/(a+r) 的值。

v=v =v=允许误差:± 0.001
6选择10 分未作答

测量噪声 rrr 变得大得多,其他不变。增益和更新后的均值会怎样?

7配对10 分未作答

把每个术语与其含义配对。

为左侧每一项选择配对项。

增益 K
新息
预测方差
更新后的方差
8用文字解释15 分未作答

用一两句话解释:为什么更新后的方差 P=(1−K)P−P=(1-K)P^{-}P=(1−K)P− 绝不会比预测方差大。

数学模式键盘或 LaTeX 输入第 1/8 题 · Ctrl ⏎ 检查

登录后即可运行你的代码。

所有讲解均可免费阅读。创建账户或登录后,即可运行代码、运行测试并提交此实验。

  • 浏览器内 IDE 与 bash 终端,无需安装
  • 即时评分、XP 以及你在学习路线中的位置
创建账户登录
或
使用 GitHub 继续