σμ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 を作らない:分解で解くほうが速く正確。
共線な資産は共分散行列を壊すほぼ同じ指数に連動する2つの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点未回答

コードで inv(A) @ b ではなく solve(A, b) を呼ぶべき理由を2文で述べてください。

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

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

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

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