σμ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 Particle filtering and non-stationary time series›Módulo 2 · Filtros de partículas›Tópico 2.3·~55 minAvançado

Esquemas de reamostragem e degeneração

Posição no roteiroMódulo 2 de 3 · Tópico 3 de 5
2.3 Esquemas de reamostragem e degeneração
—Nota do lab— / 100XP total—Sequência—

Quando os pesos colapsam

Os pesos se multiplicam de um dia para o outro. Uma partícula que explica mal algumas observações recebe um peso minúsculo, e pesos minúsculos nunca se recuperam. Depois de dias suficientes uma partícula carrega quase todo o peso e as demais são trabalho perdido. É a degeneração dos pesos, que o tamanho amostral efetivo 1/∑w~21/\sum\tilde w^{2}1/∑w~2 mede. No SPY com 100010001000 partículas e sem reamostragem, ele vale 960960960 após o primeiro dia, 368368368 após dez, 525252 após cinquenta e cerca de 1,31{,}31,3 após 250250250 dias. Um filtro reduzido a uma partícula é um único palpite, não uma nuvem.

ESS=1∑iw~i2,ESS<N2,u+iN,  i=0,…,N−1\mathrm{ESS}=\frac{1}{\sum_i\tilde w_i^{2}},\qquad \mathrm{ESS}<\frac N2,\qquad \frac{u+i}{N},\ \ i=0,\dots,N-1ESS=∑i​w~i2​1​,ESS<2N​,Nu+i​,  i=0,…,N−1(1)

A reamostragem cura a degeneração: sorteia-se uma nuvem nova de NNN partículas na qual uma partícula de peso w~i\tilde w_iw~i​ aparece em média Nw~iN\tilde w_iNw~i​ vezes, e dá-se a cada cópia o mesmo peso 1/N1/N1/N. Há várias formas de sortear. O multinomial sorteia NNN uniformes independentes. O sistemático sorteia um único uniforme uuu e usa NNN ponteiros igualmente espaçados (u+i)/N(u+i)/N(u+i)/N sobre os pesos acumulados. O estratificado sorteia um uniforme em cada uma das NNN fatias iguais. O residual dá primeiro a cada partícula suas cópias inteiras, ⌊Nw~i⌋\lfloor N\tilde w_i\rfloor⌊Nw~i​⌋, e só sorteia o resto. O sistemático é o mais barato e o mais usado.

Reamostragem sistemática à mãoTomemos quatro partículas em −1,0,1,2-1,0,1,2−1,0,1,2 com pesos 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. Os pesos acumulados são 0,058, 0,316, 0,742, 10{,}058,\ 0{,}316,\ 0{,}742,\ 10,058, 0,316, 0,742, 1. Sorteemos u=0,3u=0{,}3u=0,3. Os ponteiros são (0,3+i)/4=0,075, 0,325, 0,575, 0,825(0{,}3+i)/4=0{,}075,\ 0{,}325,\ 0{,}575,\ 0{,}825(0,3+i)/4=0,075, 0,325, 0,575, 0,825. O ponteiro 0,0750{,}0750,075 cai na segunda casa, 0,3250{,}3250,325 e 0,5750{,}5750,575 na terceira, 0,8250{,}8250,825 na quarta. A nuvem nova é 0,1,1,20,1,1,20,1,1,2: a partícula em −1-1−1 sumiu e a pesada em 111 foi copiada duas vezes, perto das 4×0,426=1,74\times0{,}426=1{,}74×0,426=1,7 cópias esperadas. Caso perigoso: com pesos 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01 o tamanho amostral efetivo é 1,061{,}061,06 e qualquer conjunto de ponteiros escolhe só a primeira partícula, então a nuvem se reduz a um valor distinto.
  • Reamostrar quando o ESS é baixo: regra comum ESS<N/2\mathrm{ESS}<N/2ESS<N/2. Reamostrar todo dia desperdiça diversidade.
  • O sistemático usa um único número aleatório: ponteiros (u+i)/N(u+i)/N(u+i)/N sobre os pesos acumulados dão um sorteio barato e de pouco ruído.
  • O preço é o empobrecimento: muitas cópias de poucas partículas deixam a nuvem com pouca variedade, que o próximo passo de propagação precisa reconstruir.

Lab: o que é avaliado

  • Cópias esperadas15 pts
  • Um ponteiro sistemático10 pts
  • Pesos colapsados15 pts
  • Sobreviventes distintas10 pts
  • Devemos reamostrar hoje10 pts
  • Reamostragem sistemática10 pts
  • Associar os esquemas15 pts
  • Por que reamostrar15 pts
0/ 8 respondidas
Tarefas
1Resposta numérica15 ptsSem resposta

Com N=4N=4N=4 partículas, uma partícula com peso normalizado 0,4260{,}4260,426 aparece em média Nw~N\tilde wNw~ vezes após a reamostragem. Quanto é esse número? Dê duas casas decimais.

n=n =n=aceito: ± 0.02
2Resposta numérica10 ptsSem resposta

A reamostragem sistemática com N=4N=4N=4 e u=0,3u=0{,}3u=0,3 usa os ponteiros (u+i)/N(u+i)/N(u+i)/N. Quanto vale o segundo ponteiro, i=1i=1i=1? Dê três casas decimais.

p=p =p=aceito: ± 0.001
3Resposta numérica15 ptsSem resposta

Os pesos normalizados são 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01. Quanto vale o tamanho amostral efetivo 1/∑w~21/\sum\tilde w^{2}1/∑w~2? Dê duas casas decimais.

ess=ess =ess=aceito: ± 0.01
4Resposta numérica10 ptsSem resposta

Os pesos são 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01 e a reamostragem sistemática usa N=4N=4N=4 ponteiros com u=0,5u=0{,}5u=0,5. Quantas partículas originais diferentes aparecem na nova nuvem?

k=k =k=aceito: ± 0.01
5Escolha10 ptsSem resposta

Há N=1000N=1000N=1000 partículas e o tamanho amostral efetivo de hoje é 620620620. Com a regra «reamostrar se ESS < N/2», o que você faz?

6Ordenar os passos10 ptsSem resposta

Coloque em ordem os passos da reamostragem sistemática.

Arraste para reordenar, ou foque um item, Espaço para levantá-lo, setas para movê-lo e Espaço para soltá-lo.

  1. 1Para cada ponteiro, tomar a primeira partícula cujo peso acumulado o alcança
  2. 2Somar os pesos em pesos acumulados
  3. 3Colocar N ponteiros igualmente espaçados (u + i) / N
  4. 4Sortear um único número uniforme u
  5. 5Dar a cada cópia o peso 1 / N
7Associar15 ptsSem resposta

Associe cada esquema de reamostragem à sua forma de sortear.

Escolha o par de cada item da esquerda.

Multinomial
Sistemático
Estratificado
Residual
8Explicar com palavras15 ptsSem resposta

Em uma ou duas frases, explique por que um filtro de partículas precisa de reamostragem e quanto ela custa.

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