324
Radioprotection et ingénierie nucléaire
Désignons par f*(X) une densité de probabilité arbitraire.
On peut écrire :
Cela montre qu’il est possible d’estimer la grandeur Y en construisant une chaîne d’événements avec une densité de probabilité f*(X) différente de la densité de probabilité naturelle
f(X) à condition de prendre comme nouvel estimateur de la grandeur Y non pas X mais
.
On dit aussi que l’on corrige l’estimateur X par le rapport de la densité naturelle à la
densité arbitraire. La densité f*(X) est souvent appelée ici densité biaisée. Ce qualificatif
« biaisée » est trompeur car cette technique de réduction de la variance dans son principe
n’introduit pas de biais.
Le problème revient donc à trouver la densité f*(X) conduisant à la meilleure qualité
possible de la simulation. On montre en théorie que si l’on prend pour f*(X) une forme
construite avec la solution du problème adjoint (résolution de l’équation du transport
Figure 10.26. Atténuation du flux neutronique d’environ 14 décades entre le cœur du réacteur et
l’extérieur de l’enceinte de confinement ou bâtiment réacteur (Baur, 1985).
E Y
( )
Xf X
( ) X
d
D
#
X
f X
( )
f * X
( )
------------- -f * X
( ) X
d
D
#
ϭ
ϭ
X
f X
( )
f * X
( )
------------- -