Power spectrum of inter-event intervals

spikespec [-F# -## -w# -i -o outfile -l# -x# -c# -V# -h] file(s)

-F maximal frequency (2*l / total time)
-# number of frequencies (F* total time /2)
-w frequency resolution (0)
-i expect intervals rather than times
-l number of values to be read (all)
-x number of values to be skipped (0)
-c column to be read (1 or file,#)
-o output file name, just -o means file_sp
-V verbosity level (0 = only fatal errors)
-h show this message

Computes a power spectrum assuming that the data are the times of singular events, e.g. heart beats. With the flag -i set, the data is taken to be inter-event intervals rather than event times. If the event times are t(n), n=1,...,l, the spectrum is defined by

            l                    2
         | --- 	                |
         | \     -i 2 pi f t(n) |
  S(f) = |  |  e	        |
   	 | /   	                |
         | --- 	                |
           n=1

that is, the signal is taken to be a sum of delta functions at the t(n). S(f) is computed for # frequencies between 0 and -F . The result is binned down to a frequency resolution -w . Output file file_ss (frequency, spectrum).

Table of Contents * TISEAN home