σμ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 AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • 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 AcademyRoadmapsRoteiros
Ctrl K
AO
+1,2 mil alunos
EntrarParticipar ↗
RoadmapsRoteiros
AO
+1,2 mil alunos
Entrar
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Roteiro de Linear Algebra›Módulo 2 · Resolver sistemas lineares›Tópico 2.3·~50 minIntermediário

Inversa e condicionamento

Posição no roteiroMódulo 2 de 4 · Tópico 3 de 3
2.3 Inversa e condicionamento
—Nota do lab— / 100XP total—Sequência—

Desfazer uma transformação

Uma matriz quadrada AAA tem inversa A−1A^{-1}A−1, com A−1A=AA−1=IA^{-1}A=AA^{-1}=IA−1A=AA−1=I, exatamente quando sua transformação não perde informação: entradas diferentes dão saídas diferentes. O determinante mede como AAA altera volumes (com sinal para inversões), então AAA é invertível exatamente quando det⁡A≠0\det A\neq0detA=0. Para uma matriz 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​). Regras úteis: det⁡(AB)=det⁡A det⁡B\det(AB)=\det A\,\det Bdet(AB)=detAdetB e 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)

Invertível não significa seguro. O número de condição κ(A)\kappa(A)κ(A) diz quanto um pequeno erro relativo em b\mathbf{b}b (ou nas entradas de AAA) pode ser ampliado na solução x\mathbf{x}x. Regra prática: perdem-se cerca de log⁡10κ\log_{10}\kappalog10​κ dos aproximadamente 16 dígitos da precisão dupla. Um determinante pequeno por si só não prova nada, porque depende da escala: 0,001 I0{,}001\,I0,001I tem determinante minúsculo, mas κ=1\kappa=1κ=1. O que torna κ\kappaκ enorme são colunas quase colineares.

Uma inversa e um quase erroPara 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, então 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​). Compare agora B=(1111,0001)B=\begin{pmatrix}1&1\\1&1{,}0001\end{pmatrix}B=(11​11,0001​): suas colunas são quase iguais, det⁡B=0,0001\det B=0{,}0001detB=0,0001 e κ(B)\kappa(B)κ(B) é cerca de 4×1044\times10^{4}4×104, de modo que um erro relativo de 10−610^{-6}10−6 nos dados pode virar um erro de cerca de 4 %4\,\%4% na solução.
  • Invertível: quadrada com det⁡A≠0\det A\neq0detA=0; nenhuma informação se perde.
  • Número de condição: quanto os erros nos dados são ampliados; um valor grande indica uma resolução frágil.
  • Nunca forme A−1A^{-1}A−1 para resolver: uma resolução por fatoração é mais rápida e precisa.
Ativos colineares quebram matrizes de covariânciaDois ETFs que acompanham quase o mesmo índice dão uma matriz de covariância com linhas quase iguais. O seu κ\kappaκ é enorme, então um otimizador de variância mínima baseado na inversa oscila muito com pequenas mudanças nos dados. Verifique κ\kappaκ (ou elimine ou combine os quase duplicados) antes de confiar nos pesos.

Lab: o que é avaliado

  • A inversa de uma matriz 2 por 220 pts
  • Um determinante 3 por 310 pts
  • Achar a matriz singular15 pts
  • Dígitos confiáveis15 pts
  • Boas práticas numéricas15 pts
  • Regras do determinante10 pts
  • Por que não inverter15 pts
0/ 7 respondidas
Tarefas
1Matriz · vetor20 ptsSem resposta

Preencha A−1A^{-1}A−1 para A=(4726)A=\begin{pmatrix}4&7\\2&6\end{pmatrix}A=(42​76​) (decimais são aceitos).

A−1=A^{-1} =A−1=
2Resposta numérica10 ptsSem resposta

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

det⁡=\det =det=aceito: ± 0.001
3Escolha15 ptsSem resposta

Qual destas matrizes é singular (sem inversa)?

4Resposta numérica15 ptsSem resposta

Um sistema tem κ(A)=106\kappa(A)=10^{6}κ(A)=106 e você o resolve em precisão dupla (cerca de 16 dígitos). Aproximadamente quantos dígitos corretos tem a solução?

[objectObject]=[object Object] =[objectObject]=aceito: ± 1
5Escolha15 ptsSem resposta

Selecione todas as boas práticas.

6Preencher as lacunas10 ptsSem resposta

Suponha det⁡A=4\det A=4detA=4 e det⁡B=3\det B=3detB=3. Preencha os espaços.

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)=
7Explicar com palavras15 ptsSem resposta

Em duas frases, diga por que o código deve chamar solve(A, b) em vez de inv(A) @ b.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 7 · Ctrl ⏎ para verificar

Entre para executar o seu código.

Todas as aulas são gratuitas para ler. Crie uma conta ou entre para executar código, rodar os testes e enviar este lab.

  • IDE no navegador com terminal bash — nada para instalar
  • Correção instantânea, XP e o seu lugar no roteiro
Criar uma contaEntrar
OU
Continuar com o GitHub