Introduccion a la Programacion en

C++
UNIVERSIDAD
NACIONAL DE
INGENIERIA
Facultad de Ciencias
Abraham Zamudio

GNUPLOT
Gnuplot es un poderoso programa
freeware para hacer graficas con
datos en 2D y en 3D. Gnuplot puede
usarse
en
muchos
ambientes
computacionales, incluyendo Linux,
IRIX, Solaris, Mac OS X, Windows y
DOS. Gnuplot requiere las mnimas
capacidades graficas y puede usarse
aun en una terminal de tipo vt100.
Tiene una amplia variedad de
opciones de salidas para que el
usuario pueda usar las graficas
resultantes como lo desee, ya sea
para ser visualizados o para incluirlos
en sus propios documentos

Graficas basicas en 2D
En sistemas UNIX, si el programa fue
correctamente instalado, simplemente debera
ejecutar: gnuplot
Esto hara que gnuplot inicie. Podra ver un
mensaje de iniciacion que entre otras cosas
menciona la version que se esta usando.

En el paquete de archivos, podra observar uno que se llama

ejem01.dat que es un archivo de texto que tiene numeros en una
sola columna, los numeros empiezan con la serie:
28.062000
52.172000
55.703000
64.281000
43.438000
6.781000
-31.281000
Para graficar estos datos, simplemente debe escribir
gnuplot> plot ejem01.dat using 1

Gnuplot pone la escala mas

conveniente para incluir
todos los datos. Si no se
especifica de antemano algo,
gnuplot dibuja poniendo
pequeos crculos (puntos)

Dibujemos de nuevo estos datos pero

poniendo mas divisiones en el eje X y
modificando la escala
tambien en el eje X.

gnuplot> set style data lines

gnuplot> set xtics 0,50,1000
gnuplot> set xrange [0:500]
gnuplot> plot ../data/ejem01.dat using 1

Hay muchas opciones para el estilo al dibujar los datos, las

opciones son:

Opciones para los estilos del dibujo de los datos. Algunas

opciones requieren listas de datos de 2 o mas columnas
La sintaxis para escribir el estilo del dibujo de los datos es:
gnuplot>set style data <style>
Donde, claro, <style> toma alguno de los valores recien listados.
El comando set xtics tiene 3 argumentos <start>,<increment>,
<end>, y el comando para modificar el rango en X, set xrange,
tiene 2 argumentos, [<start>:<end>].

Y si lo que queremos graficar es una coleccion de datos

multicolumnas como:

gnuplot> plot './ejem02.dat' using 1:2

Es posible dibujar solamente

una columna de esos datos?
Claro que s, dibujaremos
solamente la columna primera
y segunda de los datos del
archivo ejem02.dat

gnuplot>reset
gnuplot>set style data lines
gnuplot> plot ejem02.dat using 1:2

El comando using especifca la columna que sera dibujada, si

los datos estan en 3 columnas, hay muchas posibilidades de
dibujar, por ejemplo al usar... using 3:2 se dibujara la
columna 3 en el primer eje X y la columna 2 en el segundo eje
Y.

Si necesitaramos hacer operaciones con los datos de una o

algunas columnas y esos resultados graficarlos como otra
columna, es posible hacerlo en gnuplot Bueno, hay varias
maneras de hacerlo . Pero, la manera mas facil (si ya se tiene
gnuplot ) es usar esas capacidades que se incorporaron desde la
version 3.7. Vamos a suponer por ahora que el archivo
ejem02.dat contiene los datos de las componenetes x, y, z de la
velocidad que es la que queremos dibujar. La velocidad esta
determinada mediante la expresion

La manera de hacerlo en gnuplot es calcular esa velocidad como:

gnuplot> plot 'ejem02.dat' using (sqrt($1**2+$2**2+$3**2))

Tambien hay que notar que las columnas se especifican con $1,
o $2, o $3 y las potencias con ** . Note que los parentesis que
acotan la expresion son importantes.

Operadores, Constantes y Funciones

Gnuplot no solo lee datos de archivos, tambien grafica funciones
analticas. De manera que gnuplot ofrece los operadores usuales
+, -, *, /, **, etc, y funciones como sin(), cos(), log(), exp(), etc.
Grafiquemos la siguiente funcion

gnuplot> set xrange [0:250]

gnuplot> plot sin(x)*(x**2)

La grafica no parece tan suave como se supone de una funcion

sin, lo que sucede es que gnuplot> evalua la funcion con
solamente algunos puntos (muy pocos para esta funcion), para
cambiar el numero de puntos que se deben graficar en el mismo
rango de x, debemos modificar el numero de muestras:
gnuplot> set samples 1000
gnuplot> replot

