next up previous
Next: ...una parte molto opzionale Up: Simulazione dell'esperienza dello spettrometro Previous: Misura degli indici di


Determinazione di \( A\) e \( B\)

Lo scopo finale dell'esperimento è ricavare, a partire dai valori osservati di \( n\left( \lambda \right) \) i valori delle costanti \( A\) e \( B\) della parametrizzazione (1).

Essendo la parametrizzazione lineare, sarà possibile fare un fit ad una linea retta di \( n\left( \lambda \right) ^{2} \) come funzione di \( 1/\lambda ^{2} \).

A tale proposito, ricordiamo un po' di formule senza dimostrazione.

Nel caso di un fit di un insieme di punti sperimentali \( \left( x_{i},y_{i}\right) \), con errore \( \sigma _{i} \) su \( y_{i} \), ad una retta \( y=a+bx \), la migliore stima dei parametri \( a \) e \( b \) è quella data dalla minimizzazione del \( \chi ^{2} \):

\begin{displaymath}
\chi ^{2}\left( a,b\right) =\sum _{i=1}^{N}\left( \frac{y_{i}-a-bx_{i}}{\sigma _{i}}\right) ^{2}.\end{displaymath}

Definendo le quantità

\begin{displaymath}
\begin{array}{rcl}
W & = & \sum ^{N}_{i=1}\frac{1}{\sigma _{...
...}},\\
\Delta & = & W\cdot XX-\left( X\right) ^{2},
\end{array}\end{displaymath}

la migliore stima di \( a \) e \( b \) è data da
\begin{displaymath}
\begin{array}{rcl}
a & = & \frac{XX\cdot Y-X\cdot XY}{\Delta },\\
b & = & \frac{W\cdot XY-X\cdot Y}{\Delta },
\end{array}\end{displaymath} (5)

e le relative varianze sono
\begin{displaymath}
\begin{array}{rcl}
\sigma _{a}^{2} & = & \frac{XX}{\Delta },\\
\sigma _{b}^{2} & = & \frac{W}{\Delta }.
\end{array}\end{displaymath} (6)

Poiché gli errori sui valori di \( n\left( \lambda \right) \) sono noti dopo lo studio della sezione 4, l'utilizzarli nel fit è un semplice esercizio di propagazione degli errori.

Utilizzare il file scritto alla fine della sezione 4 per produrre 10000 stime dei parametri \( A\) e \( B\) dell'equazione (1) e delle relative incertezze \( \sigma _{A} \) e \( \sigma _{B} \).

Confrontare i valori calcolati di \( \sigma _{A} \) e \( \sigma _{B} \) con gli errori ``veri'' ottenuti dalle distribuzioni dei residui \( A-A_{\textrm{vero}} \) e \( B-B_{\textrm{vero}} \). Qual è la ragione della differenza?

Provare ad applicare lo stesso procedimento al caso in cui \( \alpha \) e \( \theta _{0} \) sono noti con precisione assoluta.


next up previous
Next: ...una parte molto opzionale Up: Simulazione dell'esperienza dello spettrometro Previous: Misura degli indici di
Attilio Andreazza 2000-11-15