σμ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 Calculus›Módulo 4 · Otimização e sensibilidades›Tópico 4.1·~45 minIntermediário

O jacobiano e o hessiano

Posição no roteiroMódulo 4 de 4 · Tópico 1 de 3
4.1 O jacobiano e o hessiano
—Nota do lab— / 100XP total—Sequência—

Muitas saídas, muitas entradas

Uma função FFF pode mapear várias entradas em várias saídas, por exemplo fatores de risco nos valores de várias posições. Seu jacobiano JJJ é a matriz de todas as derivadas parciais de primeira ordem, com entrada Jij=∂Fi∂xjJ_{ij}=\frac{\partial F_i}{\partial x_j}Jij​=∂xj​∂Fi​​. É a melhor aproximação linear: para um pequeno movimento Δ\DeltaΔ, F(x+Δ)≈F(x)+JΔF(\mathbf{x}+\Delta)\approx F(\mathbf{x})+J\DeltaF(x+Δ)≈F(x)+JΔ.

J=(∂F1∂x∂F1∂y∂F2∂x∂F2∂y)J=\begin{pmatrix}\frac{\partial F_1}{\partial x} & \frac{\partial F_1}{\partial y}\\[2pt] \frac{\partial F_2}{\partial x} & \frac{\partial F_2}{\partial y}\end{pmatrix}J=(∂x∂F1​​∂x∂F2​​​∂y∂F1​​∂y∂F2​​​)(1)

Para uma única saída fff, as segundas derivadas formam o hessiano HHH, com entrada Hij=∂2f∂xi∂xjH_{ij}=\frac{\partial^2 f}{\partial x_i\partial x_j}Hij​=∂xi​∂xj​∂2f​. Para funções suaves ele é simétrico. Descreve a curvatura: f(x+Δ)≈f+∇f⋅Δ+12ΔTHΔf(\mathbf{x}+\Delta)\approx f+\nabla f\cdot\Delta+\tfrac{1}{2}\Delta^{T}H\Deltaf(x+Δ)≈f+∇f⋅Δ+21​ΔTHΔ. Onde ∇f=0\nabla f=0∇f=0, um HHH definido positivo indica um mínimo local, um definido negativo um máximo local e um misto um ponto de sela. O método de Newton para otimização usa HHH para dar passos melhores que a simples descida do gradiente.

Um jacobiano e um hessiano à mãoPara F(x,y)=(x2y, x+y2)F(x,y)=(x^2y,\ x+y^2)F(x,y)=(x2y, x+y2): J=(2xyx212y)J=\begin{pmatrix}2xy & x^2\\ 1 & 2y\end{pmatrix}J=(2xy1​x22y​), então em (1,2)(1,2)(1,2) é (4114)\begin{pmatrix}4&1\\1&4\end{pmatrix}(41​14​) com determinante 15: áreas pequenas são multiplicadas por cerca de 15. Para f(x,y)=x2y+y2f(x,y)=x^2y+y^2f(x,y)=x2y+y2: fxx=2yf_{xx}=2yfxx​=2y, fxy=2xf_{xy}=2xfxy​=2x, fyy=2f_{yy}=2fyy​=2, então em (1,2)(1,2)(1,2) o hessiano é (4222)\begin{pmatrix}4&2\\2&2\end{pmatrix}(42​22​).
  • Jacobiano: a matriz das parciais de primeira ordem de uma função vetorial; seu melhor modelo linear.
  • Hessiano: a matriz simétrica das parciais de segunda ordem de uma função escalar; sua curvatura.
  • Hessiano definido positivo em um ponto crítico: um mínimo local.
Não estime um hessiano por diferenças finitas se puder evitarDiferenças segundas dividem diferenças minúsculas por h2h^2h2, o que amplifica o ruído de arredondamento. Prefira fórmulas ou diferenciação automática e, se precisar usar diferenças, verifique se o resultado é simétrico.

Lab: o que é avaliado

  • O jacobiano em um ponto20 pts
  • O hessiano em um ponto20 pts
  • O determinante do jacobiano15 pts
  • Lendo um hessiano15 pts
  • A aproximação linear15 pts
  • O que cada matriz descreve15 pts
0/ 6 respondidas
Tarefas
1Matriz · vetor20 ptsSem resposta

Seja F(x,y)=(x2y, x+y2)F(x,y)=(x^2y,\ x+y^2)F(x,y)=(x2y, x+y2). Preencha o jacobiano J(1,2)J(1,2)J(1,2).

J(1,2)=J(1,2) =J(1,2)=
2Matriz · vetor20 ptsSem resposta

Para f(x,y)=x2y+y2f(x,y)=x^2y+y^2f(x,y)=x2y+y2, preencha o hessiano H(1,2)H(1,2)H(1,2).

H(1,2)=H(1,2) =H(1,2)=
3Resposta numérica15 ptsSem resposta

Calcule det⁡J(1,2)\det J(1,2)detJ(1,2) para o jacobiano (4114)\begin{pmatrix}4&1\\1&4\end{pmatrix}(41​14​). É o fator pelo qual FFF multiplica áreas pequenas perto de (1,2)(1,2)(1,2).

det⁡J=\det J =detJ=aceito: ± 0.001
4Escolha15 ptsSem resposta

Em um ponto onde ∇f=0\nabla f=0∇f=0, o hessiano é (2006)\begin{pmatrix}2&0\\0&6\end{pmatrix}(20​06​). Que tipo de ponto é?

5Resposta numérica15 ptsSem resposta

Aqui F(1,2)=(2,5)F(1,2)=(2,5)F(1,2)=(2,5) e J(1,2)=(4114)J(1,2)=\begin{pmatrix}4&1\\1&4\end{pmatrix}J(1,2)=(41​14​) para F(x,y)=(x2y, x+y2)F(x,y)=(x^2y,\ x+y^2)F(x,y)=(x2y, x+y2). Use F(x+Δ)≈F(x)+JΔF(\mathbf{x}+\Delta)\approx F(\mathbf{x})+J\DeltaF(x+Δ)≈F(x)+JΔ para estimar a primeira componente de F(1,1; 2,0)F(1{,}1;\,2{,}0)F(1,1;2,0).

F1=F_1 =F1​=aceito: ± 0.001
6Explicar com palavras15 ptsSem resposta

Em duas ou três frases, diga o que o jacobiano e o hessiano descrevem de uma função.

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