#! /bin/bash 
# Last edited on 2006-02-28 13:46:46 by stolfi

cmd="${0##*/}"
usage="${cmd} NAME NSEQS SZSEQ"

# Plots a power spectrum which is read from file "{NAME}.pot". The
# encapsulated Postscript output is written to "{NAME}-pot.eps".
# Assumes that the spectrum was obtained by averaging the power
# spectra of {NSEQS} DNA sequences, with {SZSEQ} samples each.

if [ $# -ne 3 ]; then
  echo "usage: ${usage}"; exit 1
fi

name="$1"; shift;
nseqs="$1"; shift;
szseq="$1"; shift;
echo "plotting ${name}.pot -> ${name}-pot.ps"

gnuplot <<EOF
set terminal postscript eps color "Helvetica" 28
set output "${name}-pot.eps"
set size 1.5,1.5
szseq=${szseq}
nseqs=${nseqs}
frmax=(szseq/2)
pwmax=(2.5+0.012*szseq)
set xlabel "freq"
set xrange [(-0.05*frmax):(1.05*frmax)]
set ylabel "power"
set yrange [0:pwmax]
plot \
  "${name}.pot" using 1:2 \
  title "${name}" \
  with impulses
quit
EOF

gv ${name}-pot.eps