σμ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.1·~1小时中级

奇异值分解

路线进度模块 2/5 · 主题 1/2
2.1 奇异值分解
—实验得分— / 100总 XP—连续天数—

每个矩阵都是旋转、拉伸再旋转

特征分解需要方阵,即使是方阵也可能失败。奇异值分解对任何矩阵(无论是否为长方形)都适用:A=UΣVTA=U\Sigma V^{T}A=UΣVT,其中 UUU 和 VVV 为正交矩阵,Σ\SigmaΣ 为对角矩阵,元素非负 σ1≥σ2≥⋯≥0\sigma_1\ge\sigma_2\ge\cdots\ge0σ1​≥σ2​≥⋯≥0,即奇异值。几何上,AAA 先用 VTV^{T}VT 旋转,再把第 iii 个坐标轴拉伸 σi\sigma_iσi​ 倍,最后用 UUU 旋转。单位圆变成半轴长恰为奇异值的椭圆。

A=UΣVT,σi=λi(ATA),∥A∥F2=∑iσi2,κ(A)=σ1σmin⁡A=U\Sigma V^{T},\qquad \sigma_i=\sqrt{\lambda_i(A^{T}A)},\qquad \lVert A\rVert_F^{2}=\sum_i\sigma_i^{2},\qquad \kappa(A)=\frac{\sigma_1}{\sigma_{\min}}A=UΣVT,σi​=λi​(ATA)​,∥A∥F2​=i∑​σi2​,κ(A)=σmin​σ1​​(1)

VVV 的列是 ATAA^{T}AATA 的特征向量,UUU 的列是 AATAA^{T}AAT 的特征向量,对应的特征值为 σi2\sigma_i^2σi2​。所以 SVD 就是对 ATAA^{T}AATA 应用谱定理,但计算时从不真正构造这个乘积,因为那会把条件数平方并损失一半的有效数字。很容易读出各种事实:秩是非零奇异值的个数,∥A∥2=σ1\lVert A\rVert_2=\sigma_1∥A∥2​=σ1​ 是最大拉伸,Frobenius 范数的平方等于奇异值平方和,方阵有 ∣det⁡A∣=∏σi|\det A|=\prod\sigma_i∣detA∣=∏σi​,条件数 σ1/σmin⁡\sigma_1/\sigma_{\min}σ1​/σmin​ 衡量 AAA 能放大多少误差。对 m≥nm\ge nm≥n 的 m×nm\times nm×n 数据矩阵,完整的 UUU 是 m×mm\times mm×m,Σ\SigmaΣ 是 m×nm\times nm×n,VVV 是 n×nn\times nn×n。

2×2 的例子取 A=(3045)A=\begin{pmatrix}3&0\\4&5\end{pmatrix}A=(34​05​)。ATA=(25202025)A^{T}A=\begin{pmatrix}25&20\\20&25\end{pmatrix}ATA=(2520​2025​) 的特征值为 45 和 5,所以奇异值为 σ1=45≈6.708\sigma_1=\sqrt{45}\approx6.708σ1​=45​≈6.708,σ2=5≈2.236\sigma_2=\sqrt5\approx2.236σ2​=5​≈2.236。校验:σ12+σ22=50=9+16+25=∥A∥F2\sigma_1^2+\sigma_2^2=50=9+16+25=\lVert A\rVert_F^2σ12​+σ22​=50=9+16+25=∥A∥F2​,σ1σ2=15=∣det⁡A∣\sigma_1\sigma_2=15=|\det A|σ1​σ2​=15=∣detA∣。条件数为 σ1/σ2=3\sigma_1/\sigma_2=3σ1​/σ2​=3。若 σ2\sigma_2σ2​ 是 0.0010.0010.001,条件数约为 6700,解 Ax=bA\mathbf x=\mathbf bAx=b 会损失近四位有效数字。
  • SVD:每个矩阵都能写成 UΣVTU\Sigma V^{T}UΣVT,UUU、VVV 正交,奇异值非负且降序。
  • 事实:秩等于非零 σi\sigma_iσi​ 的个数,∥A∥F2=∑σi2\lVert A\rVert_F^2=\sum\sigma_i^2∥A∥F2​=∑σi2​,∣det⁡A∣=∏σi|\det A|=\prod\sigma_i∣detA∣=∏σi​。
  • 条件数:σ1/σmin⁡\sigma_1/\sigma_{\min}σ1​/σmin​;数值大意味着求解不稳定、逆矩阵不可靠。
从收益率到成分只需一次调用若 XXX 是去均值后的 T×NT\times NT×N 收益矩阵,其协方差为 XTX/(T−1)X^{T}X/(T-1)XTX/(T−1),所以 XXX 的右奇异向量就是主方向,σi2/(T−1)\sigma_i^2/(T-1)σi2​/(T−1) 是沿这些方向的方差。直接对 XXX 做 SVD 比对协方差矩阵做特征分解更精确,并且在 N>TN>TN>T 时也能工作。

实验:评分内容

  • 最大奇异值20 分
  • Frobenius 恒等式10 分
  • 奇异值之积10 分
  • 完整 SVD 的尺寸15 分
  • 奇异值与特征值15 分
  • 条件数15 分
  • 为什么不对 AᵀA 做特征分解15 分
0/ 已答 7 题
任务
1数值作答20 分未作答

已知 ATA=(25202025)A^{T}A=\begin{pmatrix}25&20\\20&25\end{pmatrix}ATA=(2520​2025​),求 A=(3045)A=\begin{pmatrix}3&0\\4&5\end{pmatrix}A=(34​05​) 的最大奇异值,保留三位小数。

sigma1=sigma_1 =sigma1​=允许误差:± 0.01
2数值作答10 分未作答

对同一个 AAA,所有元素的平方和等于 σ12+σ22\sigma_1^2+\sigma_2^2σ12​+σ22​。∥A∥F2\lVert A\rVert_F^{2}∥A∥F2​ 是多少?

Frobeniussquared=Frobenius squared =Frobeniussquared=允许误差:± 0.01
3数值作答10 分未作答

对同一个 AAA,乘积 σ1σ2\sigma_1\sigma_2σ1​σ2​ 等于 ∣det⁡A∣|\det A|∣detA∣。它是多少?

product=product =product=允许误差:± 0.01
4填空15 分未作答

数据矩阵 AAA 有 5 行 3 列。补全其完整 SVD A=UΣVTA=U\Sigma V^{T}A=UΣVT 各因子的尺寸。

U: ×,Σ: ×,V: ×\displaystyle U:\ \htmlData{blank=a}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}\times\htmlData{blank=a}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}},\qquad \Sigma:\ \htmlData{blank=a}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}\times\htmlData{blank=b}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}},\qquad V:\ \htmlData{blank=b}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}\times\htmlData{blank=b}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}U: ×,Σ: ×,V: ×
5选择15 分未作答

AAA 的奇异值与特征值有什么关系?

6数值作答15 分未作答

某矩阵的奇异值为 6、3、0.5,它的条件数是多少?

kappa=kappa =kappa=允许误差:± 0.01
7用文字解释15 分未作答

用一两句话说明为什么计算 AAA 的 SVD 在数值上比对 ATAA^{T}AATA 做特征分解更好。

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

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

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

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