blob: dd9c4d409796a9140d742e125b880620a17b737a [file] [log] [blame]
avm99963545e7802020-11-11 17:00:13 +01001#!/usr/bin/env gnuplot -c
2# == DEFINICIONS ==
3outputfile = "../output/finalgraph" # Nom de la imatge resultant (sense extensió)
4datafile = "../data/finaldata.dat" # Nom del fitxer de dades que es vol usar
5error_sist_b = 0.11 # Error sistemàtic en
6
7# == CONFIGURACIÓ DE L'OUTPUT PEL LATEX ==
8set terminal cairolatex size 10cm, 7cm
9set output outputfile.'.tex'
10
11# == CONFIGURACIÓ DEL PLOT ==
12set xlabel "$T \\, (\\si{\\celsius})$"
13set ylabel "$S \\, (\\si{\\per\\ampere})$"
14
15# Opcions per la llegenda:
16set key above
17set key spacing 1.5
18
19# == CONFIGURACIÓ DEL FIT ==
20set fit quiet
21
22f(x) = a*x + b # Funció a fitar
23fit f(x) datafile u 2:4 via a, b # Fem el fit de les dades
24
25# Això s'usa per obtenir el valor del coeficient de correlació "r", que estarà
26# guardat a la var. "STATS_correlation"
27stats datafile u 2:4 name "STATS" nooutput
28
29# Funció per obtenir els dígits significatius d'una expressió
30f_sd(n, i) = (int(n) == 0 ? f_sd(n*10, i+1) : (int(n) == 1 ? i+1 : i))
31significant_digits(n) = f_sd(n - 10*int(n/10), 0)
32
33significant_digits_r(n) = (n == 1 ? 0 : significant_digits(1 - n))
34
35error_a = STATS_slope_err
36sd_error_a = significant_digits(error_a)
37error_b = sqrt(STATS_slope_err + error_sist_b)
38sd_error_b = significant_digits(error_b)
39r = STATS_correlation
40sd_r = significant_digits_r(r)
41
42print("=== ".datafile." ===")
43print(sprintf('delta(a) = %.'.sd_error_a.'f, sd=%.0f', error_a, sd_error_a))
44print(sprintf('delta(b) = %.'.sd_error_b.'f, sd=%.0f', error_b, sd_error_b))
45
46# Aquesta funció escriu, a partir dels valors del paràmetres a, b i r,
47# l'equació de la regressió:
48title_f(a, b, r) = sprintf('$S(T) = %.'.sd_error_a.'f T + %.'.sd_error_b.'f$, $r = %.'.sd_r.'f$', a, b, r);
49
50plot datafile u 2:4:(0.1):5 t "Dades experimentals" w xyerr, f(x) t title_f(a, b, r)
51
52# == CONFIGURACIÓ DE L'OUTPUT PER SVG ==
53# Això ho uso per generar també una imatge de previsualització que puc carregar
54# a l'ordinador per veure més o menys com a sortit el plot sense haver
55# d'inserir-ho al LaTeX per veure-ho.
56set terminal svg dashed size 600, 600 font "Computer Modern,Tinos,Helvetica,15"
57set output outputfile.'.svg'
58
59replot