Esercitazioni
del
corso
di Laboratorio di Calcolo (I)
n. 10
fit lineare
questa esercitazione verra' valutata
- create la cartella ex10 e spostatevici
- scaricate il programma
che effettua il fit lineare e calcola il chi2 dal sito del libro di
testo
- copiatelo (con cp) in un file con un altro
nome (ex10.c) compilatelo, eseguitelo, osservatene il funzionamento
- commentate dettagliatamente il programma
(inserendo i vostri commenti tra /* e */ nel file stesso) illustrandone
il funzionamento
- riprendete il programma originale (senza i
commenti) e modificatelo in modo tale che legga i dati da un file
esterno
- eseguite il programma sui file di dati ex10.dat: la prima colonna rappresenta la pressione atmosferica misurata in mm di mercurio ad una certa quota sulle alpi, la seconda colonna la temperatura di ebollizione dell'acqua in gradi C alla stessa quota, la terza colonna una prima stima dell'errore sulla temperatura.
Inserite una stampa nella funzione chi2 per vedere come le diverse misure contribuiscono al Chi2 finale. Editate una copia del file di dati togliendo la misura peggiore (che scarta maggiormente dalla retta, ovvero che fornisce il maggiore contributo al Chi2). Girate nuovamente per verificare che i parametri della retta siano sostanzialmente immutati e che il Chi2 per grado di liberta' sia migliorato. Se il modello lineare e' corretto buone misure con errori stimati correttamente dovrebbero fornire un Chi2 per grado di liberta' prossimo a 1. Cosa potete dire degli errori assunti: sono stimati correttamente? sottostimati? sovrastimati? Che errore dovreste usare per ottenere 1 per il chi2 per grado di liberta'? Scrivete le risposte in un commento in fondo al vostro programma.