σμ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
  • 日本語
  • 中文
Linear Algebra 学习路线›模块 2 · 求解线性方程组›主题 2.3·~50分钟中级

逆矩阵与条件数

路线进度模块 2/4 · 主题 3/3
2.3 逆矩阵与条件数
—实验得分— / 100总 XP—连续天数—

撤销一个映射

方阵 AAA 有逆矩阵 A−1A^{-1}A−1(A−1A=AA−1=IA^{-1}A=AA^{-1}=IA−1A=AA−1=I),当且仅当它的映射不丢失信息,即不同输入给出不同输出。行列式衡量 AAA 对体积的缩放(翻转带符号),所以 AAA 可逆当且仅当 det⁡A≠0\det A\neq0detA=0。对 2×22\times22×2 矩阵,A−1=1ad−bc(d−b−ca)A^{-1}=\frac{1}{ad-bc}\begin{pmatrix}d&-b\\-c&a\end{pmatrix}A−1=ad−bc1​(d−c​−ba​)。有用的规则:det⁡(AB)=det⁡A det⁡B\det(AB)=\det A\,\det Bdet(AB)=detAdetB,det⁡(A−1)=1/det⁡A\det(A^{-1})=1/\det Adet(A−1)=1/detA。

κ(A)=∥A∥ ∥A−1∥=σmax⁡σmin⁡,∥δx∥∥x∥≲κ(A) ∥δb∥∥b∥\kappa(A)=\|A\|\,\|A^{-1}\|=\frac{\sigma_{\max}}{\sigma_{\min}},\qquad \frac{\|\delta\mathbf{x}\|}{\|\mathbf{x}\|}\lesssim\kappa(A)\,\frac{\|\delta\mathbf{b}\|}{\|\mathbf{b}\|}κ(A)=∥A∥∥A−1∥=σmin​σmax​​,∥x∥∥δx∥​≲κ(A)∥b∥∥δb∥​(1)

可逆并不等于安全。条件数 κ(A)\kappa(A)κ(A) 表示 b\mathbf{b}b(或 AAA 的元素)中的微小相对误差在解 x\mathbf{x}x 中最多会被放大多少。经验法则:在双精度约 16 位有效数字中,大约要损失 log⁡10κ\log_{10}\kappalog10​κ 位。仅凭行列式很小说明不了什么,因为它与尺度有关:0.001 I0.001\,I0.001I 的行列式极小,但 κ=1\kappa=1κ=1。真正让 κ\kappaκ 变得巨大的是几乎共线的列。

一个逆矩阵和一个近似奇异的例子对 A=(4726)A=\begin{pmatrix}4&7\\2&6\end{pmatrix}A=(42​76​):det⁡A=24−14=10\det A=24-14=10detA=24−14=10,所以 A−1=110(6−7−24)A^{-1}=\frac{1}{10}\begin{pmatrix}6&-7\\-2&4\end{pmatrix}A−1=101​(6−2​−74​)。再看 B=(1111.0001)B=\begin{pmatrix}1&1\\1&1.0001\end{pmatrix}B=(11​11.0001​):它的两列几乎相等,det⁡B=0.0001\det B=0.0001detB=0.0001,κ(B)\kappa(B)κ(B) 约为 4×1044\times10^{4}4×104,因此数据中 10−610^{-6}10−6 的相对误差可能变成解中约 4 %4\,\%4% 的误差。
  • 可逆:det⁡A≠0\det A\neq0detA=0 的方阵;不丢失信息。
  • 条件数:数据误差被放大的程度;很大就意味着求解不稳。
  • 不要为了求解而计算 A−1A^{-1}A−1:用分解求解更快也更准。
共线资产会破坏协方差矩阵两只几乎跟踪同一指数的 ETF 会给出各行几乎相等的协方差矩阵。它的 κ\kappaκ 极大,因此基于其逆矩阵的最小方差优化器会随数据的微小变化而剧烈摆动。在相信权重之前,先检查 κ\kappaκ(或剔除、合并近乎重复的资产)。

实验:评分内容

  • 2×2 矩阵的逆20 分
  • 3×3 行列式10 分
  • 找出奇异矩阵15 分
  • 可信的位数15 分
  • 良好的数值实践15 分
  • 行列式规则10 分
  • 为什么不求逆15 分
0/ 已答 7 题
任务
1矩阵 · 向量20 分未作答

对 A=(4726)A=\begin{pmatrix}4&7\\2&6\end{pmatrix}A=(42​76​) 填写 A−1A^{-1}A−1(可用小数)。

A−1=A^{-1} =A−1=
2数值作答10 分未作答

计算 det⁡(201132014)\det\begin{pmatrix}2&0&1\\1&3&2\\0&1&4\end{pmatrix}det​210​031​124​​。

det⁡=\det =det=允许误差:± 0.001
3选择15 分未作答

下列哪个矩阵是奇异的(没有逆)?

4数值作答15 分未作答

一个方程组的 κ(A)=106\kappa(A)=10^{6}κ(A)=106,用双精度(约 16 位)求解。解大约有几位正确数字?

[objectObject]=[object Object] =[objectObject]=允许误差:± 1
5选择15 分未作答

选出所有良好做法。

6填空10 分未作答

设 det⁡A=4\det A=4detA=4,det⁡B=3\det B=3detB=3。填空。

det⁡(AB)=,det⁡(A−1)=\displaystyle \det(AB)=\htmlData{blank=p}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}},\qquad \det(A^{-1})=\htmlData{blank=q}{\phantom{\rule[-0.35em]{2.4em}{1.35em}}}det(AB)=,det(A−1)=
7用文字解释15 分未作答

请用两句话说明为什么代码应调用 solve(A, b) 而不是 inv(A) @ b。

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

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

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

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