σμ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 AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • 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 AcademyFeuilles de routeParcours
Ctrl K
AO
+1,2 k apprenants
ConnexionRejoindre ↗
Feuilles de routeParcours
AO
+1,2 k apprenants
Connexion
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Parcours Particle filtering and non-stationary time series›Module 2 · Filtres particulaires›Sujet 2.3·~55 minAvancé

Schémas de rééchantillonnage et dégénérescence

Position dans le parcoursModule 2 sur 3 · Sujet 3 sur 5
2.3 Schémas de rééchantillonnage et dégénérescence
—Score du lab— / 100XP total—Série—

Quand les poids s'effondrent

Les poids se multiplient d'un jour à l'autre. Une particule qui explique mal quelques observations reçoit un poids minuscule, et les poids minuscules ne remontent jamais. Après assez de jours, une particule porte presque tout le poids et les autres sont du travail perdu. C'est la dégénérescence des poids, que mesure la taille d'échantillon effective 1/∑w~21/\sum\tilde w^{2}1/∑w~2. Sur SPY avec 100010001000 particules et sans rééchantillonnage, elle vaut 960960960 après le premier jour, 368368368 après dix, 525252 après cinquante et environ 1,31{,}31,3 après 250250250 jours. Un filtre réduit à une particule est une seule hypothèse, pas un nuage.

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)

Le rééchantillonnage guérit la dégénérescence : on tire un nouveau nuage de NNN particules où une particule de poids w~i\tilde w_iw~i​ doit apparaître Nw~iN\tilde w_iNw~i​ fois en moyenne, puis on donne à chaque copie le même poids 1/N1/N1/N. Il y a plusieurs façons de tirer. Le multinomial tire NNN uniformes indépendants. Le systématique tire un seul uniforme uuu et utilise NNN pointeurs régulièrement espacés (u+i)/N(u+i)/N(u+i)/N sur les poids cumulés. Le stratifié tire un uniforme dans chacune des NNN tranches égales. Le résiduel donne d'abord à chaque particule ses copies entières, ⌊Nw~i⌋\lfloor N\tilde w_i\rfloor⌊Nw~i​⌋, et ne tire que le reste. Le systématique est le moins cher et le plus utilisé.

Le rééchantillonnage systématique à la mainPrenons quatre particules en −1,0,1,2-1,0,1,2−1,0,1,2 de poids 0,058, 0,258, 0,426, 0,2580{,}058,\ 0{,}258,\ 0{,}426,\ 0{,}2580,058, 0,258, 0,426, 0,258. Les poids cumulés sont 0,058, 0,316, 0,742, 10{,}058,\ 0{,}316,\ 0{,}742,\ 10,058, 0,316, 0,742, 1. Tirons u=0,3u=0{,}3u=0,3. Les pointeurs sont (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. Le pointeur 0,0750{,}0750,075 tombe dans la deuxième case, 0,3250{,}3250,325 et 0,5750{,}5750,575 dans la troisième, 0,8250{,}8250,825 dans la quatrième. Le nouveau nuage est 0,1,1,20,1,1,20,1,1,2 : la particule en −1-1−1 a disparu et la lourde en 111 est copiée deux fois, proche de ses 4×0,426=1,74\times0{,}426=1{,}74×0,426=1,7 copies attendues. Cas dangereux : avec des poids 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01, la taille d'échantillon effective est 1,061{,}061,06, et tout jeu de pointeurs ne choisit que la première particule : le nuage se réduit à une seule valeur.
  • Rééchantillonner quand l'ESS est basse : règle courante ESS<N/2\mathrm{ESS}<N/2ESS<N/2. Rééchantillonner chaque jour gaspille de la diversité.
  • Le systématique utilise un seul nombre aléatoire : les pointeurs (u+i)/N(u+i)/N(u+i)/N sur les poids cumulés donnent un tirage économique et peu bruité.
  • Le prix est l'appauvrissement : beaucoup de copies de quelques particules laissent le nuage peu varié, ce que l'étape de propagation suivante doit reconstruire.

Lab : ce qui est noté

  • Copies attendues15 pts
  • Un pointeur systématique10 pts
  • Poids effondrés15 pts
  • Survivantes distinctes10 pts
  • Faut-il rééchantillonner aujourd'hui10 pts
  • Rééchantillonnage systématique10 pts
  • Associer les schémas15 pts
  • Pourquoi rééchantillonner15 pts
0/ 8 traités
Exercices
1Réponse numérique15 ptsSans réponse

Avec N=4N=4N=4 particules, une particule de poids normalisé 0,4260{,}4260,426 doit apparaître Nw~N\tilde wNw~ fois après rééchantillonnage. Que vaut ce nombre ? Donnez deux décimales.

n=n =n=accepté : ± 0.02
2Réponse numérique10 ptsSans réponse

Le rééchantillonnage systématique avec N=4N=4N=4 et u=0,3u=0{,}3u=0,3 utilise les pointeurs (u+i)/N(u+i)/N(u+i)/N. Que vaut le deuxième pointeur, i=1i=1i=1 ? Donnez trois décimales.

p=p =p=accepté : ± 0.001
3Réponse numérique15 ptsSans réponse

Les poids normalisés sont 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01. Que vaut la taille d'échantillon effective 1/∑w~21/\sum\tilde w^{2}1/∑w~2 ? Donnez deux décimales.

ess=ess =ess=accepté : ± 0.01
4Réponse numérique10 ptsSans réponse

Les poids sont 0,97, 0,01, 0,01, 0,010{,}97,\ 0{,}01,\ 0{,}01,\ 0{,}010,97, 0,01, 0,01, 0,01 et le rééchantillonnage systématique utilise N=4N=4N=4 pointeurs avec u=0,5u=0{,}5u=0,5. Combien de particules d'origine différentes apparaissent dans le nouveau nuage ?

k=k =k=accepté : ± 0.01
5Choix10 ptsSans réponse

Il y a N=1000N=1000N=1000 particules et la taille d'échantillon effective du jour est 620620620. Avec la règle « rééchantillonner si ESS < N/2 », que faites-vous ?

6Remettre les étapes dans l'ordre10 ptsSans réponse

Mettez dans l'ordre les étapes du rééchantillonnage systématique.

Faites glisser pour réordonner, ou placez le focus sur un élément, Espace pour le soulever, les flèches pour le déplacer et Espace pour le déposer.

  1. 1Pour chaque pointeur, prendre la première particule dont le poids cumulé l'atteint
  2. 2Additionner les poids en poids cumulés
  3. 3Placer N pointeurs régulièrement espacés (u + i) / N
  4. 4Tirer un seul nombre uniforme u
  5. 5Donner à chaque copie le poids 1 / N
7Associer15 ptsSans réponse

Associez chaque schéma de rééchantillonnage à sa façon de tirer.

Choisissez le partenaire de chaque élément de gauche.

Multinomial
Systématique
Stratifié
Résiduel
8Expliquer avec des mots15 ptsSans réponse

En une ou deux phrases, expliquez pourquoi un filtre particulaire a besoin du rééchantillonnage et ce que cela coûte.

Mode mathsSaisie au clavier ou en LaTeXExercice 1 sur 8 · Ctrl ⏎ pour vérifier

Connectez-vous pour exécuter votre code.

Tous les cours sont gratuits à lire. Créez un compte ou connectez-vous pour exécuter du code, lancer les tests et soumettre ce lab.

  • IDE dans le navigateur avec terminal bash — rien à installer
  • Correction instantanée, XP et votre place dans le parcours
Créer un compteSe connecter
OU
Continuer avec GitHub