validation

Sensibilité et budget d'erreur

Comme P(t) est un produit, chaque facteur a une élasticité unitaire — une variation d'un pour cent de n'importe quel facteur change la sortie d'un pour cent. Ce qui distingue les facteurs, c'est leur incertitude : une décomposition de variance et un budget d'erreur montrant quel terme domine.

Dernière mise à jour:

Une analyse de sensibilité pose une question tranchante : si une entrée bouge, de combien la sortie bouge-t-elle, et quelle entrée importe le plus ? Pour P(t)P(t) la réponse a une structure exceptionnellement nette, parce que l’équation maîtresse est multiplicative. Ce chapitre sépare les deux choses qu’une analyse de sensibilité confond — la réponse à une variation et la contribution à l’incertitude — et montre que pour P(t)P(t) elles se comportent très différemment. Il s’appuie directement sur l’analyse d’incertitude et partage son outillage avec le carnet de reproductibilité.

Chaque facteur a une élasticité unitaire

Écrivons la sortie agrégée schématiquement comme un produit de ses facteurs, P=∏ixiP = \prod_i x_i avec xi∈{Φ,σ,Δp,vph,ηconv}x_i \in \{\Phi, \sigma, \Delta p, v_{\text{ph}}, \eta_{\text{conv}}\}. L’élasticité de PP par rapport à n’importe quel facteur — la variation fractionnaire de la sortie par variation fractionnaire de l’entrée — vaut alors exactement un :

∂ln⁡P∂ln⁡xi=1for every i.(1)\frac{\partial \ln P}{\partial \ln x_i} = 1 \qquad \text{for every } i. \tag{1}

Une variation d’un pour cent de n’importe quel facteur seul change la sortie d’un pour cent, ni plus ni moins. La réponse au premier ordre est donc uniforme : aucun facteur n’est intrinsèquement plus « à effet de levier » qu’un autre. C’est une propriété authentique de la forme multiplicative, et cela signifie que le classement habituel par dérivée locale n’est pas informatif ici — il classerait les cinq facteurs à égalité.

Ce qui diffère réellement, c’est l’incertitude

Si les facteurs ne sont pas également réactifs, ils ne sont certainement pas également connus. La question de sensibilité utile pour P(t)P(t) est donc un budget d’erreur : quelle part de l’incertitude sur la sortie chaque facteur apporte-t-il ? En traitant les facteurs comme indépendants et distribués log-normalement — le choix naturel pour des grandeurs incertaines sur un ordre de grandeur — les log-variances s’additionnent :

Var⁡(ln⁡P)=∑iVar⁡(ln⁡xi)=∑iσi2,(2)\operatorname{Var}(\ln P) = \sum_i \operatorname{Var}(\ln x_i) = \sum_i \sigma_i^{2}, \tag{2}

de sorte que la contribution fractionnaire du facteur ii à la variance de la sortie est simplement

ci=σi2∑jσj2.(3)c_i = \frac{\sigma_i^{2}}{\sum_j \sigma_j^{2}}. \tag{3}

Le budget d’erreur est fixé entièrement par les incertitudes relatives σi\sigma_i, non par les valeurs des facteurs. C’est la grandeur qu’il vaut la peine de classer.

Budget d’erreur (illustratif)

Le tableau utilise des incertitudes log-échelle illustratives — signalées comme telles, choisies pour refléter que la section efficace du neutrino est le terme le moins contraint et la vitesse des phonons parmi les mieux contraints. Les valeurs réelles proviendraient de la mesure (mesure et falsifiabilité).

Facteurσi\sigma_i illustratif (unités ln)Part de variance cic_i
σ(E)\sigma(E) — section efficace de diffusion0,640,0 %
Φ\Phi — densité de flux de rayonnement0,527,8 %
Δp\Delta p — transfert de quantité de mouvement0,417,8 %
ηconv\eta_{\text{conv}} — rendement de conversion0,310,0 %
vphv_{\text{ph}} — vitesse des phonons0,24,4 %
Diagramme à barres horizontales de la contribution de chaque facteur à la variance de ln P : section efficace 40 pour cent, densité de flux 27,8 pour cent, transfert de quantité de mouvement 17,8 pour cent, rendement 10 pour cent, vitesse des phonons 4,4 pour cent.
Contribution de chaque facteur à la variance de ln⁡P\ln P, triée. Le terme de section efficace domine le budget d’erreur ; la vitesse des phonons est quasi négligeable.

Ce que le budget vous apprend

Deux conclusions pratiques en découlent, et ce sont elles qui rendent l’analyse digne d’être menée. Premièrement, la priorité de mesure : parce que le terme de section efficace porte la plus grande part de la variance, un programme expérimental réduit l’incertitude de sortie le plus vite en contraignant d’abord σ(E)\sigma(E) ; diviser par deux son σi\sigma_i retire plus de la variance totale que diviser par deux n’importe quel autre. Deuxièmement, où la précision est gaspillée : resserrer la vitesse des phonons, déjà le terme le mieux contraint, ne bouge presque pas le total. Combiné au résultat d’élasticité unitaire, le tableau est complet — chaque facteur importe également pour la valeur de P(t)P(t), mais ils importent très inégalement pour son incertitude, et le budget d’erreur dit exactement où regarder. Comme pour tous les nombres de ce site, les σi\sigma_i ici sont illustratifs ; la méthode — décomposer la log-variance, classer les parts, agir sur la plus grande — est le contenu transférable, et elle est reproductible depuis le carnet.

↑