σμ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 4 · Mínimos quadrados e fatorações›Tópico 4.3·~50 minIntermediário

Cholesky e matrizes definidas positivas

Posição no roteiroMódulo 4 de 4 · Tópico 3 de 4
4.3 Cholesky e matrizes definidas positivas
—Nota do lab— / 100XP total—Sequência—

Raízes quadradas para matrizes

Uma matriz simétrica AAA é definida positiva quando xTAx>0\mathbf{x}^{T}A\mathbf{x}>0xTAx>0 para todo x≠0\mathbf{x}\neq\mathbf{0}x=0; de forma equivalente, todos os seus autovalores são positivos. Matrizes de covariância são sempre semidefinidas positivas, porque wTΣw\mathbf{w}^{T}\Sigma\mathbf{w}wTΣw é a variância de uma carteira e variância não pode ser negativa. São definidas positivas, com toda carteira de variância estritamente positiva, exatamente quando nenhum ativo é redundante.

A=LLT,Ljj=Ajj−∑k<jLjk2,Lij=Aij−∑k<jLikLjkLjj (i>j)A=LL^{T},\qquad L_{jj}=\sqrt{A_{jj}-\sum_{k<j}L_{jk}^{2}},\qquad L_{ij}=\frac{A_{ij}-\sum_{k<j}L_{ik}L_{jk}}{L_{jj}}\ (i>j)A=LLT,Ljj​=Ajj​−k<j∑​Ljk2​​,Lij​=Ljj​Aij​−∑k<j​Lik​Ljk​​ (i>j)(1)

Toda matriz definida positiva admite uma fatoração de Cholesky A=LLTA=LL^{T}A=LLT com LLL triangular inferior de diagonal positiva: o análogo matricial de uma raiz quadrada. Ela aproveita a simetria, custa cerca de 13n3\tfrac13n^331​n3, metade da LU, e dispensa pivotamento. Se o algoritmo encontra um número negativo sob uma raiz, a matriz não é definida positiva, o que faz do Cholesky um teste barato dessa propriedade.

Um fator e uma simulaçãoPara A=(4225)A=\begin{pmatrix}4&2\\2&5\end{pmatrix}A=(42​25​): L11=4=2L_{11}=\sqrt4=2L11​=4​=2, L21=2/2=1L_{21}=2/2=1L21​=2/2=1, L22=5−12=2L_{22}=\sqrt{5-1^2}=2L22​=5−12​=2, então L=(2012)L=\begin{pmatrix}2&0\\1&2\end{pmatrix}L=(21​02​) e LLT=(4225)LL^{T}=\begin{pmatrix}4&2\\2&5\end{pmatrix}LLT=(42​25​). Para simular dois retornos com correlação 0,60{,}60,6 e variâncias unitárias, fatore Σ=(10,60,61)\Sigma=\begin{pmatrix}1&0{,}6\\0{,}6&1\end{pmatrix}Σ=(10,6​0,61​) em L=(100,60,8)L=\begin{pmatrix}1&0\\0{,}6&0{,}8\end{pmatrix}L=(10,6​00,8​) e tome x=Lz\mathbf{x}=L\mathbf{z}x=Lz com z\mathbf{z}z normais padrão independentes: então Cov⁡(x)=LLT=Σ\operatorname{Cov}(\mathbf{x})=LL^{T}=\SigmaCov(x)=LLT=Σ.
  • Definida positiva: xTAx>0\mathbf{x}^{T}A\mathbf{x}>0xTAx>0 para todo x≠0\mathbf{x}\neq\mathbf{0}x=0; todos os autovalores positivos.
  • Cholesky: A=LLTA=LL^{T}A=LLT, metade do custo da LU, sem pivotamento e um teste barato de definição positiva.
  • Sorteios correlacionados: x=Lz\mathbf{x}=L\mathbf{z}x=Lz transforma normais independentes z\mathbf{z}z em retornos com covariância Σ\SigmaΣ.
Uma covariância amostral pode não fatorarSe você estima uma covariância com menos observações que ativos, ou dois ativos têm retornos idênticos, a matriz amostral é apenas semidefinida e o Cholesky para num pivô nulo ou ligeiramente negativo. Os remédios: mais dados, remover duplicados, encolhimento para um alvo bem condicionado ou somar à diagonal um pequeno múltiplo da identidade.

Lab: o que é avaliado

  • Fatorar uma matriz20 pts
  • Fatorar uma matriz de correlação20 pts
  • Um sorteio correlacionado15 pts
  • Qual matriz é definida positiva15 pts
  • Quando o Cholesky falha10 pts
  • Simular retornos correlacionados20 pts
0/ 6 respondidas
Tarefas
1Matriz · vetor20 ptsSem resposta

Preencha o fator de Cholesky LLL de A=(9335)A=\begin{pmatrix}9&3\\3&5\end{pmatrix}A=(93​35​) (triangular inferior, diagonal positiva).

L=L =L=
2Matriz · vetor20 ptsSem resposta

Preencha o fator de Cholesky da matriz de correlação (10,50,51)\begin{pmatrix}1&0{,}5\\0{,}5&1\end{pmatrix}(10,5​0,51​) (decimais aceitos, três dígitos).

L=L =L=
3Resposta numérica15 ptsSem resposta

Com Σ=(10,60,61)\Sigma=\begin{pmatrix}1&0{,}6\\0{,}6&1\end{pmatrix}Σ=(10,6​0,61​) e L=(100,60,8)L=\begin{pmatrix}1&0\\0{,}6&0{,}8\end{pmatrix}L=(10,6​00,8​), o sorteio independente z=(1,2)\mathbf{z}=(1,2)z=(1,2) dá x=Lz\mathbf{x}=L\mathbf{z}x=Lz. Quanto vale x2x_2x2​?

x2=x_2 =x2​=aceito: ± 0.001
4Escolha15 ptsSem resposta

Qual destas matrizes simétricas é definida positiva?

5Escolha10 ptsSem resposta

O Cholesky falha numa matriz de covariância amostral. Selecione todas as causas prováveis.

6Explicar com palavras20 ptsSem resposta

Em duas frases, diga como o fator de Cholesky de uma matriz de covariância é usado para simular retornos correlacionados.

Modo matemáticaEntrada por teclado ou LaTeXTarefa 1 de 6 · 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