d4r.gnuplot

   1 #!/usr/bin/gnuplot
   2 # vim: set tabstop=4 softtabstop=4 shiftwidth=4 expandtab:
   3 #
   4 # Trabajo Práctico II de Análisis Numérico I
   5 # Genera gráficos con los resultados de las corridas utilizando
   6 # GNU Plot.
   7 # Copyright (C) 2002 Leandro Lucarella <leandro@lucarella.com.ar>
   8 #
   9 # Este programa es Software Libre; usted puede redistribuirlo
  10 # y/o modificarlo bajo los términos de la "GNU General Public
  11 # License" como lo publica la "FSF Free Software Foundation",
  12 # o (a su elección) de cualquier versión posterior.
  13 #
  14 # Este programa es distribuido con la esperanza de que le será
  15 # útil, pero SIN NINGUNA GARANTIA; incluso sin la garantía
  16 # implícita por el MERCADEO o EJERCICIO DE ALGUN PROPOSITO en
  17 # particular. Vea la "GNU General Public License" para más
  18 # detalles.
  19 #
  20 # Usted debe haber recibido una copia de la "GNU General Public
  21 # License" junto con este programa, si no, escriba a la "FSF
  22 # Free Software Foundation, Inc.", 59 Temple Place - Suite 330,
  23 # Boston, MA  02111-1307, USA.
  24 #
  25 # $URL: http://www.llucax.hn.org:81/svn/facultad/75.12/tp2/d4r.gnuplot $
  26 # $Date: 2002-11-25 03:26:50 -0300 (lun, 25 nov 2002) $
  27 # $Rev: 30 $
  28 # $Author: luca $
  29 #
  30
  31 # Seteo terminal para que "dibuje" en un PS.
  32 set term postscript eps enhanced color
  33 set encoding iso_8859_1
  34 set output "d4r.eps"
  35
  36 # Seteos generales.
  37 set title "Estimación del tiempo de reposo para RK4 en el caso sin depósitos y con fricción laminar"
  38 set key right top
  39
  40 # Eje X.
  41 set xlabel  "Tiempo (t)"
  42 set xrange  [0:90000]
  43 set mxtics  5
  44
  45 # Eje Y.
  46 set ylabel "Altura (z)"
  47 set mytics  5
  48
  49 # Primero aproximo por una recta para hallar valores aproximados de los parámetros.
  50 M(x) = a - b * x
  51 m(x) = -c - d * x
  52 fit M(x) 'd3r.max.txt' using 1:2 via a, b
  53 fit m(x) 'd3r.min.txt' using 1:2 via c, d
  54
  55 # Ahora ahora la verdadera aproximación con una exponencial.
  56 M(x) = a * exp( -b * x )
  57 m(x) = -c * exp( -d * x )
  58 fit M(x) 'd3r.max.txt' using 1:2 via a, b
  59 fit m(x) 'd3r.min.txt' using 1:2 via c, d
  60
  61 # Plotea
  62 plot M(x) title "Función aproximante de máximos" with lines linetype 1, \
  63      m(x) title "Función aproximante de mínimos" with lines linetype 8
  64
  65 # Se muestran las funciones aproximantes utilizadas.
  66 print "Funciones aproximantes:"
  67 print "M(x) = ", a, " * exp( -", b, " * x )"
  68 print "m(x) = -", c, " * exp( -", d, " * x )"