Mucho mejor!,
el comando
replot repite el
ultimo comando
plot usado.

Con gnuplot se pueden definir constantes y funciones,

supongamos la siguiente expresion:

Lo que queremos es graficar la funcion con (llamada eps) en el

eje x y (llamada sts) sobre el eje y. La ecuacion es valida en el
rango = [0 :600]. Primero, usaremos el comando reset para poner
los valores iniciales,luego pondremos las condiciones necesarias
en gnuplot para graficar parametricamente (dando valores a
algunos parametros), entonces cambiaremos la variable
independiente t para que sea sts. Como sabemos el rango valido
para sts, podemos cambiar eso tambien. La funcion no esta
escrictamente definida en 0, as que podemos poner un valor muy
pequeo en el lmite inferior. Entonces vamos a definir desde
gnuplot la funcion eps(sts) y las constantes E, Kp, np.

gnuplot> reset
gnuplot> set parametric
gnuplot> set dummy sts
gnuplot> set trange [1.0e-15:600]
gnuplot> eps(sts)=sts/E+(sts/Kp)**(1.0/np)
gnuplot> E = 206000.0
gnuplot> Kp = 1734.7
gnuplot> np = 0.2134
gnuplot> plot eps(sts), sts
Note que gnuplot distingue entre
letras mayusculas y minusculas: E
es diferente que e. Para las graficas
con parametros, el formato de los
comandos para graficar es plot
<x(t)>, <y(t)> donde t es la variable
independiente y x(t) y y(t) pueden
ser cualesquiera funciones de t. Las
constantes se pueden introducir
con formato decimal (600.456) o en
notacion cientfica (1.0e-15). Las
cadenas de caracteres no se
permiten (p.ej. a=tiempo).

En cualquier momento podemos saber que

funciones hemos definido, usando:
gnuplot> show functions
User-Defined Functions:
eps(sts)=sts/E+(sts/Kp)**(1.0/np)
Tambien podemos saber el valor de las funciones
almacenadas, el comando show funciona tambien
para las variables:
gnuplot> show variables
Variables:
pi = 3.14159265358979
E = 206000.0
Kp = 1734.7
np = 0.2134

Frecuentemente se deben comparar unos datos aislados con los

puntos que definen alguna funcion. Este problema se resuelve
muy bien dibujando la curva con lineas y los datos con puntos:
gnuplot> set xrange [0:0.01]
gnuplot> plot eps(sts) , sts with lines , './ejem03.dat' with point

Dando formato a las graficas

Gnuplot tiene varios parametros que se pueden usar para cambiar
la apariencia de la grafica. Esos parametros se pueden accesar
usando el comando set.
Gnuplot ubica los identificadores de las graficas de manera
automatica llamados la clave (key) en la esquina superior
derecha de la grafica. En nuestro ejemplo, la curva analtica pasa
justo encima de la clave ocultando parcialmente las etiquetas de
las lneas, haciendo ligeramente mas difcil leerlas.
Pero podemos cambiar el lugar de la clave usando: set key <x>,
<y>.

Cambiaremos tambien las etiquetas.

gnuplot> plot eps(sts), sts title 'curva analitica' with lines,
'./ejem03.dat' title 'Datos experimentales' with points
Los parametros <x> y <y> del comando set keyse refieren a la
ubicacion del sistema local de coordenadas de la grafica.El
comando title en plot debe escribirse antes del comando with.
Agregar el ttulo de la grafica y las etiquetas de los ejes tambien es
facil.
gnuplot> set title Acero 1045
gnuplot> set ylabel Esfuerzo (MPa)
gnuplot> set xlabel Tiempo (mm/mm)
gnuplot> replot

Los comandos set xlabel y set ylabel toman algunos argumentos

opcionales: set xlabel string <xoffset>, <yoffset>. Los
desplazamientos (offsets) se miden en caracteres. Agragaremos
mas divisiones en el eje X, le pondremos una cuadricula y
moveremos la etiqueta del eje Y mas cerca de la grafica actual.
gnuplot> set xtics -1, 0.001, 1
gnuplot> set yrange [1:1000]
gnuplot> set grid xtics
gnuplot> set grid ytics
gnuplot> set ylabel Stress (MPa) 2, 0
gnuplot> replot

Por default una cuadrcula sera dibujada a lo largo del eje mayor.
La opcion ytics junto con la opcion xtics le dicen a gnuplot que
dibuje lneas tanto horizontales como verticales en cada marca en
los ejes. Hay 3 parametros basicos para los ejes:
(1) set xrange [<x>:<y>], que permite especificar el rango visible
(2) set autoscale que obliga a gnuplot a usar el rango visible
(3) set logscale <x|y> que utiliza una escala logartmica..
A la grafica anterior le pondremos una escala logartmica sobre el
eje Y .
gnuplot> set logscale y
gnuplot> replot

