DAMPTRYK AF METHANOL
Øvelse i fysisk kemi, forår 2004
Vejledningen er opdateret 1.marts 2004. Den opdaterede version
kan findes her
En vejledning i at gennemføre databehandlingen ved brug af programmet
gnuplot findes her
Eksempel på brug af programmet Mathematica til at fitte modelparametre
til de eksperimentelle data.
Eksempel på et Mathematica script der læser
data fra filen med navnet dat med celciustemperaturer i første kolonne
og tryk (mbar) i anden kolonne.
xy = ReadList["dat", {Number,Number} ];
Needs["Statistics`NonlinearFit`"];
YofX[x_] = a Exp[b/(x+273.15)+c/(x+273.15)^2];
initialguesses = {{a, 10^8},{b, -5000.0},{c,0}};
ft = NonlinearFit[xy,YofX[x],x,initialguesses];
gt=NonlinearRegress[xy,YofX[x],x,initialguesses,
RegressionReport->ParameterTable];
fp = Plot[ft,{x,0,100}];
gp = ListPlot[xy];
Show[fp,gp]
Dette script forudsætter at datafilen er indeholder to kolonner
af data hvor den første er temperaturen i grader C og den anden
trykket i mbar. Datafilen data fra programmet clausius har tiden
i første kolonne, trykket i anden kolonne og temperaturen (grader
C) i tredie kolonne. En ny datafil dat som indeholder målepunkter
fra tidintervallet 2000-3000 og som placerer kolonne 3 i kolonne 1 og
kolonne 2 i kolonne 2 kan genereres med programmet awk.på
følgende måde.
awk '{if ($1>2000 && $1<3000) print $3, $2}' data
> dat
Derefter køres Mathematica scriptet med kommandoen \verb!
<<script !. Fittet kan gå galt hvis initialgættene ligger
for langt fra det endelige fit. Prøv med bedre startgæt.
Som output fås de fittede parametre med standard afvigelser
og signifikansniveauer for hypotese at den pågældende parameter
er 0 samt tegninger der viser modellem med de fittede parametre og de eksperimentelle
data.