from gzp import load import numpy as np import pylab import matplotlib matplotlib.use('Agg') """ Loads & plots data from disk """ data, times, breaktime, channels = load("experiment.dat") m, n = data.shape pylab.figure() for i in xrange(len(channels)): pylab.subplot(14, 1, i + 1) pylab.plot(data[i]) pylab.figure() for i in xrange(m): pylab.subplot(14, 1, i + 1) fd = np.hanning(n) * (data[i] - data[i].mean()) fd = np.abs(np.fft.rfft(fd))**2 freqs = np.linspace(0, 64., len(fd)) pylab.plot(freqs[10:], fd[10:]) pylab.xlim(0, 64) pylab.show()
from gzp import load import pylab """ Loads & plots data from disk """ print "this is not working yet!" data, times, breaktime = load("experiment.dat") print data pylab.figure() pylab.plot(times, data) pylab.plot([breaktime, breaktime],[-1e6, 1e6], 'k--') pylab.ylim(min(data), max(data)) pylab.xlabel("time") pylab.ylabel("GSR") pylab.show()
from gzp import load import pylab """ Loads & plots data from disk """ data, times, breaktime = load("experiment.dat") print data pylab.figure() pylab.plot(times, data) pylab.plot([breaktime, breaktime],[-1e6, 1e6], 'k--') pylab.ylim(min(data), max(data)) pylab.xlabel("time") pylab.ylabel("GSR") pylab.show()