σμ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 Optimization›Módulo 5 · Otimização na prática›Tópico 5.1·~50 minAvançado

Armadilhas do não convexo e reinícios

Posição no roteiroMódulo 5 de 5 · Tópico 1 de 2
5.1 Armadilhas do não convexo e reinícios
—Nota do lab— / 100XP total—Sequência—

Quando há mais de um vale

A maioria dos modelos interessantes não é convexa. Seus objetivos têm vários mínimos locais, pontos de sela e platôs, e a descida do gradiente simplesmente rola para o vale em que começa. Tome f(x)=x4−3x2+xf(x)=x^4-3x^2+xf(x)=x4−3x2+x. Sua derivada 4x3−6x+14x^3-6x+14x3−6x+1 se anula em x=−1.301x=-1.301x=−1.301, x=0.170x=0.170x=0.170 e x=1.131x=1.131x=1.131. O primeiro e o último são mínimos com f=−3.514f=-3.514f=−3.514 e f=−1.070f=-1.070f=−1.070, e o do meio é uma lombada entre eles. A descida com passo pequeno a partir de x0=2x_0=2x0​=2 termina em 1.1311.1311.131, apenas um mínimo local, e a partir de x0=−2x_0=-2x0​=−2 termina em −1.301-1.301−1.301, o global. Dois inícios tão próximos quanto x0=0.1x_0=0.1x0​=0.1 e x0=0.2x_0=0.2x0​=0.2, dos dois lados da lombada em 0.1700.1700.170, terminam em vales diferentes.

f(x)=x4−3x2+x,f′(x)=4x3−6x+1,f(−1.301)=−3.514,f(1.131)=−1.070f(x)=x^{4}-3x^{2}+x,\qquad f'(x)=4x^{3}-6x+1,\qquad f(-1.301)=-3.514,\quad f(1.131)=-1.070f(x)=x4−3x2+x,f′(x)=4x3−6x+1,f(−1.301)=−3.514,f(1.131)=−1.070(1)

Pontos de sela são a outra armadilha. Neles o gradiente é zero mas a hessiana tem autovalores de ambos os sinais, então a função desce em algumas direções e sobe em outras. Em altas dimensões a maioria dos pontos estacionários são selas, e a descida do gradiente simples se arrasta perto delas. O momento, o ruído dos minilotes e a informação de segunda ordem ajudam a escapar. Platôs, em que o gradiente é minúsculo numa região ampla, parecem convergência mas não são.

O remédio padrão é reiniciar. Execute o otimizador a partir de KKK pontos de partida diferentes e fique com o resultado de menor valor do objetivo. Se cada início cai independentemente no vale global com probabilidade ppp, então pelo menos um dos KKK inícios o faz com probabilidade 1−(1−p)K1-(1-p)^K1−(1−p)K. Para p=0.2p=0.2p=0.2 e K=10K=10K=10 isso dá 0.8930.8930.893, e chegar a 99%99\%99% exige K≥ln⁡0.01/ln⁡0.8=20.6K\ge\ln0.01/\ln0.8=20.6K≥ln0.01/ln0.8=20.6, ou seja 21 inícios. Melhores pontos de partida também ajudam: uma solução de um relaxamento convexo, a calibração de ontem ou uma busca em grade nos poucos parâmetros que mais importam.

  • Mínimos locais: a descida termina no vale em que começa, então o início decide a resposta.
  • Selas e platôs: gradiente zero ou minúsculo sem mínimo; momento e ruído ajudam.
  • Reinícios: KKK inícios têm sucesso com probabilidade 1−(1−p)K1-(1-p)^K1−(1−p)K; guarde o melhor valor do objetivo.
Um bom ajuste não é um ajuste estávelNuma calibração, vários conjuntos de parâmetros bem diferentes podem atingir quase o mesmo valor do objetivo. Depois dos reinícios, compare não só o melhor objetivo mas também os parâmetros entre inícios e entre dias. Se eles pulam enquanto o ajuste continua igualmente bom, o problema é mal identificado e um otimizador mais bonito não resolve: restrinja, regularize ou reduza os parâmetros.

Lab: o que é avaliado

  • O valor do mínimo global15 pts
  • Começar em x = 215 pts
  • Dois inícios próximos15 pts
  • Sucesso de reinícios aleatórios15 pts
  • Quantos reinícios15 pts
  • Perto de um ponto de sela10 pts
  • Comparar os resultados dos reinícios15 pts
0/ 7 respondidas
Tarefas
1Resposta numérica15 ptsSem resposta

Para f(x)=x4−3x2+xf(x)=x^4-3x^2+xf(x)=x4−3x2+x o mínimo global está em x=−1.3008x=-1.3008x=−1.3008. Quanto vale fff ali? Dê com duas casas decimais.

fglob=fglob =fglob=aceito: ± 0.01
2Escolha15 ptsSem resposta

A descida do gradiente com passo pequeno em f(x)=x4−3x2+xf(x)=x^4-3x^2+xf(x)=x4−3x2+x começa em x0=2x_0=2x0​=2. Onde termina?

3Escolha15 ptsSem resposta

A lombada de fff está em x=0.170x=0.170x=0.170. Onde terminam as descidas com passo pequeno a partir de x0=0.1x_0=0.1x0​=0.1 e de x0=0.2x_0=0.2x0​=0.2?

4Resposta numérica15 ptsSem resposta

Cada início aleatório alcança o vale global com probabilidade p=0.2p=0.2p=0.2, de forma independente. Com K=10K=10K=10 inícios, qual é a probabilidade de pelo menos um ter sucesso? Dê com três casas decimais.

psucc=psucc =psucc=aceito: ± 0.002
5Resposta numérica15 ptsSem resposta

Com p=0.2p=0.2p=0.2, qual é o menor número de inícios KKK que faz a probabilidade de sucesso ser pelo menos 99%99\%99%?

k=k =k=aceito: ± 0.5
6Escolha10 ptsSem resposta

Qual afirmação sobre um ponto de sela está correta?

7Explicar com palavras15 ptsSem resposta

Em uma ou duas frases, diga como escolhe entre os resultados de vários reinícios de uma calibração e o que mais verifica.

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