Con gnuplot tambien es posible personalizar la apariencia de las

marcas en los ejes. Primero, desactivaremos la escala logartmica
del eje Y para luego establecer las marcas en el eje Y para que se
lean en notacion exponencial sin dgitos despues del punto
decimal.
gnuplot> unset logscale y
gnuplot> set format y %.0e
gnuplot> replot

El comando format usa la misma sintaxis de printf() del lenguaje C,

as set format x %f wombats es un comando valido. Tambien es
posible especificar etiquetas no numericas. Los parentesis son
necesarios cuando se especifcan etiquetas no numericas.
gnuplot> set xtics ('low' 0 , 'medium' 0.005 , 'high' 0.009)
gnuplot> replot

Si queremos cambiar el tamao de las proporciones a su grafica,

puede usar el comando set size. El siguiente ejemplo crea una
grafica cuadrada.
gnuplot> set size square
gnuplot> replot

Gnuplot permite gran flexibilidad en el formato de los archivos de datos.

Por ejemplo, gnuplot permite lneas de comentarios en los archivos de datos.
Las lneas comentadas se especifcan con un #. En consecuencia los siguientes
archivos son funcionalmente los mismos

# 3 fuentes de datos desde

4.563000 5.078000 -1.719000

# 5 Julio 1997

-5.641000 3.148500 0.156000

# exx exy eyy

4.563000 5.078000 -1.719000

-0.781000 3.453500 -1.250000

-5.641000 3.148500 0.156000

6.860000 8.836000 -2.188000

-0.781000 3.453500 -1.250000

0.953000 5.586000 0.469000,

etc.

6.860000 8.836000 -2.188000

0.953000 5.586000 0.469000, etc.

Esta caracterstica hace que sea posible poner encabezados en

los archivos, para que podamos saber que es lo que esos datos
representan. Hasta ahora los datos que hemos usado estan
separados por un espacio, pero tambien pueden ser serparados
por comas.
4.563000,5.078000,-1.719000
-5.641000,3.148500,0.156000
-0.781000,3.453500,-1.250000
6.860000,8.836000,-2.188000
0.953000,5.586000,0.469000, etc.
plot './ejem04.dat' using 1:2 '%lf,%lf.%lf'
Donde %lf,%lf,%lf especifica el formato del archivo que se va a
leer. La sintaxix es identica a scanf de C.

En ocasiones, se pueden agrupar diferentes conjuntos de datos en

un mismo archivo. Gnuplot le permite graficar solo ciertos
conjuntos en un archivo, no importa que tan grande sean los
conjuntos, deberan estar separados por dos lneas en blanco,
como en el siguiente ejemplo:
1 23.4
2 24.5
3 25.8
4 27.9
1 23.1
2 24.4
3 25.5
4 28.3
1 22.9
2 24.6
3 25.9
4 27.8

Podemos usar el comando index para especificar cual conjunto de

datos queremos graficar. Si quiere graficar solamente el segundo y
el tercer conjunto:
gnuplot> reset
gnuplot> set style data linespoints
gnuplot> set xrange [0:5]
gnuplot> plot './ejem05.dat' index 1:2 using 1:2

Donde el ndice
index 1:2 denota
el conjunto 1
hasta el 2, la
enumeracion de
los
conjuntos
empieza
desde
0.

Graficas en 3D
Gnuplot tambien puede hacer graficas en 3D con
una lista de datos o con funciones analticas
como en el caso de las graficas en 2D, las
graficas son como la que resulta al interactuar
con gnuplot con la siguiente lista de comandos.
gnuplot> set style data lines
gnuplot> set parametric
gnuplot> set view 60, 60, 1, 1
gnuplot> set xlabel x
gnuplot> set ylabel y
gnuplot> set zlabel z
gnuplot> splot u,u+v,sin(0.5*(u+v))

Las graficas en 3D usan splot en lugar de plot, la s es por

superficie. La perspectiva se puede modificar usando los
comandos set view <rot x>, <rot y>,<rot z>, <scale x>,<scale y>,
<scale z>, que sirven para rotar las vistas en los ejes x, y y z y
para escalar cada uno de los ejes. Gnuplot Tambien puede graficar
con lneas ocultas
gnuplot> set hidden3d
gnuplot> splot u,u+v,sin(0.5*(u+v))

Tambien puede graficar los contornos basados en los valores del

eje z (gnuplot no puede dibijar datos en 4D solamente 3D mas
un canal de datos), se puede hacer que se grafiquen los
contornos, dibujados en la superficie, en la base, o en ambos
lugares.
gnuplot> set contour both
gnuplot> splot u,u+v,sin(0.5*(u+v))

