σμ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 の2つだけなので、フィルタはその2つを持ち運べばよいのです。これがカルマンフィルタで、この設定では厳密な答えです。毎日2つの段階を踏みます。予測では平均はそのままで分散が 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 が小さい)なときは 111 に近づき、データを信頼します。逆の場合は 000 に近づき、予測を信頼します。差 y−m−y-m^{-}y−m− は驚きで、イノベーションと呼ばれます。更新のたびに 1−K<11-K<11−K<1 なので分散は減り、予測のたびに qqq 増えます。十分に長く回すと、両者は一定の水準でつり合います。

2日間のフィルタリング信念 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)と見積もります。

  • 1日2つの数:すべてが線形で釣鐘形なら、平均と分散で足りる。
  • ゲインは信頼のつまみ:KKK が 111 に近いと新しいデータを、000 に近いと予測を信頼する。
  • 限界:ボラティリティ自体が隠れている、あるいはノイズが厚い裾をもつと、信念は釣鐘形でなくなり粒子が必要になる。

ラボ:採点される内容

  • 最初のゲイン15点
  • 更新後の平均15点
  • 更新後の分散10点
  • 2日目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) はいくつですか。小数第3位まで答えてください。

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−) はいくつですか。小数第3位まで答えてください。

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− はいくつですか。小数第3位まで答えてください。

P=P =P=許容:± 0.003
4数値で答える15点未回答

2日目の予測平均は 1.4291.4291.429、ゲインは K=0.548K=0.548K=0.548、新しい観測は y2=1y_2=1y2​=1 です。更新後の平均 m2m_2m2​ はいくつですか。小数第3位まで答えてください。

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点未回答

1〜2文で、更新後の分散 P=(1−K)P−P=(1-K)P^{-}P=(1−K)P− が予測分散より大きくなることがない理由を説明してください。

数学モードキーボードまたは LaTeX で入力問題1/8 · Ctrl ⏎ で確認

コードを実行するにはログインしてください。

レクチャーはすべて無料で読めます。コードの実行、テスト、ラボの提出にはアカウントの作成かログインが必要です。

  • ブラウザ内IDEとbashターミナル。インストール不要
  • 即時採点、XP、スキルでの現在地
アカウントを作成ログイン
または
GitHubで続ける