σμ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
  • 日本語
  • 中文
Matrix Decompositions スキル›モジュール 2 · SVDと低ランク›トピック 2.2·~50 分中級

低ランク近似

スキル上の位置モジュール 2/5 · トピック 2/2
2.2 低ランク近似
—ラボスコア— / 100合計XP—連続日数—

行列の重要な部分だけを残す

SVDは行列を、大きい層から小さい層の順に並んだランク1の層の和 A=∑iσiuiviTA=\sum_i\sigma_i\mathbf u_i\mathbf v_i^{T}A=∑i​σi​ui​viT​ として表します。最初の kkk 層だけを残すと、ランク kkk の行列 Ak=UkΣkVkTA_k=U_k\Sigma_kV_k^{T}Ak​=Uk​Σk​VkT​ が得られます。エッカート=ヤングの定理は、誤差をフロベニウスノルムで測ってもスペクトルノルムで測っても、この切り捨てが最良のランク kkk 近似であると述べます。巧妙に選んだ kkk 個の方向でも、上位 kkk 個の特異方向には勝てません。

Ak=∑i=1kσiuiviT,∥A−Ak∥F=∑i>kσi2,∥A−Ak∥2=σk+1A_k=\sum_{i=1}^{k}\sigma_i\mathbf u_i\mathbf v_i^{T},\qquad \lVert A-A_k\rVert_F=\sqrt{\sum_{i>k}\sigma_i^{2}},\qquad \lVert A-A_k\rVert_2=\sigma_{k+1}Ak​=i=1∑k​σi​ui​viT​,∥A−Ak​∥F​=i>k∑​σi2​​,∥A−Ak​∥2​=σk+1​(1)

誤差は捨てた特異値そのものなので、残る2乗サイズの割合は ∑i≤kσi2/∑iσi2\sum_{i\le k}\sigma_i^2\big/\sum_i\sigma_i^2∑i≤k​σi2​/∑i​σi2​ で、AAA が中心化データなら寄与率(説明される分散)と呼ばれます。m×nm\times nm×n 行列のランク kkk 近似は mnmnmn 個ではなく k(m+n+1)k(m+n+1)k(m+n+1) 個の数で済むので、低ランクは圧縮でもあります。金融では σi\sigma_iσi​ の小さい層は主に標本ノイズなので、切り捨てはノイズ除去になり、kkk は残す割合のしきい値や特異値のギャップで決めます。

特異値が10, 6, 3, 1行列の特異値が10, 6, 3, 1だとします。2乗は100, 36, 9, 1で合計146です。最良のランク2近似は146のうち 100+36=136100+36=136100+36=136 を残し、割合は 93.2%93.2\%93.2%、フロベニウス誤差は 9+1=10≈3.16\sqrt{9+1}=\sqrt{10}\approx3.169+1​=10​≈3.16、スペクトル誤差は σ3=3\sigma_3=3σ3​=3 です。最良のランク1近似のスペクトル誤差は σ2=6\sigma_2=6σ2​=6。1000×2501000\times2501000×250 の行列ではランク5は250,000個の代わりに 5×(1000+250+1)=62555\times(1000+250+1)=62555×(1000+250+1)=6255 個の数で済みます。
  • エッカート=ヤング:SVDを kkk 項で打ち切ると、どちらのノルムでも最良のランク kkk 近似になる。
  • 誤差:フロベニウス ∑i>kσi2\sqrt{\sum_{i>k}\sigma_i^2}∑i>k​σi2​​、スペクトル σk+1\sigma_{k+1}σk+1​。
  • 記憶量:mnmnmn 個ではなく k(m+n+1)k(m+n+1)k(m+n+1) 個の数。
kの選択は判断仮定のない規則はありません。90%90\%90% のような固定割合は簡単ですが恣意的で、特異値に明確なギャップがあればそれが自然な基準になり、このスキルの後半の共分散の場面ではマルチェンコ・パスツールの端がシグナルとノイズを分ける根拠のあるしきい値になります。どれを使うにしても標本外で選択を確かめてください。訓練データを最もよく説明するランクが、新しいデータを最もよく予測するとは限りません。

ラボ:採点される内容

  • フロベニウス誤差15点
  • 残る割合15点
  • ランク1近似15点
  • 低ランク行列の記憶量10点
  • エッカート=ヤングの定理15点
  • 行列の圧縮15点
  • なぜ切り捨てがノイズを除くのか15点
0/ 7問回答済み
問題
1数値で答える15点未回答

特異値が10, 6, 3, 1の行列について、最良のランク2近似のフロベニウス誤差を小数第2位まで求めてください。

error=error =error=許容:± 0.02
2数値で答える15点未回答

同じ行列で、ランク2近似はフロベニウスノルムの2乗の何パーセントを残しますか。小数第1位まで答えてください。

kept=kept =kept=許容:± 0.1
3行列・ベクトル15点未回答

対称行列 A=(2112)A=\begin{pmatrix}2&1\\1&2\end{pmatrix}A=(21​12​) の特異値は3と1で、主方向は u=v=(1,1)/2\mathbf u=\mathbf v=(1,1)/\sqrt2u=v=(1,1)/2​ です。最良のランク1近似 A1=3 uvTA_1=3\,\mathbf u\mathbf v^{T}A1​=3uvT を計算してください。

A1=A_1 =A1​=
4数値で答える10点未回答

1000×2501000\times2501000×250 の行列を UkU_kUk​, Σk\Sigma_kΣk​, VkV_kVk​ でランク5近似します。k(m+n+1)k(m+n+1)k(m+n+1) を使って、保存すべき数の個数を求めてください。

numbers=numbers =numbers=許容:± 0.5
5選択15点未回答

エッカート=ヤングの定理は何を述べていますか。

6手順を並べ替える15点未回答

低ランク近似の手順を正しい順に並べてください。

ドラッグして並べ替えます。または項目にフォーカスし、スペースで持ち上げ、矢印キーで移動、スペースで置きます。

  1. 1誤差または残る分散の割合を測る
  2. 2特異値から kkk を選ぶ
  3. 3行列のSVDを計算する
  4. 4上位 kkk 項から Ak=UkΣkVkTA_k=U_k\Sigma_kV_k^{T}Ak​=Uk​Σk​VkT​ を再構成する
7言葉で説明する15点未回答

リターン行列の小さな特異値を捨てるとノイズを除ける理由を1〜2文で述べてください。

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

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

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

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