Gnuplot puede leer datos desde un archivo de texto.Pero antes

vamos a quitar la opcion de lneas ocultas y dejar los ejes como al
inicio
gnuplot> unset hidden
gnuplot> set style data lines
gnuplot> set view 60,20,1,1
gnuplot> set xtics -3000,1000,3000
gnuplot> set xlabel Eje Axial
gnuplot> set ylabel Esfuerzo
gnuplot> set zlabel Esfuerzo transversal
gnuplot> unset key
gnuplot> splot ejem02.dat

Hasta ahora, hemos interactuado con gnuplot tecleando los

comandos en el prompt. Pero lo mejor y maas facil es interactuar
con gnuplot desde un archivo de texto, editandolo desde nuestro
editor de texto favorito, especialmente si estamos creando una
grafica compleja que poner y quitar parametros y definir variables.
Ahora edite el archivo en nuestro editor favorito. Ahora vamos a
cargarlo en gnuplot .
gnuplot> load script1.gp
Gnuplot ignora las lneas
en blanco y las lneas que
empiezan con #. Todo lo
demas gnuplot lo interpreta
como un comando. Si se
nos olvida algo o hemos
cometido algun error
gnuplot marca el error

Programando en C con salida para Gnuplot

Quizas una de las caractersticas mas utiles para los cientficos (o
para los que algun dia lo seran y por lo pronto deben entregar
tareas con los resultados graficados) es hacer un programa en C y
graficar los resultados. Una solucion rapida y casi inmediata es
generar un archivo de texto llamado nombreArchivo.dat y luego
llamarlo desde gnuplot de la manera ya antes estudiada, con el
comando plot nombreArchivo.dat
Pero la pregunta de cualquier usuario de computadoras sera
puedo de alguna manera hacer eso de manera automatica?, y la
respuesta es claro que s! Una solucion es usar las llamadas
tuberias .

Una tuberia es un proceso que establece un

puente donde se envan mensajes a comandos
externos (del entorno del programa). Para
estudiar de manera completa y profunda el uso
de las tuberas en C, revise cualquier manual de
C.
Enseguida transcribo un programa en ANSI C que
es completamente funcional en las teminales
UNIX compilandolo con gcc

#include <stdio.h>
#include <stdlib.h>
#define GNUPLOT_PATH "/usr/local/bin/gnuplot"
int main (int argc, const char * argv[]) {
(1) FILE *gp;
(2) gp = popen(GNUPLOT_PATH, "w");
if(gp == NULL){
fprintf(stderr, "Error.No se puede encontrar %s.", GNUPLOT_PATH);
exit(EXIT_FAILURE);

}
(3) fprintf(gp, "set term aqua title \"Function A\" \n");
(4) fprintf(gp, "set samples 2048 \n");
(5) fprintf(gp, "plot [-512:512] -abs(x*sin(sqrt(abs(x))))+400");
(6) fflush(gp); /* No hay que olvidar vaciar el buffer
getchar();
(7) pclose(gp);
return(EXIT_SUCCESS);
}

En el programa anterior se muestran algunas lneas enumeradas en el

margen izquierdo, estos numeros no se deben escribir en el programa (y
tampoco los parentesis), solamente son ndices para los siguientes
comentarios.
(1) Se crea un apuntador de tipo archivo, con el identificador gp
(2) Se abre una tuberia en la ubicacion definida por GNUPLOT PATH en
donde se puede escribir "w"
(3) Se escribe en la tuberia
(4) Se escribe en la tuberia
(5) Se escribe en la tuberia. Los tres anteriores puntos escriben la
secuencia de comandos que se pueden ejecutar desde gnuplot ,
note que los cambios de lneas se escriben en el mismo formato que
fprint
(6) Al terminar de usar la tuberia se debe limpiar para liberar la
memoria usada por los mensajes enviados.
(7) Se cierra la tuberia.

#include <stdio.h>
#define GNUPLOT_PATH "/usr/bin/gnuplot -persist"
int main()
{
FILE* gp;
gp = popen(GNUPLOT_PATH, "w");
if (gp == NULL)
return -1;
fprintf(gp, "set isosample 100\n");
fprintf(gp, "min=-1\n");
fprintf(gp, "max=1\n");
fprintf(gp, "pi=3.141592\n");
fprintf(gp, "set hidden3d\n");
fprintf(gp, "set pm3d\n");
fprintf(gp, "set contour\n");
fprintf(gp, "splot [min:max] [min:max] x*x+2*y*y-0.3*cos(3*pi*x)-0.4*cos(4*pi*y)
+0.7\n");
fprintf(gp, "pause -1\n");
fflush(gp);
pclose(gp);
return 0;
}